A Zero-Dissipative Runge-Kutta-Nyström Method with Minimal Phase-Lag

Figure 2

The error at each integration point when solving the harmonic oscillator with $\omega =8$ with initial condition $q_0=1\text{,}$$q_0^{\prime }=-2$ and $\mathrm{\Delta }t=1/20$.