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

Figure 5

The global error at each integration point when solving the “almost” periodic problem (Problem 3) with .