Cделать программу в pascal по : "найти в двумерном массиве сумму массива,сумму элементов по столбцам и среднее арифметическое по строкам"

Пусть дан массив real arr[1..n][1..n1]
заведём массивы real s1[1..n], s2[1..n1]

n - кол - во строк
n1 - кол -во столбцов

Основная часть:

...
for i:= 1 to n do
s1[i] := 0;

for i:= 1 to n1 do
s2[i] := 0;

for i:=1 to n do begin
for j:=1 to n1 do begin
s1[i] := s1[i] + arr[i][j];
s2[j] := s2[j] + arr[i][j];
end;
s1[i] := s1[i] / n
end;

...

s1[i] - среднее арифметическое для i строки
s2[i] - сумма элементов для i столбца