برنامه به زبان برنامه نویسی Matlab :

Function newton(f,x0,tol,n)

%Newton’s.method for solving the nonlinear’

%equation

%f(x)=0

iter=0

u=feval(f,x0);

v=feval(df,x0);

er=abs(u/v);

disp(‘                       ‘)

disp(‘iter x f(x) df(x) xn+1-xn’)

disp(‘                       ‘)

fprintf(‘%2.0f % 12.6f %12.6f %12.6f \n’,iter,x0,u,v)

while(err>tol)&(iter<=n)&(v=1)

x1=x0-u/v;

err=abs(x1-x0);

x0=x1;

u=feval(f,x0);

v=feval(df,x0);

iter=iter+1;

fprintf(’%2.0f %12.6f %12.6f %12.6f/n’,iter,x0,u,v,err)

end

if(v==0)

disp(‘division by zero’?)

end

if(iter>n)

disp(‘method failed to converge’)

end

خروجی مسئله به روش نیوتن :

X^4-4X^3-x^2+4X-10=0             f(x)

فرمول      Xn+1=Xn-  

f(x)َ               بادقت ماشین حساب∑=

         X^4n-4X^3n-X^2n+4Xn-10

                                   Xn+1=Xn-

         4X^3n-12X^2n-2Xn+4

N=0,1,2,…          X0=1                 

  با دقت ماشین حساب∑=

α≈X13=-1.490591294 f(X13)=0.000000000       

تهیه کننده : جعفر فتحی عزیزآباد