Python中的多维集成

Question

我想整合一个简单的函数f（r1，r2）= d / | r2-r1 |。其中r1 =（x1，y1，d），r2 =（x2，y2,0）。但以下情况永远存在：

import sympy as sp

antenna_h = 0.05 #m 
antenna_w = 0.05 #m
bb_h = 1 #m
bb_w = 1 #m
d = 0.30 #m

xa,ya,xb,yb = sp.symbols('xa,ya,xb,yb')
Aeff = sp.integrate(d/((xa-xb)**2+(ya-yb)**2+d**2)**0.5,
                    (xa,(-antenna_w/2,+antenna_w/2)),
                    (ya,(-antenna_h/2,+antenna_h/2)),
                    (xb,(-bb_w/2,+bb_w/2)),
                    (yb,(-bb_h/2,+bb_h/2)))
print(Aeff)

有人想知道如何在合理的时间内完成整合？