function  semi
n=2;     %size of matrix
s=30000;   %number of samples
d=.05;     %sample distance
e=[];     %gather up eigenvalues
im=1;     %imaginary(1) or real(0)

for i=1:s,
  a=randn(n)+im*sqrt(-1)*randn(n);a=(a+a')/(2*sqrt(2*n*(im+1)));  %The GE
  v=eig(a)';
  e=[v e];
end
hold off; [m x]=hist(e,-1.5:d:1.5); bar(x,m*pi/(2*d*n*s));
axis('square'); axis([-1.5 1.5 -1 2]); hold on;
t=-1:.01:1; plot(t,sqrt(1-t.^2),'r');
levels(n,im);
im=0;
