Boththe penetrated and the scattered fields in the presence of a non-uniform cold magnetized plasma column are investigated in this work, and the direct integral equation is formulated for describing the currents on the surrounded multiple perfectly conducting strips on the outer boundary of plasma column due to a normally incident TMwave. In order to formulate such an integral equation, two-dimensional dyadic Green's functions in multi-layered cylindrical shells are derived at first in terms of the cylindrical vector wave functions. Furthermore, by using the Galerkin's method, the surface current of the strips is expanded in the form of a series of Chebyshev polynomials of the first kind. Parameter studies are performed to show the variation of the near- and far-zone fields with the operating frequency, the effective collision frequency, as well as the angular width and the number of surrounded strips.