| quad(func, a, b[, args, full_output, ...]) | Compute a definite integral. | 
| dblquad(func, a, b, gfun, hfun[, args, ...]) | Compute the double integral of func2d(y,x) from x=a..b and y=gfun(x)..hfun(x). | 
| tplquad(func, a, b, gfun, hfun, qfun, rfun) | Compute a triple (definite) integral. | 
| fixed_quad(func, a, b[, args, n]) | Compute a definite integral using fixed-order Gaussian quadrature. | 
| quadrature(func, a, b[, args, tol, maxiter, ...]) | Compute a definite integral using fixed-tolerance Gaussian quadrature. | 
| romberg(function, a, b[, args, tol, show, ...]) | Romberg integration of a callable function or method. | 
| trapz(y[, x, dx, axis]) | Integrate along the given axis using the composite trapezoidal rule. | 
| cumtrapz(y[, x, dx, axis]) | Cumulatively integrate y(x) using samples along the given axis and the composite trapezoidal rule. | 
| simps(y[, x, dx, axis, even]) | Integrate y(x) using samples along the given axis and the composite | 
| romb(y[, dx, axis, show]) | Romberg integration using samples of a function | 
See also
scipy.special for orthogonal polynomials (special) for Gaussian quadrature roots and weights for other weighting factors and regions.