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