The Gabor expansion of a function f_(x) decomposes it into a double sum over integers m and n of a product of basis functions g(x_m X) and Fourier_series exponentials exp(2πi n/X) for given spacing X. The choice of basis function determines the coefficients a m n of the expansion. If f_(x) is band limited, the double sum can for all practical purposes be replaced by a single sum over Gaussian basis functions. This is extremely useful for expansion of multidimensional functions such as beams in phase space. Conditions of validity are given, and several examples illustrate the technique.