Discrete gradient methods for solving ODEs numerically while preserving a first integral

Abstract

We show that all direct methods for preserving a first integral during the numerical integration of an ordinary differential equation fit into the unified framework of discrete gradient methods. Using this framework we construct several new integral-preserving schemes.