How can I speed up scipy.integrate.quad?

Question

Consider the following code:

def disFunc(x,n):
        y = np.power(x, n-1)*np.exp(-x)
        return y
    

def disFunc2(t,n,beta):
    y = integrate.quad(d