Вычислить определитель матрицы методом разложения по Лапласу

const Size=8;
type TMatrix=array [1..Size, 1..Size] of Integer;
 
Procedure InitMatrix(var Matrix:TMatrix);
var i,j:Byte;
Begin
  Randomize;
  for i:=1 to Size do
    for j:=1 to Size do
      Matrix[i,j]:=Round(Random(19)-9);
End;
 
Function Input:Byte;
Begin
  Write('Введите размерность матрицы (max 8x8): ');
  Readln(Result);
End;
 
Procedure PrintMatrix(Matrix:TMatrix;N:Byte);
var i,j:Byte;
Begin
  Writeln('Matrix A:');
  for i:=1 to N do
    for j:=1 to N do
      Begin
        Write(Matrix[i,j]:2,' ');
        if j=N then Writeln;//перевод курсора на след. строку
      End;
End;
 
Function Power(a,b:Integer):Integer; // Возведение в степень (a^b)
Begin
  if b=0 then Result:=1
    else Result:=Power(a,b-1)*a;
End;
 
Procedure MakeNewMatrix(Matrix:TMatrix;N:Byte; Var NewMatrix:TMatrix); // N-столбец, который не нужно вписывать в новую матрицу
var i,j:Byte;
Begin
  for i:=1 to N-1 do    // сокращаем новую матрицу
    for j:=1 to N-1 do //  на 1 строку и столбец
      Begin
        if j<N then NewMatrix[i,j]:=Matrix[i,j];
        if j>=N then NewMatrix[i,j]:=Matrix[i,j+1];
      End;
End;
 
Function CalcDet(Matrix:TMatrix;M,N:Byte):Integer; //M-размерность исходной матрицы
var NewMatrix:TMatrix;
Begin
  Result:=0;
  if N=1 then Result:=Matrix[1,1] else 
    Begin
      MakeNewMatrix(Matrix,N,NewMatrix);
      Result:=Result+Matrix[N,N]*Power(-1,N+N)*CalcDet(NewMatrix,N-1)+CalcDet(Matrix,N-1); (*ИНДЕКСЫ МАССИВОВ*)
    End;
End;
 
Procedure Output(Det:Integer);
Begin
  Writeln('Determinant(A): ',Det);
End;
 
var Matrix:TMatrix;
    M,N:Byte;
    Det:Integer;
Begin
  InitMatrix(Matrix);
  N:=Input;//Размерность исходной матрицы
  PrintMatrix(Matrix,N);
  M:=N;
  Det:=CalcDet(Matrix,N); (*КОЛ-ВО ПАРАМЕТРОВ*)
  Output(Det);
End.