[SciPy-user] Triple integral with Gaussian quadrature

Hi all, Is there a way to compute a triple integral via Gaussian quadrature in scipy ? A small example would be appreciated. What integration method is used in tplquad(func, a, b, gfun, hfun, qfun, rfun, args=(), epsabs=1.4899999999999999e-08, epsrel=1.4899999999999999e-08) ? quadrature(func, a, b, args=(), tol=1.4899999999999999e-08, maxiter=50) Nils