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

Compute a triple (definite) integral.

Description:

Return the triple integral of func3d(z, y,x) from x=a..b, y=gfun(x)..hfun(x), and z=qfun(x,y)..rfun(x,y)

Inputs:

func3d – a Python function or method of at least three variables in the
order (z, y, x).

(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). qfun – the lower boundary surface in z. It must be a function that takes

two floats in the order (x, y) and returns a float.

rfun – the upper boundary surface in z. (Same requirements as qfun.) args – extra arguments to pass to func3d. epsabs – absolute tolerance passed directly to the innermost 1-D quadrature

integration.

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

Outputs: (y, abserr)

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