newton.m

数值分析

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