-d**4*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)/(24*(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**4*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(144*(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**4*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(64*(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**2*(-16*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)/(8*(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**2*(-16*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(32*(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**2*(-16*PI_mu*d**2*np.log(d) + 4*d**2)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(24*(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**2*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)/(16*(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**2*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(64*(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**2*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(48*(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**2*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(48*(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**2*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)**2/(8*(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) + 3*d**2*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)/(32*(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) - 3*d**2*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(128*(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**2*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)/(24*(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) + 5*d**2*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(288*(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*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)**2/(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) + (-16*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(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) - (-16*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(8*(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*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)*np.log(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) + (-16*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(8*(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*PI_mu*d**2*np.log(d) + 4*d**2)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(16*(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*PI_mu*d**2*np.log(d) + 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)/(16*(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*PI_mu*d**2*np.log(d) + 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(64*(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*PI_mu*d**2*np.log(d) + 4*d**2)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(16*(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*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)/(8*(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*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(16*(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*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)/(8*(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*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(16*(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*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)*np.log(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) + (16*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(8*(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*PI_mu*d**2*np.log(d) - 4*d**2)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(32*(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*PI_mu*d**2*np.log(d) - 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(16*(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*PI_mu*d**2*np.log(d) - 4*d**2)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(16*(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 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)**3/(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) + (PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)**2/(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) - 9*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)/(16*(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) + 9*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(32*(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 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)*np.log(d)**2/(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) + 3*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)*np.log(d)/(8*(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) - 3*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(16*(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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)/(16*(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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(32*(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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)**2/(16*(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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)/(16*(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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(32*(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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)*np.log(d)/(16*(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) + 3*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(64*(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*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(32*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*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(8*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*PI_mu*d**2*np.log(d) + 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(8*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*PI_mu*d**2*np.log(d) + 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)/(8*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*PI_mu*d**2*np.log(d) + 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(16*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*PI_mu*d**2*np.log(d) + 4*d**2)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(8*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) - 5*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(72*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) + 3*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(64*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*PI_mu*d**2*np.log(d) - 4*d**2)*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(8*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*PI_mu*d**2*np.log(d) - 4*d**2)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(16*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*PI_mu*d**2*np.log(d) - 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(16*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*PI_mu*d**2*np.log(d) - 4*d**2)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(8*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) - (PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)**2/(16*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) + 3*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(128*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) + (PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)**2/(16*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) - 9*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(32*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) - (PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)*np.log(d)**2/(8*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) + 3*(PI_mu*(-d**4 + 2*d**2 - 1)/(d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)) - 1)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(16*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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(32*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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)**2/(8*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) - 3*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)/(16*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) + 5*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(32*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*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)*np.log(d)/(8*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) - 3*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(16*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*PI_mu*d**2*np.log(d) + 4*d**2)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(48*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) + 3*(-16*PI_mu*d**2*np.log(d) + 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(64*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) - (-16*PI_mu*d**2*np.log(d) + 4*d**2)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(16*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) - 7*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(192*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) + (16*PI_mu*d**2*np.log(d) - 4*d**2)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(48*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) - (16*PI_mu*d**2*np.log(d) - 4*d**2)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(16*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) + (-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)*np.log(d)/(16*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) + 13*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 4*PI_mu*d**2 - 2*d**4 + 2*d**2)/(288*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) - (-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)*np.log(d)/(16*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) - (-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-4*PI_mu*d**6 + 8*PI_mu*d**4 - 4*PI_mu*d**2 + 4*d**6 - 4*d**4)/(8*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) + (-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)*np.log(d)/(8*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) + 9*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(4*PI_mu*d**6*np.log(d) - 4*PI_mu*d**6 + 4*PI_mu*d**4 + 4*PI_mu*d**2*np.log(d) - 4*d**6*np.log(d) + 3*d**6 - 2*d**4 - d**2)/(64*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)