Работа с матрицами (pascal).

[Ethem]

Здравствуйте!
У меня такое задание даны 2 матрицы A и B, надо вычислить матрицу C:
C=-B^3+A;

program program1;
 type Matr=array[1..4,1..4] of integer;
  var A,B,C: Matr;
Procedure Matrixini(var M:Matr);
 var i,j:integer;
begin
 randomize;
for i:=1 to 4 do
for j:=1 to 4 do
M[i,j]:=random (20)-10;
end;
Procedure Matrixtable(var M:Matr);
var i,j:integer;
begin
for i:=1 to 4 do
begin
for j:=1 to 4 do
Write(M[i,j],' ');
Writeln();
end;
end;
procedure Matrixmulti(A,B:Matr;var M:Matr);
var i,j,k:integer;
var S:integer;
begin
for i:=1 to 4 do
for j:=1 to 4 do
begin
S:=0;
for k:=1 to 4 do
S:=S+A[i,k]*B[k,j];
end;
M[i,j]:=S;
end;
begin
randomize;
Matrixini(A);
Writeln('Вывод матрицы A:');
Matrixtable(A);
Matrixini(B);
Writeln('Вывод матрицы B:');
Matrixtable(B);
Matrixini(C);
Writeln('Вывод матрицы C:');
Matrixtable(C);
Matrixmulti(A,B,C);
end.

Вот у меня есть код для перемножения матриц, я знаю, чтобы возвести в степень надо перемножить матрицу, а как организовать это в паскале не знаю.