newton.m

来自「这是一些关于MATLAB的小程序」· M 代码 · 共 14 行

function x = newton(fname,dfname,x0,e,N)
% 用途：Newton迭代法解非线性方程f(x) = 0
% 格式：x = newton(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