The singly averaged differential equations of satellite motion for e greater than or equal to 0 and less than 1

The singly-averaged differential equations of motion of a satellite are developed in terms of parameters valid for all eccentricities less than one. The perturbations included in the acceleration model are due to an aspherical central planet (zonal harmonics up to degree 20 and resonant harmonics up to degree and order 20), atmospheric drag for a time-varying atmosphere, third-body gravity (the sun and moon for an earth satellite), solar radiation pressure with shadowing, and impulsive maneuvers. Analytic averaging is used to remove short-period terms due to the aspherical central planet and third-body gravity. Numerical averaging is used to remove short-period terms due to atmospheric drag and solar radiation pressure.