The stationary probability distribution vector, x, associated with an ergodic finite Markov chain satisfies a homogeneous singular system of equations, Ax=0, where A is a real and generally unsymmetric square matrix of the form A=I-T. Here I is the identity matrix and T is the chain’s column stochastic matrix. In many applications A is very large and sparse, and in such cases it is desirable to exploit this property in computing x. This paper reviews some of the literature dealing with sparse techniques for solving the above system of equations, and in so doing attempts to present a variety of methods from a unified point of view.