clear; n=5000; a=1; clf; subplot(1,2,1); hold on; axis equal; 
S=sign(randn(n,n)); X=(rand(n,n).^(-1/a)-1).*S; E=eig(X/n^(1/a)); M=abs(E);
plot(real(E),imag(E),'ro'); title('Spectrum'); subplot(1,2,2); 
hist(M,fix(sqrt(n)),1); title('Spectrum module'); print('simu.png','-dpng')
