Abstract Many complex applications require the solution of initial-value problems where some components change fast, while others vary slowly. Multirate schemes apply different step sizes to resolve system, according their dynamics, in order achieve increased computational efficiency. The stiff fast or slow, are best discretized with implicit base methods ensure numerical stability. To this end...