The code that I used for my bounded box simulations