ar momentum, which should be rigorously preserved up to machine-ε precision.
2.4.2 Error in planetary longitudes
Since the symplectic maps preserve total energy and total angular momentum of N-body dynamical systems inherently well, the degree of their preservation may not be a good measure of the accuracy of numerical integrations, especially as a measure of the positional error of planets, i.e. the error in planetary longitudes. To estimate the numerical error in the planetary longitudes, we performed the following procedures. We compared the result of our main long-term integrations with some test integrations, which span much shorter periods but with much higher accuracy than the main integrations. For this purpose, we performed a much more accurate integration with a stepsize of 0.125 d (1/64 of the main integrations) spanning 3 × 105 yr, starting with the same initial conditions as in the N−1 integration. We consider that this test integration provides us with a ‘ps
本章未完,请点击下一页继续阅读!