In this paper, an accurate and efficient method for the analysis of coupled perfectly conducting annular rings is presented. The problem is first formulated as a couple of Integral Equation (IEs) in the Vector Hankel Transform (VHT) domain, considered as unknowns in the cylindrical harmonics of the unknown surface current density. As a second step, Galerkin’s method is applied with suitable expansion functions. The selected functions have two main properties: they reconstruct the expected physical behavior of the nth cylindrical harmonic at the edges of the annular rings, and their VHT transform is analytical and can be expressed in closed-form. Consequently, the method is effective and the problem is regularized, as testified by the truncation error. Comparisons with the commercial software CST Microwave Studio have been carried out and are presented to validate the method.