In this paper, we present the details of a new method which is a hybrid met
hod of MF-BDF and MEBDF. To obtain this new method which we call MF-MEBDF,
we compose the matrix free properties of the first method and the accuracy
of MEBDF, elaborately. Application of this new method to some important sti
ff problems show that MF-MEBDF is generally faster than MEBDF and more accu
rate than MF-BDF and MEBDF. Since in MEBDF, the LU factorization is used, w
e expect MF-MEBDF to be more efficient than the other two methods for large
stiff systems of ODEs. (C) 1999 Elsevier Science Ltd. All rights reserved.