Найти алгебраические дополнения. Треугольная матрица

{Треугольная мартица по методу Гаусса с выбором максимума в столбце}
procedure Per(k,n:integer;var a:Tmatrix);
var j,i,p:integer;z:Real;
begin
z:=a[k,k];i:=k; //текущая строка
for j:=k+1 to n-1 do //ниже ее
 begin
  if abs(a[j,k])>z then //ищем в этом столбце  максимальный элемент
   begin
    z:=a[j,k];
    i:=j;
   end;
 end;
if i>k then  //если нашли
 begin
  for p:=k to n-1 do  //переставляем строки
    begin
     z:=a[i,p];a[i,p]:=a[k,p];a[k,p]:=z;
    end;
 end;
end;
 
 
procedure TForm1.Button7Click(Sender: TObject);
var i,j,k:integer;
    r:real;
begin
  Form1.ClientWidth:=630;
  Form1.Label2.Visible:=True;
  Form1.StringGrid4.Visible:=true;
//создание матрицы
Setlength(a,n,n);
for i:=0 to n-1 do
for j:=0 to n-1 do
a[i,j]:=StrToFloat(StringGrid1.Cells[j,i]);
//преобразование
for k:=0 to n-1 do
 begin
  Per(k,n,a);  //перестановка
  for j:=k+1 to n-1 do  //прямой ход Гаусса
   begin
    r:=a[j,k]/a[k,k];
    for i:=k to n-1 do
    a[j,i]:=a[j,i]-r*a[k,i];
    a[j,n]:=a[j,n]-r*a[k,n];
   end;
 end;
//вывод результата
for i:=0 to n-1 do
for j:=0 to n-1 do
StringGrid4.Cells[j,i]:=FloatToStrF(a[i,j],ffFixed,6,2);
end;