An integral representation technique, which was developed for calculating the modal properties of a single homogeneous weakly guided dielectric waveguide of arbitrary cross-sectional shape, is extended to the calculation of the modes of an ensemble of two or more such guides arrayed in arbitrary spatial configurations. Detailed numerical results are presented for two identical circular, two identical square, and two identical rectangular guides and are in excellent agreement with almost all published data. Similar calculations are made for two dissimilar guides. The technique is also applied to a weakly coupled linear array of N identical, circular, uniformly spaced guides, where N equals 3,4,5 and 6.