Gaussian basis functons, routinely employed in molecular electronic structure calculations, can be combined with numerical grid-based functions in a discrete variable representation to provide an efficient method for computing molecular continuum wave functions. This approach, combined with exterior complex scaling, obviates the need for slowly convergent single-center expansions, and allows one to study a variety of electron-molecule collision problems. The method is illustrated by computation of various bound and continuum properties of H + 2 .