newton.m

来自「数值分析」· M 代码 · 共 12 行

function y=Newton(x0,a)
y=x0-f(x0)/df(x0);
n=1;
while(abs(y-x0)>=a)&(n<=100000)
    x0=y;
    y=x0-f(x0)/df(x0);
    n=n+1;
    n
    y
end