Stability of numerical method as applied to motion in a circular orbit

Abstract

In this paper we discuss the numerical stability of Cowell's method when applied to the Keplerian circular motion. The critical stepsize h_m is given for the PE, PECE and CE algorithms of various orders. A comparison of our results with others' shows our method to be more precise.