A closed queueing system has Q service channels and a waiting line. There are Ni customers of type i in the system, i = 1, ···, m, ∑mi=1Ni = N > Q. Q customers are served and Q0 = N – Q stay in the waiting line. Q channels are partitioned into n groups with Qj channels in the jth group, j = 1, ···, n. The service time of the ith type customer by a channel of the jth group is τij ~ Exp (λij). When a customer leaves the channel, it is immediately replaced by another one picked up randomly from the waiting time. The customer which has cleared service joins the waiting line without delay. Let Xij be the number of ith type customers served by jth group channels in equilibrium. An explicit formula for P(Xij = kij, i = 1, ···, m; j = 1, ···, n) is found. It is shown in a form of a local limit theorem that the asymptotic distribution of {Xij} is a multidimensional normal, if Qj/N and Qj/N have positive limits as N → ∞. Formulas for mean values and covariances are given. It turns out that the means of Xij and covariances between Xij and Xrk can be found, using an efficient iterative algorithm, from the deterministic version of the system. A numerical example demonstrates that the normal approximation is rather accurate.