求助大侠：非线性方程组求解出现问题。

东南小忽悠

Ｍ文件　如下：
function x=mzzzfun1
format long
x0 = [0;5;0.3;4000];
options=optimset('MaxFunEvals',330,'MaxIter',7000);
[x,fval] = fsolve(@f,x0,options) ;
function F=f(x)
k=3000
ga=10322.87553
ei=178892900.6
p=60000
c=0
h=173
ae=0
ea=95600000
a=0.0478
F =[k*x(2)-ea*(((h-ae-x(2))^2+(x(1))^2)^0.5-h)*cos(x(3))/h;
x(1)-(h-ae-x(2))*tan(x(3));
x(4)*(1-cos(a*h)-c*h)/p+(p*(c-k*x(2)*tan(x(3)))+ga*c)*sin(a*h)/(ga*a*p)-x(1);
(ga+p)*f*cos(a*h)/ga/a/p+(k*x(2)*tan(x(3))*(1-cos(a*h))*p-(ga+p)*f+a*ga*x(4)*sin(a*h))/(ga+p)/p]


end

end

运行出现一些错误，（原先能运行的，把五元方程改为四元反而不行了，其他语句都一样），哪位高手帮忙检查一下，谢谢啦！
Input argument "x" is undefined.

Error in ==> mzzzfun1>f at 16
F =[k*x(2)-ea*(((h-ae-x(2))^2+(x(1))^2)^0.5-h)*cos(x(3))/h;

Error in ==> mzzzfun1>f at 16
F =[k*x(2)-ea*(((h-ae-x(2))^2+(x(1))^2)^0.5-h)*cos(x(3))/h;

Error in ==> fsolve at 180
        fuser = feval(funfcn{3},x,varargin{:});

Error in ==> mzzzfun1 at 5
[x,fval] = fsolve(@f,x0,options) ;
　　

ChaChing

回复 1 # 东南小忽悠 的帖子

Ref 4F, 常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html
(from http://forum.vibunion.com/home-s ... -blog-id-18250.html)