In order to determine the stationary distribution for discrete time quasi-birthdeath
Markov chains, it is necessary to find the minimal nonnegative solution of a
quadratic matrix equation. The Newton-Shamanskii method is applied to solve this
equation, and the sequence of matrices produced is monotonically increasing and converges to its minimal nonnegative solution. Numerical results illustrate the effectiveness of this procedure.