the density function is where the magic is... that and turbulence
the povray manual states that the built-in 'spherical' function is 1.0-min(1, sqrt((x*x) + (y*y) + (z*z)))
i love how a simple blind random change can create such beauty