By assuming series solutions for the total field, the SMM method transforms the domain into a cylindrical problem. In this domain total field is written in terms of Bessel and Hankel function solutions to the cylindrical Helmholtz equation. SMM method formulation, finally helps compute these coefficients of the cylindrical harmonic functions within the cylinder and outside it, at the same time satisfying EM boundary conditions.
Finally, SMM accuracy can be increased by adding (removing) cylindrical harmonic terms used to model the scattered fields.
SMM, eventually leads to a matrix formalism, and the coefficients are calculated through matrix inversion. For N-cylinders, each scattered field modeled using 2M+1 harmonic terms, SMM requires to solve a N(2M + 1) system of equations.