我找到了RK-Dopri(4,5)的玩具实现，用 Python 编写。
但是，我担心第 118 行：
y = y + h * (b1*K1+b3*K3+b4*K4+b5*K5+b6*K6)
b2*K2是否故意省略了该术语？在 Hairer、Wanner 和 Norsett 中给出（1.8），我有一种感觉，它不是：
你可否确认？
RK-Dopri(4,5) 中的最终结果是如何计算的？
计算科学
龙格库塔
教育
2021-12-19 20:38:38
1个回答
在Dormand-Prince的 Butcher 画面中为零，因此您可以跳过该术语。
请注意，您链接到的代码还包含行b2 = 0.0, b2p = 0.0，它们只是作为提醒，因为这些变量从未使用过。
其它你可能感兴趣的问题