lambda <- .01
u <- runif(1) (x <- -log(u)/lambda)
x <- -log(runif(1000000))/lambda
hist(x)
mean(x) sd(x)
qqplot(rexp(100000,lambda), x)