clear;
clf;
N=input('Nombre de realisations:');
n=input('Premier parametre de la binomiale:');
p=input('Deuxieme parametre de la binomiale:');
X=rbinom(N,n,p);
plot(cumsum(X)'./[1:length(X)],'b')
title('Loi des grands nombres')
xlabel('Nombre de realisations')
ylabel('Moyennes empiriques')
hold on
plot(n*p*ones(length(X),1),'r--')
legend('Empirique','Theorique')
hold off