The scattering of a Gaussian beam by an array of parallel circular cylinders is investigated using an iterative scattering procedure. The cylinders are made of homogeneous, isotropic dielectric material, perfectly conducting material, or a combination thereof. The iterative scattering technique is developed to evaluate both the near and far fields due to any arbitrary number of interactions between the cylinders. The technique starts with the exact series representation of the initial scattered fields from each cylinder due to the incident Gaussian beam. Then, the sum of the initial scattered fields from all cylinders is considered to be the field of an incident wave on each individual cylinder. The additional scattered field contribution is then obtained and successive scattering takes place until the boundary conditions on the surface of every cylinder are satisfied. An error criterion is developed to predict the number of iterations required to satisfy certain predetermined percentage error in the near field computations. Transverse magnetic waves (TM) are considered in detail and the extension to transverse electric waves (TE) is outlined. The effects of the incident beam parameters on the far scattered field pattern are investigated and numerical results are presented for a few examples.