Discrete hill's equations

In the work, based on the formulas of discrete mechanics in a rotating frame, a discretization of classical Hill's equations of the moon motion is worked out. It is proved that the proposed discretization conserves Jacobi's constant of motion. For computational purposes an algorithm to the solution of obtained discrete Hill's equations is given.