龙格库塔（RK4）僵硬的IVP

Question

我需要一个解决僵硬的初始值问题（IVP）的蟒蛇利用RK4最好明确。过去几天我一直在寻找，但找不到它。以下是我的查询：

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

Answer 1

SciPy的似乎是能够做到这一点（http://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.ode.html），如果设置了set_integrator到dopri5你得到一个RK4求解

Answer 2

RK4或经典的龙格 - 库塔方法一个具体的积分方法。作为一个明确的方法，对于僵硬的问题来说，它是非常不适合的。作为单一方法的方法，它没有步长控制的内在特征。

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