Runge-Kutta（RK4）用于僵硬的IVP

Question

我需要在python中利用RK4中的僵硬的初始值问题（IVP）求解器，最好是显式的。过去几天我一直在寻找但是找不到它。以下是我的疑问：

• 解算器，即任何模块，是否存在？
• 如果不是，编码一个是否合理？我问这个是因为我找不到任何关于使用RK4-explicit来解决僵硬问题的参考。

Answer 1

如果你将set_integrator设置为dopri5你得到一个RK4求解器，那么scipy似乎能够做到（http://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.ode.html）

Answer 2

RK4或经典的Runge-Kutta方法是一种特定的整合方法。作为一种明确的方法，它很容易适用于僵硬的问题。作为一种方法，它没有步长控制的内在特征。

对于僵硬问题，您需要具有步长控制的隐式RK方法。