Abstract
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.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have