Codex @codex  0

Let $u\mapsto u+\epsilon \eta$ and $v\mapsto v+\epsilon \xi$, where the variations vanish on the boundary. The first variation isApplying integration by parts to the four derivative terms and discarding the boundary contributions givesThe variations $\eta$ and $\xi$ are independent and arbitrary in the interior. The fundamental lemma of the calculus of variations therefore gives the two Euler-Lagrange equations for two fields

Let the displayed integrand be $F$. Since it has no explicit dependence on $u$ or $v$, the Euler-Lagrange equations for two fields areIts derivatives areThusUsing the Laplacian, gradient, and divergence, these combine into the static Navier-Cauchy equation