Вычислить матрицу по формуле

type Tmatr=array[1..nmax,1..nmax] of real;
procedure Vvod(var mt:Tmatr;k:byte;c:string);
{параметры - матрица, размер и имя}
var i,j:byte;
begin
writeln('Матрица ',c);
writeln('Введите элементы матрицы, вещественные числа,');
writeln('ниже главной диагонали и на ней');
for i:=1 to k do
for j:=1 to i do
 begin
  write('el[',i,',',j,']=');
  readln(mt[i,j]);
 end;
{симметрично заполняем верх матрицы}
for i:=1 to k-1 do
for j:=i+1 to k do
mt[i,j]:=mt[j,i];
clrscr;
end;
 
procedure Vyvod(var mt:Tmatr;k,p:byte;c:string);
{параметры - матрица, размер, количество позиций для вывода чисел, имя}
var i,j:byte;
begin
writeln(c,':');
for i:=1 to k do
  begin
    for j:=1 to k do
    write(mt[i,j]:p:2);
    writeln;
  end;
writeln;
end;
 
procedure UmnMatr(m1,m2:Tmatr;k:byte;var m:Tmatr);{умножение 2х матриц}
{параметры - матрицы, которые перемножаем, размер и полученная матрица}
var i,j,p:byte;
begin
for p:=1 to k do
for j:=1 to k do
  begin
   m[p,j]:=0;
   for i:=1 to k do
   m[p,j]:=m[p,j]+m1[p,i]*m2[i,j];
  end;
end;
 
var a,b,c,a2,b2,b3,c2,c4,x,y:Tmatr;
    k:byte;
begin
clrscr;
randomize;
repeat
write('Размер матриц до ',nmax,' k=');
readln(k);
until k in [1..nmax];
Vvod(a,k,'A');
Vvod(b,k,'B');
Vvod(c,k,'C');
Vyvod(a,k,6,'Matrica A');
Vyvod(b,k,6,'Matrica B');
Vyvod(c,k,6,'Matrica C');
write('Press Enter...');
readln;
UmnMatr(a,a,k,a2);
Vyvod(a2,k,7,'Matrica A2=A^2');
UmnMatr(b,b,k,b2);
UmnMatr(b,b2,k,b3);
Vyvod(b3,k,8,'Matrica B3=B^3');
UmnMatr(c,c,k,c2);
UmnMatr(c2,c2,k,c4);
Vyvod(c4,k,9,'Matrica C4=C^4');
write('Press Enter...');
UmnMatr(a2,b3,k,x);
UmnMatr(x,c4,k,y);
Vyvod(y,k,12,'Matrica Y=A^2*B^3*C^4');
write('Press Enter...');
readln
end.