An efficient procedure for constructing quasi-exactly solvable matrix
models is suggested. It is based on the fact that the representation s
paces of representations of the algebra sl(2,R) within the class of fi
rst-order matrix differential operators contain finite dimensional inv
ariant subspaces. (C) 1997 Elsevier Science B.V.