Матрица. Если максимальные элементы столбцов находятся на главной диагонали, необходимо найти среднее арифметическое

const n=3;
 
var i,j,max,minI,minJ:integer;
    average,count:double;
    Matrix:array[1..n,1..n]of double;
    temp:boolean;
begin
  for i:=1 to n do begin
    for j:=1 to n do begin
      Read(Matrix[i,j]);
    end;
  end;
  
  for i:=1 to n do begin
    for j:=1 to n do begin
      Write(Matrix[i,j]+'  ');
    end;
    writeln();
  end;
  
  temp:=true;  
  i:=1;
  while (temp=true) and(i<=n) do begin
    max:=1;
    j:=1;    
    while(temp=true) and (j<=n) do begin
      if (Matrix[i,j]>Matrix[i,max]) then max:=j;
      j+=1;
    end;
    if (max<>i)then temp:=false;
    inc(i);
  end;
  
  if (temp=true) then begin
      minI:=1;
      minJ:=1;
      for i:=1 to n do
        for j:=1 to n do
          if(Matrix[i,j]<Matrix[minI,minJ]) then begin
            minI:=i;
            minJ:=j;
          end;
      
      count:=0;
      average:=0;
      for i:=2 to n do begin
        for j:=1 to i-1 do begin
          count:=count+1;
          average:=average+Matrix[i,j];
        end;
      end;
      
      Matrix[minI,minJ]:=Matrix[minI,minJ]+(average/count);
      
      Writeln('Результрующяя матрица');
      for i:=1 to n do begin
        for j:=1 to n do begin
          Write(Matrix[i,j]:2:1,'  ');
        end;
        writeln();
      end;
      
  end
  else Writeln('На главной диагонале не максимальные элементы');
end.

begin
  var (n, count) := (ReadlnInteger('N ='), 0);
  var a := MatrGen(n, n, (i, j) -> ReadlnReal(Format('a[{0},{1}] =', i, j))); a.Println;
  var maxs := a.Cols.Select(v -> v.Max).ToArray;
  for var j := 0 to a.GetLength(1) - 1 do
    if maxs[j] = a[j, j] then count += 1;
  if count = n then begin
    var average := a.ElementsWithIndexes.Where(v -> v[1] > v[2]).Average(v -> v[0]);
    var min := a.ElementsByCol.Min;
    a.Transform(v -> v = min ? v + average : v);
    writeln('Новая матрица:'); a.Println;
  end else writeln('Матрица остаётся без изменений');
end.

begin
  var n := ReadlnInteger('n = ');
  var mx := MatrGen(n, n, (i, j) -> ReadlnReal(Format('mx[{0}, {1}] = ', i, j))).Println(5);
  
  if mx.Cols.Numerate(0).All(v -> v.Item2.ToArray.IndexesOf(it -> it = v.Item2.Max).Contains(v.Item1)) then
  begin
    var avg := mx.ElementsWithIndexes.Where(v -> v[1] > v[2]).Average(v -> v[0]);
    var min_value := mx.ElementsByRow.Min;
    mx.Transform(v -> v = min_value ? v + avg : v);
    Println('Измененная матрица:'); mx.Println;
  end
  else writeln('Матрица остается без изменений');
  
end.