Il s'agit d'une simulation du pendule sphérique, ou pendule à deux dimensions, dont le lagrangien est donné par
$$
\mathcal{L} = \dfrac{1}{2}mL^2\dot{\theta}^2 + \dfrac{1}{2}mL^2\dot{\phi}^2\sin^2\ thêta +
$$
$$
+ mgL\cos\thêta,
$$
où $\theta$ et $\phi$ sont les coordonnées sphériques traditionnelles.
Les équations du mouvement sont intégrées à l'aide de la méthode Runge-Kutta du quatrième ordre avec pas de temps $h$.
Pour les grands angles initiaux (supérieurs à 120º),
il faut diminuer $h$ pour garantir la convergence de la solution.
L'accélération de la gravité est définie comme $g = 9,8$ m/s².
En savoir plus: