newton.m

来自「数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序」· M 代码 · 共 15 行

function x=Newton(fname,dfname,x0,e,N)
%用途：Neton迭代法解非线形方程f(x)=0;
%fname和dfname分别表示f(x)及其导函数的M函数句柄或内嵌函数表达式
%x0为迭代初值，e为精度
%x为返回数值解，并显示计算结果，设置迭代次数的上限
if nargin<5,N=500;end
if nargin<4,e=1e-4;end
x=x0;x0=x+2*e;k=0;
fprintf('It.no=%2d      x[%2d]=%12.9f\n',k,k,x)
while abs(x0-x)>e&k<N,
    k=k+1;
    x0=x;x=x0-feval(fname,x0)/feval(dfname,x0);
    fprintf('It.no=%2d     x[%2d]=%12.9f\n',k,k,x)
end
if k==N,fprintf('已达到迭代次数上限');end