import random
import numpy as np

numbers = 10
tRes = 100
mass = 0.1
tStep = 1/tRes

simTime = 10

px = []
py = []
pz = []
vx = []
vy = []
vz = []

for i in range(numbers):
    vx.append()

for t in range(simTime*100):


    for i in range(numbers):
        for j in range(numbers):
            if(j > i):
                if(dist(px[i], py[i], pz[i], px[j], py[j], pz[j]) < 0.5):
                    temp = px[i]
                    px[i] = px[j]
                    px[j] = temp

                    temp = py[i]
                    py[i] = py[j]
                    py[j] = temp

                    temp = py[i]
                    py[i] = py[j]
                    py[j] = temp


    for i in range(numbers):
        px[i] = px[i] + vx[i]*tStep
        py[i] = py[i] + vy[i]*tStep
        pz[i] = pz[i] + vz[i]*tStep