[shw1220]

m文件为：
function dydt=vdt1(t,y)
dydt=[(7.7794e-14)*360000*sinh(y(2)/21.7455);
y(2)*(7.7794e-14)*360000*sinh(y(2)/21.7455)];
计算命令为：
[t,y]=ode45(@vdt1,[0 100000],[0 150]);
plot(t,y(:,1))
图像显示计算到t为9000多的时候就中断了，改成ode23s和ode15s都显示Unable to meet integration tolerances without reducing the step size below the smallest value allowed at time t.
本人新手，最近需要计算这个方程组才翻开matlab的，诚请各位帮助，谢谢！