Euler–Lagrange equation explained

In the calculus of variations and classical mechanics, the Euler–Lagrange equations^[1] are a system of second-order ordinary differential equations whose solutions are stationary points of the given action functional. The equations were discovered in the 1750s by Swiss mathematician Leonhard Euler and Italian mathematician Joseph-Louis Lagrange.

Because a differentiable functional is stationary at its local extrema, the Euler–Lagrange equation is useful for solving optimization problems in which, given some functional, one seeks the function minimizing or maximizing it. This is analogous to Fermat's theorem in calculus, stating that at any point where a differentiable function attains a local extremum its derivative is zero. In Lagrangian mechanics, according to Hamilton's principle of stationary action, the evolution of a physical system is described by the solutions to the Euler equation for the action of the system. In this context Euler equations are usually called Lagrange equations. In classical mechanics,^[2] it is equivalent to Newton's laws of motion; indeed, the Euler-Lagrange equations will produce the same equations as Newton's Laws. This is particularly useful when analyzing systems whose force vectors are particularly complicated. It has the advantage that it takes the same form in any system of generalized coordinates, and it is better suited to generalizations. In classical field theory there is an analogous equation to calculate the dynamics of a field.

History

The Euler–Lagrange equation was developed in the 1750s by Euler and Lagrange in connection with their studies of the tautochrone problem. This is the problem of determining a curve on which a weighted particle will fall to a fixed point in a fixed amount of time, independent of the starting point.

Lagrange solved this problem in 1755 and sent the solution to Euler. Both further developed Lagrange's method and applied it to mechanics, which led to the formulation of Lagrangian mechanics. Their correspondence ultimately led to the calculus of variations, a term coined by Euler himself in 1766.^[3]

Statement

Let

(X,L) be a real dynamical system with n degrees of freedom. Here X is the configuration space and L = L(t,{\boldsymbol q}(t), {\boldsymbol v}(t)) the Lagrangian, i.e. a smooth real-valued function such that {\boldsymbol q}(t) \in X, and {\boldsymbol v}(t) is an n -dimensional "vector of speed". (For those familiar with differential geometry, X is a smooth manifold, and L : { R} _{t} x TX \to { R}, where TX is the tangent bundle of X).

Let

{\cal P}(a,b,\boldsymbol x _{a,\boldsymbol} x _{b)} be the set of smooth paths \boldsymbol q: [a,b] \to X for which \boldsymbol q(a) = \boldsymbol x _{a} and \boldsymbol q(b) = \boldsymbol x _{b.} S : {\cal P}(a,b,\boldsymbol x _{a,\boldsymbol} x _{b)} \to R is defined via$$S[\backslash boldsymbol\; q]\; =\; \backslash int\_a^b\; L(t,\backslash boldsymbol\; q(t),\backslash dot(t))\backslash ,\; dt.$$

A path

\boldsymbol q \in {\cal P}(a,b,\boldsymbol x _{a,\boldsymbol} x _{b)} is a stationary point of S if and only if

Here,

• \boldsymbol q (t) is the time derivative of \boldsymbol q(t). When we say stationary point, we mean a stationary point of S with respect to any small perturbation in \boldsymbol q . See proofs below for more rigorous detail.

Example

A standard example is finding the real-valued function y(x) on the interval [''a'', ''b''], such that y(a) = c and y(b) = d, for which the path length along the curve traced by y is as short as possible.

s = \int ^{b} _{a} \sqrt{ d x ^{2+dy} ^{2}} = \int ^{b} _{a} \sqrt{ 1+ y' ^{2}dx,} the integrand function being $L(x,y,\; y\text{'})\; =\; \backslash sqrt$.

The partial derivatives of L are:

(\partial L(x, y, y'))/(\partial y') = (y')/(\sqrt{ 1+ y' ^{2}) } \quad \text \quad\frac = 0.By substituting these into the Euler–Lagrange equation, we obtain\begin{align} (d)/(d x) (y'(x))/(\sqrt{ 1+ (y'(x)) ^{2}) } &= 0 \\ \frac &= C = \text \\\Rightarrow y'(x)&= \frac =: A \\\Rightarrow y(x) &= Ax + B\endthat is, the function must have a constant first derivative, and thus its graph is a straight line.

Generalizations

Single function of single variable with higher derivatives

The stationary values of the functional

I[f] = \int ^{x _{1}} _{x _{0}} l{L}(x, f, f', f'', ... , f ^{(k)} )~ d x ~;~~ f' : = \cfrac{ d f}{ d x}, ~f'' : = \cfrac{ d ^{2f}{dx} ^{2},} ~ f ^{(k)} : = \cfrac{ d ^{kf}{dx} ^{k}} can be obtained from the Euler–Lagrange equation^[4] \cfrac{\partial l{L}}{\partial f} - \cfrac{ d }\left(\cfrac\right) + \cfrac\left(\cfrac\right) - \dots + (-1)^k \cfrac\left(\cfrac\right) = 0 under fixed boundary conditions for the function itself as well as for the first k -1 derivatives (i.e. for all f ^{(i)} , i \in \{ 0 , ..., k -1 \} ). The endpoint values of the highest derivative f ^{(k)} remain flexible.

Several functions of single variable with single derivative

If the problem involves finding several functions (

f _{1,} f _{2,} ... , f _{m} ) of a single independent variable (x ) that define an extremum of the functionalI[f _{1,f} _{2,} ... , f _{m]} = \int ^{x _{1}} _{x _{0}} l{L}(x, f _{1,} f _{2,} ... , f _{m,} f _{1',} f _{2',} ... , f _{m')~dx} ~;~~ f _{i'} : = \cfrac{ d f _{i}{dx}} then the corresponding Euler–Lagrange equations are^[5] \begin{align} (\partial l{L)/() } - \frac\left(\frac\right) = 0 ; \quad i = 1, 2, ..., m \end

Single function of several variables with single derivative

A multi-dimensional generalization comes from considering a function on n variables. If

\Omega is some surface, thenI[f] = \int _{\Omega} l{L}(x _{1,} ... , x _{n,} f, f _{1} , ... , f _{n} ) d x ~;~~ f _{j} : = \cfrac{\partial f}{\partial x _{j}}

is extremized only if f satisfies the partial differential equation

(\partial l{L)/() } - \sum_^ \frac\left(\frac\right) = 0.

When n = 2 and functional

l I is the energy functional, this leads to the soap-film minimal surface problem.

Several functions of several variables with single derivative

If there are several unknown functions to be determined and several variables such that

I[f _{1,f} _{2,...,f} _{m]} = \int _{\Omega} l{L}(x _{1,} ... , x _{n,} f _{1,} ... , f _{m,} f _{1,1} , ... , f _{1,n} , ... , f _{m,1} , ... , f _{m,n} ) d x ~;~~ f _{i,j} : = \cfrac{\partial f _{i}{\partial} x _{j}} the system of Euler–Lagrange equations is^[4] \begin{align} (\partial l{L)/() } - \sum_^ \frac\left(\frac\right) &= 0_1 \\ \frac - \sum_^ \frac\left(\frac\right) &= 0_2 \\ \vdots \qquad \vdots \qquad &\quad \vdots \\ \frac - \sum_^ \frac\left(\frac\right) &= 0_m. \end

Single function of two variables with higher derivatives

If there is a single unknown function f to be determined that is dependent on two variables x₁ and x₂ and if the functional depends on higher derivatives of f up to n-th order such that

\begin{align} I[f] & = \int _{\Omega} l{L}(x _{1,} x _{2,} f, f _{1} , f _{2} , f _{11} , f _{12} , f _{22} , ... , f _{22...} ) d x \\ & f _{i} : = \cfrac{\partial f}{\partial x _{i}} , f _{ij} : = \cfrac{\partial ^{2} f}{\partial x _{i\partial} x _{j}} , ... \end{align} then the Euler–Lagrange equation is^[4] \begin{align} (\partial l{L)/() } & - \frac\left(\frac\right) - \frac\left(\frac\right) + \frac\left(\frac\right) + \frac\left(\frac\right) + \frac\left(\frac\right) \\ & - \dots + (-1)^n \frac\left(\frac\right) = 0 \end which can be represented shortly as:(\partial l{L)/() } +\sum_^n \sum_ (-1)^j \frac \left(\frac\right)=0 wherein \mu _{1} ... \mu _{j} are indices that span the number of variables, that is, here they go from 1 to 2. Here summation over the \mu _{1} ... \mu _{j} indices is only over \mu _{1} \leq \mu _{2} \leq \ldots \leq \mu _{j} in order to avoid counting the same partial derivative multiple times, for example f _{12} = f _{21} appears only once in the previous equation.

Several functions of several variables with higher derivatives

If there are p unknown functions f_i to be determined that are dependent on m variables x₁ ... x_m and if the functional depends on higher derivatives of the f_i up to n-th order such that

\begin{align} I[f _{1,\ldots,f} _{p]} & = \int _{\Omega} l{L}(x _{1,} \ldots, x _{m;} f _{1,\ldots,f} _{p;} f _{1,1} ,\ldots, f _{p,m} ; f _{1,11} ,\ldots, f _{p,mm} ;\ldots; f _{p,1\ldots} , \ldots, f _{p,m\ldots} ) d x \\ & f _{i,\mu} : = \cfrac{\partial f _{i}{\partial} x _{\mu}} , f _{i,\mu _{1\mu} _{2}} : = \cfrac{\partial ^{2} f _{i}{\partial} x _{\mu _{1}} \partial x _{\mu _{2}} } \;, \;\; \dots \end

where

\mu _{1} ... \mu _{j} are indices that span the number of variables, that is they go from 1 to m. Then the Euler–Lagrange equation is(\partial l{L)/() } +\sum_^n \sum_ (-1)^j \frac \left(\frac\right)=0

where the summation over the

\mu _{1} ... \mu _{j} is avoiding counting the same derivative f _{i,\mu _{1\mu} _{2}} = f _{i,\mu _{2\mu} _{1}} several times, just as in the previous subsection. This can be expressed more compactly as\sum ^{n} _{j =0} \sum _{\mu _{1} \leq \ldots \leq \mu _{j}} ( -1 ) ^{j} \partial ^{j} _{\mu _{1} \ldots \mu _{j}} \left( (\partial l{L)/(}{\partial) f _{i,\mu _{1...\mu} _{j}} }\right)=0

Generalization to manifolds

Let

M be a smooth manifold, and let C ^{infty([a,b])} denote the space of smooth functions f\colon [a,b]\to M . Then, for functionals S\colon C ^{infty} ([a,b])\to R of the formS[f] = \int ^{b} _{a} (L\circ • f )(t) d t where L\colon TM\to R is the Lagrangian, the statement d S _{f=0} is equivalent to the statement that, for all t\in [a,b] , each coordinate frame trivialization (x ^{i,X} ^{i)} of a neighborhood of • f (t) yields the following \dim M equations:\forall i: (d)/(d t) (\partial L)/(\partial X ^{i}) | _{• f (t)} = (\partial L)/(\partial x ^{i}) | _{• f (t)} . Euler-Lagrange equations can also be written in a coordinate-free form as ^[6] l{L} _{\Delta} \theta _{L=dL} where \theta _{L} is the canonical momenta 1-form corresponding to the Lagrangian L . The vector field generating time translations is denoted by \Delta and the Lie derivative is denoted by l{L} . One can use local charts (q ^{• \alpha, q} ^{\alpha)} in which \theta _{L = (\partial L)/(\partial • q ^{\alpha})} dq ^{\alpha} and \Delta: = (d)/(= dt) • q ^{\alpha (\partial)/(\partial q ^{\alpha})} + \ddot{q} ^{\alpha (\partial)/(\partial • q ^{\alpha})} and use coordinate expressions for the Lie derivative to see equivalence with coordinate expressions of the Euler Lagrange equation. The coordinate free form is particularly suitable for geometrical interpretation of the Euler Lagrange equations.

See also

• Lagrangian mechanics
• Hamiltonian mechanics
• Analytical mechanics
• Beltrami identity
• Functional derivative

References

• • • • Book: Gelfand, Izrail Moiseevich . Israel Gelfand . Calculus of Variations . Dover . 1963 . 0-486-41448-5.
• Roubicek, T.: Calculus of variations. Chap.17 in: Mathematical Tools for Physicists. (Ed. M. Grinfeld) J. Wiley, Weinheim, 2014,, pp. 551–588.

Notes and References

1. Book: Fox, Charles. An introduction to the calculus of variations. Courier Dover Publications. 1987. 978-0-486-65499-7.
2. Book: Herbert Goldstein. Charles P. Poole. H.. Goldstein. C.P.. Poole. J.. Safko. Classical Mechanics. Addison Wesley. 2014. 3rd.
3. http://numericalmethods.eng.usf.edu/anecdotes/lagrange.pdf A short biography of Lagrange
4. Book: Courant . R . Richard Courant . Hilbert . D . David Hilbert . Methods of Mathematical Physics . I . First English . Interscience Publishers, Inc . 1953 . New York . 978-0471504474.
5. Book: Weinstock, R. . 1952 . Calculus of Variations with Applications to Physics and Engineering . registration . McGraw-Hill . New York .
6. Web site: José . Saletan . 1998 . Classical Dynamics: A contemporary approach . 2023-09-12 . Cambridge University Press . en . 9780521636360.