# funcao_gva_exp
# Autor
# 13/03/2020
# Função para gerar valores com distribuição exponencial 
#  com parâmetro lambda
rm(list = ls())

# Semente
set.seed(23)

#função
expo <- function(lambda,n){
  y <- c(rep(0,n))
  for(i in 1:n){
    u <- runif(1)
    y[i] <- (-log(1-u))/lambda
  }
  return(y)
}

# 100 Valores gerados de Exp(4)
x <- expo(4,100)

# Para verificar se a amostra gerada segue a distribuição 
# Função de distribuição empírica
plot(ecdf(x))
# Função distribuição exponencial 
curve(pexp(x, 4), add = TRUE, type = 'l', col = 'red')

# ou
# Histograma
hist(x, freq = FALSE)
# Função densidade exponencial
curve(dexp(x, 4), add = TRUE, type = 'l', col = 'red')