clf; hold on;
title('Simulation d''une loi normale');
ylabel('Effectifs'); xlabel('Valeurs');
[E,C]=histo(sqrt(-2*log(rand(5000,1))).*cos(2*pi*rand(5000,1)),100,0,1);
plot(C,(2*pi)^(-1/2)*exp(-C.^2/2),'k-*');
legend('Empirique','Theorique');