We consider the efficient numerical solution of coupled dynamical systems, consisting a low dimensional nonlinear part and high linear time invariant part, e.g., stemming from spatial discretization an underlying partial differential equation. The subsystem can be eliminated in frequency domain for resulting integro-differential algebraic equations, we propose combination Runge–Kutta or multist...