The basic problem in approximate integration is to compute an approximate solution to a definite integral $$\int_a^b f(x)\;dx$$ to a given degree of accuracy. If $f$ is a smooth function integrated over a small number of dimensions, and the domain of integration is bounded, there are many methods for approximating the integral to the desired accuracy.