Abstract
Based on the results of the article, perturbations of the first and second order are determined in rectangular coordinates and components of the speed of the body. The constructed system of differential equations of the perturbed motion of the body makes it possible to create a unified algorithm for determining perturbations of the second and higher orders in the form of finite polynomials in powers of regularizing variables, which are selected at each stage of the approximation and with the help of approximating polynomials it is possible to calculate any intermediate point of the trajectory of the motion of the body. It was proved that the introduction of regularizing variables provides the representation of the right-hand sides of the system of differential equations of perturbed motion in the form of polynomials with respect to this regularizing variable. A procedure for reducing the degree of approximating polynomials using Chebyshev polynomials was also performed.