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