 Solution

Let $P^{\mu }=\partial \mathcal{L}/\partial ({\partial }_{\mu }h)$ and $Q^{\mu \nu }=\partial \mathcal{L}/\partial ({\partial }_{\mu }{\partial }_{\nu }h)$. Define the second coefficient with ordered repeated indices and $Q^{\mu \nu }=Q^{\nu \mu }$, so thatThis convention avoids double-counting a mixed derivative such as $h_{,12}=h_{,21}$. The first derivative term needs one integration by parts. The second derivative term needs two:Combining the terms givesFor compactly supported variations, or fixed $\delta h$ and its first derivatives at the boundary, the last integral vanishes. The fundamental lemma of the calculus of variations then yields the Euler-Lagrange field equation with second derivatives:The plus sign of the final term comes from the two integrations by parts. In general this gives fourth-order equations, although the order in a particular coordinate depends on which second derivatives occur in the density.