//[]
clf;
N=10000;
n=3000;
pas_histogramme=1000/N;
x_max=4;
x_min=-x_max;

U=[];

for k=1:N;

S=sqrt(12*n)*(sum(rand(1,n))/n-0.5);

U=[U,S];

end

function y=densite(x); y=(exp(-x^2/2))/sqrt(2*%pi); endfunction

xdens=x_min:0.01:x_max;
plot2d(xdens,densite(xdens),5)

rect_histo=x_min:pas_histogramme:x_max;
histplot(rect_histo,U)

xtitle(justify(["VISUALISATION DU TCL :"; 
"histogramme pour un echantillon de "+string(N);
 "copies independantes de sqrt(n)(S_n-m)/s"; 
"pour n = "+string(n)+" et des v.a. de loi uniforme sur [0,1]"],'c'))

legend("y=(exp(-x^2/2))/sqrt(2*%pi)", "histogramme")
xselect()