IMEX extensions of linear multistep methods with general monotonicity and boundedness properties

For solving hyperbolic systems with stiff sources or relaxation terms, time stepping methods should combine favorable monotonicity properties for shocks and steep solution gradients with good stability properties for stiff terms. In this paper we consider implicit–explicit (IMEX) multistep methods. Suitable methods will be constructed, based on explicit methods with general monotonicity and boundedness properties for hyperbolic equations. Numerical comparisons are made with several implicit–explicit Runge–Kutta methods.