什么ODE求解器在步进函数中使用计算进行插值？

Question

我对具有不同初始条件的多个ODE解决方案求平均值，因此所有解决方案在同一时间具有值是非常重要的;例如，增量为0.01。

我一直在使用数字配方3（nr3）中的ODE例程。它们执行自适应大小步骤并使用计算值执行相同的插值顺序。我不能使用它们因为它们与提升冲突。有没有其他类似的例程？

我看了GSL，它非常好但它没有内置插值。我能做到的一种方法是用自适应大小求解ODE，然后运行Akima插值。但似乎nr3解决方案会更快，更准确。

Answer 1

我使用了来自http://www.unige.ch/~hairer/software.html的DOPRI5，密集输出=插值。我觉得它很可靠。我使用的是原始版本（在Fortran中）;同一个网页上还有一个我没有用过的C版本，但我似乎记得人们对它很满意。

Answer 2

您可以使用odeint。它有Dopri5，Rosenbrock4和Burlish-Stoer用于密集输出。