On the numerical integration of Hill's problem

Abstract

We present explicit high order composition methods based on a second order symmetric method to numerically integrate Hill's lunar problem. The linear/nonlinear splitting of the non-separable Hamiltonian allows us to build a class of integrators that are simple to use and efficient in comparison with other standard symplectic methods. Our numerical results show that the methods preserve the energy very well in long time integration.