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")