Post-Newtonian Lagrangian planetary equations

We present a method to study the time variation of the orbital parameters of a Post-Keplerian binary system undergoing a generic external perturbation. The method is the relativistic extension of the planetary Lagrangian equations. The theory only assumes the smallness of the external perturbation while relativistic effects are already included in the unperturbed problem. This is the major advantage of this novel approach over classical Lagrangian methods.