ProofComplete

Proof of the Euler-Lagrange Equation

The Euler-Lagrange equation is the foundational result of the calculus of variations, providing necessary conditions for a function to be an extremal of a functional. Its derivation combines integration by parts with the fundamental lemma of the calculus of variations.


Statement

Theorem7.3Euler-Lagrange Equation

Let L(x,y,p)∈C2L(x,y,p) \in C^2 and suppose yβˆ—βˆˆC2([a,b])y^* \in C^2([a,b]) is a stationary point of the functional J[y]=∫abL(x,y(x),yβ€²(x)) dxJ[y] = \int_a^b L(x,y(x),y'(x))\,dx subject to fixed boundary conditions y(a)=Ay(a) = A, y(b)=By(b) = B. Then yβˆ—y^* satisfies the Euler-Lagrange equation: βˆ‚Lβˆ‚y(x,yβˆ—,(yβˆ—)β€²)βˆ’ddxβˆ‚Lβˆ‚p(x,yβˆ—,(yβˆ—)β€²)=0\frac{\partial L}{\partial y}(x,y^*,(y^*)') - \frac{d}{dx}\frac{\partial L}{\partial p}(x,y^*,(y^*)') = 0 for all x∈[a,b]x \in [a,b].


Proof

Proof

Step 1: Setup the variation. Let η∈C2([a,b])\eta \in C^2([a,b]) satisfy Ξ·(a)=Ξ·(b)=0\eta(a) = \eta(b) = 0 (admissible variation). For Ρ∈R\varepsilon \in \mathbb{R}, define Ξ¦(Ξ΅)=J[yβˆ—+Ρη]\Phi(\varepsilon) = J[y^* + \varepsilon\eta]. Since yβˆ—y^* is a stationary point, Ξ¦β€²(0)=0\Phi'(0) = 0.

Step 2: Compute the derivative. Differentiating under the integral: Ξ¦β€²(Ξ΅)=∫ab[βˆ‚Lβˆ‚y(x,yβˆ—+Ρη,(yβˆ—)β€²+Ρη′)β‹…Ξ·+βˆ‚Lβˆ‚p(x,yβˆ—+Ρη,(yβˆ—)β€²+Ρη′)β‹…Ξ·β€²]dx\Phi'(\varepsilon) = \int_a^b \left[\frac{\partial L}{\partial y}(x, y^*+\varepsilon\eta, (y^*)'+\varepsilon\eta')\cdot\eta + \frac{\partial L}{\partial p}(x, y^*+\varepsilon\eta, (y^*)'+\varepsilon\eta')\cdot\eta'\right]dx

Evaluating at Ξ΅=0\varepsilon = 0: 0=Ξ¦β€²(0)=∫ab[Ly(x,yβˆ—,(yβˆ—)β€²)Ξ·+Lp(x,yβˆ—,(yβˆ—)β€²)Ξ·β€²]dx0 = \Phi'(0) = \int_a^b\left[L_y(x,y^*,(y^*)')\eta + L_p(x,y^*,(y^*)')\eta'\right]dx

where we write Ly=βˆ‚L/βˆ‚yL_y = \partial L/\partial y and Lp=βˆ‚L/βˆ‚pL_p = \partial L/\partial p for brevity.

Step 3: Integration by parts. Apply integration by parts to the second term: ∫abLp η′ dx=[Lp η]abβˆ’βˆ«abddx(Lp) η dx\int_a^b L_p\,\eta'\,dx = \left[L_p\,\eta\right]_a^b - \int_a^b\frac{d}{dx}(L_p)\,\eta\,dx

The boundary term [Lp η]ab=Lp(b)Ξ·(b)βˆ’Lp(a)Ξ·(a)=0[L_p\,\eta]_a^b = L_p(b)\eta(b) - L_p(a)\eta(a) = 0 since Ξ·(a)=Ξ·(b)=0\eta(a) = \eta(b) = 0. Substituting: 0=∫ab[Lyβˆ’ddxLp]η dx0 = \int_a^b\left[L_y - \frac{d}{dx}L_p\right]\eta\,dx

Step 4: Fundamental lemma. We invoke the fundamental lemma of the calculus of variations: if f∈C([a,b])f \in C([a,b]) and ∫abf(x)Ξ·(x) dx=0\int_a^b f(x)\eta(x)\,dx = 0 for every η∈C2([a,b])\eta \in C^2([a,b]) with Ξ·(a)=Ξ·(b)=0\eta(a) = \eta(b) = 0, then f(x)=0f(x) = 0 for all x∈[a,b]x \in [a,b].

Proof of the lemma: Suppose f(x0)>0f(x_0) > 0 for some x0∈(a,b)x_0 \in (a,b). By continuity, f>0f > 0 on (x0βˆ’Ξ΄,x0+Ξ΄)(x_0-\delta, x_0+\delta) for some Ξ΄>0\delta > 0. Choose Ξ·(x)=[(xβˆ’x0+Ξ΄)(x0+Ξ΄βˆ’x)]3\eta(x) = [(x-x_0+\delta)(x_0+\delta-x)]^3 on (x0βˆ’Ξ΄,x0+Ξ΄)(x_0-\delta,x_0+\delta) and Ξ·=0\eta = 0 outside. Then Ξ·β‰₯0\eta \geq 0, η∈C2\eta \in C^2, and ∫fη dx>0\int f\eta\,dx > 0, contradicting the hypothesis. Similarly for f(x0)<0f(x_0) < 0. Therefore f≑0f \equiv 0.

Step 5: Conclusion. Since Lyβˆ’ddxLpL_y - \frac{d}{dx}L_p is continuous (by the C2C^2 hypotheses on LL and yβˆ—y^*), the fundamental lemma gives: βˆ‚Lβˆ‚y(x,yβˆ—,(yβˆ—)β€²)βˆ’ddxβˆ‚Lβˆ‚p(x,yβˆ—,(yβˆ—)β€²)=0forΒ allΒ x∈[a,b]\frac{\partial L}{\partial y}(x,y^*,(y^*)') - \frac{d}{dx}\frac{\partial L}{\partial p}(x,y^*,(y^*)') = 0 \quad \text{for all } x \in [a,b]

which is the Euler-Lagrange equation. β–‘\square

β– 

ExampleExpanded Form of the Euler-Lagrange Equation

Expanding ddxLp\frac{d}{dx}L_p by the chain rule: Lpx+Lpy yβ€²+Lpp yβ€²β€²=LyL_{px} + L_{py}\,y' + L_{pp}\,y'' = L_y. When Lppβ‰ 0L_{pp} \neq 0 (the Legendre condition, also called regularity), this is a second-order ODE: yβ€²β€²=1Lpp[Lyβˆ’Lpxβˆ’Lpy yβ€²]y'' = \frac{1}{L_{pp}}[L_y - L_{px} - L_{py}\,y']. The condition Lpp>0L_{pp} > 0 (strict convexity in pp) ensures that the extremal is a local minimum (Jacobi's strengthened condition).

RemarkWeak and Strong Extremals

The proof above yields a necessary condition for weak extremals (variations small in the C1C^1 norm). For strong extremals (small in C0C^0 norm), additional conditions are needed: the Weierstrass excess function E(x,y,p,q)=L(x,y,q)βˆ’L(x,y,p)βˆ’(qβˆ’p)Lp(x,y,p)β‰₯0\mathcal{E}(x,y,p,q) = L(x,y,q) - L(x,y,p) - (q-p)L_p(x,y,p) \geq 0 for all qq (Weierstrass necessary condition), and the Jacobi condition (no conjugate points). Together, the Euler-Lagrange equation, Legendre condition, Weierstrass condition, and Jacobi condition form the complete set of sufficient conditions for a strong minimum.