A new method of computation of wave fields using Gaussian beams