Посчитать определитель матрицы вещественных чисел

const N=3;
type matr=array[1..N,1..N] of real;
var
  a:matr;
  i,j:integer;
function det(N:integer;a:matr):double;{взял в книге А.Е. Мудрова}
var
  i,j,k:integer;
  p,r,s:real;
begin
  p:=1.0;  s:=0;
  for k:=1 to n-1 do begin
    s:=a[k,k];
    j:=k;
    for i:=k+1 to n do begin
      r:=a[i,k];
      if abs(r)>abs(s) then begin 
        s:=r;
        j:=i;
      end;
    end;
    if s=0.0 then exit;
    if j<>k then begin
      p:=-p;
      for i:=k to n do begin
        r:=a[k,i];
        a[k,i]:=a[j,i];
        a[j,i]:=r;
     end;
   end;
   for j:=k+1 to n do a[k,j]:=a[k,j]/s;
   for i:=k+1 to n do begin
     r:=a[i,k];
     for j:=k+1 to n do a[i,j]:=a[i,j]-a[k,j]*r;
   end;
   p:=p*s;
 end;
 det:=p*a[n,n];
end; 
 
procedure printm(N:integer;a:matr);
var
  i,j:integer;
begin
  for i:=1 to N do begin
    for j:=1 to N do write(a[i,j]:5);
    writeln;
  end;  
end;
 
begin
  randomize;
  for i:=1 to N do 
    for j:=1 to N do 
      a[i,j]:=random(10)-5;
  printm(N,a);
  writeln(' Det(A)=',det(N,a):0:6);
end.