Equatorial density cuts for τ₀ = 0.018 (initial conditions, highest and narrowest peak), τ₀ = 0.04 and τ₀ = 0.09 (widest distribution). In each case the diamond’s curve represents the distribution inferred in the Z = Z₀ plane from the numerical simulations, using a bilinear interpolation, while the solid curve represents the exact solution of Lynden-Bell & Pringle (1974). In addition, the dotted and dashed curves represent the situation obtained for τ₀ = 0.09 in numerical calculations in which we set and to zero, respectively.