nadja in wbo

This article is a program of calculating Wigner's semi-cyclic rule.

J is symmetry matrix (dimensionality N) with on diagonal 0. Therefore off diagonals of the symmetry matrix are distributed on N(0,1/N);

The followings are programs on Wigner semicircle rule utilizing Matlab.

%%%%%%%%%%%%%%%%%%%%%%%

N=1000;
FidJ=fopen('JJ.dat','wt');
A=randn(N);
A=A-diag(diag(A));
normal_J=1.0/sqrt(2.0*N);
J=(A+A')*normal_J;
JJ=eig(J);
fprintf(FidJ,'%f\n',JJ);
fclose(FidJ);
exit

%%%%%%%%%%%%%%%%%%%%%%%

or we describe histogram.

%%%%%%%%%%%%%%%%%%%%%%%

N=1000;
A=randn(N);
A=A-diag(diag(A));
normal_J=1.0/sqrt(2.0*N);
J=(A+A')*normal_J;
JJ=eig(J);
hist(JJ,100);
pause
exit

%%%%%%%%%%%%%%%%%%%%%%%

Yesterday, the first time is that I have made the program by Matlab.