使用`scipy`固定$ n $维积分的阶高斯正交

Question

我希望使用scipy.integrate.nquad方法来集成n变量（为简单起见，使用双变量）函数。但是，我也想像scipy.integrate.fixed_quad中那样使用定序正交规则，但是我找不到找到这种方法的方法。

这是一个已知问题吗，有什么聪明的方法吗？

提前谢谢！