Simulation, Physics, Numerical Methods, Classical Mechanics, Computation.
-l1*np.cos(theta1_list) - l2*np.cos(theta2_list) #Pendulum 2 y return x1,y1,x2,y2 #Define system parameters g = 9.8 #m/s^2 m1 = 1 #kg m2 = 1 #kg l1 = 1 #m l2 = 1 #m theta1 = np.radians(90) theta2 = np.radians(45) theta1_velocity = 0 #m/s theta2_velocity = 0 #m/s theta1_list = [theta1] theta2_list =