各位大佬这是我的代码，但是计算后报错，写的是The integral is probably divergent, or slowly convergent.(积分可能是发散的，或者是缓慢收敛的),有大佬知道这个该怎么更改吗？
from scipy.integrate import dblquad
import math as ma
a=c=5
# 定义被积函数，这里以 x*y 的形式为例
def f(x, y):
return ma.sqrt(2*ma.sqrt((x-ma.sqrt(x**2*a**2/(x**2+y**2)))**2+(y-ma.sqrt(y**2*a**2/(x**2+y**2)))**2))/(ma.pi**1.5*(x-a)**2+(y)**2)
h = lambda x :ma.sqrt(a**2-x**2)
# 计算二重积分，注意顺序是先内后外
result, error = dblquad(f, -5,5,0,h)
print('二重积分的结果是:', result)