微分方程组求解遇到问题，求帮助
以下时.m文件function fde=def(t,x)
fde=zeros(12,1);
m0=1560;
J=m0/3*;
c=0.2*sqrt(300*11*10^6);
k=;
xdamper=[-0.6;-0.6;0;0;0.6;0.6];
ydamper=[-0.3;0.3;-0.3;0.3;-0.3;0.3];
F=;
x0=;
m=0;n=0;
for i=1:6
m=m+xdamper(i)^2;
n=n+ydamper(i)^2;
end
fde=
-6*k(1)*x(1)/m0
x(4)
(F(1)-6*k(2)*x(3))/m0
x(6)
(F(2)-6*k(3)*x(5)-6*c*x(6))/m0
x(8)
1/J(1)*(-c*x(8)*n-k*x(7)*n-F(1)*(x0(3)+x(7)*x0(2)-x(9)*x0(1))+F(2)*(x0(2)+x0(1)*x(11)-x(7)*x0(3)))
x(10)
1/J(2)*(-c*x(10)*m-k*x(9)*m-F(2)*(x0(1)-x0(2)*x(11)+x0(3)*x(9)))
x(12)
1/J(3)*(F(1)*(x0(1)-x0(2)*x(11)+x0(3)*x(9))-k*x(11)*(m+n))];
end
M文件似乎没显示有错误。但是调用的时候出现问题，调用程序如下：
tsp=;
xt0=;
=ode45(@def,tsp,xt0);
plot(t,x(:,1),'r-')
目的是求振动模型的振动过程，六个自由度。错误显示如下：
??? Error using ==> odearguments at 116
Solving DEF requires an initial condition vector of length 12.
Error in ==> ode45 at 173
本人是在是看不出错在哪里，昨天编了个简单一点的通过了，这是今天改编的复杂一点的。
已自行解决问题，.m文件中有k，写错了，应该是k(3)、k(2)、k(1)写法
页:
[1]