An Encke-type special perturbation method

In this paper, a combination analytical-numerical integration method for solving the differential equations of a modified set of Lagrange's planetary equations is described. The integration method is an Encke-type method because it involves integrating the deviations between the actual trajectory and a reference trajectory. The reference trajectory is obtained from an analytical solution containing the dominant secular and periodic effects of the gravitational field of the primary body. A set of nonsingular elements is used so that the method will be valid for all circular and elliptical motions. It is shown that the method is an accurate and efficient means of satellite ephemeris generation.This paper was presented at the AIAA/AAS Meeting, Princeton University, August 1969.