s， variation of delaunay elements and angular momentum deficit， and results of our simple time–frequency analysis on all of our integrations.
in section 2 we briefly explain our dynamical model， numerical method and initial conditions used in our integrations. section 3 is devoted to a description of the quick results of the numerical integrations. very long-term stability of solar system planetary motion is apparent both in planetary positions and orbital elements. a rough estimation of