This paper is concerned with the asymptotic expansion and numerical solution of systems of linear delay differential equations with highly oscillatory forcing terms. The computation of such problems using standard numerical methods is exceedingly slow and inefficient, indeed standard software is practically useless for this purpose. We propose an alternative, consisting of an asymptotic expansion of the solution, where each term can be derived either by recursion or by solving a non-oscillatory problem. This leads to methods which, counter-intuitively to those developed according to standard numerical reasoning, exhibit improved performance with growing frequency of oscillation.
