Se trata de una simulación del péndulo esférico, o péndulo en dos dimensiones, cuyo lagrangiano viene 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,
$$
donde $\theta$ y $\phi$ son las coordenadas esféricas tradicionales.
Las ecuaciones de movimiento se integran utilizando el método de Runge-Kutta de cuarto orden con paso de tiempo $h$.
Para ángulos iniciales grandes (mayores a 120º),
es necesario disminuir $h$ para garantizar la convergencia de la solución.
La aceleración de la gravedad se define como $g = 9,8$ m/s².
Leer más: