working molecular dynamics code - vel verlet