set.seed(424)
a = 3;  b = 8;  v = runif(10^6, a,b);  cutp=seq(a, b, by=.5)
par(mfrow=c(1,2)) # enables 2 panels per plot
 hist(v, prob=T, br=cutp, col=rainbow(11))
 r = (3*v/(4*pi))^(1/3); cutpr = (3*cutp/(4*pi))^{1/3}
 hist(r, prob=T, br=cutpr, col=rainbow(11))
  curve(4*pi*x^2/(b-a), add=T, lwd=2)
par(mfrow=c(1,1))