The axi-symmetric satellite problem including radiation pressure and drag is treated. The equations of motion of the satellite are derived. An energy-like is given for a general drag force function of the polar angle θ, and then it is used to find a relation for the orbit equation of the satellite with initial conditions satisfying the vanishing of arbitrarily choosing higher derivatives of the velocity.
Jezwski, D. and Mittleman, D. (1983) Integrals of Motion for the Classical Two-Body Problem with Drag. International Journal of Non-Linear Mechanics, 18, 119-124. http://dx.doi.org/10.1016/0020-7462(83)90039-2
Leach, P.G.L. (1987) The First Integrals and Orbit Equation for the Kepler Problem with Drag. Journal of Physics A: Mathematical and General, 20, 1997-2004. http://dx.doi.org/10.1088/0305-4470/20/8/019
McMahon, J. and Scheeres, D. (2010) Secular Orbit Variation due to solar Radiation Effects: A Detailed Model for BYORP. Celestial Mechanics and Dynamical Astronomy, 106, 261-300. http://dx.doi.org/10.1007/s10569-009-9247-9
El-Shaboury, S.M. and Mostafa, A. (2014) The Motion of Axisymmetric Satellite with Drag and Radiation Pressure. Astrophysics and Space Science, 352, 515-519. http://dx.doi.org/10.1007/s10509-014-1975-y