An orbit-preserving discretization of the classical Kepler problem

We present a remarkable discretization of the 3-dimensional classical Kepler problem which preserves its trajectories and all integrals of motion. The points of any discrete orbit belong to some exact continuous trajectory.