Написать подпрограмму умножение двух матриц (с пояснением)

program Matrix_Umnoj; uses Crt;
  const n = 5;//задает максимальное количество строк и столбцов для всех матриц в проге
type
   nm = 1..n;
   my = array[nm,nm]of integer;
   mmy = array[nm,nm]of integer;
   myy = array[nm,nm]of integer;
var m: my;mm: mmy;//m,mm переменные для умножаемых матриц
   nn: myy;// nn переменная для матрицы с результатами умножения(результирующей)
   n1,m1,n2,m2 : byte;//переменные для ввода кол-ва строк и столбцов матриц m и mm(тех которые будем умножать)
   f: Boolean;
//procedure LongMatrix проверка ввода строк и столбов
//матриц m и mm чтобы в пределах const n =
procedure LongMatrix(var n1,m1,n2,m2: byte);
begin
  repeat
    write('Введите число строк первой матрицы не более ',n,' -> ');
    readln(n1);
  until (0 < n1) and (n1 <= 10);
    writeln;//пустая строка на экране между правильно введенными данными
  repeat
    write('Введите число столбцов первой матрицы не более ',n,' -> ');
    readln(m1);
  until (0 < m1) and (m1 <= n);
     writeln;
   repeat
     write('Введите число строк второй матрицы не более ',n,' -> ');
     readln(n2);
   until (0 < n2) and (n2 <= 10);
      writeln;
   repeat
     write('Введите число столбцов второй матрицы не более ',n,' -> ');
     readln(m2);
   until (0 < m2) and (m2 <= n);
     writeln;
end;//LongMatrix
//function Simile равно ли кол-во столбов m
//кол-ву строк mm если нет,то умножать нельзя и f = false
function Simile(m1,n2: byte):Boolean;
begin
  Simile:= m1 = n2;
end;//Simile
//заполняем матрицы m и mm случайными числами
procedure zap(var m: my;var mm: mmy);
  var i,j: byte;
begin randomize;
  for i:=1 to n1 do
   for j:=1 to m1 do m[i,j]:= random(10);
  for i:=1 to n2 do
   for j:=1 to m2 do mm[i,j]:= random(10);
end;//zap
//выводим обе заполненные матрицы(m и mm) на экран
procedure output(m: my;c,d: byte);
  var i,j: byte;
begin
  for i:=1 to c do begin
   for j:=1 to d do  write(m[i,j],' ':3);
    writeln;end;//for i
end;//output
procedure output(mm: mmy;c,d: byte);
  var i,j: byte;
begin
  for i:=1 to c do begin
   for j:=1 to d do  write(mm[i,j],' ':3);
    writeln;end;//for i
end;//output
//procedure Action умножаем m * mm и
//результат заносим в матрицу nn
procedure Action(var m: my;var mm: mmy;
var nn: myy);
 var i,j,l: byte;
begin
  for i:=1 to n1 do //в nn строки = n1
   for j:=1 to m2 do //в nn кол-во столбцов = m2
    for l:=1 to m2 do
{вначале умножаем первое число в строке 1 матрицы
на каждое значение в 1-й строке 2-й матрицы и заносим в nn}
     if j = 1 then nn[i,l]:= m[i,j]*mm[j,l] else
{2 и т.д.числа в строке 1-й матрицы на числа во 2-й и т.д.
строках 2-й матрицы и плюсуем их к обрабатываемой строке nn[i,}
      nn[i,l]:= nn[i,l] + m[i,j]*mm[j,l];//проверить nn[j,l +
end;//Action
//выводим матрицу nn на экран
procedure output(nn: myy;c,d: byte);
 var i, j : byte;
begin
  for i:=1 to c do begin
   for j:=1 to d do
     write(nn[i,j],' ':3);
     writeln;
  end;//for i
end;//output
begin ClrScr;//ОСНОВНОЕ ТЕЛО ПРОГРАММЫ
//procedure LongMatrix ввод и проверка строк и столбов умножаемых массивов
LongMatrix(n1,m1,n2,m2);
f:= Simile(m1,n2); {кол-во столбов Массива1(m
  должно быть = кол-ву строк Массива2 mm}
if f then begin
 zap(m,mm);//если f = true заполняем массивы
 output(m,n1,m1);writeln;//выводим обе
 output(mm,n2,m2);writeln;//матрицы
   FillChar(nn,SizeOf(nn),0);//инициация массива nn нулями
 Action(m,mm,nn);//здесь умножаем m * mm
 //результаты заносим в массив nn
  writeln('вывод матрицы nn');
  writeln;
 output(nn,n1,m2);end//вывод результата умножения
//матриц m * mm
  else
 writeln('Матрицы умножать нельзя');
 writeln('нажмите любую клавишу');
 readkey;
end.//begin
{НАГЛЯДНОЕ ПОСОБИЕ УМНОЖЕНИЯ МАТРИЦ
 1 2             А B C
 3 4  уножаем на D E F
 5 6
 вначале умножаем 1 на каждое значение
 первой строки
 второй матрицы и заносим эти значения
 в результирующую матрицу
 
 1*A 1*B 1*C
 
 затем умножаем 2 на каждое значение
 второй строки второй матрицы и плюсуем
 этот результат в первой строке
 
 1*A+2*D 1*B+2*E 1*C+2*F -это полученные
 
 значения первой строки в результирующей
 матрице.Первая строка результирующей
 матрицы заполнена.Внешний цикл становится
 равным цифре 2 - это переход ко 2 строке
 первой матрицы.
 Теперь первое значение во 2-й строке
 первой матрицы 3 ,умножаем на каждое значение
 в первой строке 2-й матрицы то есть ,те же
 действия что и вначале ,только теперь работает
 2-я строка 1-й матрицы и снова все значения
 2-й матрицы...итак
 
 3*A 3*B 3*C
 
 далее
 
 3*A+4*D 3*B+4*E 3*C+4*F
 
 получились 2 строки результирующей матрицы
 
  1*A+2*D 1*B+2*E 1*C+2*F
  3*A+4*D 3*B+4*E 3*C+4*F
 
  те же действия с третьей строкой 1-й матрицы и
  результат
 
  1*A+2*D 1*B+2*E 1*C+2*F
  3*A+4*D 3*B+4*E 3*C+4*F
  5*A+6*D 5*B+6*E 5*C+6*F}