Вычисление определителя матрицы

det A= sum (-1)^x*a[1,s1]*a[2,s2]...a[n,sn].

function TMatrix.Det:extended;
var Size,C,R:integer; K:extended; Check:boolean;
begin
  if fCols <> fRows then begin Result:=0.0; Exit end;
  Result:=1.0; Size:=fCols;
  while Size > 1 do begin
    if abs(M[Size-1,Size-1]) <= 1e-100 then begin
      Check:=false;
      for R:=0 to Size-2 do if abs(M[Size-1,R]) > 1e-100 then begin
        Check:=true; for C:=0 to Size-1 do M[C,Size-1]:=M[C,Size-1]+M[C,R]
      end;
      if not Check then begin Result:=0.0; Break end
    end;
    for R:=0 to Size-2 do begin
      K:=-M[Size-1,R]/M[Size-1,Size-1];
      for C:=0 to Size-1 do M[C,R]:=M[C,R]+K*M[C,Size-1]
    end;
    Result:=Result*M[Size-1,Size-1]; Dec(Size)
  end;
  Result:=Result*M[0,0]
end;