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