Вычислить среднее арифметическое каждого столбца над главной диагональю (диагональ включать) матрицы B[1..n,1..n] и запишите данные значения в последнюю строку данной матрицы — Pascal(Паскаль)

uses crt;
const n=4;
var b:array[1..n,1..n] of real;
    a:array[1..n] of real;
    i,j,k,x:integer;
begin
ClrScr;
Randomize;
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    b[i,j]:=-10+random(20);
    Write(b[i,j]:5:1);
   end;
  Writeln;
 end;
Writeln;
for j:=1 to n do
 begin
  a[j]:=0;
  k:=0;
   for i:=1 to n do
    if i<=j then
     begin
      k:=k+1;
      a[j]:=a[j]+b[i,j];
     end;
    a[j]:=a[j]/k;
 end;
Writeln;
Writeln;
for i:=1 to n do
 begin
  for j:=1 to n do
   b[n,j]:=a[j];
 end;
for i:=1 to n do
 begin
  for j:=1 to n do
   Write(b[i,j]:5:1);
  Writeln;
 end;
Readln;
end.