Esta é uma simulação do pêndulo esférico, ou pêndulo em duas dimensões, cujo lagrangeano é dado por
$$
\mathcal{L} = \dfrac{1}{2}mL^2\dot{\theta}^2 + \dfrac{1}{2}mL^2\dot{\phi}^2\sin^2\theta +
$$
$$
+ mgL\cos\theta,
$$
onde $\theta$ e $\phi$ são as coordenadas esféricas tradicionais.
As equações do movimento são integradas usando o método de Runge-Kutta de quarta ordem com passo temporal $h$.
Para ângulos iniciais grandes (maiores do que 120º),
é necessário diminuir $h$ para garantir a convergência da solução.
A aceleração da gravidade está definida como $g = 9.8$ m/s².
Leia mais: