# Lecture 3A demonstration: Sum of exponential waiting times
# (This is R code)

rate    <- 1/100
repeats <- 1000000
n       <- 1        # Try changing this value
X       <- matrix(rexp(repeats*n, rate), repeats, n);
S       <- rowSums(X)
hist(S,100)