Normal approximation in an urn model with indistinguishable balls

For the Bose-Einstein Statistics, where n indistinguishable balls are distributed in m urns such that all the arrangements are equally likely, define the random variables

M_k = number of urns containing exactly k balls each;

N_k = number of urns containing at least k balls each.

We consider the approximation of the distributions of M_k and N_k by suitable normal distributions, for large but finite m. Estimates are found for the error in the approximation to both the probability mass function and the distribution function in each case. These results apply also to the alternative model where no urn is allowed to be empty. The results are illustrated by some numerical examples.