Выполнение операций сложения, вычитания и умножения над матрицами целых чисел

uses crt;
type matrix=array[1..10,1..10]of integer;
var a,b,c:matrix;
    choose,x,y,z:integer;
 
procedure Zapolnenie(var m:matrix; x,y:integer);
var i,j:integer;
begin
 for i:=1 to x do
 for j:=1 to y do
 m[i,j]:=random(10);
end;
 
procedure Vyvod(m:matrix; x,y:integer);
var i,j:integer;
begin
 for i:=1 to x do
 begin
  for j:=1 to y do 
  write(m[i,j]:4);
  writeln;
 end;
end;
 
procedure Peremnozhenie(a:matrix; x,y:integer; b:matrix; z:integer; var c:matrix);
var i,j,k:integer;
    s:integer;
begin
 for i:=1 to x do
 for j:=1 to z do
 begin
  S:=0;
  for k:=1 to y do
  S:=S+a[i,k]*b[k,j];
  c[i,j]:=S;
 end;
end;
 
procedure slozhenie(x,y:integer; a,b:matrix; var c:matrix);
var i,j:integer;
begin
 for i:=1 to x do for j:=1 to y do c[i,j]:=a[i,j]+b[i,j];
end;
 
procedure vychitanie(x,y:integer; a,b:matrix; var c:matrix);
var i,j:integer;
begin
 for i:=1 to x do for j:=1 to y do c[i,j]:=a[i,j]-b[i,j];
end;
 
begin
 clrscr;
 write('Операция: 1-сложение, 2-вычитание, 3-умножение ');
 readln(choose);
 case choose of
 1:begin
   write('Введите кол-во строк матриц:');
   readln(x);
   write('Введите кол-во столбцов матриц:');
   readln(y);
   zapolnenie(a,x,y);
   zapolnenie(b,x,y);
   writeln('Матрица A');
   vyvod(a,x,y);
   writeln('Матрица B');
   vyvod(b,x,y);
   slozhenie(x,y,a,b,c);
   writeln('Матрица C');
   vyvod(c,x,y);
   end;
  2:begin
   write('Введите кол-во строк матриц:');
   readln(x);
   write('Введите кол-во столбцов матриц:');
   readln(y);
   zapolnenie(a,x,y);
   zapolnenie(b,x,y);
   writeln('Матрица A');
   vyvod(a,x,y);
   writeln('Матрица B');
   vyvod(b,x,y);
   vychitanie(x,y,a,b,c);
   writeln('Матрица C');
   vyvod(c,x,y);
   end;    
  3:begin
   write('Введите кол-во строк матрицы A:');
   readln(x);
   write('Введите кол-во столбцов матрицы A:');
   readln(y);
   Zapolnenie(a,x,y);
   Vyvod(a,x,y);
   write('Введите количество столбцов матрицы B:');
   readln(z);
   Zapolnenie(b,y,z);
   Vyvod(b,y,z);
   Peremnozhenie(a,x,y,b,z,c);
   writeln('Матрица C ');
   Vyvod(c,x,z);
   end;
  end;
  readln;
end.

uses crt;
const nmax=20;
var n:byte;
type Tmatr=array[1..nmax,1..nmax] of integer;
procedure Vvod(var mt:Tmatr;c:string);
var i,j:byte;
begin
writeln('Введите ',n*n,' элементов матрицы ',c);
for i:=1 to n do
for j:=1 to n do
 begin
  write('el[',i,',',j,']=');
  readln(mt[i,j]);
 end;
clrscr;
end;
 
procedure Vyvod(var mt:Tmatr;c:string);
var i,j:byte;
begin
writeln(c,':');
for i:=1 to n do
  begin
    for j:=1 to n do
    write(mt[i,j]:5);
    writeln;
  end;
writeln;
end;
 
procedure UmnMatr(m1,m2:Tmatr;var m:Tmatr);
var i,j,k:byte;
begin
for k:=1 to n do
for j:=1 to n do
  begin
   m[k,j]:=0;
   for i:=1 to n do
   m[k,j]:=m[k,j]+m1[k,i]*m2[i,j];
  end;
end;
 
procedure SumMatr(m1,m2:Tmatr;var m:Tmatr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
m[i,j]:=m1[i,j]+m2[i,j];
end;
 
procedure RazMatr(m1,m2:Tmatr;var m:Tmatr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
m[i,j]:=m1[i,j]-m2[i,j];
end;
 
var a,b,c:Tmatr;
    w:char;
begin
clrscr;
randomize;
write('Введите размер матриц n=');
readln(n);
Vvod(a,'A');
Vvod(b,'B');
Vyvod(a,'Matrica A');
Vyvod(b,'Matrica B');
write('Press Enter...');
readln;
repeat
clrscr;
Vyvod(a,'Matrica A');
Vyvod(b,'Matrica B');
writeln('Выберите действие');
writeln('1-сложение');
writeln('2-вычитание');
writeln('3-умножениe');
writeln('другое-выход');
write('w=');readln(w);
case w of
'1':begin
    SumMatr(a,b,c);
    Vyvod(c,'Матрица C=A+B');
    write('Press Enter...');
    readln
    end;
'2':begin
    RazMatr(a,b,c);
    Vyvod(c,'Матрица C=A-B');
    write('Press Enter...');
    readln
    end;
'3':begin
    UmnMatr(a,b,c);
    Vyvod(c,'Матрица C=A*B');
    write('Press Enter...');
    readln
    end;
else exit;
end;
until not(w in ['1'..'3']);
end.