Problem matched basis functions are proposed for the method of moments analysis of printed slot coupled microstrips. The appropriate equivalent currents of the integral equation kernel are represented in terms of two sets of entire domain basis functions. These functions synthesize on one hand the resonant behavior of slots, microstrips or dipoles and on the other hand the field in proximity of the feeding source and of the discontinuities. In order to define these basis functions, canonical geometries are identified, whose Green's functions have been found in semi-analytical form. The accuracy and the effectiveness of the method in terms of convergence rate and number of unknowns is demonstrated by comparison with a standard fine meshing full-wave analysis. The method is extremely convenient for large arrays, where the subwavelength details should be treated together with large global dimensions. Since the proposed solution is independent of the dimensions of these details, it provides dramatic reduction of the number of unknowns and improvement of condition number.