Solution (source code)

= Solution

Assume $\rho=\lambda/[\mu(1-p)]<1$. Detailed balance gives
$$
\pi_{n+1}\mu(1-p)=\pi_n\lambda,
$$
so normalization yields
$$
\boxed{\pi_n=(1-\rho)\rho^n,\qquad n\geq0.}
$$
This is a geometric distribution on $\mathbb Z_+$, and therefore
$$
\boxed{\mathbb E_\pi L=\frac{\rho}{1-\rho}
=\frac{\lambda}{\mu(1-p)-\lambda}.}
$$