A MATLAB eig () függvénye nagy vonalakban milyen algoritmus alapján dolgozik?

Itt van némi infó az algoritmusáról:

[link]

Először a karakterisztikus polinom alapján ki kell számítani a mátrixod sajátértékeit. (III. fokú polinomig létezik megoldóképlet, a fölött vannak közleítő megoldások). A sajátértékek ismeretében visszahelyettesítesz (Av_i=lambda_i v_i), és lineáris egyenletrendszereket kapsz, azt pedig átrendezed, és megoldod Gauss eliminációval. Így lesznek sajátvektoraid.

[link]