-6*d**6*ep**2*np.log(d)**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 8*d**6*ep**2*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 77*d**6*ep**2/(24*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - 6*d**4*ep**2*np.log(d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 11*d**4*ep**2/(2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - 9*d**2*ep**2/(4*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - ep**2/(6*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + ep**2/(8*d**2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1))
