Several nonstationary extensions of the generalized Langevin equation describing chemical reaction dynamics have involved a multiplicative noise term. This includes the cases of space dependent friction, time dependent friction, and nontrivial combinations thereof through the so-called irreversible generalized Langevin equation (iGLE). In the present work, a fourth-order numerical integration scheme that leads to dramatically reduced computation times is shown to be applicable in all of these cases.