In this paper we present a high order method for the evaluation of
integrals of highly oscillatory fuctions. The integration problem is transformed
into an initial value problem for ordinary differential equations and
solved by Runge-Kutta methods. After the comparison with different
Runge-Kutta schemes, we propose an algorithm based on Gaussian method
of order $p=8$. A lot of examples are given.
