Heavy symmetric top reduction (source code)

= Heavy symmetric top reduction

Fixing the cyclic momenta $p_\psi=\lambda$ and $p_\phi=\mu$ leaves the canonical pair $(\theta,p_\theta)$ and a one-degree-of-freedom Hamiltonian with effective potential
$$
V_{\lambda,\mu}(\theta)
=\frac{(\mu-\lambda\cos\theta)^2}{2I\sin^2\theta}
+\frac{\lambda^2}{2J}+Mga\cos\theta.
$$