Для рассчитанной матрицы C расположить (сортировать) ее элементы в порядке возрастания по столбцам и вывести на печать

const n=5;
var a,a2,b,e:array[1..n,1..n] of integer;
    i,j,k,max,imx,s:integer;
begin
writeln('Матрица А');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    if j=i then a[i,j]:=i
    else a[i,j]:=n-i-j;
    write(a[i,j]:5);
   end;
  writeln;
 end;
//матрица A^2
for i:=1 to n do
for j:=1 to n do
 begin
  a2[i,j]:=0;
  for k:=1 to n do
  a2[i,j]:=a2[i,j]+a[i,k]*a[k,j];
 end;
writeln('Матрица В=5А^2-7A-5E');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    b[i,j]:=5*a2[i,j]-7*a[i,j];
    if j=i then b[i,j]:=b[i,j]-5;
    write(b[i,j]:5);
   end;
  writeln;
 end;
max:=-maxint-1;
for i:=1 to n do
 begin
  s:=0;
  for j:=1 to n do
  s:=s+sqr(b[i,j]);
  if s>max then
   begin
    max:=s;
    imx:=i;
   end;
 end;
writeln('Максимальная сумма квадратов в строке ',imx);
for j:=1 to n do
write(b[imx,j]:6);
end.

const nmax=10;
//процедура обмена 2х элементов
procedure Swap(var a,b:integer);
var tmp:integer;
begin
tmp:=a;
a:=b;
b:=tmp;
end;
var a:array[1..nmax,1..nmax] of integer;
    m,n,i,j,k:byte;
    x:integer;
begin
repeat
randomize;
write('Количество строк до ',nmax,' m=');
readln(m);
until m in [1..nmax];
repeat
write('Количество столбцов до ',nmax,' n=');
readln(n);
until n in [1..nmax];
writeln('Исходная матрица:');
for i:=1 to m do
 begin
  for j:=1 to n do
   begin
    a[i,j]:=-10+random(21);
    write(a[i,j]:4);
   end;
  writeln;
 end;
for i:=1 to m*n-1 do
for j:=i+1 to m*n do
if a[(i-1)mod m+1,(i-1)div m+1]>a[(j-1)mod m+1,(j-1)div m+1]
then Swap(a[(i-1)mod m+1,(i-1)div m+1],a[(j-1)mod m+1,(j-1)div m+1]);
writeln('Отсортированная матрица:');
for i:=1 to m do
 begin
  for j:=1 to n do
  write(a[i,j]:4);
  writeln;
 end;
end.