Den Fall a0 ≠ 0 hast du ja.
Für den Fall ao = 0 kannst du doch dein Verfahren iterieren bis zu ersten Koeffizient der nicht 0 ist. Das muss spätestens bei an-1 der Fall sein, da anAn nicht einfach = 0 sein kann, wegen dem Rang der Matrix. Sagen wir bei ak
g(A)=anAn+an-1An-1+...+a1A+a0 = 0
(anAn-1+an-1An-2+...+a1)A = - a0
|Falls a0 =0
= anAn-1+an-1An-2+...+a1 = 0
= (anAn-2+an-1An-3+...+a2)A = -a1
…
Allgemein: Sek k das kleinste k, für das ak ≠ 0
g(A)=anAn+an-1An-1+...+akAk = 0
(anAn-k+an-1An-k-1+...+ak+1)Ak+1 + akAk= 0
(anAn-k+an-1An-k-1+...+ak+1)Ak+1 = - akAk
(anAn-k+an-1An-k-1+...+ak+1)Ak+1 A^{-k-1} = - akAk * A^{-k-1}
(anAn-k+an-1An-k-1+...+ak+1) = - akA-1 |:(-ak)
usw, wie du's gemacht hast.
Sollte mE genügen.
Nenne aber dein erstes Polynom nicht schon g, sondern zB. p, damit du das resultierende Polynom für A^{-1} mit g bezeichnen kannst.