请教解微分方程组的问题
各位好！我现在想请教的问题是，解2阶微分方程组的问题，我的方程是这样的，
M*D2X+C*DX+K*X=F;其中M是一个3*3的对角阵，C和K是一个3*3的对角阵，F是一个3*1的一列，
例如：M=；C=；K=；F=
这样的一个2阶方程组应该怎么解呢，
我也想把它划成一阶微分方程组，可划完之后还是不会往下做了，希望各位指点迷津 应该有不少方法吧，子空间迭代法呢？ 非常感谢您！
子空间迭代法,不太懂这个方法，能说的具体点不？
我在matlab中自己这样做的，可是出问题了function y=vdp_eq1(t,x,flag,k,c)
y=[eye(3)*x(2);
-k*x(1)-c*x(2)];
这是我编写的函数，x0=;k=;c=;=ode45('vdp_eq1',,x0,h_opt,k,c)，在matlab空间下运行这句话，就出现了这样的错误，??? Error using ==> funfun\private\odearguments
VDP_EQ1 must return a column vector.
Error in ==> ode45 at 171
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...
能给指点一下这是怎么回事啊？
页:
[1]