We develop a Fourier–Galerkin spectral technique for computing the stationary solutions of 2D generalized wave equations. To this end a special complete orthonormal system of functions in L2(−∞,∞) is used for which product formula is available. The exponential rate of convergence is...