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

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls;
 
type
  TForm1 = class(TForm)
    Memo1: TMemo;
    Button1: TButton;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
matr = array [1..5, 1..5] of real;
rslt = array[1..5] of real;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
procedure matr_init(var matrix: matr);
begin
//Инициализация определителя матрицы
 matrix[1, 1] := 0;
 matrix[1, 2] := 8;
 matrix[1, 3] := 1;
 matrix[1, 4] := 2;
 matrix[1, 5] := -3;
 
 matrix[2, 1] := -8;
 matrix[2, 2] := -9;
 matrix[2, 3] := 9;
 matrix[2, 4] := 0;
 matrix[2, 5] := 9;
 
 matrix[3, 1] := -5;
 matrix[3, 2] := -2;
 matrix[3, 3] := 3;
 matrix[3, 4] := -1;
 matrix[3, 5] := 0;
 
 matrix[4, 1] := -9;
 matrix[4, 2] := 5;
 matrix[4, 3] := -8;
 matrix[4, 4] := 7;
 matrix[4, 5] := -6;
 
 matrix[5, 1] := 1;
 matrix[5, 2] := -7;
 matrix[5, 3] := 2;
 matrix[5, 4] := 2;
 matrix[5, 5] := 0;
end;
 
function DivZero(var A: matr; rozm: integer):matr;
var temp1, temp2: matr;
i, j: byte;
begin
 //Перестановка 1 и 2 строк местами
  for j:=1 to rozm do
   begin
    temp1[1, j] := A[1, j];
    temp2[2, j] := A[2, j];
   end;
 
  for j:=1 to rozm do
   begin
    A[1, j] := temp2[2, j];
    A[2, j] := temp1[1, j];
   end;
 
 Result := A;
end;
 
function DetN(var A: matr; rozm: integer):real; //Функция нахождения определителя
var i, j:integer;
temp: Matr;
del: real;
begin
 
while rozm <> 3 do //Понижение порядка от rozm до 3-го
begin
 
 if A[1, 1] = 0 then a := DivZero(A, rozm); //Если первый элемент равен 0, то переставляем 1 и 2 строки
   del := a[1, 1];
 
for i:=2 to rozm do
 begin
 for j:=2 to rozm do
  begin
   temp[i-1, j-1] := (a[1, 1] * a[i, j] - a[i, 1] * a[1, j]) / del; //Понижение порядка до 3-го
  end;
 end;
 
 a := temp;
 
 rozm := rozm - 1;
 
end; {while}
 
 if rozm = 3 then
 begin
 //Расчёт определителя 3 на 3
 Result := del * (a[1,1]*a[2,2]*a[3,3]+a[1,2]*a[2,3]*a[3,1]
          +a[2,1]*a[3,2]*a[1,3]-a[1,3]*a[2,2]*a[3,1]
          -a[2,1]*a[1,2]*a[3,3]-a[1,1]*a[2,3]*a[3,2]); 
 end;
end;
 
procedure TForm1.Button1Click(Sender: TObject);
var matrica: matr;
begin
 matr_init(matrica);
 Memo1.Lines.Append(FloatToStr(DetN(matrica, 5)));
end;
 
end.