def _j_roots_newton(n, alpha, beta, releps=3e-14, max_iter=10):