Abstract-In the present paper, a combined method of auxiliary sources (MAS)-reaction matching (RM) approach is presented for the analysis of arrays of arbitrarily located cylindrical dipoles. It is shown that the addition of auxiliary monopole terminal sources to each array element results in a superior solution with regard to the numerical stability of the computed quantities, the behavior of the current distributions of the array elements and the resulting errors of the electric field boundary condition. Numerical results are presented for various representative array configurations, in order to illustrate the features of the proposed method and exhibit its advantages over conventional Method of Moments (MoM) schemes, especially in cases of moderately large-scale arrays. Finally, a few concluding remarks are discussed.