龍格庫塔（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方法。