Perturbed Kepler-orbit conservation laws

If the relative acceleration is $\ddot{\mathbf{r}}=-GM\mathbf{r}/r^3+\mathbf{f}$, with constant masses, then the specific orbital energy $\epsilon$, specific angular momentum $\mathbf{h}$ and eccentricity vector obeyThe inverse-square-force terms cancel in each derivative. These identities express perturbing work and torque, and track the instantaneous change of the ellipse's shape and orientation.