A new numerical integral algorithm is presented according to the idea of multiresolution sampling. Its precision is improved by the forcing function. It's easy to use in high-dimensions integral by tensor extending method. The results show that it is not only in good agreement with traditional Simpson algorithm but also more precise to appreciate the high-vibrating integral.