Numerical Integration in PyTorch

Hi,

Is it possible to perform multi-dimensional numerical integration using PyTorch? I have a particular interest in Gauss Quadrature schemes, but other approaches are also welcomed.

Thanks!