###### Abstract
The rigorous coupled wave analysis in bipolar coordinates has been used to investigate scattering from eccentric multi-cylinder configurations [24] and the T-matrix method has been employed to study scattering by a cylinder with eccentric cylindrical inclusions [25]. The more complicated case of bianisotropic cylindrical inclusions in a host cylinder which is embedded in a bianisotropic space, at oblique incidence, has also been treated [26]. The scattering of electromagnetic waves by structures that involve several parallel cylindrical rods has been considered by many researchers; a common geometric configuration comprises cylinders with a circular [27] or arbitrary cross section [28; 29] at normal incidence, while obliquely incident plane waves have been also examined [30]. More complicated problems/geometries include scattering from multiple dielectric and metallic objects studied by using the T-matrix method [31], oblique incidence at multilayered crossed arrays of circular cylinders [32], an S-matrix solution to scattering from periodic arrays of metallic cylinders with an arbitrary cross section [33], an iterative scattering procedure applied to parallel anisotropic chiral cylinders [34] and scattering from various periodic arrays of cylinders [35; 36; 37]. Over the last few years, interest in problems of scattering by cylindrical structures remains vivid, albeit the geometries/problems treated become even more sophisticated [38]. A series solution based on continuity has been presented for scattering from multilayer cylinders of arbitrary shapes [39; 40], a Voltera integral equation formulation has been proposed for scattering by radially inhomogeneous cylinders in the case of oblique incidence [41], the spectral integral method (SIM) has been employed in the case of scattering by multilayer magnetodielectric cylinders [42], a hybrid integral equation method has been investigated in the case of scattering from doubly connected cylinders at oblique incidence [43] and a generalized vector cylinder harmonics (VCH) expansion has been applied for the scattering of an inhomogeneous plane wave by various cylindrical structures, such as a multilayer cylinder and an ensemble of stratified cylinders [44; 45; 46]. Recently, multilayered cylindrical geometries have been adopted to simulate photonic crystals [47; 48]. The absorption and reflection of EM-waves by such crystals has been examined by using the transfer matrix method [47], while a multiphysics sensor consisting of magnetized plasma and isotropic dielectric layers with a periodic cylindrical structure has been proposed by [PERSON] et al. [48]. In this paper, we extend the method applied in [26] as a means to address and efficiently solve the problem of electromagnetic scattering of plane waves with arbitrary polarization, as they are obliquely incident upon a complex cylindrical structure of infinite length consisting of a group of parallel stratified circular lossy dielectric cylinders, embedded in a dielectric circular cylindrical housing surrounded by dielectric space. To the best of our knowledge, except [49] where scattering from a similar configuration is addressed as an engineering electromagnetics exercise from an educational point of view, this cylindrical structure is one of the most complicated configurations considered so far at oblique incidence. Admittedly, even more complex geometries may be found in the literature [23; 24]; however, the EM wave is assumed to impinge normally on these structures. Thus, the main contribution of this paper is a solution to EM-wave scattering by combining a rather complicated cylindrical configuration with EM waves at obliquely incidence, which, in addition, are arbitrarily polarized. Furthermore, a semi-analytical method is used to address the problem, the only approximation being the truncation of the infinite series and the (numerical) matrix inversion.