有没有高阶辛方法是的'= f(是)y′=f(y)?

计算科学 保护
2021-12-05 19:04:46

是否有高阶能量守恒或辛方法求解y=f(y)?

1个回答

Ernst Hairer & co 编写了几何积分:

  • E. Hairer、C. Lubich 和 G. Wanner (2002):几何数值积分。常微分方程的结构保持算法。计算机中的 Springer 系列。数学,卷。31

  • E. Hairer 和 M. Hairer (2002):GniCodes - 用于几何数值积分的 Matlab 程序

您可以从Hairer 教授的网站下载 FORTRAN77、Matlab 或 C++ 中的例程他实施的隐式龙格-库塔 (IRK) 方案的顺序为 8。