Programming a Double Pendulum

How do I model and solve for theta1 and theta2 for double pendulum problem? I am a complete beginner using Matlab and any programming language. I have solved for the Lagrangian by hand. I'm just not sure what needs to be in the 'for loop'. I also need to plot time theta1 and theta2. Please help me!

Answer 1

[Jan]

There are mnany different solutions. Perhaps you simply want to integrate the Lagrange function with ODE45.

Note that the simulation is more efficient using depending coordinates and a constraint matrix (see odeset->mass matrix and DAE options), that solving the equations to independent coordinates, because the trigonometrix stuff requires a lot of time.