F and G Taylor series solutions to the Stark and Kepler problems with Sundman transformations

The classic $F$ and $G$ Taylor series of Keplerian motion are extended to solve the Stark problem and to use the generalized Sundman transformation. Exact recursion formulas for the series coefficients are derived, and the method is implemented to high order via a symbolic manipulator. The results lead to fast and accurate propagation models with efficient discretizations. The new $F$ and $G$ Stark series solutions are compared to the Modern Taylor Series (MTS) and 8th order Runge–Kutta–Fehlberg (RKF8) solutions. In terms of runtime, the $F$ and $G$ approach is shown to compare favorably to the MTS method up to order 20, and both Taylor series methods enjoy approximate order of magnitude speedups compared to RKF8 implementations. Actual runtime is shown to vary with eccentricity, perturbation size, prescribed accuracy, and the Sundman power law. The method and results are valid for both the Stark and the Kepler problems. The effects of the generalized Sundman transformation on the accuracy of the propagation are analyzed. The Taylor series solutions are shown to be exceptionally efficient when the unity power law from the classic Sundman transformation is applied. An example low-thrust trajectory propagation demonstrates the utility of the $F$ and $G$ Stark series solutions.