d**2*(-d**4 + 2*d**2 - 1)*(2*d**6 - 8*d**4*np.log(d) - 2*d**2)*((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)))/(64*(-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)) + d**2*(-d**4 + 2*d**2 - 1)*((2*d**6 - 8*d**4*np.log(d) - 2*d**2)*(-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) + (12*d**5 - 32*d**3*np.log(d) - 8*d**3 - 4*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*(2*d**6 - 8*d**4*np.log(d) - 2*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)))*(16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)/(64*(-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)) + d**2*(2*d**6 - 8*d**4*np.log(d) - 2*d**2)*((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 + 2*d**4*np.log(d) - 2*d**4 + 3*d**2 - 1)/(12*(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)) + d**2*((2*d**6 - 8*d**4*np.log(d) - 2*d**2)*(-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) + (12*d**5 - 32*d**3*np.log(d) - 8*d**3 - 4*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*(2*d**6 - 8*d**4*np.log(d) - 2*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)))*(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 + 2*d**4*np.log(d) - 2*d**4 + 3*d**2 - 1)/(12*(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)) + (-d**4 + 2*d**2 - 1)*(-d**8 + 4*d**6*np.log(d) + d**4)*((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)))/(24*(-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)) + (-d**4 + 2*d**2 - 1)*((2*d**6 - 8*d**4*np.log(d) - 2*d**2)*(-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) + (12*d**5 - 32*d**3*np.log(d) - 8*d**3 - 4*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*(2*d**6 - 8*d**4*np.log(d) - 2*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*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*(-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)) + (-d**8 + 4*d**6*np.log(d) + d**4)*((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 + 2*d**4*np.log(d) - 2*d**4 + 3*d**2 - 1)/(4*(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*d**6 - 8*d**4*np.log(d) - 2*d**2)*(-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) + (12*d**5 - 32*d**3*np.log(d) - 8*d**3 - 4*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*(2*d**6 - 8*d**4*np.log(d) - 2*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 + 2*d**4*np.log(d) - 2*d**4 + 3*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*(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)) + (-d**4 + 2*d**2 - 1)*(-d**8 + 4*d**6*np.log(d) + d**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)*(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)))/(16*d**2*(-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)) + (-d**4 + 2*d**2 - 1)*((-d**8 + 4*d**6*np.log(d) + d**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) + (-8*d**7 + 24*d**5*np.log(d) + 4*d**5 + 4*d**3)/(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*(-d**8 + 4*d**6*np.log(d) + d**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*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)/(16*d**2*(-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)) + (-d**8 + 4*d**6*np.log(d) + d**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)*(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 + 2*d**4*np.log(d) - 2*d**4 + 3*d**2 - 1)/(2*d**2*(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)) + ((-d**8 + 4*d**6*np.log(d) + d**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) + (-8*d**7 + 24*d**5*np.log(d) + 4*d**5 + 4*d**3)/(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*(-d**8 + 4*d**6*np.log(d) + d**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 + 2*d**4*np.log(d) - 2*d**4 + 3*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*d**2*(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))