  Given a non-variational system of differential equations, the simplest way of
turning it into a variational one is by adding a correction term. In the paper,
we propose a way of obtaining such a correction term, based on the so-called
Vainberg-Tonti Lagrangian, and present several applications in general
relativity and classical mechanics.
