newton（计算方法）

最新推荐文章于 2022-10-12 18:34:46 发布

EX_NET

最新推荐文章于 2022-10-12 18:34:46 发布

阅读量2k

点赞数

本文链接：https://blog.csdn.net/zjianbo/article/details/595398

版权

读书期间作业留念专栏收录该内容

30 篇文章 0 订阅

订阅专栏

% 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)