Найти произведение положительных элементов каждого столбца матрицы - Turbo Pascal

type
  TIntArr = array [0..32764] of Integer;
  PIntArr = ^TIntArr;
var
  a: PIntArr;
  m, n, mn, i, j: Integer;
  p: Real;
begin
  Randomize;
  repeat Write('m n: '); Read(m,n); until (m>1) and (n>1);
  mn:=m*n;
  GetMem(a,SizeOf(Integer)*mn);
  for j:=0 to mn-1 do begin
    if j mod n=0 then WriteLn;
    a^[j]:=-9+Random(19); Write(a^[j]:8);
  end; WriteLn;
  WriteLn('Произведения:');
  for j:=0 to n-1 do begin
    p:=1; i:=j;
    while i<mn do begin
      if a^[i]>0 then p:=p*a^[i];
      Inc(i,n);
    end; Write(' ',p:7:0);
  end; WriteLn;
  FreeMem(a,SizeOf(Integer)*mn);
end.