Составить подрограмму вычислений k-ой степени матрицы из (n x n) элементов, используя подпрограмму умножения двух матриц

uses crt;
const na=3;
      nb=4;
type matr=array[1..nb,1..nb] of real;
{умножение матрицы на матрицу}
procedure Umn(a1,a2:matr;p:integer;var a3:matr);
var k,j,i:byte;
begin
for k:=1 to p  do
for j:=1 to p do
 begin
  a3[k,j]:=0;
  for i:=1 to p do
  a3[k,j]:=a3[k,j]+a1[k,i]*a2[i,j];
 end;
end;
{вывод на экран}
procedure Vyvod(a:matr;p:integer;st:string);
var i,j:byte;
begin
writeln(st);
for i:=1 to p do
 begin
  for j:=1 to p do
  write(a[i,j]:9:1);
  writeln;
 end;
end;
procedure stepen(a:matr;p,k:integer;st,s:string);
var a1,a2:matr;
    i,j:integer;
begin
for i:=1 to p do
for j:=1 to p do
a[i,j]:=-2+5*random;
vyvod(a,p,'Исходная матрица '+st);
a1:=a;
for i:=1 to k-1 do{k-1 раз}
 begin
  Umn(a1,a,p,a2);{умножаем A*A}
  a1:=a2;{запомним результат}
 end;
Vyvod(a2,p,'Матрица '+st+'^'+s+':');
end;
var a,b:matr;
    k:byte;
    s:string;
begin
clrscr;
randomize;
repeat
write('Введите натуральное число от 2 до 7 - степень матриц k=');
readln(k);
until k in[2..7];
str(k,s);
stepen(a,na,k,'A',s);
stepen(b,nb,k,'B',s);
readln
end.