Normal modes of a double pendulum

For equal lengths and mass ratio $\mu =m/M$, the squared frequencies are ${\omega }_0^2(1+\mu \pm \sqrt{\mu (1+\mu )})$.