Abstract. The paper discusses both theoretical properties and practical implementation of product integration rules of the formwhere / is continuous, k is absolutely integrable, the nodes {*",} are roots of the Hermite polynomials H"(x), and the weights {%,•} are chosen so that the rule is exact if/ is any polynomial of degree < n. Convergence of the rule to the exact integral as n -» oo is proved for a wide class of functions/and k (including singular or oscillatory functions k), and rates of convergence are estimated. The rules are shown to have the property of asymptotic positivity, and as a consequence exhibit good numerical stability. Numerical calculations for some practical cases are presented, which show the method to be computationally effective for integrands (including highly oscillatory ones) that decay suitably at infinity. Applications of the method to integration over [ 0, oo) are also discussed.