This paper introduces an extension of the classical BDF formulae, which int
egrate, without local truncation error, homogeneous linear ODEs with consta
nt coefficients. In particular, they integrate exactly scaler linear equati
ons whose characteristic roots have a large stiffness ratio. (C) 2001 Elsev
ier Science Ltd. All rights reserved.