I am attempting to simulate a double spherical pendulum, i.e. a combination of the spherical pendulum and the double pendulum.

I understand that the equations of motion can be derived via the Lagrangian and the Euler-Lagrange equations. However this method rapidly becomes very messy.

Is there an alternate method that could be used to simplify the calculations required?