Simulation, Physics, Numerical Methods, Classical Mechanics, Computation.
np.radians(15) #degrees to radians g = 9.8 #m/s^2 l = 1.0 #m omega = np.sqrt(g/l) phi = 0 #for small angle time_span = np.linspace(0,20,300) #simulate for 20s split into 300 time intervals theta = [] for t in time_span: theta.append(simple_pendulum(theta_0, omega, t, phi)) x = l*np.sin(theta) y =