$$\theta_1 ''=\frac{-g(2m_1+m_2)sin\theta_1-m_2gsin(\theta_1-2\theta_2)-2sin(\theta_1-\theta_2)m_2(\theta_2'^2l_2+\theta_1'^2l_1cos(\theta_1-\theta_2)}{l_1(2m_1+m_2-m_2cos(2cos(2\theta_1-2\theta_2)))}$$
$$\theta_2 ''=\frac{2sin(\theta_1-\theta_2)(\theta_1'l_1(m_1+m_2)+g(m_1+m_2)cos\theta_1+\theta_2'^2l_2m_2cos(\theta_1-\theta_2))}{l_2(2m_1+m_2-m_2cos(2cos(2\theta_1-2\theta_2)))}$$
These are the equations, and How can I possibly solve this simultaneoouslt in RK4 method?
$m_1,m_2=$masses of pendulum 1 and 2, $\theta_1,\theta_2=$ angles formed by the pendulums, $\theta_1'=\omega_1,\theta_2'=\omega_2$