#random numbers set.seed(42) r <- runif(10000, 0.0, 1.0) x<-ifelse(r<=0.5,6*r-3,4*sqrt(2*r-1))
## Warning in sqrt(2 * r - 1): NaNs produced
hist(x)
#QED :)