Research of efficiency of algorithms of method Everhart with high order of approximating formulas

The modified algorithm for the numerical integration of the equations of celestial motion by the Everhart's method is developed. The study of the effectiveness of the algorithm for approximating high-order formulas is carried. High efficiency of the method is shown on the example of joint integration of the equations of motion of major planets, the Moon, the Sun and the small bodies of Solar system.