updated crystal.cpp <- changed the random generator to initialize
velocities + also added printing velocities when printing the whole
crystal (although not in SI units!!)