A Trigonometrically Fitted Block Method for Solving Oscillatory Second-Order Initial Value Problems and Hamiltonian Systems

Figure 7

Perturbed Kepler problem: the logarithm of the global error of the Hamiltonian against for and the momentum are presented in (a) and (b), respectively. In (c) we compare the efficiency curves for the BHTRKNM and N4. Timing comparison is provided in (d). It is clear from the timing curves that BHTRKNM is very efficient.

(a) ;

(b) ;

(c) Efficiency curves for BHTRKNM and N4

(d) Timing curves for BHTRKNM and N4

(e) BHTRKNM solution,

(f) N4 solution,