(n, mode="numpy")
return E1r2Scheme(f"Gauss-Hermite ({n})", weights, points, 2 * n - 1)
Gauss-Hermite quadrature for integrals of the form

        int_{-inf}^{+inf} exp(-x^2) f(x) dx.
Gauss-Hermite quadrature for integrals of the form
Gauss - Hermite quadrature for integrals of the form
def gauss_hermite(n, mode="numpy"):
    """
    Gauss-Hermite quadrature for integrals of the form

        int_{-inf}^{+inf} exp(-x^2) f(x) dx.
    """
    if mode == "numpy":
        points, weights = np.polynomial.hermite.hermgauss(n)
    else:
        rc = orthopy.e1r2.RecurrenceCoefficients("physicists", "monic", symbolic=True)
        _, alpha, beta = np.array([rc[k] for k in range(n)]).T
        points, weights = scheme_from_rc(alpha, beta, rc.int_1, mode=mode)
    return E1r2Scheme(f"Gauss-Hermite ({n})", weights, points, 2 * n - 1)
def gauss_hermite ( n , mode = "numpy" ) : if mode == "numpy" : points , weights = np . polynomial . hermite . hermgauss ( n ) else : rc = orthopy . e1r2 . RecurrenceCoefficients ( "physicists" , "monic" , symbolic = True ) _ , alpha , beta = np . array ( [ rc [ k ] for k in range ( n ) ] ) . T points , weights = scheme_from_rc ( alpha , beta , rc . int_1 , mode = mode ) return E1r2Scheme ( f"Gauss-Hermite ({n})" , weights , points , 2 * n - 1 )