function y=f(x) // simulation d'une fonction
    y=1/sqrt(2*%pi)*exp({-x^{2}/2})
endfunction

M=grand(10000,200,"exp",1); // on veut illustrer le tcl pour la loi exponentielle: on considère 10000 échantillons de longueur 200 
u=[-5:0.1:5]; // on part de -5 et on va jusqu'à 5 de 0,1 en 0,1
a=zeros(1,10000);
for i=1:10000
a(i)=((sum(M(i,:))/200)-1)*sqrt(200); // on applique le tcl pour chaque échantillon
end
histplot(100,a)
plot2d(u,f(u)) 
