clear();
N=1000;
Esp=0;
Var=0;
alpha=3;
valeur=0;
Variance=0;
r=0;
for i=1:N
cumul=0;
ajout = 0;
valeur=0;
r=rand();
J=1;
cumul=exp(-1*alpha);
ajout=cumul;
while cumul<r do
valeur=valeur+1;
ajout=ajout*alpha/valeur;
cumul=cumul+ajout; 
end;
Esp=Esp+valeur/N;
Var=Var+ valeur*valeur /N;
end;
Variance = Var - (Esp*Esp);
disp(Esp);
disp(Variance);
