Adaptive symplectic conservative numerical methods for the Kepler problem

We suggest and substantiate a unified form of a family of adaptive conservative numerical methods for the Kepler problem. The family contains methods of the second, fourth, and sixth approximation order as well as an exact method. The methods preserve all the global properties of the exact solution of the problem. The variable time step is chosen automatically depending on the properties of the solution.