我想使用辛普森的规则来评估以下双积分：

$$\int_{a}^{b}\left|\int_{0}^{z}x\cdot \mathrm{erf}(x-10)\cdot J_{0}(x) \mathrm{dx}\right|^{2}\exp(-0.5*(z-40)^2)\mathrm{dz}$$

通常，如果两个限制都是数字的，我只需申请两次就可以逃脱simps，scipy.integrate但在这里我不知道如何继续。