We propose extension of the algorithm for numerical modelling of Bose-Einst
ein correlations (BEC), which was presented some time ago in the literature
. It is formulated on quantum statistical level for a single event and uses
the fact that identical particles subjected to Bose statistics do bunch th
emselves, in a maximal possible way, in the same cells in phase-space. The
bunching effect is in our case obtained in novel way allowing for broad app
lications and fast numerical calculations. First comparison with e(+)e(-) a
nnihilations data performed by using simple cascade hadronization model is
very encouraging. (C) 2001 Elsevier Science B.V. All rights reserved.