Сформировать матрицу D из трех других матриц

Procedure Form_D(A,B,C:matrix;Var D:matrix; size:integer);
 Var j,i,k:integer;
 Begin
 For i:=1 to size do
  For j:=1 to size do
   Begin
    D[i,j]:=A[i,j];
    For k:=1 to size do
     D[i,j]:=D[i,j]*C[i,k]-B[k,j];
   End;
End;

Procedure Form_D(A,B,C:matrix;Var D:matrix; size:integer);
Var 
  j,i,k:integer;
Begin
  For i:=1 to size do
    For j:=1 to size do Begin
      D[i,j]:=0;
      For k:=1 to size do D[i,j]:= D[i,j]+A[i,k]*C[k,j];
      D[i,j]:=D[i,j]-B[i,j];
   End;
End;

Procedure MatrixAddWithCoeff(A,B: tMatrix; k: double; var C: tMatrix);
{ Performs C=A+k*B }
var
  i,j: integer;
begin
  for i:=1 to MatrixSize do
    for j:=1 to MatrixSize do
      C[i,j]:= A[i,j]+B[i,j]*k
end;
 
Procedure MatrixMult(A,B: tMatrix; var C: tMatrix);
{ Performs C=A*B }
var
  i,j,k: integer;
  s: tMatrixElement;   // переменная s того же типа, что и элементы матрицы
begin
  for i:=1 to MatrixSize do
    for j:=1 to MatrixSize do begin
      s:=0;
      for k:=1 to MatrixSize do C[i,j]:= s+A[i,k]*B[k,j]
    end
end;
 
// ...
 
begin
  // ... подготавливаем A, B и C
  MatrixMult(A,B,D);
  MatrixAddWIthCoeff(D,B,-1,D);
  // вывод D
  // ...
end.

Procedure Vvod_Form(Var A,B,C:matrix; size:integer);
 Var i,j:integer;
 Begin
   For i:=1 to size do
    For j:=1 to size do
     Begin
      A[i,j]:=(-i/4+1)*j/2;
      B[i,j]:=i*(A[i,j]-j);
      C[i,j]:=-2*A[i,j]+B[i,j];
     End;