n<-1000
U<-runif(n,min=0,max=1)
X<- -log(1-U)
Y<-exp(-X)
Z<-X*X*X*X*X*X*X
hist(Y, breaks = 50, col = "skyblue", freq = FALSE,
     main = "Distribución de Y = exp(-X), X~Exp(1)",
     xlab = "Y")

hist(Z, breaks = 3, col = "skyblue", freq = FALSE,
     main = "Distribución de Z = X^7, X~Exp(1)",
     xlab = "Z")