получить матрицу

uses crt;
{$N+}//директива для работы с extended
const max=4;//макс. разм. матриц
type matr=array[1..max,1..max] of extended;
var m,n:byte;
procedure UmnMt(m1,m2:matr;var m3:matr);/умножение матрицы на матрицу
var k,j,i:byte;
begin
for k:=1 to m do
for j:=1 to m do
   begin
     m3[k,j]:=0;
     for i:=1 to m do
     m3[k,j]:=m3[k,j]+m1[k,i]*m2[i,j];
   end;
end;
procedure Power(var a:matr;k:byte);//возведение матрицы в степень
var b,c:matr;
    i:byte;
begin
if k=1 then b:=a
else
 begin
  c:=a;//запомним исходную
  for i:=1 to k-1 do
   begin
    UmnMt(c,a,b);//умножаем ее на себя, запоминаем в b
    c:=b;//снова она исходная
   end;
 end;
a:=b;//запомним последнюю
end;
 
procedure Summ(var m1:matr;m2:matr);//сумма матриц
var i,j:byte;
begin
for i:=1 to m do
for j:=1 to m do
m1[i,j]:=m1[i,j]+m2[i,j];//к первой добавим вторую, первую запомним
end;
var a,b,e:matr;
    i,j:byte;
begin
clrscr;
randomize;
repeat
write('Введите размер матрицы до ',max,' m=');
readln(m);
until m in [1..max];
writeln('Матрица A:');
for i:=1 to m do
 begin
  for j:=1 to m do
   begin
    a[i,j]:=2*random-1;//числа от -1 до +1
    write(a[i,j]:5:1);
   end;
  writeln;
 end;
writeln;
repeat
write('Введите степень матрицы до ',max+1,' n=');//до 5
readln(n);
until n in [1..max+1];
for i:=1 to m do
for j:=1 to m do
if i=j then e[i,j]:=1 else e[i,j]:=0;//создадим единичную
Summ(e,a);//прибавим к ней первую степень
writeln('Матрица E+A:');
for i:=1 to m do
 begin
  for j:=1 to m do
  write(e[i,j]:5:1);
  writeln;
 end;
writeln;
for i:=2 to n do
 begin
  Power(a,i);//возводим в очередную степень
  Summ(e,a);//суммируем
 end;
writeln('Результат:');
for i:=1 to m do
 begin
  for j:=1 to m do
  write(e[i,j]:16:1);
  writeln;
 end;
readln
end.