расчет алгебраического дополнения, минор каждого элемента, транспортирование матрицы

program opred;
uses crt;
const c=10;
type Tmatr=array [1..c, 1..c] of real;
var a:Tmatr;
    i, j, n:integer;
    det:real; {Определитель}
 
{Процедура переустановки строк, чтобы главный элемент не оказался 0 или
 близким к 0 значением}
procedure Per(k,n:integer;var a:Tmatr; var p:integer);
var i, j:integer;
       z:real;
begin
   z:=a[k,k]; {После...}
   i:=k; {каждого...}
   p:=0; {преобразования...}
   for j:=k+1 to n do  {ищем по оствшимся строкам...}
     begin
       if abs(a[j,k])>z then {максимальный по модулю элемент}
          begin
            z:=abs(a[j,k]);i:=j; {Запоминаем номер строки}
            p:=p+1; {Считаем кол-во переустановок, т.к. в каждой...}
                    {переустановке меняется знак определителя}
          end;
     end;
   if i>k then  {Если эта строка ниже данной}
   for j:=k to n do
     begin
       z:=a[i,j];a[i,j]:=a[k,j];a[k,j]:=z; {тогда делаем переустановку}
     end;
end;
 
{Функция знака определителя}
function Znak(p:integer):integer;
 begin
  if p mod 2=0 then {Если четное кол-во переустановок...}
  znak:=1 else znak:=-1; {"+", если нет "-"}
 end;
 
{Процедура вычисления определителя}
procedure opr(n:integer;var a:Tmatr;var det:real);
var k, i, j, p:integer;
    delenie:real;
begin
 det:=1;
 for k:=1 to n do
  begin
   if a[k,k]=0 then {Если главный элемент = 0,}
      Per(k,n,a,p); {делаем переустановку}
    det:=Znak(p)*det*a[k,k]; {Меняем знак определителя}
    for j:=k+1 to n do  {Ниже делаем преобразования}
     begin
      delenie:=a[j,k]/a[k,k];
       for i:=k to n do
        begin
         a[j,i]:=a[j,i] - delenie * a[k,i];
        end;
     end;
   end;
end;
 
{Основная программа}
begin
 clrscr;
 writeln('Вычислить определитель матрицы.');
 writeln('(Определитель можно вычислить только для квадратной матрицы).');
 writeln;
 writeln('Введите кол-во элементов в строке матрицы и нажмите ENTER');
 writeln('(Число элементов в строке будет равно числу элементов в столбце!):');
 readln(n);
 writeln;
 writeln('Вводите коэфф-ты матpицы A по стpокам нажимая ENTER:');
  for i:=1 to n do
   begin
    for j:=1 to n do
        read(a[i,j]);
   end;
  writeln;
  writeln;
  writeln('Матрица:');
  for i:=1 to n do
   begin
    for j:=1 to n do
        write (a[i, j]:5:3,' ':2); {Вывод матрицы с отступами}
        writeln;
   end;
  writeln;
 
 opr(n,a,det); {Процедура определителя}
 
 write('Определитель = ',det:2:0, '.');
 
 readkey;
 
end.