Circular law states the eigenvalues of a matrix with random entries of mean 0 and variance 1/n are approximately uniformly distributed in the unit disk of the complex plane. To see this, create a random matrix:
n = 1000
x = randomu(seed, n, n) - 0.5
x *= sqrt(12.0 / n)
Find the eigenvalues:
eigenvalues = la_eigenproblem(x, eigenvectors=eigenvectors)
Plot them:
plot, real_part(eigenvalues), imaginary(eigenvalues), $
psym=3, $
xstyle=1, xrange=[-1.5, 1.5], $
ystyle=1, yrange=[-1.5, 1.5]
This gives a plot like below:
Via the excellent John D. Cook blog. I recommend reading his site if you are interested in a combination of mathematics and Python.