Necessary Condition for Integral Functional to have Extremum for given function/Dependent on n Variables

Theorem
Let $ \mathbf x $ be an n-dimensional vector.

Let $ u \left ( { \mathbf x } \right ) $ be a real function.

Let $ J $ be a functional such that


 * $ \displaystyle J \left [ { u } \right ] = \idotsint_R F \left ( { \mathbf x, u, u_{ \mathbf x } } \right ) \mathrm d x_1 \dots \mathrm d x_n $

Then a necessary condition for $ J \left [ { u } \right ] $ to have an extremum (strong or weak) for a given function $ u \left ( { \mathbf x } \right )$ is that $ u \left ( { \mathbf x } \right )$ satisfies Euler's equation:


 * $ \displaystyle F_u - \frac{ \partial }{ \partial \mathbf x } F_{ u_{ \mathbf x } } = 0 $

Proof

 * $ \displaystyle F \left [ { \mathbf x, u + h, u_{ \mathbf x } + h_ { \mathbf x } } \right ] = F \left [ { \mathbf x, u, u_{ \mathbf x } } \right ] + \frac{ \partial F \left [ { \mathbf x, u, u_{ \mathbf x } } \right ] }{ \partial u } h + \frac{ \partial F \left [ { \mathbf x, u, u_{ \mathbf x } } \right ] }{ \partial u_{ \mathbf x} } h_{ \mathbf x } + \mathcal O \left ( { h^2, hh_{ \mathbf x }, h_{ \mathbf x}^2  } \right ) $


 * $ \displaystyle \Delta J \left [ { u, h } \right ] = \idotsint_R \left ( { \frac{ \partial F \left [ { \mathbf x, u, u_{ \mathbf x } } \right ] }{ \partial u } h + \frac{ \partial F \left [ { \mathbf x, u, u_{ \mathbf x } } \right ] }{ \partial u_{ \mathbf x } } h_{ \mathbf x } + \mathcal O \left ( { h^2, hh_{ \mathbf x }, h_{ \mathbf x}^2 } \right ) } \right ) \mathrm d x_1 \dots \mathrm d x_n $


 * $ \displaystyle \idotsint_R \frac{ \partial }{ \partial \mathbf x } \left [ { \frac{ \partial F \left [ { \mathbf x, u, u_{ \mathbf x } } \right ] }{ \partial u_{ \mathbf x } } h \left ( { \mathbf x } \right )  } \right ] \mathrm d x_1 \dots \mathrm d x_n = \idotsint_\Gamma h \left ( { \mathbf x } \right ) F_{ u_{ \mathbf x } } \left [ { \mathbf x, u, u_{ \mathbf x } } \right ]  \boldsymbol \nu \mathrm d \sigma$


 * $ \displaystyle \delta J = \idotsint_R \left ( { \frac{ \partial F \left [ { \mathbf x, u, u_{ \mathbf x } } \right ] }{ \partial u } - \frac{ \partial F_{ u_{ \mathbf x } } \left [ { \mathbf x, u, u_{ \mathbf x } } \right ]  }{ \partial \mathbf x }  } \right ) h \left ( { \mathbf x } \right ) \mathrm d x_1 \dots \mathrm d x_n$