from scipy.integrate import odeint import numpy as np %pylab inline
g= -9.8 def f(y, t): return g * t t0=0. t1=3. y0=200 t=np.arange(t0, t1, 0.01) y= odeint(f, y0, t)
plt.plot(t,y) plt.xlabel('$t$') plt.ylabel('$y(t)$')