-d**3*((64*d**3*np.log(d) - 32*d**3 + 32*d)/(d**4*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)) + (16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)*(-16*PI_mu*d**3*np.log(d) + 12*PI_mu*d**3 - 16*PI_mu*d + 4*PI_mu/d + 16*d**3*np.log(d) - 8*d**3 + 8*d)/(d**4*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) - 4*(16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)/(d**5*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)))*(8*d**4*np.log(d) - 5*d**4 + 6*d**2 - 1)/(24*(8*PI_mu*d**2*np.log(d) - 2*PI_mu*d**2 + 2*PI_mu - d**2)) - d**3*((64*d**3*np.log(d) - 32*d**3 + 32*d)/(d**4*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)) + (16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)*(-16*PI_mu*d**3*np.log(d) + 12*PI_mu*d**3 - 16*PI_mu*d + 4*PI_mu/d + 16*d**3*np.log(d) - 8*d**3 + 8*d)/(d**4*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) - 4*(16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)/(d**5*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)))*(2*PI_mu*d**4*np.log(d) - 2*PI_mu*d**4 + 4*PI_mu*d**2 - 2*PI_mu*np.log(d) - 2*PI_mu - 6*d**4*np.log(d) + 4*d**4 - 5*d**2 + 1)/(4*(32*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu - 4*d**2)) - d*((-16*PI_mu*d**3*np.log(d) + 12*PI_mu*d**3 - 16*PI_mu*d + 4*PI_mu/d + 16*d**3*np.log(d) - 8*d**3 + 8*d)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(d**4*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (48*PI_mu*d**5*np.log(d) - 40*PI_mu*d**5 + 64*PI_mu*d**3 - 16*PI_mu*d*np.log(d) - 24*PI_mu*d - 96*d**5*np.log(d) + 56*d**5 - 64*d**3 + 8*d)/(d**4*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)) - 4*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(d**5*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)))*(8*d**4*np.log(d) - 5*d**4 + 6*d**2 - 1)/(16*(8*PI_mu*d**2*np.log(d) - 2*PI_mu*d**2 + 2*PI_mu - d**2)) - d*((-16*PI_mu*d**3*np.log(d) + 12*PI_mu*d**3 - 16*PI_mu*d + 4*PI_mu/d + 16*d**3*np.log(d) - 8*d**3 + 8*d)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(d**4*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (48*PI_mu*d**5*np.log(d) - 40*PI_mu*d**5 + 64*PI_mu*d**3 - 16*PI_mu*d*np.log(d) - 24*PI_mu*d - 96*d**5*np.log(d) + 56*d**5 - 64*d**3 + 8*d)/(d**4*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)) - 4*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(d**5*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)))*(2*PI_mu*d**4*np.log(d) - 2*PI_mu*d**4 + 4*PI_mu*d**2 - 2*PI_mu*np.log(d) - 2*PI_mu - 6*d**4*np.log(d) + 4*d**4 - 5*d**2 + 1)/(2*(32*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu - 4*d**2)) + (8*d**4*np.log(d) - 5*d**4 + 6*d**2 - 1)*(16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)**2/(64*d**2*(8*PI_mu*d**2*np.log(d) - 2*PI_mu*d**2 + 2*PI_mu - d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)**2*(2*PI_mu*d**4*np.log(d) - 2*PI_mu*d**4 + 4*PI_mu*d**2 - 2*PI_mu*np.log(d) - 2*PI_mu - 6*d**4*np.log(d) + 4*d**4 - 5*d**2 + 1)/(12*d**2*(32*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu - 4*d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (8*d**4*np.log(d) - 5*d**4 + 6*d**2 - 1)*(16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(24*d**4*(8*PI_mu*d**2*np.log(d) - 2*PI_mu*d**2 + 2*PI_mu - d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)*(2*PI_mu*d**4*np.log(d) - 2*PI_mu*d**4 + 4*PI_mu*d**2 - 2*PI_mu*np.log(d) - 2*PI_mu - 6*d**4*np.log(d) + 4*d**4 - 5*d**2 + 1)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(4*d**4*(32*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu - 4*d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (8*d**4*np.log(d) - 5*d**4 + 6*d**2 - 1)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)**2/(16*d**6*(8*PI_mu*d**2*np.log(d) - 2*PI_mu*d**2 + 2*PI_mu - d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (2*PI_mu*d**4*np.log(d) - 2*PI_mu*d**4 + 4*PI_mu*d**2 - 2*PI_mu*np.log(d) - 2*PI_mu - 6*d**4*np.log(d) + 4*d**4 - 5*d**2 + 1)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)**2/(2*d**6*(32*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu - 4*d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2)