Метод использует для нахождения собственных значений характеристический полином, который по созданному Леверрье и усовершенствованному Фаддеевым алгоритму получают следующим образом:
строят последовательность матриц
A1=A; Sp A1=p1; B1=A1–p1E;
A2=AB1; Sp A2=p2; B2=A2–p2E;
………………………………………………..
An-1=ABn-2; (1/(n-1))Sp An-1=pn-1; Bn-1=An-1–pn-1E;
An=ABn-1; (1/n)Sp An=pn; Bn=An–pnE;
В результате получают уравнение
λn–p1λn-1–p2λn-2–…–pn=0.
Попутно можно получить обратную матрицу
A-1=Bn-1/pn
Собственные векторы находим по формулам
x0=e; 
i=1, 2, …, n-1
где e – столбец единичной матрицы, 
– одноименный столбец матрицы Вk, а собственный вектор 
соответствует λk.
