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