hill <- function(x) {
  x <- log(rev(sort(x[x > 0])))
  n <- length(x)
  cumsum(x)[ - n]/(1:(n - 1)) - x[-1]
}

sim <- function(n) {
  i <- 1
  while (i <= 10000) {
    a <- hill (1.0 / runif (n))
    i <- i + 1
  }
}