Зображену діаграму побудовано за таким діапазоном клітинок​

что-то типа этого:

const n = 8, m = 8;
var
    A: array[1..n,1..m] of integer;
    i, j, sum, product: integer;
    isnotnull: boolean;
begin
    randomize;
    writeln('Случайная матрица:');    for i:=1 to n do begin
        for j:=1 to m do begin
            A[i,j] := random(51) - 25;
            write(A[i,j]:5);
        end;
        writeln;
    end;    sum := 0;
    for i:=1 to n do
        if A[i,n-i+1] < 0 then
            sum := sum + A[i,n-i+1];
    writeln('Сумма отрицательных элементов побочной диагонали = ', sum);    product := 1;
    isnotnull := False;
    for i:=1 to n-1 do
        for j:=2 to n do
            if (j > i) and (A[i,j] <> 0) then begin
                isnotnull := True;
                product := product * A[i,j];
            end;
    if isnotnull
        writeln('Произведение ненулевых элементов в области выше главной диагонали = ', product)
    else
        writeln('Ненулевых элементов в области выше главной диагонали нет.', product);    readln;
end.

Sub z()
    Range(Cells(1, 1), Cells(100, 100)).Clear
    n1=5
    n2=3
    n3=4

    k = 1
    A = InitMatrix(n1, k, 1)
    k = k + n1 + 2
    B = InitMatrix(n2, k, 1)
    k = k + n2 + 2
    C = InitMatrix(n3, k, 1)

End Sub

Function InitMatrix(n, cx, cy)
    ReDim A(n, n)
    For i = LBound(A, 1) To UBound(A, 1)
        For j = LBound(A, 2) To UBound(A, 2)
            A(i, j) = Rnd * 200 - 100
            Cells(cx + i, cy + j) = A(i, j)
        Next
    Next

    Cells(cx, cy + n + 1) = "PositiveAverage ="
    Cells(cx, cy + n + 2) = PositiveAverage(A)

    InitMatrix = A
End Function

Function PositiveAverage(A)
    For i = LBound(A, 1) To UBound(A, 1)
        For j = LBound(A, 2) To UBound(A, 2)
            If A(i, j) > 0 Then
                s = s + A(i, j)
                k = k + 1
            End If
        Next
    Next

    PositiveAverage = s / k
End Function