Abstract
In this paper we introduce and describe a new scheme for the numerical integration of smooth functions. The scheme is based on the modified Taylor expansion and is suitable for functions that exhibit near-sinusoidal or repetitive behaviour. We discuss the method and its rate of convergence, then implement it for the approximation of certain integrals. Examples include integrands involving Airy wave, Bessel, Gamma, and elliptic functions. The results, from the data in the tables, demonstrate that the method converges rapidly and approximates the integral as well as some well-known numerical integration methods used with sufficiently small step sizes.