Вычислть сумму элементов двумерного массива по столбцам

uses crt;
var A:Array [1..10,1..10] of real;
    N,M,i,j:integer;
    S:array[1..10] of Real;
procedure InputOutput;
begin {описание процедуры ввода-вывода исходных данных}
 clrscr;
 write('Количество строк: ');
 readln(N);
 write('Количество столбцов: ');
 readln(M);
 for i:=1 to N do {Ввод элементов матрицы}
  for j := 1 to M do
   begin 
    write('A[' , i , ', ' , j , ']= ? ');
    readln(A[i, j]) {ввод значения}
   end; WriteLn;
 clrscr;
 writeln('Матрица А: ');
 for i:=1 to N do {Вывод матрицы по строкам}
  begin
   for j := 1 to M do
   write(A[i, j] : 5 : 1); {вывод i-ой строки}
   writeln {перенос курсора на начало следующей строки}
  end;
 writeln
end;
 
procedure SumCol;
begin {описание процедуры вычисления сумм элементов столбцов}
 for j:=1 to M do {цикл по столбцам матрицы}
  begin
   S[j]:=0; {обнуление суммы элементов j-го столбца}
   for i:=1 to n do 
    S[j]:=S[j]+A[i, j] {накопление суммы}
  end;
end;
 
procedure OutResult; {описание процедуры вывода результатов}
begin
 write( 'О т в е т : Суммы элементов столбцов равны ');
 for j:=1 to M do
 write(S[j]:5:1);
 writeLn; 
 readLn
end;
 
begin
  InputOutput; {вызов процедуры ввода-вывода исходных данных }
  SumCol; {вызов процедуры вычисления сумм }
  OutResult; {вызов процедуры вывода результатов}
end.

uses crt;
const n=5; m=4;
var a:array[1..n,1..m]of integer;
i,j,s:integer;
 
begin
 clrscr;
 randomize;
 writeln('Исходный массив');
 for i:=1 to n do
 begin 
  for j:=1 to m do
  begin
    a[i,j]:=random(100);
    write(a[i,j]:3);
  end;
  writeln;
 end;
 
 for j:=1 to m do
 begin
  s:=0;
  for i:=1 to n do s:=s+a[i,j];
  writeln('Столбец ',i,' = ',s);
 end;
 readln;
end.