Slow cubature

[stevengj]

@e3c6

but it seems to support only 1-dimensional quadrature, there is no method for cubature.

For cubature you can just 1d quadrature rules along each dimension (a “tensor product” of 1d rules). See e.g. Gauss quadrature in 3D with change of variables - #2 by stevengj