Найти матрицу k*A(в квадрате)+mE

uses crt;
const nmax=20;
type matr=array[1..nmax,1..nmax] of integer;
procedure Print(var mt:matr;n:byte;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]:6);
  writeln;
 end;
writeln;
end;
 
var a,b,e:matr;
    n,i,j,p:byte;
    k,m:integer;
begin
clrscr;
randomize;
repeat
write('Размер матрицы=');
readln(n);
until n in [1..nmax];
for i:=1 to n do
for j:=1 to n do
a[i,j]:=random(10)+1;
Print(a,n,'A:');
write('k=');readln(k);
write('m=');readln(m);
writeln;
{произведение А*А}
for p:=1 to n do
for j:=1 to n do
 begin
  b[p,j]:=0;
  for i:=1 to n do
  b[p,j]:=b[p,j]+a[p,i]*a[i,j];
 end;
Print(b,n,'B=A^2:');
{умножение на число}
for i:=1 to n do
for j:=1 to n do
b[i,j]:=b[i,j]*k;
Print(b,n,'B=k*A^2:');
writeln('Press Enter...');
readln;
{создание и умножение на число единичной матрицы}
for i:=1 to n do
for j:=1 to n do
if j=i then e[i,j]:=m
else e[i,j]:=0;
Print(e,n,'E=m*E:');
{сумма матриц}
for i:=1 to n do
for j:=1 to n do
b[i,j]:=b[i,j]+e[i,j];
Print(b,n,'B=k*A^2+m*E:');
writeln('Press Enter...');
readln
end.