An experiment on random matrices

The distribution of the bulk of the spectrum is an example of the circular law. For the model you selected (where each entry is uniformly chosen at random from an interval), the law was first proven by Bai (at least in the case where the entries are normalised to have mean zero), building upon previous work of Girko; the non-central case (non-zero mean) was recently established by Chafai.

The non-central case is a rank one perturbation of the central case (by the matrix whose entries are all equal to the mean) and should therefore cause one exceptional eigenvalue.
(In the central case it can be shown that the spectral radius is close to the radius of the disk, so there are basically no exceptional eigenvalues.) If the mean of each entry is $\mu$, then the rank one perturbation has an eigenvalue at $\mu n$, so one expects the exceptional eigenvalue to linger near this number. In your case, $\mu = 0.05$ and $n = 301$, so the exceptional eigenvalue should linger near $15.05$.

There is a paper of Silverstein which makes the above heuristics precise; see also the earlier work of Andrew.

There is quite a bit of recent literature on the circular law, see for instance this survey by Van Vu and myself, or my lecture notes on this topic.

It looks like Joseph O'Rourke has beaten me, but I've generated a bit more data than he has:

Let $l$ be the lower bound so that the non-diagonal entries are uniformly distributed between $l$ and 1 and the diagonal entries are distributed between $l+1$ and 2. Note that Piero's original characterization of his code was inaccurate, but I just used his code to generate the figures below.

Here is a grid of eigenvalues of such $N$ by $N$ matrices with $l=-.9$ (plots are labeled by $N$):

Here's the picture with $l=-.87$:

Here's the picture with $l=-.93$:

And here's a picture of the absolute value of the largest eigenvalue of the matrix as a function of $N$ (now the label is $l$):

Edit based on Helge's comment: Here's a picture of $\max_{1\leq j \leq N} \sum_{n=1}^{N} u_j^N(n)$ as a function of $N$ (the $u_j^N$ are normalized eigenvectors) with $l=-0.9$:

And here's a picture of the further normalized version $\max_{1\leq j \leq N} \frac{1}{\sqrt{N}}\sum_{n=1}^{N} u_j^N(n)$:

Edit 2: In response to Helge's comment asking for a plot of $\max_{1\leq j \leq N} \sum_{n=1}^{N} u_j^N(n)$ when $l=-1$:

Here is the output from Piero's posted code: