Найти определитель квадратной матрицы

program gauss;
uses crt;
const y=10;
type t=array [1..y,1..y] of real;
var i,j,k,n,u,m,cif:integer;
    z,opr:real;
    a:t;
begin
     clrscr;
     writeln('vvdite rang');
     readln (n);
     writeln('vvedite elementi');
     for i:=1 to n do
         for k:=1 to n do readln (a[i,k]);
         u:=0;
     cif:=0;
     opr:=1;
     repeat
           u:=u+1;
           k:=u;
           while (a[k,u]=0) and (k<=n) do
                 begin
                 k:=k+1;
                 if k>=n then begin
                              writeln (' opr=0 ');
                              halt (1)
                              end;
                 end;
           cif:=cif+1;
           opr:=opr*a[u,u];
           if k<>u then for m:=u to n do begin
                                           z:=a[u,m];
                                           a[u,m]:=a[k,m];
                                           a[k,m]:=z
                                           end;
           for j:=n downto u do a[u,j]:=a[u,j]/a[u,u];
           for i:=u+1 to n do
               for j:=u+1 to n do
                   a[i,j]:=a[i,j]-a[i,u]*a[u,j]
     until u=n;
     if (cif mod 2)=0 then opr:=opr
     else opr:=-opr;
     writeln(opr:6:1);
     readln
end.