-d**4*(16*PI_mu*d**2*np.log(d) - 4*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)*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*(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)/(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)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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*(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)*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*(16*PI_mu*d**2*np.log(d) - 4*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)*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*(16*PI_mu*d**2*np.log(d) - 4*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)/(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)*(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)*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*(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)/(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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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) - 4*d**2)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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) - 4*d**2)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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*(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**4 + 16*PI_mu*d**2 + 8*d**4)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)*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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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) - (-16*PI_mu*d**2*np.log(d) + 4*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)*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*(-16*PI_mu*d**2*np.log(d) + 4*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)*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*(-16*PI_mu*d**2*np.log(d) + 4*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)/(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*(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)*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*(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)/(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*(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)*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*(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)/(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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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) + 4*d**2)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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) - 8*PI_mu*d**2 + 8*PI_mu)/(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*(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)*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*(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)/(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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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) - 8*PI_mu*d**2 + 8*PI_mu)/(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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)*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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)*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**4 + 16*PI_mu*d**2 + 8*d**4)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)*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**4 + 16*PI_mu*d**2 + 8*d**4)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)*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)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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*(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)*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**2*np.log(d) + 4*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)/(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*(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)/(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*(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)/(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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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) + 4*d**2)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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*(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)*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*(16*PI_mu*d**2*np.log(d) - 4*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)/(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) - (16*PI_mu*d**2*np.log(d) - 4*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)/(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) - 5*(16*PI_mu*d**2*np.log(d) - 4*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)/(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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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) - 8*PI_mu*d**2 + 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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)*(-16*PI_mu*d**4 + 16*PI_mu*d**2 + 8*d**4)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)*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**4 + 16*PI_mu*d**2 + 8*d**4)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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) + 4*d**2)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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) - 8*PI_mu*d**2 + 8*PI_mu)/(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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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) - 4*d**2)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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) - 4*d**2)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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) - 7*(16*PI_mu*d**2*np.log(d) - 4*d**2)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)*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)*(-8*PI_mu*d**4 + 16*PI_mu*d**2 - 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(-16*PI_mu*d**2*np.log(d) + 8*PI_mu*d**2 - 8*PI_mu)/(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**4 + 16*PI_mu*d**2 + 8*d**4)*(16*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 + 8*PI_mu)/(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)