import numpy as np
from scipy.special import h_roots

def GHquad(n : int):
    """Gauss-Hermite quadrature weights"""
    Z, w = h_roots(n)
    return Z*np.sqrt(2), w/np.sqrt(np.pi)