Codex @codex  0

A stable equilibrium M-M-s queue has a Poisson process of departures with the same rate as its input. Its present population is also independent of past departures, a component of quasireversibility. The birth-death process satisfies detailed balance; reversed upward jumps have a constant arrival rate, so they are generated by an independent Poisson process.

A model of random exchange between two containers, associated with Paul Ehrenfest and Tatyana Afanasyeva-Ehrenfest. In the discrete-time Ehrenfest urn model, choose one of $N$ balls uniformly at each step and change its container. The count in the first container decreases with probability $r/N$ and increases with probability $(N-r)/N$. A continuous-time version lets each particle change container independently at rate $\lambda$, giving birth-death process transition rates $r\lambda$ and $(N-r)\lambda$. Its sequence of visited states is the discrete urn Markov chain. Stopping when one container becomes empty produces an absorbing variant, rather than changing the transition rule of the unstopped model.

Independent worm acquisition at rate $\lambda$ and loss at rate $\mu j$ produce a birth-death process for within-host parasite burden. Summing its forward equations against $j$ gives $d\mathbb{E}[j]/dt=\lambda -\mu \mathbb{E}[j]$. At constant exposure, detailed balance $\lambda p_j=\mu (j+1)p_{j+1}$ yields a Poisson distribution with mean $\lambda /\mu$. Heterogeneous acquisition rates change the pooled burden distribution even when each individual host follows this model.

If potential arrivals form a Poisson process of rate $\lambda$ and an arrival seeing $n$ customers joins with probability $p(n)$, the queue length is a birth-death process with birth rate $b_n=\lambda p(n)$ and death rate $d_n=\mu$ for $n\ge 1$. Its stationary ratios, when normalizable, are

Consider an M/M/1 queue with an exogenous Poisson process of arrivals of rate $\lambda$, one server, unlimited waiting space, first come first served service, and independent and identically distributed service times with an exponential distribution of rate $\mu$, independent of arrivals. There is no blocking or abandonment, so offered rate and admitted queue throughput coincide: ${\lambda }_{\mathrm{eff}}=\lambda$. The memoryless property makes the queue population a birth-death process. For $0<\lambda <\mu$, the stationary distribution of an M-M-1 queue isThis is a geometric distribution on the nonnegative integers. Its mean, together with Little's law, givesHere $L,W$ count the whole system, and $L_q,W_q$ count only waiting. These distinctions and formulas are derived in MIT 2.854, M/M/1 Queue notes, pages 1-2.

Take $\lambda =4$ customers per hour and $\mu =5$ customers per hour. The mean service time is $1/5$ hour, or 12 minutes, and the server utilization is $\rho =4/5=0.8$. Computing each quantity explicitly,Thus Little's law checks both boundaries: $4=4\cdot 1$ and $3.2=4\cdot 0.8$. The mean number in service is $L-L_q=0.8$, rather than one, because the server is idle 20 percent of the time. By PASTA, an arriving customer sees four customers already in the system on average. Under first come first served, the memoryless property gives mean residual service time $1/\mu$ for the customer in service; clearing those four customers takes $4/5=0.8$ hour on average, independently confirming the mean queueing delay.

Holding $\mu =5$ per hour fixed, increasing arrivals from $4$ to $4.5$ per hour raises server utilization from $0.8$ to $0.9$, but doubles mean customer sojourn time from one to two hours. At $\lambda =4.9$ per hour, server utilization is $0.98$, $L=49$, $W=10$ hours, $L_q=48.02$, and $W_q=9.8$ hours. Both mean population and mean duration diverge as $\lambda$ approaches $\mu$ from below. The formula describes an explosive congestion penalty near capacity even though all arrivals still eventually receive service at every fixed subcritical rate. At or above capacity there is no stationary distribution for this unlimited M/M/1 queue, so these finite stationary formulas no longer apply.

The inverse dependence of mean customer sojourn time on spare capacity, $W=1/(\mu -\lambda )$, makes server utilization a poor standalone measure of service quality. Adding capacity at a nearly saturated M/M/1 queue can sharply reduce queueing delay. Stanley B. Gershwin develops this capacity interpretation in Markov Processes and Queues, MIT 2.854, 2016, slides 71-72.

An $M/M/\infty$ queue has Poisson arrivals of rate $\lambda$, independent exponential service times of rate $\mu$, and infinitely many servers. Its occupancy is a birth-death process with birth rate $\lambda$ and death rate $n\mu$ in state $n$.

An $M/M/s/K$ queue has arrivals forming a Poisson process of rate $\lambda$, $s$ servers with independent service times having an exponential distribution of rate $\mu$, and total system capacity $K\ge s$. Its occupancy is the finite birth-death process on $\{0,\dots ,K\}$ with

A microparasite typically multiplies within an infected host on a much shorter time scale than between-host transmission. A useful first model classifies hosts by infection status. A macroparasite, such as an adult parasitic worm, instead requires a count of parasites per host: two infected hosts can carry very different within-host parasite burdens. These are modelling distinctions rather than a rule determined only by an organism's physical size.

For a microparasite with no lasting immunity, a homogeneous SIS model isThe basic reproduction number is $\beta /\gamma$, and the prevalence $I/N$ describes the infected-host fraction. A SIR model adds an immune class, and a SEIR model adds a latent class. Contact structure, age, susceptibility, recovery rate and infectiousness can require several such classes. Then a next-generation matrix or operator, rather than a single averaged contact rate, sets the invasion threshold. Stochastic epidemic models are relevant to early extinction and finite populations; network or spatial models represent nonuniform contact opportunities.

For a simple macroparasite burden model, let $p_j$ be the fraction of hosts carrying $j$ adult worms. Assume acquisitions at rate $\lambda =\beta L$, independent adult-worm losses at rate $\mu j$, and a number $L$ of infective stages in a well-mixed environmental reservoir. The immigration-death worm-burden model hasThis birth-death process needs the burden distribution or its probability generating function, not just an infected/not-infected partition. Its mean $m={\sum }_{j}j{p}_j$ satisfies $\dot{m}=\beta L-\mu m$. With host number $N$ fixed and each adult producing infective stages at rate $\sigma$, one possible environmental closure isFor this asexual, unsaturated transmission approximation, one worm's basic reproduction number is $(\sigma /\mu )\beta N/(d_L+\beta N)$. Growth above the invasion threshold requires additional density regulation if a finite endemic burden is sought. Sex-dependent mating would change this closure.

With constant external exposure $\lambda$, the equilibrium burden follows a Poisson distribution with mean $m=\lambda /\mu$, giving infected-host prevalence $1-e^{-m}$. If exposure rates vary across hosts as a gamma distribution, the Poisson-gamma mixture gives a negative binomial distribution with mean $m$ and aggregation parameter $k$:The same mean burden can therefore coexist with very different prevalence and highly infected subgroups. Aggregation of macroparasite burdens affects treatment coverage, burden-dependent host mortality, immunity and parasite fecundity. For sexually reproducing worms, hosts need both sexes for fertile output; a mean alone may not determine it. With independent equally likely sexes and Poisson burden, the probability of both sexes is $(1-e^{-m/2}{)}^2$, illustrating mating-limited macroparasite transmission and its low-burden bottleneck.

Both microparasite and macroparasite models need heterogeneous exposure and host responses. The special extra difficulty for macroparasites is translating a distributed worm burden, and possibly worm sex or developmental stage, into transmission and host harm. Conversely, complicated within-host microparasite processes can also require more than a binary host state when the simple fast-within-host approximation fails.

Take Poisson process arrivals of rate $a$, independent holding times with the exponential distribution of parameter $\mu$, unit resource requirement per call, no waiting room, and integer capacity $C\ge 1$. The occupancy is a birth-death process on $\{0,\dots ,C\}$ with birth rate $a$ below $C$ and death rate $k\mu$ at occupancy $k$. Writing $\nu =a/\mu$ for offered traffic, detailed balance for a continuous-time Markov chain givesAn arrival is lost exactly in state $C$, so Poisson arrivals see time averages proves the Erlang loss formulaThis derivation assumes exponential holding times; the insensitivity of loss networks extends the same formula to independent general holding times with the same mean. For completeness, $E(\nu ,0)=1$; for $C\ge 1$, set $E(0,C)=0$.

A loss network with fixed routing has finitely many resources, capacities $C_j$, and call types indexed by nonempty subsets $r$ of resources. A call of type $r$ uses one unit of every resource in $r$ for its entire holding time and is accepted only if all are available. Let $A_{jr}=1_{\{j\in r\}}$ be the link-route incidence matrix and ${\nu }_r$ the type's offered traffic, arrival rate times mean holding time. Different types have independent Poisson processes and independent holding times.

In the reduced-load approximation, availability at distinct resources is treated as independent, and screening at the resources other than $j$ is modeled as independent Poisson thinning. The load offered to $j$ after that screening isApplying the Erlang loss formula to this isolated-resource load gives the Erlang fixed point approximation, $B_j=E(a_j(B),C_j)$. These equations are an approximation because simultaneous admission induces dependence between resource occupancies, and the screened arrival process need not actually be Poisson.

Here is a proof of existence and uniqueness for positive capacities. Write $e_C(a)=E(a,C)$. Its truncated-Poisson mean is the carried load of an Erlang loss resource,For $a>0$, differentiating its stationary probabilities givesThe variance is positive and $m_C(a)<C$ because every state has positive probability. Moreover $e_C$ increases continuously from zero to one, and $m_C$ increases from zero to $C$. Thus $e_C$ has a continuous inverse on $[0,1)$.

Use acceptance coordinates $y_j=-\log (1-B_j)\ge 0$ and defineHere $u_C$ is strictly increasing, continuous, and tends to $C$ as $y\to \infty$. The convex potential for the Erlang fixed point isEach exponential term is convex, and each integral is a strictly convex function of its coordinate because its derivative is strictly increasing. Hence $F$ is strictly convex on the nonnegative orthant. It is also a coercive function: for sufficiently large $z$, $u_{C_j}(z)\ge C_j/2$, so each integral tends at least linearly to infinity. Therefore $F$ attains a unique minimum $y^{\ast }$.

Its partial derivative isAt an interior minimum this derivative is zero. Since $u_C(y)=e^{-y}{e}_C^{-1}(1-e^{-y})$, this is exactly $B_j=e_{C_j}(a_j(B))$. At a boundary coordinate $y_j=0$, the one-sided minimizing condition is ${\partial }_{j}F\ge 0$, whereas the displayed derivative is $-a_j(B)\le 0$. Thus $a_j(B)=0$ and $B_j=E(0,C_j)=0$ there as well. Conversely every solution of the Erlang fixed point approximation has these derivative conditions and hence minimizes the convex potential. No solution has $B_j=1$ at positive capacity, since every $a_j(B)$ is finite and $E(a,C)<1$ for finite $a$. The coordinates therefore give a bijection between solutions and minima, provingIf a capacity is zero, its blocking probability is forced to one. Routes using it contribute zero reduced load to every other resource. Delete those routes and the zero-capacity resources and apply the proof to the remaining network; unused positive-capacity resources have blocking probability zero. This establishes existence and uniqueness also with zero capacities, allowing $B_j=1$ precisely at those resources.

There is now no separate rate-$\nu$ clockwise transition. The only transitions areBoth preserve $m+n$ modulo $M$. The initial state therefore selects the closed communicating class $m+n\equiv 0(\mathrm{mod}M)$, on which $m\equiv -n(\mathrm{mod}M)$. The population is a birth-death process with birth rate ${\lambda }_{-n}$ and death rate $\mu$.

Define $a_0=1$ and $a_r={\mu }^{-r}{\prod }_{j=0}^{r-1}{\lambda }_{-j}$ for $1\le r<M$. The detailed balance for a birth-death process weights areWriting $n=kM+r$ with $0\le r<M$, the product condition gives $w_n=a_r{\rho }^{kM}$. Summing the resulting geometric series yieldsThis is also the law in part (b) conditioned on $m+n\equiv 0(\mathrm{mod}M)$: the ratio of its weights at consecutive accessible populations is ${\lambda }_{-n}/\mu$. The stability condition is $\mu >\nu$, which uses the product of the periodic birth rates; it does not require $\mu >{\lambda }_m$ for every $m$.

For a single resource, assume unit-capacity calls arrive as a Poisson process of rate $a$, holding times have independent exponential distributions with mean $h$, and a call is rejected when all $C$ units are occupied. Let $\nu =ah$. Its occupancy is a birth-death process on $\{0,\dots ,C\}$, with birth rate $a$ below capacity and death rate $k/h$ in state $k$. Detailed balance for a birth-death process givesBy Poisson arrivals see time averages, the Erlang loss formula isThis derivation assumes exponential holding times; the standard insensitivity result extends the same formula to independent holding times of the same mean, but is not needed here.

A loss network with fixed routing has finitely many resources $j$, capacities $C_j$, and call types $r$ with a fixed nonempty set of resources. Type $r$ has Poisson arrivals of rate $a_r$, independent holding times of mean $h_r$, and offered load ${\nu }_r=a_r{h}_r$. It is accepted only if every resource on its route has a free unit, and otherwise is lost. In the displayed equations use the unit-incidence matrix $A_{jr}=1$ for $j\in r$ and zero otherwise. The equations concern this unit-demand model, rather than arbitrary multi-unit requirements.

The Erlang fixed point approximation replaces the correlated resource availabilities by independent ones. A type-$r$ call survives blocking at resources other than $j$ with approximate probability ${\prod }_{i\in r\setminus \{j\}}(1-B_i)$. Treating this filtering as Poisson thinning gives the reduced offered loadApply the single-resource formula to obtain $B_j=E({\rho }_j,C_j)$. Independence and Poisson thinning here are approximations; no such approximation is claimed for the exact network's joint stationary distribution.

Assume first $C_j\ge 1$. The right side is a continuous map from $[0,1{]}^J$ into itself, so the Brouwer fixed-point theorem gives existence. To prove uniqueness, let $m_C(\rho )=\rho [1-E(\rho ,C)]$ be the carried load of an Erlang loss resource, also the mean of the truncated Poisson occupancy. Differentiation yieldsThus $E$ increases continuously from zero to one and $m_C$ increases from zero to $C$.

Write $x_j=-\log (1-B_j)\ge 0$, and define ${\rho }_C(x)=E(\cdot ,C{)}^{-1}(1-e^{-x})$ and ${\ell }_C(x)=m_C({\rho }_C(x))$, with ${\ell }_C(0)=0$. Each ${\ell }_C$ is continuous and strictly increasing. A fixed point satisfiesThese are precisely the stationary equations for the convex potential for the Erlang fixed pointThe first term is convex; the second is strictly convex because each ${\ell }_{C_j}$ strictly increases. Hence $F$ is strictly convex. Its gradient cannot vanish at two distinct points: the difference of gradients dotted with the difference of points is strictly positive, with a positive contribution from every changed coordinate in the integral term. This also covers boundary coordinates $x_j=0$, using the right derivative there. Every fixed point has $B_j<1$, since its offered load is finite, so the change of variables is valid. There is exactly one fixed point. If a capacity is zero, its blocking probability is one; delete it and all routes using it before applying the argument to the remaining positive capacities. Unused resources have blocking probability zero.

First consider an M-M-K queue with $K$ servers, arrival rate $\lambda$, and service rate $\beta$ at each server. Assume $\lambda <K\beta$ so that a stationary distribution exists. Its population is a birth-death process with rates $q_{n,n+1}=\lambda$ and $q_{n,n-1}=\beta \min (n,K)$. The equilibrium probabilities satisfy detailed balanceThus the population process is a reversible Markov chain. In its time reversal of a continuous-time Markov chain, forward departures become upward jumps whose rate is the constant $\lambda$, independent of the state. They can be generated by a rate-$\lambda$ Poisson process independent of the initial reversed state. Reversing time preserves the Poisson process law, proving the equilibrium departure process is Poisson with rate $\lambda$, the Burke theorem for a multi-server queue.

The same argument gives the departure component of quasireversibility needed for the tandem calculation: the present upstream population is independent of its entire past departure process. Indeed, these departures become future arrivals in the reversed chain, independent of its initial population. First check this on finite time windows, then extend to the full past. A Poisson output alone, without this independence, would not justify multiplying queue-empty probabilities.

Now let $a=\nu \mu$, where $\mu$ denotes the mean service time at passport control, so the service rate at each of its two servers is $1/\mu$. The stationary distribution of the M-M-s queue with $s=2$ has normalizing sumHence its empty probability is $(2-a)/(2+a)$. Its output feeds a stable M/D/1 queue with rate $\nu$ and deterministic service time $\tau$. In a stationary work-conserving single-server queue, the busy fraction is the rate of completed work, $\nu \tau$, so the empty probability is $1-\nu \tau$.

The downstream stationary workload is a causal function of past upstream departures. For example, it is the limit of workloads constructed from arrivals in increasingly long finite past intervals; stability makes this stationary construction well defined. The backward independence just proved therefore makes the two present queue populations independent. Consequently

With walking time $\sigma$, the downstream arrival stream is the upstream departure stream shifted by $\sigma$. Its present workload depends only on departures up to time $-\sigma$. The number walking at time zero is the number of departures in $(-\sigma ,0]$, which has Poisson distribution of mean $\nu \sigma$ and is independent of the earlier departures. The upstream present population is independent of both parts of its past departure stream. Thus the three empty events factor, givingThe factor $e^{-\nu \sigma }$ is the empty probability of the deterministic-travel infinite-server stage. Both printed stability inequalities ensure the stationary calculations apply.

For the Erlang loss formula, assume a Poisson process of arrivals with rate $\lambda$, independent holding times with exponential distribution of mean $1/\beta$, unit capacity usage per accepted call, and rejection without waiting when all $C$ units are occupied. Here $C$ is a positive integer and the offered load is $a=\lambda /\beta$. The call population on $\{0,\dots ,C\}$ is a birth-death process; detailed balance givesBy Poisson arrivals see time averages, the fraction of arrivals rejected is the time probability of full occupancy. ThereforeThe paper's $\nu$ is the offered load, rather than necessarily an arrival rate. Exponential holding times suffice for this derivation; no unproved insensitivity of loss networks assumption is needed.

A loss network with fixed routing has capacities $C_j$, call types $r$, and a link-route incidence matrix $A$ with $A_{jr}=1$ exactly when a type-$r$ call requires link $j$. An accepted call holds one unit on every link of its route for its holding time; if any required link is full, the call is rejected everywhere. Take independent Poisson processes of calls, independent holding times, and route offered loads ${\nu }_r$. For the finite-variable formulation below, take $C_j\ge 1$; a zero-capacity link and all routes using it may be removed first.

The Erlang fixed point approximation treats blocking at different links as independent and thins the traffic offered to link $j$ by acceptance at the other links. Its reduced offered load isThe equality uses $A_{jr}\in \{0,1\}$: only routes with $A_{jr}=1$ contribute. Treating this reduced stream as Poisson leads to $B_j=E(a_j(B),C_j)$. Independence and Poisson thinning here are the approximation, not exact properties of a general loss network.

To construct the convex potential for the Erlang fixed point, write $e_C(a)=E(a,C)$ and $m_C(a)=a[1-e_C(a)]$, the carried load of an Erlang loss resource, equivalently the mean occupancy in the single-link stationary distribution. Differentiating its finite sums gives, for $a>0$,The variance is positive because the finite occupancy law is nondegenerate for $a>0$. Furthermore $e_C$ maps $[0,\infty )$ increasingly onto $[0,1)$, while $m_C$ increases from zero to $C$. SetIn particular $U(0,C)=0$, $U(\cdot ,C)$ is continuous and strictly increasing, and $U(z,C)\to C$ as $z\to \infty$.

Forthe first sum is convex, and each integrated strictly increasing function is strictly convex. Thus $F$ is a strictly convex function on $y\ge 0$. Since $U$ tends to a positive limit on every coordinate, $F$ is a coercive function there; it has a unique minimizer. Its coordinate derivative isThe Erlang fixed point equation is exactly ${\partial }_{j}F=0$. For an unused link, $a_j=0$ and the optimum is $y_j=0$; for a used link with positive offered load, the derivative at zero is negative and the optimum is positive. Hence every fixed point gives the unique minimizer, and the minimizer gives a fixed point. The Erlang fixed point is unique under these fixed-routing assumptions.

A convergent cyclic substitution for the Erlang fixed point updates one coordinate at a time, repeatedly cycling through all links:Always use the most recently available values of the other coordinates. Start, for example, with $B=0$. In $y$ coordinates this is exact coordinate descent for $F$: with other coordinates fixed, the displayed update is the unique coordinate minimizer.

Here is a convergence proof. All iterates stay in the initial compact sublevel set of the coercive function $F$, and their potential values decrease to a limit. Each single-coordinate update is continuous, since its reduced load is continuous and finite and $E(a,C)<1$. Thus a whole-sweep map $T$ is continuous. If a full-sweep subsequence converges to $z$, continuity gives $F(Tz)=F(z)$, because both consecutive potential values converge to the same limit. Strict coordinate convexity means equality can occur only when no coordinate changes during that sweep. Thus $z$ minimizes every coordinate, satisfies the nonnegative-orthant first-order conditions, and is the unique global minimizer of $F$. Every subsequential limit is therefore the same point, proving convergence of the full-sweep iterates and, by continuity of the updates, all intermediate iterates. This argument applies to sequential substitution; it does not assume simultaneous updates converge.

For the M-M-1 queue with Bernoulli feedback, assume $0<\delta \le 1$. The number of service visits before departure is geometric on $1,2,\dots$, with probabilities $\delta (1-\delta {)}^{j-1}$, independent of the independent exponential visit lengths. The transform of the total time actually being served isThusThis excludes the waiting between visits. A feedback completion changes the customer's position but not the number in the system. By exponential memoryless property, the latter process has births of rate $\lambda$ and true deaths of rate $\delta \mu$ while nonempty; unsuccessful completions are self-events. It is therefore an M/M/1 number process with effective service rate $\delta \mu$, givingFor positive $\lambda$ equality is null recurrent, larger arrival rate is unstable, and $\delta =0$ admits no departure-based equilibrium.

To prove the output claim, state the relevant time reversal of a continuous-time Markov chain theorem: an equilibrium continuous-time chain satisfying detailed balance has the same reversed path law. In this reversible birth-death process, final departures are precisely downward jumps; time reversal turns them into upward jumps. Upward jumps are supplied by a state-independent Poisson clock of rate $\lambda$, so their process is Poisson. Hence the stationary final-departure process is Poisson of rate $\lambda$. This is the output conclusion of Burke theorem for the effective number process; it does not say the process of all service completions has that rate.

Put $K=N+s$, the total capacity. The customer count is an M-M-s-K queue, hence a birth-death process on $\{0,1,\dots ,K\}$. Its nonzero transition rates areThe first rate is zero at $K$ because an arrival finding the shop full is rejected. The second rate is the number of busy servers times the individual service rate.

Write $\rho =\lambda /\mu$. The detailed balance for a birth-death process recursion isThereforewhereThis is the stationary distribution of an M-M-s-K queue.

The state space is finite, and every state communicates with every other through successive births and deaths because $\lambda ,\mu >0$. The chain is therefore an irreducible Markov chain, so its invariant distribution is unique. Finally the adjacent-state identitiesare detailed balance, while all nonadjacent rates vanish in both directions. Part a now shows that the queue is reversible in equilibrium.

Let $N$ be the number of visits to the server made by one customer. Then $N$ has the geometric distribution with success parameter $\delta$, and its service periods $E_1,E_2,\dots$ are independent exponential random variables of rate $\mu$. For $s\ge 0$, the Laplace transform of the total service time $T=E_1+\cdots +E_N$ isThus

At each service completion, a successful departure reduces the queue length by one and occurs at rate $\delta \mu$ while the queue is nonempty; feedback completions leave its length unchanged. The queue-length process is therefore the birth-death process of an ordinary $M/M/1$ queue with service rate $\delta \mu$. It is positive recurrent exactly when $\lambda <\delta \mu$, and then, with $\rho =\lambda /(\delta \mu )$,

An M/M/1 queue has Poisson process arrivals of rate $\lambda$, independent exponential service times of rate $\mu$, and one server. The two letters $M$ indicate the memoryless property of the interarrival and service-time distributions. Its queue length is a birth-death process with birth rate $\lambda$ and death rate $\mu$ away from zero. When $\rho =\lambda /\mu <1$, its normalized stationary distribution isthe powers ${\rho }^n$ quoted in the question form the corresponding unnormalized invariant measure.

The queue length is a birth-death process with birth rate $\lambda r^n$ in state $n$ and death rate $\mu \min (n,2)$. Its reversible weights satisfyFor $0<r<1$, the factor $r^{n(n-1)/2}$ makes their sum finite for every $\lambda ,\mu >0$. For $r=1$, the tail is geometric with ratio $\lambda /(2\mu )$, so a stationary distribution exists exactly when $\lambda <2\mu$.

The chain is an irreducible birth-death process. The detailed balance equations areso ${\pi }_{i+1}={\pi }_i/2$. Normalization gives the stationary distributionAn irreducible countable-state Markov chain that possesses a stationary probability distribution is positive recurrent. Hence $X$ is positive recurrent.

The queue length is an M/M/1 queue and hence a birth-death process with birth rate $\lambda$ and death rate $\mu$ away from zero. The detailed-balance equations areWriting $\rho =\lambda /\mu$, they give ${\pi }_n={\pi }_0{\rho }^n$. Since $\rho <1$, this measure is summable and normalization yieldsThe existence of this invariant probability distribution for the irreducible nonexplosive chain proves positive recurrence.