Simulation, Physics, Numerical Methods, Classical Mechanics, Computation.
theta2)) normalization = l1*(2*m1 + m2 - m2*np.cos(2*theta1 - 2*theta2)) theta2_ddot = system/normalization return theta2_ddot #Update theta1 def theta1_update(m1,m2,l1,l2,theta1,theta2,theta1_velocity,theta2_velocity,g,time_step): #Numerical Integration theta1_velocity +=