N=1000000
M=0
lambda=2
tirages<-rpois(N,lambda);
for (i in 1:N) {
  if (tirages[i]%%2==0) {
    M=M+1;
  }
}
h<-hist(tirages,breaks=200)
cat('proba de trouver un nombre pair est de',M/N);