from pylab import *

dt=0.001
N=550

x=empty(N+1)
z=empty(N+1)
vx=empty(N+1)
vz=empty(N+1)

m=0.057
g=-9.81

x[0]=0.
z[0]=1.5
vx[0] = 35.
vz[0] = 0.

for n in range(0,N):
    x[n+1]=x[n]+vx[n]*dt
    z[n+1]=z[n]+vz[n]*dt
    vx[n+1]=vx[n]         # Fx=0
    vz[n+1]=vz[n]+g*dt    # Fz=m*g

plot(x,z)

print(z[547])
print(x[340])
print(z[340])