Numerical orbit integration efficiency of the Delaunay-Similar elements

Orbit equations with a set of conservative and a set of nonconservative perturbing potentials were considered. Scheifele's DS formulation of these equations has dependent variables similar to Delaunay's orbital elements with the true anomaly as the independent variable. Efficiency curves of computing cost v.s. accuracy were constructed for Adams integrators of order of 2 through 15 with several correcting algorithms and for a Runga-Kutta integrator. Considering stability regions, choices were made for the optimally efficient integration modes for the DS elements. Integrating in these modes reduces computing costs for a specified accuracy