d**6*ep**2*(2*d**2 - 2)*np.log(d)/(2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - 29*d**6*ep**2*(2*d**2 - 2)/(96*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + 5*d**4*ep**2*(2*d**2 - 2)/(16*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + 3*d**4*ep**2*(-4*d**4 + 4*d**2)*np.log(d)**2/(4*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - 3*d**4*ep**2*(-4*d**4 + 4*d**2)*np.log(d)/(4*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + d**4*ep**2*(-4*d**4 + 4*d**2)/(4*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + 3*d**4*ep**2*(4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)*np.log(d)/(4*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - d**4*ep**2*(4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)/(2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + 3*d**2*ep**2*(-4*d**4 + 4*d**2)*np.log(d)/(8*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - 9*d**2*ep**2*(-4*d**4 + 4*d**2)/(32*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + 9*d**2*ep**2*(4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)/(16*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + ep**2*(2*d**2 - 2)/(48*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - ep**2*(2*d**2 - 2)/(32*d**2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + ep**2*(-4*d**4 + 4*d**2)/(32*d**2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) - ep**2*(4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)/(16*d**2*(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1))