global cta1 L1 L2 L3 L4 
omega1=1;L1=20;L2=100;L3=60;L4=120;
cta=[];omega=[];alpha=[];xo=[0.6423,1.2143];
t=0:0.1:2*pi;
for cta1=omega1*t
    xo=fsolve('position',xo,optimset('Display','off'));
    s1=sin(cta1);s2=sin(xo(1));s3=sin(xo(2));
    c1=cos(cta1);c2=cos(xo(1));c3=cos(xo(2));
    A=[-L2*s2,L3*s3,0,0
        L2*c2 -L3*c3 0 0];
    A1=[-L2*c2 L3*c3 0 0
         -L2*c2 L3*s3 0 0];
     b=[s1;-c1];
     b1=[c1;s1];
     w=A\(omega1*L1*b);
     e=A\(-A1*w.^2+omega1^2*L1*b1);
     cta=[cta;xo];
     omega=[omega;w'];
     alpha=[alpha;e'];   
end

subplot(3,1,1);plot(omega1*t*180/pi,cta);axis([0 360 0 3]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角位移曲线');
legend('\fontsize{6}\bf2杆角位移','\fontsize{6}\bf3杆角位移');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角位移\sl{\theta}');

subplot(3,1,2);plot(omega1*t*180/pi,omega);axis([0 360 -1 1]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角速度曲线');
legend('\fontsize{6}\bf2杆角速度','\fontsize{6}\bf3杆角速度');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角速度\sl{\omega}');

subplot(3,1,3);plot(omega1*t*180/pi,alpha);axis([0 360 -2 2]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角加速度曲线');
legend('\fontsize{6}\bf2杆角加速度','\fontsize{6}\bf3杆角加速度');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角加速度\sl{\epsilon}');
function F=position(x)
global cta1 L1 L2 L3 L4
F=[L1*cos(cta1)+L2*cos(x(1))-L3*cos(x(2))-L4;...
    L1*sin(cta1)+L2*sin(x(1))-L3*sin(x(2))];





[size=5][size=2]
为什么会出现下面的错误提示？谢谢帮忙解答！！！[/size][/size]

[color=FF0000]??? Error using ==> feval
Undefined command/function 'position'.

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

请好心的朋友帮忙解答一下
谢谢了