erri = 0.5 * (abserr / (np.abs(fy0) + np.abs(fy) + _TINY + tol))