m=1 # [kg] mass of the object
  h=5;xNum=h # initial position
  v=0 # initial velocity
  g=9.776 #[m/s^2] gravitational constant Bogota
  F=-m*g
  dt=0.03 # [s] time advance
  t=0 # [s] initial time

Time = []
PositionNum=[]
PositionAnal=[]
while xNum>0:  
  t=t+dt # time evolution
  Time.append(t)
  a=F/m # acceleration "evolution"
  v=v+a*dt # velocity evolution
  xNum=xNum+v*dt # position evolution
  xAnal=h-g*t**2/2 # position evolution
  PositionNum.append(xNum)
  PositionAnal.append(xAnal)

from matplotlib import pyplot
pyplot.plot(Time, PositionNum)
pyplot.plot(Time, PositionAnal)
pyplot.xlabel('t [s]');
pyplot.ylabel('x [m]');