【求助】龙格－库塔法求常微分方程组问题？

suck268110

<P>形式如:MX''+CX'+Kx=F,如何利用龙格库塔求解,其中M,C,K都是四维的矩阵，<BR>变形成8个状态方程后，用ode45解方程组，程序运行时间太长，如何减少运行时间？请各位高手们相助。。。<BR>ps：我的机器是p41.5，128m内存。<BR><BR></P>[em06]