[R-br] ajuda para um scrip

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