Computational efficiency of solving the dynamics of highly oscillatory systems is an important issue due to the requirement of small step size of explicit numerical integration algorithms. A system is considered to have an oscillatory solution if it contains a fast solution that varies regularly about a slow solution. This paper investigates the use of the so-called Local Linearization Method (LLM) in the integration of multibody equations of motion that exhibit oscillatory behavior. The LLM is an exponential method that is based on the piecewise linear approximation of the equations through a firstorder Taylor expansion at each time step, where the solution at the next time step is determined by the analytic solution of the approximated linear system. In this paper the LLM is applied to simple examples. The results show that the LLM can improve computational efficiency, without jeopardizing the accuracy, when the multibody system is highly oscillatory.