To predict the motion of spacecrafts, a numerical-analytical method for integrating the differential equation of the orbital motion of a spacecraft stabilized by the Baumgart differential method is proposed. The stabilization of the differential equation of motion by the Baumgart method is carried out according to the energy of the spacecraft. Stabilization is carried out to reduce the influence of the Lyapunov instability on the accumulation of numerical errors in the integration of the differential equation, which is effective when conducting a long-term numerical prediction of the motion of spacecraft. Integration of the stabilized equation is based on differential Taylor transformations. Computational schemes with a constant step and an integration order are considered, as well as schemes with adaptation by an integration step and order. For adaptive schemes, the results of forecasting the motion of spacecraft according to the criterion “accuracy-computational complexity» for a given relative error of integration with respect to integration phase variables and spacecraft energy are presented. It is shown that both options require setting various internal adaptation parameters, but they have comparable efficiency. Recommendations are proposed on the use of the developed method for integrating energy-stabilized equations for predicting the motion of spacecraft in the near space in the Greenwich rectangular coordinate system.