function part1gsim2(N,L,numavg)
  
  alleigs = [];

  for (l=1:numavg) 
    l
    G = randn(L,N);
    alleigs = [ alleigs ; eig((1/N)*G*G') ];
  end
  
  [ n , xout ] = hist(alleigs,200);
   n = n / ( numavg*L*(xout(2) - xout(1)) );
   bar(xout,n)
  axis([0 10 0 1])