 % Illustration du TLC avec un échantillon de lois uniformes
 % Fait appel à StixBox pour dnorm & histo
 clf     % Nettoyage des graphiques
 clear   % Nettoyage des variables
 i=0;    
 p=input('Nombre de v.a. uniformes (par exple 12) = ? ');
 n=input('Taille de l''echantillon (par exple 1000) = ? ');
 % On simule n fois la variable Z:=U_1+...+U_p, 
 X=rand(p,n);
 Z=sum(X)-p/2; % Centrage
 Z=Z/sqrt(p/12); % Réduction
 % Z est proche de la N(0,1) d'après le TCL
 % Nombre de classes pour histogramme
 nc=sqrt(n); 
 % On trace l'histogramme normalisé en nc classes
 % E contient les effectifs pour nos classes
 % C contient les extrémités des classes
 [E,C]=histo(Z,nc,0,1);
 hold on;
 % On trace la densité normale
 plot(C,dnorm(C),'r-');
 legend('Histogramme','Loi normale');