如何用MATLAB模拟傅科摆运动轨迹

642712189

其中地球重力加速度g=9.8。绳长67。w为地球自转角速度
程序如下

a=input(‘输入纬度=‘）；
q=input(‘按次格式依次输入[x0,vx0,y0,vy0]=')
c=a*pi/;
[t.x]=ode45(‘fkbfun’，[0:0.02:100],q,[].c);
xlabel(‘x’)
ylabel(‘y’);
comet(x(:,1),x(:,3))
函数文件是一个独立的文件，文件名为 fkbfun.m.
function  tt=fkb(t,x,flag,c)]
a=(2*pi*sin(c))/100;  (放大了地球的角速度）
b=9.8/67;
tt=[x(2);
2*a*x(4)-b*x(1);
x(4);
-2*a*x(2)-b*x(3)];

为什么运行错误呢   求高手讲解！