Estimating a complicated integral numerically is a problem which is ubiquitous in applied mathematics. Often we have to deal with functions for which the integral has no simple closed form. Even worse, often we do not even have a formula for the function we are integrating, and instead can only query its value at specific points, for example by performing an experiment. We might, however, have ...