Quasi-Exactly Solvable Matrix Schrödinger Operators

Two families of quasi-exactly solvable 2×2 matrix Schrödinger operators are constructed. The first one is based on a polynomial matrix potential and depends on three parameters. The second is a one-parameter generalization of the scalar Lamé equation. The relationship between these operators and QES Hamiltonians already considered in the literature is pointed out.