为什么MATLAB的ode45函数和我按照龙格库塔法解出来的结果不一样？

Question

如图，画出来的图不一样。求大佬指教呜呜

Answer 1

为什么MATLAB的ode45函数和我按照BK4（龙格库塔法）解出来的结果不一样？其根本的原因是时间区间取了太大而造成，因为该微分方程超过【0，1.6】后，ode45法和BK4法得到的图像就失稳。由于ode45法和BK4法虽然都是龙格库塔法的思想，但ode45法是采用变步长进行计算，而BK4法是采用定步长进行计算，所以随着时间区间范围越来越大，其就不完全一样了。

下图为时间区间【0，1.6】时的图像。