Solution (source code)

= Solution

The differential operators are self-adjoint under integration by parts, so the <Euler-Lagrange field equations> are
$$
\widetilde a\Box V^\nu
+\widetilde b\partial^\nu(\partial\cdot V)
+m^2V^\nu=0.
$$