 % Illustration du TLC avec un échantillon de lois uniformes
 % Fait appel à StixBox pour pnorm
 clear; 
 i=0;    
 p=input('Nombre de v.a. uniformes (par exemple 10) = ? ');
 n=input('Taille de l''echantillon (par exemple 50) = ? ');
 d=input('Discrétisation (par exemple 100) = ? ');
 % On simule n fois la variable Z:=(U_1+...+U_p-p/2)/(p/12)
 X=rand(p,n);
 Z=sum(X)-p/2;     % Centrage
 Z=Z/sqrt(p/12);   % Réduction
 dz=range(Z)/d;
 % Fonction de Répartition Simulée
 [Effectifs,Classes]=hist(Z,d);
 FRS=cumsum(Effectifs)/n;
 FRS(d+1)=1;
 % Tracé de la fonction de répartition simulée
 % et de celle de la loi normale
 Abscisses=[min(Z)+dz/2:dz:max(Z)+dz/2]; 
 figure; clf; hold on;
 plot(Abscisses,FRS,'b+-');
 FRN=pnorm(Abscisses);
 plot(Abscisses,FRN,'r-');
 legend('Fonc. de rep. simulee','Fonc. de rep. Normale',4);
