R functions to simulate the Petersburg Paradox.

petersburg.r

pburg <- function() {
  p <- 1
  while (sample(1:2, 1)[1] == 1) {
    p <- p * 2
  }
  return(p)
}

ipburg <- function(n) {
  r <- sapply(1:n, function(x) { return(pburg()) })
  return(r)
}

cumpburg <- function(n, v=1) {
  f <- seq(1,n,v)
  r <- cumsum(ipburg(n))[f]/f
  return(r)
}