scipy.integrate.dblquad

scipy.integrate.dblquad(func, a, b, gfun, hfun, args=(), epsabs=1.4899999999999999e-08, epsrel=1.4899999999999999e-08)

Compute a double (definite) integral.

Description:

Return the double integral of func2d(y,x) from x=a..b and y=gfun(x)..hfun(x).

Inputs:

func2d – a Python function or method of at least two variables: y must be

the first argument and x the second argument.

(a,b) – the limits of integration in x: a < b gfun – the lower boundary curve in y which is a function taking a single

floating point argument (x) and returning a floating point result: a lambda function can be useful here.

hfun – the upper boundary curve in y (same requirements as gfun). args – extra arguments to pass to func2d. epsabs – absolute tolerance passed directly to the inner 1-D quadrature

integration.

epsrel – relative tolerance of the inner 1-D integrals.

Outputs: (y, abserr)

y – the resultant integral. abserr – an estimate of the error.

See also:

quad - single integral tplquad - triple integral fixed_quad - fixed-order Gaussian quadrature quadrature - adaptive Gaussian quadrature odeint, ode - ODE integrators simps, trapz, romb - integrators for sampled data scipy.special - for coefficients and roots of orthogonal polynomials