Дан вещественный массив А[4][3]. Найти среднее арифметическое значение положительных элементов массива А - Pascal

uses crt;
 
type T2DArray = array [1..4, 1..3] of real;
 
function matrix: T2DArray; //считываем матрицу
var i, j: integer;
begin
    for i:=1 to 4 do
        for j:=1 to 3 do
            readln(matrix[i, j]);
    clrscr;
end;
 
function ArithmeticMean(a: T2DArray): real; //находим среднее арифметическое положительных элементов
var
    i, j, k: integer;
begin
    ArithmeticMean:=0;
    k:=0;
    for i:=1 to 4 do
        for j:=1 to 3 do
            if a[i, j]>0 then
            begin
                ArithmeticMean:=ArithmeticMean+a[i, j];
                inc(k);
            end;
    ArithmeticMean:=ArithmeticMean/k;
end;
 
begin
    randomize;
    write(ArithmeticMean(matrix));
    readln;
end.