newton.m

用newton迭代法求方程根的源代码 给出了newton迭代法的实现

function [x_star,index,it] = Newton(fun,x,ep,it_max)
it_max=1000;ep=1e-5;
index= 0; k=1;
while k <= it_max
    x1 = x; 
    f = feval(fun,x);% f(1)时，调用第一个方程；f(2)时，调用第二个方程
    if abs(f(2))< ep break; end
    x = x - f(1)/f(2);
    if abs (x-x1) < ep index=1; break ; end
    k= k + 1;
end
x_star = x; it = k;