(8*d**4*np.log(d) - 5*d**4 + 6*d**2 - 1)/(d**2*(8*PI_mu*d**2*np.log(d) - 2*PI_mu*d**2 + 2*PI_mu - d**2))