A numerical technique for obtaining the definite integral of a function whose smoothness, or lack thereof, is unknown, to a desired degree of accuracy, while doing only as much work as necessary on each subinterval of the interval in question.