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)

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