# lambda es el parametro de la Poisson que se desea generar.

lambda<-5

x<-vector(mode="numeric",length=100)

#
# Aqui empieza el ciclo para generar las variables.
# Mientras la suma de los tiempos exponenciales generados
# sea menor o igual que lambda, sigue generando tiempos.
#
# La expresion que genera tiempo ~ Exp(1) podria reemplazarse 
# por una expresion que calculase la exponencial a partir
# de la Uniforme(0,1) (runif).
#

for (i in 1:100)
{
 suma<-0
   k<-0
   while (suma<=lambda)
   {
      tiempo<-rexp(1,1)
      k<-k+1 
      suma<-suma+tiempo
   }
 x[i]<-(k-1)
}