请教关于matlab中求解微分方程数值解的ode函数问题！

sky8754180

本人最近在做柔性体的动力学方程求解，最终得到的动力学方程形式为：M(t,y)y'=f(t,y)，其中，质量矩阵M(t,y)不是常量，是时间t和广义坐标y的函数（当然y也是t的函数），刚开始用ode45计算，算了一天一夜没结果，但也没有报错，不知道什么问题？难道是计算量太大的原因，我这儿的M为8*8的方阵，并且右上方的4*4子方阵和左下方的4*4方阵为零矩阵，论坛里有没有哪位朋友做过类似的计算，望指教！

另外，我也考虑用ode15s试试，不知道情况是否有所变化？

顺便问一下，用matlab中的ode求解微分方程数值解和用C语言或Fortran等语言编写计算程序，调用IMSL现成的微分方程求解函数，这另种方法的计算效率差别大吗？（假设所用的算法均为Gear's method）

谢谢！