N-fold Supersymmetry in Quantum Mechanical Matrix Models

We formulate N-fold supersymmetry in quantum mechanical matrix models. As an example, we construct general two-by-two Hermitian matrix 2-fold supersymmetric quantum mechanical systems. We find that there are two inequivalent such systems, both of which are characterized by two arbitrary scalar functions, and one of which does not reduce to the scalar system. The obtained systems are all weakly quasi-solvable.