Составить программу, которая вычисляет сумму тех элементов двумерного массива, номера строк и столбцов которых принадлежат, соответственно, непустым множествам S1 и S2 — Pascal(Паскаль)

uses crt;
const nmax=20;//максимальный размер матрицы
var a:array[1..nmax,1..nmax] of integer;//тип матрицы
    n,m,i,j:byte;//переменные - размеры матрицы, счетчики циклов
    sm,b:integer;//сумма, переменная для чтения чисел
    s1,s2:set of byte;//множества чисел от 0 до 255
begin
clrscr;//очищаем экран
randomize;//генерация псевдо-случайных чисел
repeat
write('Количество строк=');
readln(m);
until m in [1..nmax];//проверяем правильность ввода
repeat
write('Количество столбцов=');
readln(n);
until n in [1..nmax];
writeln('Исходная матрица:');
for i:=1 to m do//выводим матрицу по строкам
 begin
   for j:=1 to n do
    begin
     a[i,j]:=random(20);//случайные числа от 0 до 19
     write(a[i,j]:4);
    end;
   writeln;
 end;
writeln;
writeln('Множество S1:');
s1:=[];//пустое множество
for i:=1 to m do
 begin
  b:=random(m)+1;//наполняем случайными числами от одного до количества строк
  s1:=s1+[b];//записываем в множество
 end;
for i:=1 to m do
if i in s1 then write(i,' ');//выводим элементы множества на экран
writeln;
writeln('Множество S2:');
s2:=[];
for i:=1 to n do
 begin
  b:=random(n)+1;
  s2:=s2+[b];
 end;
for i:=1 to n do
if i in s2 then write(i,' ');
writeln;
writeln;
sm:=0;//сумма пока ноль
for i:=1 to m do//идем по матрице
for j:=1 to n do
if(i in s1)and(j in s2) then sm:=sm+a[i,j];{если номер строки и номер столбца есть 
в соответствующих им множествам, суммируем}
write('Сумма=',sm);//выводим сумму на экран
readln
end.