Molecular Dynamics algorithm

start with configuration 

compute forces 

integrate acceleration numerically to obtain new configuration

go back to 1.