Не могу найти ошибки!

[nikolairzn]

В двумерном массиве A[N*N] проанализировать эллемент на главной диагонали.Если он >0 то найти среднее арифметическое этого столбца затем заменить это число на главной диагонали на это среднее арифметическое.Если <=0 то столбцы обнулить.
Паскаль изучаю недавно , но уже появилась необходимость решить задачу ,в общем то код написал но есть ошибки поправьте пожалуйста буду очень благодарен!


program dvummas;
uses crt;
const N=5;
var A:array [1..N,1..N] of real;
i,j:integer;
sr,sum:real;

BEGIN
clrscr;randomize;
for i:=1 to N do
for j:=1 to N do
A[i,j]:=random(10)-4;
readkey;
writeln;
for i:=1 to N do begin
for j:=1 to N do
write(A[i,j]:4);
writeln;
end;
readkey;
for i := 1 to N do
for j := 1 to N do
begin
if (i = j) then
if A[i, j] > 0 then
for j:=1 to N do begin
sum:=sum+A[i,j];
sr:=sum/N;
A[i,j]:=sr;
if A[i, j] <= 0 then
for j:=1 to n do
A[i,j]:=0;
end;
end;
for i:=1 to N do begin
for j:=1 to N do
write(A[i,j]:4);
end;
readkey;
END.

[BDA]

program dvummas;
uses crt;
const 
    N = 5;
var
    A: array [1..N, 1..N] of real;
    i, j: integer;
    sum: real;
 
BEGIN
    clrscr;
    randomize;
    for i := 1 to N do
        for j := 1 to N do
            A[i, j] := random(21)-10; //чила из [-10, 10]
    writeln('Source matrix');
    for i := 1 to N do
    begin
        for j := 1 to N do
            write(A[i, j]:4:2, ' ');
        writeln;
    end;
    for j := 1 to N do
        if A[j, j] > 0 then
        begin
            sum := 0;
            for i:=1 to N do
                sum := sum + A[i, j];
            A[j, j] := sum / N;
        end
        else
            for i := 1 to N do
                A[i, j] := 0;
    writeln('Edited matrix');
    for i := 1 to N do
    begin
        for j := 1 to N do
            write(A[i, j]:4:2, ' ');
        writeln;
    end;
    readln;
END.