1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
from pylab import *

t = 0.
dt = .1
x = 0
v = 2

ta = []    # list for storing time
xa = []    # and instantaneous number of atoms

while t < 5:
  ta.append(t)
  xa.append(x)
  x = x + v * dt
  t = t + dt

plot(ta,xa)
show()