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.
© 2001-2024 Fundación Dialnet · Todos los derechos reservados