Reducing round-off errors in symmetric multistep methods

Certain symmetric linear multistep methods have an excellent long-time behavior when applied to second order Hamiltonian systems with or without constraints. For high accuracy computations round-off can be the dominating source of errors. This article shows how symmetric multistep methods should be implemented such that round-off errors are minimized and propagate like a random walk.