(commands used to generate figure above)

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

>