Numerical methods that preserve a Lyapunov function for Ordinary Differential Equations | Publicación