%%% Test de Kolmogorov-Smirnov pour générateurs gaussiens
clear; clf
%
n=1000; r=1000; X=zeros(1,r);
%
h = waitbar(0,'Patience...');
for i=1:r,
 U=rand(1,n); V=rand(1,n);
 G=sqrt(-2*log(U)).*cos(2*pi*V);
 %G=randn(1,n);
 GG=sort(G); 
 % Ici, les GG(i) sont tous differents presque surement
 % car la loi gaussienne est diffuse. C'est par ailleurs nécessaire 
 % à l'application du test de K-S. 
 X(i)=sqrt(n)*max(abs([1:n]/n - pnorm(GG)));
 waitbar(i/r)
end
close(h)
X=sort(X);
hold on
plot(X,pks(X),'g-')
stairs(X,[1:r]/r)
legend('PKS','Femp')
