I need to implement a root-searching algorithm for a function, that can possibly have some (non-removable) divergence points (DPs), e.g. f(x) = x**2 * np.tan(x).