Abstract:
A family of numerical methods for solving the Keplerian problem is proposed. All the methods in this family are symplectic. They preserve the angular momentum, the total energy, the components of the Laplace–Runge–Lenz vector, and the phase volume. The underlying idea is an exact linearization of the problem based on the Levi–Civita transformation and two-stage symmetricsymplectic Runge–Kutta methods.