This paper deals with numerical methods for the discretization of highly oscillatory systems. We approach the problem

by writing the solution in terms of the Magnus expansion based on the Cayley map. The global error, obtained when the

method is applied to the linear oscillator, is investigated. Moreover, we provide numerical experiments in order to validate

our theoretical results