t the initial and final parts of the integrationsN±1. The axes units are au. The xy -plane is set to the invariant plane of Solar system total angular momentum.(a) The initial part ofN+1 ( t = 0 to 0.0547 × 10 9 yr).(b) The final part ofN+1 ( t = 4.9339 × 10 8 to 4.9886 × 10 9 yr).(c) The initial part of N−1 (t= 0 to −0.0547 × 109 yr).(d) The final part ofN−1 ( t =−3.9180 × 10 9 to −3.9727 × 10 9 yr). In each panel， a total of 23 684 points are plotted with an interval of about 2190 yr over 5.47 × 107 yr . Solid lines in each panel denote the present orbits of the four terrestrial planets (taken from DE245).
The variation of eccentricities and orbital inclinations for the inner four planets in the initial and final part of the integration N+1 is shown in Fig. 4. As expected， the character of the variation of planetary orbital elements does not differ significantly between the initial and final part of each integration， at least for Venus， Earth and Mars. The elements of Mercury，
本章未完，请点击下一页继续阅读！