Вычислить след матрицы

uses crt;
 
const
  nmax = 4;
 
type
  TMatrix = array [1..nmax, 1..nmax] of real;
  
procedure EnterMatrix(var A : TMatrix; n : integer; NameMatrix : char);
var
  i, j : integer;
begin
  Writeln('Введите элементы массива ', NameMatrix, ': ');
  for i := 1 to n do
    for j := 1 to n do
      begin
        Write(NameMatrix, '[', i, ';', j, '] = ');
        Readln(A[i, j]);
      end;
end;
 
function  trMatrix(A : TMatrix; n : integer) : real;
var
  i : integer;
  tr : real;
begin
  tr := 0;
  for i := 1 to n do
    tr := tr + A[i, i];
  trMatrix := tr;  
end;
 
procedure TabFunction (xn, xk, h, c, d : real);
var
  x, y : real;
begin
  x := xn;
  while x <= xk do
    begin
      y := c*SQR(x) + d;
      Writeln('x = ', x:0:4, ' y = ', y:0:4);
      x := x + h;
    end;
end;
 
 
var
  A, B : TMatrix;
  xn, xk, h, c, d : real;
begin
  EnterMatrix(A, 3, 'A');
  EnterMatrix(B, 4, 'B');
  c := trMatrix(A, 3);
  d := trMatrix(B, 4);
  Writeln('След матрица А = ', c:0:4);
  Writeln('След матрицы B = ', d:0:4);
  Write('X начальное = ');
  Readln(xn);
  Write('X конечное = ');
  Readln(xk);
  Write('Шаг = ');
  Readln(h);
  TabFunction(xn, xk, h, c, d);
  Readln;
end.