def calc_pts_hg(npts=20):
    """Returns Hermite-Gauss quadrature points for even functions"""
    pts_hg, wts_hg = np.polynomial.hermite.hermgauss(npts*2)
    pts_hg = pts_hg[npts:]
    wts_hg = wts_hg[npts:] * np.exp(pts_hg*pts_hg)
    return pts_hg, wts_hg
def calc_pts_hg(npts=20):
    """Returns Hermite-Gauss quadrature points for even functions"""
    pts_hg, wts_hg = np.polynomial.hermite.hermgauss(npts*2)
    pts_hg = pts_hg[npts:]
    wts_hg = wts_hg[npts:] * np.exp(pts_hg*pts_hg)
    return pts_hg, wts_hg
def calc_pts_hg ( npts = 20 ) : pts_hg , wts_hg = np . polynomial . hermite . hermgauss ( npts * 2 ) pts_hg = pts_hg [ npts : ] wts_hg = wts_hg [ npts : ] * np . exp ( pts_hg * pts_hg ) return pts_hg , wts_hg
Returns Hermite-Gauss quadrature points for even functions
Returns Hermite - Gauss quadrature points for even functions