Simulation, Physics, Numerical Methods, Classical Mechanics, Computation.
theta2_acceleration(m1,m2,l1,l2,theta1,theta2,theta1_velocity,theta2_velocity,g): system = 2*np.sin(theta1 - theta2)*(theta1_velocity**2*l1*(m1 + m2) + g*(m1 + m2)*np.cos(theta1) + theta2_velocity**2*l2*m2*np.cos(theta1 - theta2)) normalization = l1*(2*m1 + m2 - m2*np.cos(2*theta1 - 2*theta2))