A previously reported algorithm for solving the problem of coupling electromagnetic energy in the LHRF from a phased array of identical rectangular waveguides to a plane-stratified, magnetized cold plasma is numerically implemented. The resulting computer codes are sufficiently general to allow for an arbitrary number of waveguides with finite dimensions in both poloidal and toroidal directions and are thus capable of computing coupling to both slow and fast waves in the plasma. Some of the details of the implementation and the extension of the algorithm to allow study of the Fourier spectrum of slow and fast waves launched by the array are discussed. Good agreement is found with previously reported, less general work for the slow wave launching case. The effect of phasing multirow arrays in the poloidal direction is studied, and an asymmetric between phasing 'up' and 'down' is found that persists in the case where the plasma adjacent to the array is uniform. A 4 × 3 array designed to launch fast waves of high phase velocity is studied. By using the optimal poloidal phasing, low reflection coefficients (|R|2 ⪅ 20%) are found under some not unrealistic edge plasma conditions, but most of the input power is trapped in the outermost layer of the plasma. Implications of our results for fast wave current drive experiments are discussed.