Simulation, Physics, Numerical Methods, Classical Mechanics, Computation.
#Convert back to cartesian coordinates x = l*np.sin(theta) y = -l*np.cos(theta) #Use same function from simple pendulum makeGif(x,y,"pendulum.gif") Double Pendulum def theta1_acceleration(m1,m2,l1,l2,theta1,theta2,theta1_velocity,theta2_velocity,g): mass1 = -g*(2*m1 + m2)*np.sin(theta1) mass2 =