newton.m

用于求积分

%该算法是用来求解非线性方程的
%fx是由syms定义的一元函数
function Newton(fx,x0,E,N)
%digits(10)
x=0
for n=1:N
    x=x0-subs(fx,x0)/subs(diff(fx),x0);
    fprintf('%f\n',x);
    if abs(det(x-x0))<E
        break;
    end
    if n<N
        n=n+1;
        x0=x;
    else
        break;
    end
end
end