“…For this reason, among others, we focus our attention on an alternative formulation of (2) called Rosenbrock-Krylov [14] or exponential-Krylov [13] methods, in which the coefficients are chosen to permit a specific, Krylov-based, approximation A n of the Jacobian matrix J n . These methods have been shown to be at least efficient for a broad spectrum of problems through comparisons of relative performance against an assortment of standard integrators [6,[12][13][14]16], and have been effectively used to study engineering problems in [5,11,17,18].…”