A partir de la Lagrangiana he obtenido las siguientes ecuaciones de movimiento para el doble péndulo en 2D. (Las masas son diferentes pero las longitudes de las dos péndulas son iguales.) Sea $m_2$ sea la masa más baja.
$$(m_1+m_2)\ddot{\theta_1}+2m_2\ddot\theta_2\cos(\theta_2-\theta_1)=\\ -2m_2\dot\theta_1\dot\theta_2\sin(\theta_1-\theta_2)-(m_1+m_2)g/l\sin(\theta_1)$$
y
$$m_2\ddot{\theta_1}+2m_2\ddot\theta_2\cos(\theta_2-\theta_1)=\\ 2m_2\dot\theta_1\dot\theta_2\sin(\theta_1-\theta_2)-m_2g/l\sin(\theta_1)$$
En la aproximación de ángulos pequeños se convierten, respectivamente, en
$$(m_1+m_2)\ddot{\theta_1}+2m_2\ddot\theta_2= -2m_2\dot\theta_1\dot\theta_2(\theta_1-\theta_2)-\theta_1(m_1+m_2)g/l$$
y
$$m_2\ddot{\theta_1}+2m_2\ddot\theta_2= 2m_2\dot\theta_1\dot\theta_2(\theta_1-\theta_2)-\theta_1m_2g/l$$ .
La mayoría de las fuentes no tienen los términos de la orden $\dot\theta$ . Esto se debe a que aplican la aproximación de ángulo pequeño a la lagrangiana antes de tomar las derivadas, ignorando así los términos de orden $\theta.$ ¿Qué justificación tenemos para deshacernos de estos términos?