Numerical path integration with Coulomb potential

A simple and efficient method for quantum Monte Carlo simulation is presented, based on discretization of the action in the path integral, and a Gaussian averaging of the potential, which works well e.g. with the Coulomb potential.