newton（计算方法）

% newton.m
%x=newton(inline('x^3-2*x-5'),inline('3*x^2-2'),0.5,100,1e-6)

%function [y,z]=newton(fv,df,x0,n,tol)
function [y]=newton(fv,df,x0,tol)
x(1)=x0; b=1; k=1;
while or(k==1,abs(b)>tol*abs(x(k)))
    x(k+1)=x(k)-feval(fv,x(k))/feval(df,x(k));
    b=x(k+1)-x(k);
    k=k+1;
    %if(k>n)
     %   error('Error: Reached maximum iteration times');
     %end
end
y=x(k-1);
%if nargout>1
%    z=k-1;
end

//   newton_run.m  /
clear all
%x=newton(inline('x^3-2*x-5'),inline('3*x^2-2'),0.5,100,1e-6)
%[xx,k]=newton(inline('x^3-2*x-5'),inline('3*x^2-2'),0.5,100,1e-6)

x=newton(inline('x^3-2*x-5'),inline('3*x^2-2'),0.5,1e-6)