  We consider a stationary Markov process that models certain queues with a
bulk service of a fixed number $m$ of admitted customers. We find an integral
expression of its transition probability function in terms of certain
multi-orthogonal polynomials with respect to a system of distributions that
contain measures supported on starlike subsets of the complex plane. We also
study the corresponding steady state.
