This method is a particular implementation of a general technique that is often referred to as envelope following [25] and consists in replacing the derivatives of the ????1 slow aperiodic time scale by finite-difference approximations (e.g., the Backward Euler rule), to then obtain a set of successive boundary value problems with periodic boundary conditions