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