caida libre

import matplotlib.pyplot as plt

g  =  9.8
y0  =  100
v0  =  0

ti  =  0
tf  =  5
Nt  =  98
dt  =  (tf-ti)/Nt

for i in range(98):
  t  = ti+i*dt
  y  = y0+v0*t-0.5*g*t*t
  plt.ylim(-10, 110)
  plt.plot(0, y, 'o', color='black')
  if i<  10: plt.savefig("caida0{}.png".format(i))
  if i>= 10: plt.savefig( "caida{}.png".format(i))
  plt.close()

Este programa se crea con la instruccion

$convert c*.png caida.gif