num.rep<-100000
tra<-sqrt(2)-sqrt(2-2*runif(num.rep,0,1))
plot(density(tra),xlim=c(0,2),ylim=c(0,2),type="l")