High-Order Analytic Expansion of Disturbing Function for Doubly Averaged Circular Restricted Three-Body Problem

Figure 7

Variation of , , , and of representative bodies in Figures 6(a)–6(f). Initial values of and for each body are as follows: (a) , (b) , (c) , (d) , (e) , and (f) . Note that the unit of and is degree, not radian, for consistency with the unit of . All the elements are plotted with a time interval of 1500 years and are calculated on the heliocentric coordinate centered on the central mass.

(a)

(b)

(c)

(d)

(e)

(f)