Asymptotic-numerical solvers for linear neutral delay differential equations with high- frequency extrinsic oscillations

We present a method to compute efficiently and easily solutions of systems of linear neutral delay differential equations with highly oscillatory forcing terms. This method is based on asymptotic expansions in inverse powers of a perturbed oscillatory parameter. Each term of the asymptotic expansion is derived by recursion. The cost of the computation is essentially independent of the oscillatory parameter. Numerical examples are provided and show that with few terms of the asymptotic expansion,  the solutions are approximated with high accuracy.