%%
%% TP1 Exercice 3
%% Simulation de la loi Binomiale
%% a partir de la loi de Bernoulli
%%

N=input('Entrez le nombre de réalisations N : ');
n=input('Précisez la valeur du paramètre n : ');
p=input('Précisez la valeur du paramètre p : ');
X=[];
for i=1:N
    U=rand(n,1);
    P=p*ones(n,1);
    V=(U<=P);
    X=[X;sum(V)];
end;

plot(cumsum(X)'./[1:N])
title('Loi des grands nombres')
xlabel('Nombre de realisations')
ylabel('Moyennes empiriques successives')
hold on 
plot(n*p*ones(N,1),'r--')
legend('Empirique','Theorique')
hold off

