The boundary-value problem of electromagnetic wave scattering from multiple circular apertures in a thick conducting plane is rigorously solved. The eigenfunction expansion, integral transform, and superposition principle are utilized to represent the scattered field in the discrete and continuous modes. The boundary conditions are enforced to obtain a set of simultaneous equations for the discrete modal coefficients. The transmission coefficient is represented in a fast convergent series. Computation is performed to illustrate the behavior of transmission and coupling through multiple circular apertures in terms of the aperture geometry.