Fully symmetric interpolatory rules for multiple integrals over infinite regions with Gaussian weight

Fully symmetric interpolatory integration rules are constructed for multidimensional integrals over infinite integration regions with a Gaussian weight function. The points for these rules are determined by successive extensions of the one-dimensional three-point Gauss-Hermite rule. The new rules are shown to be efficient and only moderately unstable.