Tenho este algoritmo

p=1, N=0, c=e^-lambda

Repita

N=N+1

gere Ui

atualize p<c.e^-lambda

Faça X=N-1

Estou fazendo esta rotina abaixo, mas o resultado não está consistente! Qual o melhor jeito?

p=1
n=0
lambda=2
c=exp(-lambda)
for(n in 1:100){
u=runif(n)
p=p*u
p<c
}
x=n-1