(1) Department of Mathematics, and Centre for Mathematics and Statistics of Complex Systems, La Trobe University, Bundoora, Melbourne, 3086, Australia
Abstract:
We consider ordinary differential equations (ODEs) with a known Lyapunov function V. To ensure that a numerical integrator reflects the correct dynamical behaviour of the system, the numerical integrator should have V as a discrete Lyapunov function. Only second-order geometric integrators of this type are known for arbitrary Lyapunov functions. In this paper we describe projection-based methods of arbitrary order that preserve any given Lyapunov function. AMS subject classification (2000) 65L05, 65L06, 65L20, 65P40