重现原始问题的一部分以简化问题。
我需要计算这个函数的积分
def f(z):
return ((1-2*z)*np.exp(-d/z))/(((1-z)**(2+d))*(z**(2-d)))
这里d是一个常数。我正在使用这个功能:
import scipy.integrate as integrate
return integrate.quad(f,0, 5)
但是,我需要为许多不同的d值计算。我尝试将所有d值放在一个numpy数组中,但它没有用。知道如何对此进行矢量化吗?