tic;
nn=100;
for i=1:1000
 a=rand(nn); b=rand(nn);
 c=a*b;
endfor
toc
pause
tic;
n=1000;
i=1:n+1;
dx=2.0/n;
x=(i-1)*dx;
m=10000; p=3;
eve=zeros(1,n+1);
for k=1:m
 d=abs(det(rand(p)));
 j=floor(d/dx);
 if (j<=n)
 eve(j+1)=eve(j+1)+1;
 endif
endfor
toc
figure 1
plot(x,eve,"linewidth",2)
title("Eventos entre x e x + dx");
xlabel("x"); ylabel("#eve");
P(1,1)=0;
for k=2:n+1
 P(1,k)=P(1,k-1)+eve(k);
endfor
P = P/m;
figure 2
plot(x,P,"linewidth",2)
title("Probabilidade de |det(A)|<x");
xlabel("x"); ylabel("P");