Найти все комбинации последовательностей элементов квадратной матрицы, не лежащих в одной строке и столбце

...
const
  n=4;
  Mas: Arrya [1..n] of Arra [1..n] of Byte = ((1,2,3,4),
                                              (4,1,2,3),
                                              (3,4,1,2),
                                              (2,3,4,1));
...
procedure TForm1.Button1Click(Sender: TObject);
var Ar: Array [1..n] of Byte;
  i,j: Byte;
procedure MyInc(var n:Byte);
begin
If n=4 then n:=1 else Inc(n);
end;
procedure ShowMas(mas: Array of Byte);
var i:Byte;
  Buf: String;
begin
For i:=Low(mas) to High(mas) do Buf:=Buf+IntToStr(mas[i])+'; ';
Form1.Memo1.Lines.Add(Buf);
end;
begin
i:=1;
While i<n+1 do
begin
  For j:=1 to n do
  begin
    Ar[j]:=Mas[j,i];
    MyInc(i);
  end;
  Inc(i);
  ShowMas(Ar);
end;
end;
...

...
type
  TMas = Array of Array of Byte;
...
const
  n = 4;
  Mas: Array [0..n-1] of Array [0..n-1] of Byte = ((10,11,12,13),
                                                   (14,15,16,17),
                                                   (18,19,20,21),
                                                   (22,23,24,25));
...
Procedure ShowMas(Mas: TMas);
var i,j: Byte;
  Buf: String;
begin
For i:=Low(Mas) to High(Mas) do
begin
  Buf:='';
  For j:=Low(Mas[i]) to High(Mas[i]) do Buf:=Buf+IntToStr(Mas[i,j])+'; ';
  Form1.Memo1.Lines.Add(Buf);
end;
end;
 
Procedure ShowAr(Mas: TMas);
var i,j: Byte;
  Buf: String;
begin
For i:=Low(Mas) to High(Mas) do
begin
  Buf:='';
  For j:=Low(Mas[i]) to High(Mas[i]) do Buf:=Buf+IntToStr(Mas[i,j])+'; ';
  Form1.Memo1.Lines.Add(Buf);
end;
end;
 
function Variant(Mas: TMas): TMas;                                              //Возвращает перечень главных и побочных диагоналей
var BufMas: TMas;
  i,j,k,l: Byte;
begin
If Length(Mas)=2 then                                                           //массив 2х2 (меньше уже смешно, а больше будет сложнее)
begin
  SetLength(Result,2);
  SetLength(Result[0],2);
  SetLength(Result[1],2);
  Result[0,0]:=Mas[0,0];  //Главная
  Result[0,1]:=Mas[1,1];  //диагональ;
  Result[1,0]:=Mas[0,1];  //Побочная
  Result[1,1]:=Mas[1,0];  //диагональ;
end else
begin
  For i:=High(Mas) downto Low(Mas) do                                           //Проходим по первой строчки
  begin
    j:=0;
    k:=0;
    SetLength(BufMas,Length(Mas)-1);
    While j<n do
    begin
      If i<>j then
      begin
        For l:=Low(BufMas) to High(BufMas) do
        begin
          SetLength(BufMas[l],k+1);
          BufMas[l,k]:=Mas[l+1,j];                                              //Генерируем массив N-1xN-1 убирая 1-ую строчку и i-ый столбец
        end;
        Inc(k);
      end;
      Inc(j);
    end;
    BufMas:=Variant(BufMas);                                                    //Передаем сгенерированный массив снова в эту процедуру пока не будет массив 2х2, ну а при 2х2 всего 2 возможные комбинации
    For j:=Low(BufMas) to High(BufMas) do
    begin
      SetLength(Result,Length(Result)+1);
      SetLength(Result[Length(Result)-1],Length(BufMas[j])+1);
      For k:=Low(BufMas[j]) to High(BufMas[j]) do 
        Result[Length(Result)-1,k]:=BufMas[j,k];                                //Добавляем возможные результаты в нашу функцию (функции будут запускаться много раз, самая первая и будет иметь конечный результат, остальные будут иметь результат на усеченные массивы)
      Result[Length(Result)-1,Length(BufMas[j])]:=Mas[0,i];                     //Добавляем в конец результатов предыдущих функций, нашу переменную в i-ом столбце
    end;
  end;
end;
end;
 
procedure TForm1.Button3Click(Sender: TObject);
var BufMas: TMas;
  i,j: Byte;
begin
SetLength(BufMas,Length(Mas));
For i:=Low(Mas) to High(Mas) do
begin
  SetLength(BufMas[i],Length(Mas[i]));
  For j:=Low(Mas[i]) to High(Mas[i]) do BufMas[i,j]:=Mas[i,j];
end;
BufMas:=Variant(BufMas);
ShowMas(BufMas);
end;
...

function Variant(Mas: TMas): TMas;                                              //Возвращяет перечень главных и побочных диагоналей
var BufMas: TMas;
  i,j,k,l: Byte;
begin
If Length(Mas)=2 then                                                           //массив 2х2 (меньше уже смешно, а больше будет сложнее)
begin
  SetLength(Result,2);
  SetLength(Result[0],2);
  SetLength(Result[1],2);
  Result[0,0]:=Mas[0,0];  //Главная
  Result[0,1]:=Mas[1,1];  //диагональ;
  Result[1,0]:=Mas[0,1];  //Побочная
  Result[1,1]:=Mas[1,0];  //диагональ;
end else
begin
  Setlength(Result,0);  //Я тут забыл очистить память...
  For i:=High(Mas) downto Low(Mas) do                                           //Проходим по первой строчки
...

type TMas=Array of Array string;
Mas:TMas;
...
SetLength(Mas,d,u)
...
function Variant(Mas: TMas): TMas;
var BufMas: TMas;
  v,c,s,l: integer;
begin
  If Length(Mas)=2 then 
  begin
    SetLength(Result,2);
    SetLength(Result[0],2);
    SetLength(Result[1],2);
    Result[0,0]:=Mas[0,0];
    Result[0,1]:=Mas[1,1];
    Result[1,0]:=Mas[0,1];
    Result[1,1]:=Mas[1,0];
  end else
  begin
    Setlength(Result,0);
    For v:=High(Mas) downto Low(Mas) do begin
      c:=0;
      s:=0;
      SetLength(BufMas,Length(Mas)-1);
      While c<d do begin
        If v<>c then begin
          For l:=Low(BufMas) to High(BufMas) do begin
            SetLength(BufMas[l],s+1);
            BufMas[l,s]:=Mas[l+1,c]; 
          end;
          Inc(s);
        end;
        Inc(c);
      end;
    BufMas:=Variant(BufMas); \\ На эту строку выдаёт ошибку Access violation
    For c:=Low(BufMas) to High(BufMas) do
    begin
      SetLength(Result,Length(Result)+1);
      SetLength(Result[Length(Result)-1],Length(BufMas[c])+1);
      For s:=Low(BufMas[c]) to High(BufMas[c]) do Result[Length(Result)-1,s]:=BufMas[c,s];                                
      Result[Length(Result)-1,Length(BufMas[c])]:=Mas[0,v];
    end;
  end;
end;
end;