Starting point of an MD simulation of non-interacting particles.

exercise_1_start.py

import random
import pylab as plt

npart = 100
nsteps = 10000
dt = 0.001

X = [random.random() for i in range(npart)]
Y = [2*(random.random()-0.5) for i in range(npart)]

plt.plot(X,Y,'ro')
plt.axis((-1,1,-1,1))
plt.savefig('startcoord.png')