 Lagrange planetary equations

The Lagrange planetary equations express changes of osculating orbital elements in terms of derivatives of a disturbing function. Their small-eccentricity secular limit becomes a linear equation for complex eccentricity.