Hier ich erklär dir mal, wie man das richtig macht.   Die lassen dich doch dumm sterben.
   Die wichtigste Erkenntnis;  Diktat für Formelsammlung, Regelheft und Spickzettel:
   "  Jede Matrix löst ihre eigene Säkulardeterminante ( SD )  "
     Fürdiagonalisierbare  ===>  halbeinfache Matrizen lässt sich das ja trivial  einsehen; aber es gilt eben allgemein.
    Und wie findet man die  SD?  Komisch; in den BNüchern ist das immer so Mega  kompliziert erklärt.  Du machst erst mal den  quadratischen   Ansatz
      p_M  (  x  )  =  x  ²  -  p  x  +  q       (  1  )
     p  und  q  folgen über Vieta das geschmähte  Stiefkind.
       p  =  E1  +  E2  =  Sp  (  M  )         (  2a  )
       q  =  E1  E2  =  det  (  M  )        (  2b  )
   p_M  ( x )  =  x  ²  -  x  Sp  (  M  )  +  det  (  M  )    (  2c  )
    Jetzt tust du die Matrix  M einsetzen in ihr Eigenpolynom  (  2c  )
        M  ²  -  M  Sp  (  M  )  +  det  (  M  )  *  1|  =  0  |   *  M  ^ -  1    (  3a  )
     M  -  1|  *  Sp  (  M  )  +  M  ^ -  1  det  (  M  )  =  0    (  3b  )
   M  ^ -  1  =   [  1 / det ( M ) ]  [  Sp ( M ) * 1| - M ]   (  3c  )  (  3c  )