Написать программу, которая вычисляет определитель квадратной матрицы любого порядка(универсальная, для порядка n)

uses crt;
type Tmatr=array [1..20,1..20] of real;
{перестановка строк матрицы, если элемент на главной диагонали=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;
    r:real;
begin
det:=1;
for k:=1 to n do
   begin
     if a[k,k]=0 then per(k,n,a,p);
     det:=znak(p)*det*a[k,k];
     for j:=k+1 to n do
       begin
         r:=a[j,k]/a[k,k];
         for i:=k to n do
           begin
             a[j,i]:=a[j,i]-r*a[k,i];
           end;
       end;
   end;
end;
var n:integer;
    a:Tmatr;
    det:real;
    i,j:integer;
begin
clrscr;
write('n=');readln(n);
writeln('Ischodnaia matritsha');
randomize;
for i:=1 to n do
  begin
    for j:=1 to n do
      begin
        a[i,j]:=10*random-2;
        write(a[i,j]:6:2);
      end;
  writeln;
  end;
opr(n,a,det);
write('opr=',det:0:0);
readln
end.