An adaptive quadrature method for the automatic computation of integrals with strongly oscillating integrand is presented. The integration method is based on a truncated Chebyshev series approximation. The algorithm uses a global subinterval division strategy. There is a protection against the influence of round-off errors. A Fortran implementation of the algorithm is given.
