studied in averaging methods [1]. Here f and g are assumed to be periodic in φ and bounded as ε → 0. In the standard terminology of multiscale analysis, we would call x and φ the fast variables of these systems, y and I the slow variables. But we also have a particular interest on systems for which the fast and slow variables exist but cannot be explicitly identified beforehand. At the present time, there does not exist a unified strategy for dealing with both problems of type (1.1) and (1.2). There is, however, a large literature on stiff systems of the type (1.1) and oscillatory systems of the type (1.2) separately. Starting from the pioneering work of Dahlquist and Gear [5], there has been extensive work on designing efficient numerical methods for stiff ODEs [7], such as the backward differentiation formula, implicit Runge-Kutta methods [5], extrapolation methods and Rosenbrock methods [7, 8, 6]. There is also extensive work on oscillatory systems, some of which are analytical [1], and some are numerical [9]. We will apply the framework of the heterogeneous multiscale method (HMM) [3], which is a general methodology for problems with multiscales. In the case of two scales, a macroscale and a microscale, HMM consists of two components: selection of a conventional macroscale solver, here a standard ODE solver, and estimating the effective forces used in the macroscale solver by performing numerical experiments using the microscale model and processing the data obtained. The procedure can be iterated if the system has more than two separated scales. There are a number of related numerical methods for systems with multiple time scales, in particular for stiff ODEs. Even though the most popular numerical methods for stiff ODEs stem from implicit methods such as the backward differentiation formula, a number of explicit methods have also been proposed [8, 2, 7]. In the simplest version, these methods are Runge-Kutta in nature, each stage of which is a forward
[1]
J. M. Watt.
Numerical Initial Value Problems in Ordinary Differential Equations
,
1972
.
[2]
V. Arnold.
Mathematical Methods of Classical Mechanics
,
1974
.
[3]
F. Krogh,et al.
Solving Ordinary Differential Equations
,
2019,
Programming for Computations - Python.
[4]
Linda R. Petzold,et al.
Numerical solution of highly oscillatory ordinary differential equations
,
1997,
Acta Numerica.
[5]
Assyr Abdulle,et al.
Fourth Order Chebyshev Methods with Recurrence Relation
,
2001,
SIAM J. Sci. Comput..
[6]
Eric Vanden-Eijnden,et al.
NUMERICAL TECHNIQUES FOR MULTI-SCALE DYNAMICAL SYSTEMS WITH STOCHASTIC EFFECTS ⁄
,
2003
.
[7]
Ioannis G. Kevrekidis,et al.
Projective Methods for Stiff Differential Equations: Problems with Gaps in Their Eigenvalue Spectrum
,
2002,
SIAM J. Sci. Comput..
[8]
Björn Engquist,et al.
Heterogeneous multiscale methods for stiff ordinary differential equations
,
2005,
Math. Comput..