Velocity Verlet Algorithm
Roundoff advantage of leapfrog, but better treatment of velocities
Implemented in stages
- evaluate current force
- compute r at new time
- add current-force term to velocity (gives v at half-time step)
- compute new force
- add new-force term to velocity
Also mathematically equivalent to Verlet algorithm (in giving values of r)