Abstract
Abstract A nonstationary two-body problem is considered such that one of the bodies has a spherically symmetric density distribution and is central, while the other one is a satellite with axisymmetric dynamical structure, shape, and variable oblateness. Newton’s interaction force is characterized by an approximate expression of the force function up to the second harmonic. The body masses vary isotropically at different rates. Equations of motion of the satellite in a relative system of coordinates are derived. The problem is studied by the methods of perturbation theory. Equations of secular perturbations of the translational-rotational motion of the satellite in analogues of Delaunay–Andoyer osculating elements are deduced. All necessary symbolic computations are performed using the Wolfram Mathematica computer algebra system.