Быстрая рекурсивная сортировка для матриц

const n=4;
var
i,j,k,z,x,c:integer;
a:array[1..n, 1..n] of integer;
 
procedure return(k: integer; var i,j: integer);  
begin
case k of
  1:begin  i:=1;  j:=1; end;
  2:begin  i:=1;  j:=2; end;
  3:begin  i:=1;  j:=3; end;
  4:begin  i:=1;  j:=4; end;
  5:begin  i:=2;  j:=3; end;
  6:begin  i:=2;  j:=2; end;
  7:begin  i:=2;  j:=1; end;
  8:begin  i:=3;  j:=1; end;
  9:begin  i:=3;  j:=2; end;
  10:begin  i:=4;  j:=1; end;
  11:begin  i:=4;  j:=2; end;
  12:begin  i:=4;  j:=3; end;
  13:begin  i:=4;  j:=4; end;
  14:begin  i:=3;  j:=4; end;
  15:begin  i:=3;  j:=3; end;
  16:begin  i:=2;  j:=4; end;
 end;
    
end;
 
begin
  for k:=1 to n*n do 
  begin
    return(k,z,x);//заполнение 
    a[z,x]:=k;
  end;
  for i:=1 to n do
  begin
    for j:=1 to n do 
    begin
      write(a[i,j]:3);  //вывод
    end;
    writeln();
  end;
 
 
end.

program MatrixSort;
 
type
  TMatrix = array of array of integer;
 
  procedure SortMatrix(var a: TMatrix; m, n: integer);
  var
    i, j: integer;
    t: integer;
  begin
 
    for i := 0 to (m * n - 1) do
      for j := 0 to (m * n - 1) - i - 1 do
        if a[j div n, j mod n] > a[(j + 1) div n, (j + 1) mod n] then
        begin
          t := a[j div n, j mod n];
          a[j div n, j mod n] := a[(j + 1) div n, (j + 1) mod n];
          a[(j + 1) div n, (j + 1) mod n] := t;
        end;
  end;
 
var
  m, n: integer;
  i, j: integer;
  a: TMatrix;
begin
  {ввод размерности и заполнение матрицы}
  m := 2;
  n := 8;
  SetLength(a, m, n);
  for i := 0 to m - 1 do
    for j := 0 to n - 1 do
      a[i, j] := 10 * (m - i) + (n - j);
  {матрица до сортировки}
  writeln('Before:');
  for i := 0 to m - 1 do
  begin
    for j := 0 to n - 1 do
      Write(a[i, j]: 4);
    writeln;
  end;
  {вызов процедуры сортировки}
  SortMatrix(a, m, n);
  {матрица после сортировки}
  writeln('After');
  for i := 0 to m - 1 do
  begin
    for j := 0 to n - 1 do
      Write(a[i, j]: 4);
    writeln;
  end;
end.