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

Table 2

Stability test for PFAFTDRK4() using coefficient of with for variable .

		Global error
4.55	2.3772841250
4.35	1.5452837990
4.15	0.9838244901
1.00	0.3678739259
0.15	0.8607079765
0.01	0.9900498337