Phase-Fitted and Amplification-Fitted Higher Order Two-Derivative Runge-Kutta Method for the Numerical Solution of Orbital and Related Periodical IVPs

Figure 3

The error at each integration point when solving the harmonic oscillator (Problem 1) with with initial condition , , and with .