MATLAB程序设计：牛顿迭代法

function x=newton(x0,e,N,fx)
%输入x0,误差限e,迭代次数N和函数Fx
k=1;
while k<N
    if subs(diff(fx),x0)==0
        disp("输出奇异标志");
        break;
    end
    x1=x0-subs(fx,x0)/subs(diff(fx),x0);
    if abs(x1-x0)<e
        break;
    end
    x0=x1;
    k=k+1;
end
if k<N
    disp(x1);
else
    disp("迭代失败");
end
end