#---------------------------------------------------------
# Simula esponenziale mediante inversione della cumulativa
#
# GdA, PED 25/11/2010
#---------------------------------------------------------
# function pausa
pausa <- function() {
  cat ("\n >> Guarda il plot e dai enter per continuare\n")
  scan()
}

tau = 2.5
n = 100000
t <- -tau * log(runif(n))
hist(t, nc=100)

media = mean(t)
stdev = sd(t)
cat(sprintf("media = %.2f; stdev = %.2f\n", media, stdev))

pausa()
# modo alternativo,usando una funzione (misteriosa...) di R:
t <- rgamma(n, 1, 1/2.5)
hist(t, nc=100)
media = mean(t)
stdev = sd(t)
cat(sprintf("media = %.2f; stdev = %.2f\n", media, stdev))