多様体を$M$,曲線を$c: (\alpha, \beta) \to M,\, t \mapsto c(t)$,ラグランジアンを接束$TM$上の関数$L: TM \to \mathbb{R},\, (x, v) \mapsto L(x,v)$,作用を$S: C^{\infty}((\alpha, \beta), M) \to \mathbb{R},\, c \mapsto S[c]$とする。$C^{\infty}((\alpha, \beta), M)$は区間$(\alpha, \beta)$から多様体$M$への$C^{\infty}$級写像全体の集合,つまり曲線全体の集合である。閉区間$[a,b]$を$[a,b]\subset(\alpha, \beta)$となるようにとる。作用の値$S[c]$は曲線$c$の速度ベクトルを$c'$として,
\begin{equation}
S[c] = \int_{a}^{b}L(c(t), c'(t))dt
\end{equation}
で与えられる。今は時間に陽に依存しないラグランジアンを考えているが依存していても以降の議論は変わらない。
曲線の変形を$\Phi: (-\epsilon, \epsilon)\times M\to M,\,(h,q) \mapsto\Phi_{h}(q)$とする。これは$h\in(-\epsilon,\epsilon)$をひとつ決めると$M$上の点$q$を動かした点$\Phi_{h}(q)$を与える。$\Phi_{h}(c)$は変形された曲線を表す。この定義にはユークリッド空間の座標を使っていないので多様体上の変分法を表現できる。
変形の条件として以下を課す。
変形された曲線に対する作用の値は次のようになる。
\begin{equation}
S[\Phi_{h}(c)] = \int_{a}^{b}L(\Phi_{h}(c(t)), \frac{\partial}{\partial t}\Phi_{h}(c(t)))dt
\end{equation}
曲線$c$を固定すると$S[\Phi_{(-)}(c)]$は$h$の関数である。状況としてはひとつの決まった曲線をパラメータに応じてぐにゃぐにゃ曲げている。曲線$c$が作用の停留を与えているならば,$h=0$における微分係数が$0$である($h$の$1$変数関数だから例えば2次関数を思い浮かべてもらえればよい)。よって停留作用の原理は以下で与えられる。
\begin{equation} \left.\frac{d}{dh}S[\Phi_{h}(c)] \right|_{h=0} =0 \end{equation}
多様体$M$のチャートを$(U,\phi)$,局所座標を$x^{\>i}$とする。接束$TM$に自然に誘導されるチャートの局所座標を$(x^{\>i},v^{\>i})$とする。速度ベクトルの第$i$成分$\dot{c}^{\>i}$には$\displaystyle \dot{c}^{\>i} = \frac{d}{dt}(x^{\>i} \circ c)$の関係がある。合成関数の微分を行うと,
\begin{align}
&\left.\frac{d}{dh}S[\Phi_{h}(c)] \right|_{h=0} \notag \\
&=
\int_{a}^{b}
\sum_{i}
\left.\frac{\partial L}{\partial x^{\>i}} \right|_{\substack{x = c(t) \\v = c'(t)}}
\left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0} \notag \\
&\quad+ \left.\frac{\partial L}{\partial v^{\>i}} \right|_{\substack{x = c(t) \\v = c'(t)}}
\left.\frac{\partial}{\partial h}\frac{\partial}{\partial t}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0} dt
\end{align}
となる。$x^{\>i}, v^{\>i}$の微分係数をとっている点$x =c(t), v=c'(t)$は,$\Phi_{0}(c(t)) = c(t)$から来ている。計算の途中ではこの表示を省く。$\sum$の中の第2項は変形$\Phi$が$C^{\infty}$級であることから$h$の微分と$t$の微分を交換することができる[1]。
$t$の微分を前に持ってきて部分積分を行う。もうこうなれば普通のEL方程式の導出と同じである。
\begin{align}
&\int_{a}^{b}\frac{\partial L}{\partial v^{\>i}} \frac{\partial}{\partial t}
\left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0}
dt \notag \\
&= \left[ \frac{\partial L}{\partial v^{\>i}}
\left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0}
\right]_{a}^{b} \notag \\
&\quad- \int_{a}^{b}\left(\frac{d}{dt}\frac{\partial L}{\partial v^{\>i}}\right)
\left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0}dt \notag \\
&=\int_{a}^{b} - \left(\frac{d}{dt}\frac{\partial L}{\partial v^{\>i}}\right)
\left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0}dt
\end{align}
境界項は端点においては任意の$h$で変形がない(定値写像である)ので微分係数がゼロになり消える。ここで,$V^{\>i} =\displaystyle \left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0} $とおく。$V^{\>i}$は変分ベクトル場というベクトル場の第$i$成分を定める(これは点$c(t)$における接ベクトルの成分である)。部分積分で得た第2項を元の積分の中に戻すと,
\begin{equation}
\int_{a}^{b}\sum_{i}
\left(
\frac{\partial L}{\partial x^{\>i}} - \frac{d}{dt}\frac{\partial L}{\partial v^{\>i}}
\right)V^{\>i}dt = 0
\end{equation}
変分法の基本補題[2]を適用することで,多様体上のEL方程式を得る。
\begin{align} &\left.\frac{d}{dh}S[\Phi_{h}(c)] \right|_{h=0} =0 \notag \\ &\iff \left. \frac{d}{dt}\frac{\partial L}{\partial v^{\>i}} \right|_{\substack{x = c(t) \\v = c'(t)}} = \left. \frac{\partial L}{\partial x^{\>i}} \right|_{\substack{x = c(t) \\v = c'(t)}} \end{align}