A method-of-lines solution algorithm for reacting how problems modelled by
hyperbolic partial differential equations (PDEs) with stiff source terms is
presented. Monotonicity preserving advection schemes are combined with spa
ce/time error balancing and a Gauss-Seidel iteration to provide an efficien
t solver. Numerical experiments on two challenging examples are presented t
o illustrate the performance of the method. (C) 2001 Published by Elsevier
Science B.V. on behalf of IMACS.