We can write

The perturbing force causes the planet's orbital energy per unit mass to evolve in time as

(I.16) |

where use has been made of Equations (I.3) and (I.15). Thus, it follows from Equations (I.5), (I.7), (I.9), and (I.10) that

The perturbing force causes the planet's orbital angular momentum per unit mass to evolve in time as

where Equations (I.2), (I.6), (I.14), and (I.15) have been employed. Consequently,

Furthermore,

where use has been made of Equation (I.6).

The perturbing force causes the planet's eccentricity vector to evolve in time as

(I.23) |

where Equations (I.7), (I.8), (I.14), (I.15), and (I.18) have been employed. It follows that

where use has been made of Equation (I.11). Furthermore,

where Equation (I.8) has been employed.

Differentiation of Equation (I.9) with respect to time yields

(I.28) |

Making use of Equations (I.4), (I.7), and (I.10), we can rearrange the previous expression to give

(I.29) |

Differentiation of Equation (I.13) with respect to time yields

(I.30) |

The previous two equations can be combined to give

(I.31) |

where use has been made of Equations (I.9), (I.11), and (I.12). Thus, it follows from Equations (I.11), (I.17), and (I.27) that

(I.32) |