sim <- function(p) {
   # n sind Nächte die Glühbirne überlebt
   n <- 1

   x <- rbinom(1, 1, p)
   while( x == 0 ) {
      n <- n + 1
      x <- rbinom(1, 1, p)
   }

   return(n)
}

sim2 <- function(p) {
   x <- runif(1)
   n <- ceiling(log(1-x) / log(1-p))
   return(n)
}