from scipy import integrate
import numpy as np
from mpmath import coth
DvDc = 6.5
dens = 5.65
Vs = 4.6e3
a = 3.3e-9
w0 = np.sqrt(2)*Vs/a
T = 50
kb = 1.38064852e-5 # in eV`
j0 = (DvDc)**2 / ((2*np.pi)**2 *dens*Vs**5)
def func(x):
return x*np.exp(-(x/w0)**2)*coth(x/(2*kb*T))
S = j0*integrate.quad(func, 0, np.inf)[0]
print(S)
嗨，所以我没有使用数值积分的经验，但我想知道是否有人可以提供帮助。在这里，我定义了一些变量，并有一个要集成的函数。但是我得到了“积分可能是发散的，或者慢慢收敛的。” 尽管我在 coth 的分母中添加了内容，但几乎得到了相同的结果。
有谁知道我是否犯了一个错误，或者我只是不能做到这一点。
谢谢！