We present an algorithmic approach to find the stationary probability distribution of M/G/1-type Markov chains which arise frequently in performance analysis of computer and communication networks. The approach unifies finite- and infinite-level Markov chains of this type through a generalized state-space representation for the probability generating function of the stationary solution. When the underlying probability generating matrices are rational, the solution vector for level k, x/sub k/, is shown to be in the matrix-geometric form x/sub k+1/=gF/sup k/H, k/spl ges/0, for the infinite-level case, whereas it takes the modified form x/sub k+1/=g/sub 1/F/sup k//sub 1/H/sub 1/+g/sub 2/F/sup K-k-1//sub 2/H/sub 2/, 0/spl les/k/spl les/K, for the finite-level case. The matrix parameters in the above two expressions can be obtained by decomposing the generalized system into forward and backward subsystems, or, equivalently, by finding bases for certain generalized invariant subspaces of a regular pencil /spl lambda/E-A. We note that the computation of such bases can efficiently be carried out using advanced numerical linear algebra techniques including matrix-sign function iterations with quadratic convergence rates or ordered generalized Schur decomposition. The simplicity of the matrix-geometric form of the solution allows one to obtain various performance measures of interest easily, e.g., overflow probabilities and the moments of the level distribution, which is a significant advantage over conventional recursive methods.