Умножение, сложение и вычитание матриц - Turbo Pascal

uses crt;
const n=3;
type mtr=array[1..n,1..n] of integer;
const a:mtr=(( 1, 2, 4),
             ( 1, 0, 5),
             (-1, 3, 2));
      b:mtr=(( 1, 3, 1),
             ( 2, 1, 0),
             ( 1,-1, 1));
      c:mtr=(( 1, 1, 2),
             ( 0, 1, 1),
             ( 2,-1, 0));
var d,f:mtr;
    i,j,k:byte;
begin
clrscr;
{умножим А*В}
for i:=1 to n do
for j:=1 to n do
 begin
  d[i,j]:=0;
  for k:=1 to n do
  d[i,j]:=d[i,j]+a[i,k]*b[k,j];
 end;
{A*B+5*C-A}
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    f[i,j]:=d[i,j]+5*c[i,j]-a[i,j];
    write(f[i,j]:3);
   end;
  writeln;
 end;
readln
end.