Computing Highly Oscillatory Integrals