龙格-库塔(RK4)治疗硬IVP

Question

我需要一个刚性信息值问题(IVP)的求解者在python中开发RK4，最好是显式的。我已经找了几天了，但找不到。以下是我的问题：

• 求解器，即任何模块，是否存在？
• 如果没有，那麽编码是否合理呢？我问这个问题是因为我找不到任何使用RK4-显式的刚性问题的参考。

Answer 1

(http://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.ode.html)如果您将set_integrator设置为dopri5，则会得到一个RK4求解器

Answer 2

RK4或经典的Runge方法是-one专用积分法.作为一种显式方法，它不适用于刚性问题。作为一种方法，它没有固有的特点来控制步长.

对于刚性问题，您希望隐式RK方法具有步骤大小控制。