%%% Test de Kolmogorov-Smirnov pour générateurs gaussiens
clear
%
n=1000; alpha=0.05;
%
U=rand(1,n); V=rand(1,n);
G=sqrt(-2*log(U)).*cos(2*pi*V);
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=sqrt(n)*max(abs([1:n]/n - pnorm(GG)));
p=pks(x);
if p<1.-alpha,
   sprintf('On accepte H0. Niveau = %4.2f, n=%d, KS=%f.',alpha,n,p)
else	
   sprintf('On rejette H0. Niveau = %4.2f, n=%d, KS=%f.',alpha,n,p)
end
