δq(t1)=δq(t2)=0
L(q+δq,˙q+δ˙q)≃L(q,˙q)+∂L(q,˙q)∂qδq+∂L(q,˙q)∂˙qδ˙q
δS[q]≃∫t2t1dt(∂L(q,˙q)∂qδq+∂L(q,˙q)∂˙qδ˙q)=∫t2t1dt(∂L(q,˙q)∂qδq+∂L(q,˙q)∂˙qddtδq)=[∂L(q,˙q)∂˙qδq]t2t1+∫t2t1dt(∂L(q,˙q)∂q−ddt∂L(q,˙q)∂˙q)δq=∫t2t1dt(∂L(q,˙q)∂q−ddt∂L(q,˙q)∂˙q)δq 任意の微小な仮想変位 δq に対してこれが 0 であるとして、 オイラー-ラグランジュ方程式 ∂L(q,˙q)∂q−ddt∂L(q,˙q)∂˙q=0 を得る。
L(q,˙q)=12m˙q2−12kq2
∂L∂q=−kqddt∂L∂q=ddtm˙q=m¨q なので、オイラー-ラグランジュ方程式は、 m¨q=−kq となる。
E(t)=m˙q2−(12m˙q2−12kq2)=12m˙q2+12kq2ddtE(t)=m˙q¨q+kq˙q=m˙q⋅(−kmq)+kq˙q=0
ω=√k/m として、 q(t)=Asinωt+Bcosωt ただし、 A,B は積分定数である。
x, y, z 軸の正の方向の単位ベクトルをそれぞれ \boldsymbol{i}, \boldsymbol{j}, \boldsymbol{k} とする。
I \omega_1 \boldsymbol{k}
\boldsymbol{k} \times \boldsymbol{i} = \boldsymbol{j} なので、y 軸の正の方向である。
I \omega_1 \boldsymbol{k} + lF \Delta t \boldsymbol{j}
\begin{align} \cos \theta &= \frac{ \boldsymbol{k} \cdot \left( I \omega_1 \boldsymbol{k} + lF \Delta t \boldsymbol{j} \right)} { \left| I \omega_1 \boldsymbol{k} + lF \Delta t \boldsymbol{j} \right|} \\ &= \frac{I \omega_1}{\sqrt{I^2 \omega_1^2 + l^2 F^2 \Delta t^2}} \end{align}
\begin{align} \left( X \boldsymbol{i} + Y \boldsymbol{j} + Z \boldsymbol{k} \right) \times \left( - mg \boldsymbol{k} \right) = mg \left( - Y \boldsymbol{i} + X \boldsymbol{j} \right) \end{align}
O から重心に向かう単位ベクトルは、 \begin{align} \frac{X \boldsymbol{i} + Y \boldsymbol{j} + Z \boldsymbol{k}}{R} \end{align} と書けるから、求める運動方程式は、時間微分を \dot{} で表して、 \begin{align} \frac{I \omega_2}{R} \left( \dot{X} \boldsymbol{i} + \dot{Y} \boldsymbol{j} + \dot{Z} \boldsymbol{k} \right) = mg \left( - Y \boldsymbol{i} + X \boldsymbol{j} \right) \end{align} あるいは成分で書けば \begin{align} \dot{X} &= - \frac{mgR}{I \omega_2} Y \\ \dot{Y} &= \frac{mgR}{I \omega_2} X \\ \dot{Z} &= 0 \end{align} である。
まず、 Z = R \cos \alpha はすぐにわかる。 次に、コマの軸をxy平面に射影したときのx軸からの角度を \varphi とすると、 \begin{align} X &= R \sin \alpha \cos \varphi \\ Y &= R \sin \alpha \sin \varphi \end{align} であるから、これを (6) の運動方程式に代入して、次を得る: \begin{align} \dot{\varphi} = \frac{mg}{I \omega_2} \end{align} t=0 のとき \varphi = 0 であるから、次がわかる: \begin{align} \varphi = \frac{mgt}{I \omega_2} \end{align} まとめると、求める解は、 \begin{align} X &= R \sin \alpha \cos \frac{mgt}{I \omega_2} \\ Y &= R \sin \alpha \sin \frac{mgt}{I \omega_2} \\ Z &= R \cos \alpha \end{align} である。
(7) より \begin{align} \Omega &= \dot{\varphi} \\ &= \frac{mg}{I \omega_2} \end{align}