  A quantum statistical model of nuclear multifragmentation is proposed. The
recurrence equation method used within the canonical ensemble makes the model
solvable and transparent to physical assumptions and allows to get results
without involving the Monte Carlo technique. The model exhibits the first order
phase transition. Quantum statistics effects are clearly seen on the
microscopic level of occupation numbers but are almost washed out for global
thermodynamic variables and the averaged observables studied. In the latter
case, the recurrence relations for multiplicity distributions of both
intermediate-mass and all fragments are derived and the specific changes in the
shape of multiplicity distributions in the narrow region of the transition
temperature is stressed. The temperature domain favorable to search for the HBT
effect is noted.
