Программа, для обработки двухмерных матриц, Язык - Паскаль

Напишите программу, реализующую обработку двухмерных массивов. Дана матрица размера M x N. В каждом ее столбце найти количество элементов, больших среднего арифметического всех элементов этого столбца. Моя попытка -

const
    M = 5;
    N = 4;
 
var
    a: array[1..M, 1..N] of byte;
    i, j: byte;
    sum: integer;
 
begin
    randomize;
 
    sum := 0;

    for i := 1 to N do begin
        for j := 1 to M do begin
            a[i,j] := random(10);
            write(a[i,j]:3);

            sum := sum + a[i,j]

        end;
        writeln
    end;
 
    writeln ('------------');
 
    writeln('sum = ', sum);
 
end.