Burke's theorem says that the output process of an $\rm M/M/1$, an $\rm M/M/C$, and a $\rm M/M/\infty$ queue with arrival rate $\lambda$ and service rate $\mu$ follows a Poisson with parameter $\lambda$.
I was able to derive the proof for the $\rm M/M/1$ easily using Laplace transforms:
For the busy period, the time spent in between departures is given by an exponential distribution with parameter $\mu$. For the non-busy period the inter-departure is given by a sum of an inter-arrival and one service time.
\begin{align}LD(s)&=\frac{\lambda}{\mu}\cdot\frac{\mu}{\mu+s}+(1−\frac{\lambda}{\mu})\cdot\frac{\lambda}{\lambda+s}\cdot\frac{ \mu}{\mu+s}=\frac{\lambda}{\lambda+s}\\D(t)&=\lambda e^{-\lambda t}\end{align}
However, I am struggling immensely to derive the proof for the $\rm M/M/\infty$ case. Can someone point out where I am going wrong?
So far I have for the busy period: $\sum\limits_{i=1}^\infty e^{-\lambda/μ}\cdot \frac{i\mu}{ i\mu+s}$. I am not even sure if this is correct. Do we need to account for the possibility of new arrivals before the departure too? How one would go with that?
Also, how would I set it up for the non-busy period? And how do I deal with this $\frac{i\mu}{ i\mu+s}$ term in the busy one?