-6*d**6*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 + 8*d**6*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 - 77*d**6*ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)/(24*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2) - 12*d**5*ep**2*np.log(d)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 12*d**5*ep**2*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 4*d**5*ep**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 6*d**4*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 + 11*d**4*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) - 6*d**3*ep**2*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 9*d**3*ep**2/(2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - 9*d**2*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) - ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)/(6*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2) - ep**2/(2*d*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + ep**2*(-16*d**3*np.log(d) + 8*d**3 - 8*d)/(8*d**2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2)