R=1/2 +X/6; X/6=R-1/2

X=6(R-1/2) — Equation 1 for 1/2<R<0

R=1/2 +X^2/32; X^2=32*(R-1/2)

X=4(SQRT(2R-1)) — Equation 2 for 1/2>R>1

set.seed(4000)
r<-runif(10000,0.5,1)
x<-ifelse(r<=0.5,6*(r-1/2),4*sqrt(2*r-1))
 hist(x)