我需要一个在 python 中利用 RK4 的僵硬初始值问题 (IVP) 的求解器，最好是显式的。我一直在寻找过去几天，但找不到它。以下是我的查询：
- 求解器（即任何模块）是否存在？
- 如果不是，编码一个是否合理？我问这个是因为我找不到任何使用 RK4-explicit 来解决棘手问题的参考资料。
RK4 用于刚性 IVP
计算科学
刚性
2021-11-30 04:13:25
2个回答
您可以在Mathworld上查看显式 RK4 的公式。鉴于表达式的简单性，为自己编写一个代码是非常合理的，这就是您可能找不到模块的原因。
但是，如果如您所说，您正在求解一个刚性 IVP，那么显式 Runge-Kutta 方法是不合适的（绝对稳定区域是有界的）。也许您应该研究隐式 Runge-Kutta 方法，因为存在任意高阶的 A 稳定方法。
您可以在此页面中找到许多具有多种语言的 RK4 代码
其它你可能感兴趣的问题