用牛顿迭代法求解非线性方程组

用牛顿迭代法求解非线性方程组newton.mclearclc[F,v]=wuzu;x1=input('请输入初值:');eps=input('请输入误差精度:');b=jacobian(F,v);fori=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-b\F';i=1;whilemax(abs(x2-x1))=epsx1=x2;F=wuzu;b=jacobian(F,v);fori=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-b\F';i=i+1;ifi==1000break;endendifi=1000disp('newton迭代法很可能不收敛!')elsedisp('方程的根是：');x=x2Endwuzu.mfunction[F,v]=wuzusymsxyf1=x^2-10*x+y^2+8;f2=x*y^2+x-10*y+8;F=[f1,f2];v=[x,y];End