Modified multirevolution integration methods for satellite orbit computation

Multirevolution methods allow for the computation of satellite orbits in steps spanning many revolutions. The methods previously discussed in the literature are based on polynomial approximations, and as a result they will integrate exactly (excluding round-off errors) polynomial functions of a discrete independent variable. Modified methods are derived that will integrate exactly products of linear and periodic functions. Numerical examples are given that show that these new methods provide better accuracy for certain satellite problems. It is also shown that information obtained from an approximate analytical solution of the satellite equations of motion, may be used to increase the accuracy and/or efficiency of the multirevolution integration.