如何实现辛普森规则的双积分（无内部积分的数值限制）

Question

并希望每个人都安全。

我想使用Simpson规则来评估以下双重积分：

[通常，如果两个限制都是数字，那么我只需从simps两次应用scipy.integrate就可以避免，但是在这里我不知道如何进行。任何帮助将非常感激。我知道我可以使用dblquad来执行此操作，但这将非常缓慢。在这里，我更关注速度而不是准确性（在某种程度上是必须的）。如果还有其他方法可以解决，请指导我。

Answer 1

只需在一系列z值上对x进行积分，对z进行插值和积分。

为了提高速度，我怀疑内置Quad比Simp好。