d**6*ep**2*(4*d/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (2*d**2 - 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)*np.log(d)/2 - 29*d**6*ep**2*(4*d/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (2*d**2 - 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)/96 + 5*d**4*ep**2*(4*d/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (2*d**2 - 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)/16 + 3*d**4*ep**2*((-16*d**3 + 8*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (-4*d**4 + 4*d**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)*np.log(d)**2/4 - 3*d**4*ep**2*((-16*d**3 + 8*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (-4*d**4 + 4*d**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)*np.log(d)/4 + d**4*ep**2*((-16*d**3 + 8*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (-4*d**4 + 4*d**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)/4 + 3*d**4*ep**2*((-16*d**3*np.log(d) + 8*d**3 - 8*d)*(4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 + (16*d**3*np.log(d) - 8*d**3 + 4*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1))*np.log(d)/4 - d**4*ep**2*((-16*d**3*np.log(d) + 8*d**3 - 8*d)*(4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 + (16*d**3*np.log(d) - 8*d**3 + 4*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1))/2 + 3*d**2*ep**2*((-16*d**3 + 8*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (-4*d**4 + 4*d**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)*np.log(d)/8 - 9*d**2*ep**2*((-16*d**3 + 8*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (-4*d**4 + 4*d**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)/32 + 9*d**2*ep**2*((-16*d**3*np.log(d) + 8*d**3 - 8*d)*(4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 + (16*d**3*np.log(d) - 8*d**3 + 4*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1))/16 + ep**2*(4*d/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (2*d**2 - 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)/48 - ep**2*(4*d/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (2*d**2 - 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)/(32*d**2) + ep**2*((-16*d**3 + 8*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + (-4*d**4 + 4*d**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)/(32*d**2) - ep**2*((-16*d**3*np.log(d) + 8*d**3 - 8*d)*(4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)**2 + (16*d**3*np.log(d) - 8*d**3 + 4*d)/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1))/(16*d**2)
