  The equilibrium distributions of the different pasta geometries and their
linear sizes are calculated from the mean field Gibbs energy functional in
symmetric nuclear matter at finite temperature. The average sizes and shapes
coincide approximately with the ones predicted by a standard pasta calculation
in the coexisting phase approximation, but fluctuations are additionally
calculated and seen to increase with temperature and baryonic density. The
different pasta shapes are shown to coexist in a wide domain of density and
temperature, in qualitative agreement with the findings of large scale
molecular dynamics simulations, but with a much less expensive computational
cost.
