earthquake <- function(lambda=1, n_sim =10, n=100){ meanls <- c() for (i in 1:n){ meanls <- c(meanls,round(mean(rexp(n_sim,1/lambda)),2)) } retu