I'm trying to solve and understand a provided solution to a problem in Rovelli - Quantum Loop Gravity on page 56, eq. number (2.86).
In the problem among other things I'm asked to compute the equations of motion for $a$ and $b$ using following covariant Lagrangian: $$ {\cal{L}} \left( a, b, \dot{a}, \dot{b} \right) = \sqrt{\left(2E - a^2 - b^2\right) \left( \dot{a}^2 + \dot{b}^2 \right)}. \tag{2.86} $$ The hint says that:
- The action is minimized for a geodesic of $$ds^2 = \left(2E - a^2 - b^2 \right) \left( da^2 + db^2 \right).$$
- Gauge fixing can be applied to set $\dot{a}^2 + \dot{b}^2 = 1$.
- Variables can be separated.
Furthermore, I'm supposed to show that the system is formally similar to two harmonic oscillators with total energy $E$.
Hamiltonian associated to this Lagrangian is trivially zero (as it should be for covariant theory). Applying gauge-fixing leads to trivial equations of motion $a = 0$, $b = 0$ so I assume this should not be done as the first step. Otherwise the equations of motion derived from:
$$ p_a = \frac{\partial{\cal{L}}}{\partial \dot{a}} = \dot{a} \sqrt{\frac{2E - a^2 - b^2}{\dot{a}^2 + \dot{b}^2}}, \\ \frac{\partial{\cal{L}}}{\partial{a}} = -a \sqrt{\frac{\dot{a}^2 + \dot{b}^2}{2E - a^2 - b^2}} $$ are inseparable and cumbersome with no solution in sight.
There is a Hamiltonian constraint provided in the text that can be easily verified without clear steps to derivation: $$ C = \frac{1}{2} \left( p_a^2 + a^2 + p_b^2 + b^2 \right) - E = 0.\tag{2.87} $$
Technically, it should play the role of equations of motion, but substituting $p_a$ from above and similarly $p_b$ gives trivial zero yet again.
What's even more confusing is that the answers that Rovelli provides: $$\begin{align} a(\tau) =& A \sin(\tau + \phi_a),\tag{2.90} \cr b(\tau) =& B \sin(\tau + \phi_b) \tag{2.91} \end{align}$$ does not seem to be a solution either (unless I'm mistaken).
Does anyone know how this could be solved?