# UMA DEMONSTRAÇÃO DAS ALTERNATIVAS DA FÓRMULA DE STURGES...
set.seed(1908)
histo.exp <- function(n){
  # Esta funcao constroi 4 histrogramas para amostras simuladas da 
  #    distribuicao Exponencial.
  x <- rexp(n,rate=10)
  op <- par(mfrow=c(2,2))
  
  hist(x, main="Fórmula de Sturges")
  hist(x,breaks="Scott", main="Fórmula de Scott")
  hist(x,breaks="FD", main="Fórmula de Freedman-Diaconis")
  hist(x,breaks=sqrt(n), main="Raiz Quadrada de n")
  
  par(op)
  cat("Estatisticas Descritivas: \n")
  print(summary(x))
  cat("Distancia Interquartilica =",IQR(x),"\n")
  cat("Desvio Padrao =",sd(x), "\n\n\n")
  
  rm(x)
}

# Histogramas com tamanhos de amostras crescentes
histo.exp(30)
histo.exp(300)
histo.exp(3000)
histo.exp(30000)
histo.exp(3000000)

#Garbage Collection...
gc(TRUE)


# FIM DE PROGRAMA