  We propose new algorithm for numerical modelling of Bose-Einstein
correlations (BEC). It is based on the fact that identical particles subjected
to BEC do, by definition, bunch themselves in a maximal possible way,
restricted only by conservation laws, in the same cells in phase-space.
