where S is an oscillatory function and r is a large parameter and a and b are real and finite. Approximating I by usual numerical integration algorithms requires many function evaluations of g and S. Sometimes the evaluation of g is very expensive, and the computation becomes highly time consuming for very large values of r. The method in [2] is an extension of the collocation method presented in [1] for integrals of the type