-48*d**6*ep**2*(16*d**2*(2*d**2*np.log(d) - d**2 + 1)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**2*np.log(d) + d**2 - 1)*np.log(d)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 + 64*d**6*ep**2*(16*d**2*(2*d**2*np.log(d) - d**2 + 1)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**2*np.log(d) + d**2 - 1)*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 - 77*d**6*ep**2*(16*d**2*(2*d**2*np.log(d) - d**2 + 1)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**2*np.log(d) + d**2 - 1)/(3*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2) - 36*d**5*ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)*np.log(d)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 + 36*d**5*ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 - 45*d**5*ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)/(4*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2) - 36*d**4*ep**2*np.log(d)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 20*d**4*ep**2*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 2*d**4*ep**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 48*d**4*ep**2*(16*d**2*(2*d**2*np.log(d) - d**2 + 1)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**2*np.log(d) + d**2 - 1)*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 + 44*d**4*ep**2*(16*d**2*(2*d**2*np.log(d) - d**2 + 1)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**2*np.log(d) + d**2 - 1)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 - 24*d**3*ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 + 16*d**3*ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 - 18*d**2*ep**2*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 7*d**2*ep**2/(2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - 18*d**2*ep**2*(16*d**2*(2*d**2*np.log(d) - d**2 + 1)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**2*np.log(d) + d**2 - 1)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 - 9*d*ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)/(2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2) - 2*ep**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 4*ep**2*(16*d**2*(2*d**2*np.log(d) - d**2 + 1)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**2*np.log(d) + d**2 - 1)/(3*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2) + ep**2/(2*d**2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + ep**2*(16*d**2*(2*d**2*np.log(d) - d**2 + 1)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**2*np.log(d) + d**2 - 1)/(d**2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2) - ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)/(4*d**3*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2)
