  We demonstrate algorithm for numerical modelling of Bose-Einstein
correlations (BEC) formulated on quantum statistical level for a single event
and exploring the property that identical particles subjected to Bose
statistics do bunch themselves in a maximal possible way in the same cells in
phase-space.
