def test(N, rtol=1e-15, atol=1e-14): x, w = orth._roots_hermite_asy(N) H = hermite_recursion(N+1, x) assert_allclose(H[-1,:], np.zeros(N), rtol, atol) assert_allclose(sum(w), sqrt(np.pi), rtol, atol)