Получить матрицу, образованную по соответствующему закону

const
 n=6;
 m=6;
var
 i,j,k,ij : integer;
 a,b,c,d,e : array[1..n,1..m] of real;
 r,sum, min: real;
 
procedure vivod(tmp: array[1..n,1..m] of real; tx: string; k1,k2: integer);
 var
  ii, jj: integer;
 
begin
 writeln('   Матрица ', tx);
 for ii:=1 to n do
  for jj:=1 to m do
   begin
    write(tmp[ii,jj]:k1:k2);
    if jj=m then writeln;
   end;
  writeln;
 
end;
 
 
 begin
  {заполнение и вывод матрицы A}
 writeln('   Матрица А');
 for i:=1 to n do
  for j:=1 to m do
   begin
    a[i,j]:=(i+j)*sin(j)/cos(i);
    write(a[i,j]:9:3);
    if j=m then writeln;
   end;
 
  writeln;
  
  b:=a;
  for i:=1 to n do
   for k:=1 to n-1 do
    for j:=k to m do
     if b[i,k]>b[i,j] then
      begin
       r:=b[i,j];
       b[i,j]:=b[i,k];
       b[i,k]:=r;
      end;
 
  {вывод матрицы B}
  vivod(b,'B',9,3);
  
  {заполнение матрицы C}
 for i:=1 to n do
  for j:=1 to m do
   begin
    c[i,j]:=1;
    if (i>j) and (i+j<n+1) then c[i,j]:=0;
    if (i<j) and (i+j>n+1) then c[i,j]:=0;
   end;
   
 {вывод матрицы C}
  vivod(c,'C',3,0);
 
  {вычисление и вывод матрицы D}
 writeln('   Матрица D');
  for i:=1 to n do
   for j:=1 to m do
    begin
     sum:=0;
     for k:=1 to n do
      begin
       sum:=sum+b[i,k]*c[k,j];
       d[i,j]:= sum;
      end;
     write(d[i,j]:9:3);
     if j=m then writeln;
    end;
   writeln;
 
 {вычисление матрицы E} {конкретно для матрицы 6х6}
  e:=d;
  for i:=4 to n do
   for j:=1 to 3 do
    begin
     r:=e[i,j];
     e[i,j]:=e[i-3,j+3];
     e[i-3,j+3]:=r;
    end;
    
  {вывод матрицы E}
    vivod(e,'E',9,3);
 
 {вычисление минимального элемента главной диагонали}
  min:=e[1,1];
  for i:=2 to n do
   if e[i,i]<min then
    begin
     min:=e[i,i];
     ij:=i;
    end;
  writeln(' Минимальный элемент: E[',ij,',',ij,'] = ',min:8:3)
 end.