Simulation, Physics, Numerical Methods, Classical Mechanics, Computation.
l = 1.0 #m theta = [np.radians(90)] #theta_0 theta_velocity = 0 time_step = 20/300 time_span = np.linspace(0,20,300) #simulate for 20s split into 300 time intervals for t in time_span: theta_new, theta_velocity = full_pendulum(g,l,theta[-1], theta_velocity, time_step) theta.append(theta_new)