lambda=0.5; clf; hold on;
title('Simulation d''une loi exponentielle')
ylabel('Densite'); xlabel('Valeurs');
[E,C]=histo(-log(rand(500,1))/lambda,100,0,1);
plot(C,lambda*exp(-lambda*C),'k-*')
legend('Empirique','Theorique')
