This study presents a full‐wave solution for arrays of cylindrical‐rectangular microstrip antennas, with superstrate. It uses the method of moments (MoM), with full‐domain basis functions. The same set of functions is used as testing‐functions resulting in a Galerkin procedure. A numerical procedure is presented allowing an efficient and numerically stable evaluation of the elements of the MoM linear system. Results for input impedance are compared to those using HFSS for a flush‐mounted geometry. Results for radiation pattern, mutual impedance, and coupling coefficient S12 are also presented. And the effects of the superstrate on those parameters are evaluated.