This is an implementation of the model described in
Vicsek T and AS Szalay. 1987. Fractal distribution of galaxies modeled by a cellular-automaton-type stochastic process. Phys. Rev. Lett. 58, 2818-21.