Asymptotic solvers for ordinary differential equations with multiple frequencies

We construct asymptotic expansions for ordinary differential equations with highly oscillatory forcing terms, focussing on the case of multiple, non-commensurate frequencies. We derive an asymptotic expansion in inverse powers of the oscillatory parameter and use its truncation as an exceedingly effective means to discretize the differential equation in question. Numerical examples illustrate the effectiveness of the method.