N <- 10000
x <- seq(0,1,len=N+1)
X <- rbeta(x,1,2)
bornes <- seq(0,1,0.05)
hist(X,breaks=bornes,col='lightcyan')
points(x,N*0.05*dbeta(x,1,2),type="l",col="red")