newton_1.m

追赶法

function y=newton_1(a,n,x0,nn,eps1)
x(1)=x0；
b=1；
i=1；
while(abs(b)>eps1*x(i))
    i=i+1；
    x(i)=x(i-1)-n_f(a,n,x(i-1))/n_df(a,n,x(i-1))；
    b=x(i)-x(i-1)；
    if(i>nn)error(ˊnn is fullˊ)；
        return；
    end
end
y=x(i)；
i