[R-br] ajuda para um scrip

[Mauro Sznelwar]

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
 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110823/72444718/attachment.html>