Secular Variations in the Restricted Problem: A Numerical Technique for Studying Stability of Equilibrium Solutions

In this paper we develop and implement an explicit numerical technique for studying stability of equilibrium solutions concerning the secular variations in the restricted problem. In our implementation, we use fourth-order expansions for the secular terms, but the method itself is independent of the particular order used as upper limit in the required expansions.