A working code for many particles in three dimensions with periodic boundary conditions. Tests should include two- and three-particle Lennard-Jones gas. Reference calculations with the velocity Verlet integrators and performance tests should be included.