MATLAB牛顿法求解非线性方程组 部分源码 function Newton() x0=[0.1;0.5]; x1=x0-inv(myJacobi(x0))*myfun(x0); while norm(x1-x0)>1e-3 x0=x1; x1=x0-inv(myJacobi(x0))*myfun(x0); end x1 end