  We consider random Hermitian matrices made of complex or real $M\times N$
rectangular blocks, where the blocks are drawn from various ensembles. These
matrices have $N$ pairs of opposite real nonvanishing eigenvalues, as well as
$M-N$ zero eigenvalues (for $M>N$.) These zero eigenvalues are ``kinematical"
in the sense that they are independent of randomness. We study the eigenvalue
distribution of these matrices to leading order in the large $N,M$ limit, in
which the ``rectangularity" $r={M\over N}$ is held fixed. We apply a variety of
methods in our study. We study Gaussian ensembles by a simple diagrammatic
method, by the Dyson gas approach, and by a generalization of the Kazakov
method. These methods make use of the invariance of such ensembles under the
action of symmetry groups. The more complicated Wigner ensemble, which does not
enjoy such symmetry properties, is studied by large $N$ renormalization
techniques. In addition to the kinematical $\delta$-function spike in the
eigenvalue density which corresponds to zero eigenvalues, we find for both
types of ensembles that if $|r-1|$ is held fixed as $N\rightarrow\infty$, the
$N$ non-zero eigenvalues give rise to two separated lobes that are located
symmetrically with respect to the origin. This separation arises because the
non-zero eigenvalues are repelled macroscopically from the origin. Finally, we
study the oscillatory behavior of the eigenvalue distribution near the
endpoints of the lobes, a behavior governed by Airy functions. As $r\rightarrow
1$ the lobes come closer, and the Airy oscillatory behavior near the endpoints
that are close to zero breaks down. We interpret this breakdown as a signal
that $r\rightarrow 1$ drives a cross over to the oscillation governed by Bessel
functions near the origin for matrices made of square blocks.
