nt.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 11 行

function [y,n]=NT(x)
x1=x-ff(x)/dff(x);     %第一次切线选代计算
n=1;
while (abs(x1-x)>=1e-6)&(n<100000)
    x=x1;
    x1=x-ff(x)/dff(x);      %切线选代求解
    n=n+1;
end
y=x;