用以下程序解非线性方程组,如何查看解的迭代计算过程。
function x=myfun1
format long
x0 = [30;5;0;2000000];
options=optimset('MaxFunEvals',80000,'MaxIter',85000);
[x,fval] = fsolve(@f,x0,options) ;
function F=f(x)
k=5000
ga=10322.87553
ei=178892900.6
p=60000
c=7000
h=173
ae=0
ea=0
a=0.0478
F =[k*x(2)-ea*(((h-ae-x(2))^2+(x(1))^2)^0.5-h)*abs(cos(x(3)))/h;
(abs(x(1))-(h-ae-x(2))*abs(tan(x(3))))*100;
x(4)*(1-cos(a*h)-c*h)/p+(p*(c-k*x(2)*abs(tan(x(3))))+ga*c)*sin(a*h)/(ga*a*p)-abs(x(1));
x(4)-c*h-p*abs(x(1))]
end