A stochastic algorithm for high-dimensional integrals over unbounded regions with Gaussian weight

Details are given for a Fortran implementation of an algorithm that uses stochastic spherical–radial rules for the numerical computation of multiple integrals over unbounded regions with Gaussian weight. The implemented rules are suitable for high-dimensional problems. A high-dimensional example from a computational finance application is used to illustrate the use of the rules.