if wfun == 4:  # Gauss-Jacobi
        return j_roots(n, alpha, beta)