Возведение матрицы в степень.

const
  n = 5;
var
  Matr: Array[1..n, 1..n] of Real;
  Arr: Array[1..n] of Real;
  Sled: Real;
  m, i, j, k: Byte;
begin
  Randomize;
  k:=1;
  writeln('Матрица: ');
  for i:=1 to n do
  begin
    for j:=1 to n do
    begin
      Matr[i, j]:=Random(3)+1;
      write(Matr[i, j]:4);
      if (i = j) then  
        Arr[i]:=Matr[i, j];
    end;
    writeln;
  end;
  write('N: ');
  readln(m);
  for k:=1 to m do
  begin
    Sled:=0;
    write('След матрицы в ', k, ' степени: ');
    for i:=1 to n do
    begin  
      Sled:=Sled + Power(Arr[i], k);
      write(Power(Arr[i], k), ' ');
    end;
    writeln('(', Sled:0:0, ');');
    writeln;
    end;
end.

uses crt;
const nmax=10;
type matr=array[1..nmax,1..nmax] of real;
var m:byte;
procedure Umn(a1,a2:matr;var a3:matr);
var k,j,i:byte;
begin
for k:=1 to m do
for j:=1 to m do
   begin
     a3[k,j]:=0;
     for i:=1 to m do
     a3[k,j]:=a3[k,j]+a1[k,i]*a2[i,j];
   end;
end;
function Sled(a:matr):real;
var i:byte;
    s:real;
begin
s:=0;
for i:=1 to m do
s:=s+a[i,i];
Sled:=s;
end;
procedure Vyvod(a:matr;s:string);
var i,j:byte;
begin
writeln(s);
for i:=1 to m do
 begin
  for j:=1 to m do
  write(a[i,j]:6:2);
  writeln;
 end;
writeln('След матрицы=',sled(a):0:2);
writeln;
write('Нажмите Enter');
readln;
clrscr
end;
 
var a,b,c:matr;
    n,i,j:byte;
    s:string;
begin
randomize;
repeat
write('Размер матрицы до ',nmax,' m=');
read(m);
until m in [1..nmax];
for i:=1 to m do
for j:=1 to m do
a[i,j]:=-2+5*random;
repeat
write('Введите натуральное число от 2 до 7 - степень матрицы n=');
read(n);
until n in[2..7];
Vyvod(a,'Матрица A:');
c:=a;
for i:=1 to n-1 do
 begin
    Umn(c,a,b);
    str(i+1,s);
    Vyvod(b,'Матрица A^'+s+':');
    c:=b;
   end;
end.