Метод Леверье(Matlab)

[taki]

Всем привет.Я тут новенький.Надеюсь на вашу отзывчивость
Суть задания:Написать программу и решить контрольный пример (с проверкой стандартной функции MATLAB).Алгоритмы построения передаточной матрицы и характеристического полинома. Метод Леверье.
Сам метод я знаю и могу записать его для матрицы с известным числом строк и столбцов.Однако мне надо нужно написать программу для квадратной матрицы (n,n).

[taki]

Сам додумался только до такого,правда не работает
A=input('matrix A [n x n]: ');
[n,m]=size(A);
a(1)=-trace(A);
b=0;
k=n-1;
for i=2: (n+1)
k=i;
for j=1:i
b1=b+a(j)*trace(A^(k));
k1=k-1;
b=b1;
k=k1;
end;
a(i)=-1/i*b;
b1=0;
b=0;
end;
disp(a);

[taki]

Спасибо за отзычивость.