Вычислить среднее арифметическое положение элементов двумерного массива, попадающих в интервал (b,c) — Pascal(Паскаль)

program PascalGuru;

var
  A: array [1 .. 50, 1 .. 50] of integer; { объявление переменных }
  n, m, i, j, b, c, summa, k: integer;

begin
  write('N= ');
  readln(n);
  write('M= ');
  readln(m);

  for i := 1 to n do
    for j := 1 to m do
    begin
      write('A[', i, ',', j, ']=');
      readln(A[i, j]);
    end;

  writeln('Vvedennyi massiv:');
  for i := 1 to n do
  begin
    writeln;
    for j := 1 to m do
      write(A[i, j]:8);
  end;
  writeln;

  write('b= ');
  readln(b);
  write('c= ');
  readln(c);

  summa := 0;
  k := 0;
  for i := 1 to n do
    for j := 1 to m do
      if (A[i, j] >= b) and (A[i, j] <= c) then
      begin
        inc(k);
        summa := summa + A[i, j];
      end;

  writeln;
  writeln('Sr.Ar= ', summa / k:0:2);

  readln;

end.