nanewton.m

Matlab数值解算法实现代码

function x=nanewton(fname,dfname,x0,e,N)
%用途：牛顿法解非线性方程f(x)=0
%格式：x=nanewton(fname,dfname,x0,e,N)fname和dfname分别为表示f(x)及其导函数的M函数句柄或内嵌函数，
%      x0为迭代初值，e为精度要求(默认1e-4)，x返回数值解,设置迭代次数上限N以防发散(默认500)。
if nargin<5, N=500;end
if nargin<4, e=1e-4;end
x=x0;x0=x+2*e;k=0;
while abs(x0-x)>e&k<N,
   k=k+1;
   x0=x;x=x0-feval(fname,x0)/feval(dfname,x0);
   disp(x)
end
if k==N, warning('已达迭代次数上限');end