Точное сохранение инвариантов движения в численном решении задачи Кеплера

A correction algorithm has been developed to preserve required motion invariants at each step of numerical integration of gravitational N-body problem, independent of the underlying numerical integration method. The effectiveness of the proposed algorithm is demonstrated using the Kepler problem as an example. Specifically, it is shown that the algorithm extends the applicability range of classical integrators of lower order