Diagonally Implicit Symplectic Runge-Kutta Methods with High Algebraic and Dispersion Order

Figure 1

Stability region of the M968 method.