Обратная матрица - Delphi - Ответ 4615361

{$R *.dfm}
type matrice=array of array of extended;
function Determinant(A:matrice):extended;
var det:extended;
    n,k,i,j:integer;
    M:matrice;
begin
n:=high(A);
setlength(M,n,n);
if n=0 then Determinant:=A[0,0] {caz elementar}
  else
    begin
    det:=0;
    for k:=0 to n do
      begin
        for i:=0 to n-1 do
        for j:=0 to k-1 do M[i,j]:=A[i+1,j];
        for i:=0 to n-1 do
        for j:=k to n-1 do M[i,j]:=A[i+1,j+1];
        if odd(k) then det:=det-A[0,k]*Determinant(M)
          else det:=det+A[0,k]*Determinant(M);
        end;
        Determinant:=det;
    end;
end;
procedure Inversa(A:matrice;var Inv:matrice);
var minor:matrice;
    detmin,det:Extended;
    n,i,j,kx,ky:Integer;
begin
  n:=High(A);
  SetLength(minor,n,n);
  for i:=0 to n do
  for j:=0 to n do
  begin
    for kx:=0 to i-1 do
    begin
    for ky:=0 to j-1 do
    minor [kx,ky]:=A[i,j];
    for ky:=j to n-1 do
    minor[kx,ky]:=A[i,j+1];
  end;
  for kx:=i to n-1 do
  begin
    for ky:=0 to j-1 do minor[kx,ky]:=A[i+1,j] ;
      for ky:=0 to n do minor[kx,ky]:=A[i+1,j+1];
  end;
  det:=Determinant(A);
  detmin:=Determinant(minor);
  if Odd(i+1) then Inv[j,i]:=-detmin/det
  else
  Inv[j,i]:=detmin/det;
end;
end;