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