《用牛顿迭代法求解非线性方程组.doc》由会员分享,可在线阅读,更多相关《用牛顿迭代法求解非线性方程组.doc(2页珍藏版)》请在三一文库上搜索。
用牛顿迭代法求解非线性方程组newton.mclearclcF,v=wuzu;x1=input(请输入初值:);eps=input(请输入误差精度:);b=jacobian(F,v);for i=1:1:length(x1)b=subs(b,v(i),x1(i);F=subs(F,v(i),x1(i);endb=eval(b);F=eval(F);x2=x1-bF;i=1;while max(abs(x2-x1)=eps x1=x2; F=wuzu;b=jacobian(F,v);for i=1:1:length(x1)b=subs(b,v(i),x1(i);F=subs(F,v(i),x1(i);endb=eval(b);F=eval(F);x2=x1-bF; i=i+1; if i=1000 break; endendif i=1000 disp(newton迭代法很可能不收敛!)else disp(方程的根是:);x=x2Endwuzu.mfunction F,v=wuzusyms x y f1=x2-10*x+y2+8;f2=x*y2+x-10*y+8;F=f1,f2;v=x,y;End