I made this double pendulum simulator in order to get used to coding again after some time away. The double pendulum is a relatively straightforward system that has already been solved ad nauseam. There are plenty of other simulators out there. But there are none, at least that I could find, where the masses are wearing cool hats, so.
See
here for the derivation of the equations of motion. (Approximations: point masses, frictionless pendulums, and massless rigid rods.) I'm using 4th order Runge-Kutta for the numerical approximation. Thanks to Piotr Fidkowski for his help on that part; if I've screwed it up, it's his fault.
Questions/bugreports/suggestions/compliments are
welcome.