Asymptotic solvers for oscillatory systems of differential equations