Написать программу нахождения матрицы, равной (А-В)^Т(в степени Т)

for i:=1 to n do
for j:=1 to n do
a[i,j]:=a[i,j]-b[i,j];

for i:=1 to n do
for j:=1 to n do
at[i,j]:=a[j,i];

uses crt;
var a,b:array[1..50,1..50] of integer;{2 матрицы}
    n,i,j,x:integer;{размер, 2 счетчика циклов, буфер для обмена элементов}
begin
clrscr;
randomize;
repeat
write('Введите размер матриц n=');
readln(n);
until n in [1..50];{проверяем чтоб n было в диапазоне}
writeln('Матрица А:');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    a[i,j]:=random(20);
    write(a[i,j]:4);
   end;
  writeln;
 end;
writeln('Матрица B:');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    b[i,j]:=random(20);
    write(b[i,j]:4);
   end;
  writeln;
 end;
writeln('Матрица А=А-В:');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    a[i,j]:=a[i,j]-b[i,j];{вычитаем из элементов матрицы А элементы матрицы В}
    write(a[i,j]:4);
   end;
  writeln;
 end;
writeln('Транспонированная матрица At=(A-B)^T:');
for i:=1 to n do
for j:=1 to n do
if j<i then {меняем местами элементы ниже главной диагонали с элементами выше главной диагонали}
 begin
  x:=a[i,j];
  a[i,j]:=a[j,i];
  a[j,i]:=x;
 end;
for i:=1 to n do
 begin
  for j:=1 to n do
  write(a[i,j]:4);
  writeln;
 end;
readln
end.