Gaussian beam (GB) representation methods are used to analyze the electromagnetic coupling into and the scattering by a large nonuniform cavity. The aperture field in the cavity is decomposed into beams using the Gabor expansion, and shooting techniques are then employed. The method is illustrated only for the two-dimensional (2-D) case. The GBs are tracked axially using the rules of beam optics which ignore any beam distortion upon reflection at the walls. The effects of beam distortion are not significant for relatively slowly varying waveguide cavities. The field scattered into the exterior by the termination within the cavity is found using a reciprocity integral formulation which requires a knowledge of the beam fields near the termination. Numerical results based on this GB approach are presented and compared with results based on an independent reference solution.>