# Transformacao Integral da Probabilidade set.seed(24/11/2021) n <- 1e3 theta <- 2 u <- runif(n,0,1) y <- theta*u^(1/n)
# a ser feito