(-d**4 + 2*d**2 - 1)*(2*d**6 - 8*d**4*np.log(d) - 2*d**2)*(16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)/(64*d**2*(-2*PI_mu*d**2 + 2*PI_mu + d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (2*d**6 - 8*d**4*np.log(d) - 2*d**2)*(16*d**4*np.log(d) - 12*d**4 + 16*d**2 - 4)*(-2*PI_mu*d**4*np.log(d) + 2*PI_mu*d**4 - 4*PI_mu*d**2 + 2*PI_mu*np.log(d) + 2*PI_mu + 2*d**4*np.log(d) - 2*d**4 + 3*d**2 - 1)/(12*d**2*(8*PI_mu*d**2 - 8*PI_mu - 4*d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (-d**4 + 2*d**2 - 1)*(2*d**6 - 8*d**4*np.log(d) - 2*d**2)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(24*d**4*(-2*PI_mu*d**2 + 2*PI_mu + d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (2*d**6 - 8*d**4*np.log(d) - 2*d**2)*(-2*PI_mu*d**4*np.log(d) + 2*PI_mu*d**4 - 4*PI_mu*d**2 + 2*PI_mu*np.log(d) + 2*PI_mu + 2*d**4*np.log(d) - 2*d**4 + 3*d**2 - 1)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(4*d**4*(8*PI_mu*d**2 - 8*PI_mu - 4*d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (-d**4 + 2*d**2 - 1)*(-d**8 + 4*d**6*np.log(d) + d**4)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(16*d**6*(-2*PI_mu*d**2 + 2*PI_mu + d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2) + (-d**8 + 4*d**6*np.log(d) + d**4)*(-2*PI_mu*d**4*np.log(d) + 2*PI_mu*d**4 - 4*PI_mu*d**2 + 2*PI_mu*np.log(d) + 2*PI_mu + 2*d**4*np.log(d) - 2*d**4 + 3*d**2 - 1)*(8*PI_mu*d**6*np.log(d) - 8*PI_mu*d**6 + 16*PI_mu*d**4 - 8*PI_mu*d**2*np.log(d) - 8*PI_mu*d**2 - 16*d**6*np.log(d) + 12*d**6 - 16*d**4 + 4*d**2)/(2*d**6*(8*PI_mu*d**2 - 8*PI_mu - 4*d**2)*(4*PI_mu*d**4*np.log(d) - 4*PI_mu*d**4 + 8*PI_mu*d**2 - 4*PI_mu*np.log(d) - 4*PI_mu - 4*d**4*np.log(d) + 3*d**4 - 4*d**2 + 1)**2)