The generalized Sundman transformation for propagation of high-eccentricity elliptical orbits

A generalized Sundman transformation dt = cr^n ds for exponent n >= 1 may be used to accelerate the numerical computation of high-eccentricity orbits, by transforming time t to a new independent variable s. Once transformed, the integration in uniform steps of s effectively gives analytic step variation in t with larger time steps at apogee than at perigee, making errors at each point roughly comparable. In this paper, we develop techniques for assessing accuracy of s-integration in the presence of perturbations, and analyze the effectiveness of regularizing the transformed equations. A computational speed comparison is provided.