N <- 100000 U <- runif(N,0,1) X <- 1-sqrt(U) bornes <- seq(0,1,0.05) hist(X,breaks=bornes)