Векторизация интеграла в scipy

Question

Воспроизведение части исходного вопроса для упрощения проблемы.

Мне нужно рассчитать интеграл этой функции

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, но это не сработало. Есть идеи, как это векторизовать?