Numerical solution of modified-Kepler problem using a splitting method

An efficient geometric integrator is proposed for numerically solving the Kepler problem modified by an additional potential energy term. This method is stable over long integration times, uses relatively large timesteps and gives a good estimate of the global error. Two case studies are presented to demonstrate these advantages.