as well as for many other one and muti-dimensional integrals of functions with rapid irregular oscillations. the analysis provides a convergence rate and it shows that the relative error of the method is even decreasing as the frequency of the oscillations increases