Conservative discretizations of the Kepler motion

Modified vector fields are used to construct high-order conservative discretizations of the three-dimensional Kepler motion. The numerical integrators preserve the Hamiltonian function, the angular momentum and Runge-Lenz vector. In particular, the exact integrator of the Kepler motion is found. The proposed numerical schemes permit explicit implementation