Табулирование функции

program tabuliruemaya_funkciya;
 
const
  N = 10;
 
type
  arr1 = array[1..10] of real;
  arr3 = array[1..10, 1..N] of real;
 
var
  A: arr1;
  Y: arr3;
  b, Da, Xn, Xk, Dx, An, Ak: real;
  err: byte;
  m: integer;
 
 
 
 
procedure InputVariable(var An, Ak, Da, Xn, Xk, Dx: real); //ввод значений 
begin
  writeln('Введите начальное значение аргумента х:');
  write('Xn = ');readln(Xn);
  writeln('Введите конечное значение аргумента х:');
  write('Xk = ');readln(Xk);
  writeln('Введите шаг изменения аргумента х:');
  write('Dx = ');readln(Dx);
  writeln('Введите начальное значение параметра A:');
  write('An = ');readln(An);
  writeln('Введите конечное значение параметра A:');
  write('Ak = ');readln(Ak);
  writeln('Введите шаг изменения параметра A:');
  write('Da = ');readln(Da);
end;
 
 
procedure InA(An, Ak, Da: real; m: integer;  var A: arr1); //ввод массива А
var
  i: integer;
begin
  m := round((Ak - An) / Da);
  for i := 1 to m + 1  do
    A[i] := An + (i - 1) * Da;
end;
 
procedure InB(var B: real); //значение интеграла 
begin
  B := 1;
end;
 
procedure Tab(Ak, An, Dx, Xn, Da: real; m, N: integer; A: arr1; B: real; var Y: arr3); //табулирование функции
var
  i, j: integer;
  X: real;
begin
  if Dx <> 0 then
  begin
    for i := 1 to m + 1  do
      for j := 1 to N do
        if (x = 5 / 2) or (x = 1 / 3) then err := 1 else 
        begin
          X := Xn + (j - 1) * Dx;
          Y[i, j] := A[i] * B * (X * tan(pi * X / 5) + 2 * A[i] * cos(X * X / 2));
        end;
  end else err := 1;
end;
 
 
procedure outY(Y: arr3; A: arr1;  N: integer; Xn, Xk, Dx: real); //вывод двумерного массива значения функции Y
var
  m, i, j: integer;
begin
  m := round((Ak - An) / Da);
  for i := 1 to m do
  begin
    write(A[i]:6:4, ' ');
    for j := 1 to N do
      write(Y[i, j]:8:3);
    writeln();
  end;
end;
 
begin
  InputVariable(An, Ak, Da, Xn, Xk, Dx);
  m := round((Ak - An) / Da);
  InA(An, Ak, Da, m, A);
  InB(B);
  Tab(Ak,  An, Dx, Xn, Da, m, N, A, B, Y);
  outY(Y, A, N, Xn, Xk, Dx); 
  if Err = 1 then
    writeln('Несколько или все значения табулируемой функции не могут быть вычисленны');
end.

program tabuliruemaya_funkciya;