An dein Ding geh ich mal mit meinem überlegenen wissen heran;  es wäre nicht verkehrt, wenn du dich mal mit ===> Elementarteilern beschäftigst ( obwohl das zugegebener Maßen nicht nach jedermanns Geschmack ist. )  Eine der wichtigsten Erkenntnisse
    "  Jede Matrix löst ihre eigene Säkulardeterminante ( SD )  "
   Bei diagonalisierbaren Matrizen wäre das ja auch trivial; aber es gilt eben allgemein.
   Kannst du schon Eigenwerte? Die Koeffizienten der  SD  bestimmen natürlich eindeutig über die Eigenwerte;  sämtliche Lösungsmatrizen deiner  SD  haben daher die selben Eigenwerte.
    Würden sie einen Vektorraum bilden.  Angenommen Matrix  A hat Eigenwert  E ;  dann müsste ja Matrix  k  A Eigenwert  k  E  haben.  Demnach kann k A keine Lösung deines Matrixpolynoms sein.