求解微分方程初值问题

Question

用matlab仿真，ode45命令调用m文件，matlab计算了一个小时还没有计算出来，这是怎么回事，方程只有六个。大家有没有遇到过这样的问题？是怎么解决的？

Answer 1

m文件如下：%%%%%%%%%%%%%%%%%%%%%%%%%%%%function xdot=xiaxing2(t,x)c2=5e-11;r4=2.25e+10;m1=0.067;c7=0.22e-13;m2=0.068;c12=4e-12;r20=2.781e+10;c18=2.5e-10;xdot=[1.97e-3-x(1)/(r4*c2)+x(2)/(r4*m1*c7)+x(3)/(r4*c12*m1*m2) x(1)/(r4*c2*m1)-x(2)/(r4*c7*m1^2)-x(3)/(r4*c12*m2*m1^2) x(1)/(r4*c2*m1*m2)-x(2)/(r4*c7*m2*m1^2)-x(3)*(1/(r4*m1^2*m2^2)+1/r20)/c12+x(4)/(r20*c18) x(3)/(r20*c12)-x(4)/(r20*c18)];%%%%%%%%%%%%%%%%%%%%%%%%%%%%%用ode45命令调用，仿真时间：0～100，初始值为：［0，0，0，0］，这个还不是很复杂，计算机算了近一个小时才出现结果。请问这是怎么回事？

Answer 2

有木有高手指教啊～～