  The probability distribution of density in isothermal, supersonic, turbulent
gas is approximately lognormal. This behaviour can be traced back to the shock
waves travelling through the medium, which randomly adjust the density by a
random factor of the local sonic Mach number squared. Provided a certain parcel
of gas experiences a large number of shocks, due to the central limit theorem,
the resulting distribution for density is lognormal. We explore a model in
which parcels of gas undergo finite number of shocks before relaxing to the
ambient density, causing the distribution for density to deviate from a
lognormal. We confront this model with numerical simulations with various
r.m.s. Mach numbers ranging from subsonic as low as 0.1 to supersonic at 25. We
find that the fits to the finite formula are an order of magnitude better than
a lognormal. The model naturally extends even to subsonic flows, where no
shocks exist.
