(commands used to generate figure above)

>	chargerand:= rand(1..100)/100: for i from 1 to 10000 by 1 do    ppp[i]:=[1-chargerand()^2,1-chargerand()^2,2*chargerand()]: end do: pointplot3d({seq(ppp[i],i=1..10000)},symbol=point,scaling=constrained,axes=boxed,color=blue):

>