(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2))