Протабулировать функцию - Pascal (81280)

Формулировка задачи:

Народ, помогите сделать к завтрашнему утру. Вычислить значение функции, если X меняется от X(начальное) к X(конечное) с шагом ΔX. Результат вычислений вывести в виде таблицы. В случае если вычисления функции невозможно, выдать соответствующее сообщение. Вычислить наименьшее значение функции на этом интервале. Входные данные задать самостоятельно. Написать программу, привести результаты работы программы.

Решение задачи: «Протабулировать функцию»

Листинг программы

const goriz = '+------------------------+------------------------+';
var x, xn, xk, dx, a, xi, y, min: double;
    fmin: boolean;
begin
  write('Xнач = ');
  readln(xn);
  write('Xкон = ');
  readln(xk);
  if xn > xk
    then begin
      x := xn;
      xn := xk;
      xk := x
    end;
  repeat
    write('|dX| > 0;  dX = ');
    readln(dx);
    if dx < 0 then dx := -dx
  until dx > 0;
  write('A = ');
  readln(a);
  write('Ksi = ');
  readln(xi);
  fmin := false;
  x := xn;
  writeln(goriz);
  writeln('|', 'x':12, ' ':12, '|', 'y':12, ' ':12, '|');
  writeln(goriz);
  while x < xk + dx / 2 do
    begin
      write('| ', x, ' | ');
      if (x = 0) or (xi <= 0)
        then write('Функция не определена':22)
        else begin
          y := exp(a * x) + sqrt(sqrt(abs(a * x + xi * xi * xi))) + ln(xi) / 2 / x;
          if not fmin or (y < min) then min := y;
          fmin := true;
          write(y)
        end;
      writeln(' |');
      x := x + dx
    end;
  writeln(goriz);
  if fmin
    then writeln('Минимальное значение: ', min)
    else writeln('Минимальное значение не найдено: на всём интервале функция не определена.');
  readln
end.

Объяснение кода листинга программы

Создается константа goriz для форматирования вывода.
Объявляются переменные: x, xn, xk, dx, a, xi, y, min — все с типом double.
Переменная fmin устанавливается в значение false.
Значение переменной xn присваивается переменной x.
Выводится форматированная строка с константой goriz.
Читаются значения переменных xn и xk с помощью функции readln.
Если xn больше xk, то значения переменных x, xn и xk обновляются.
Запускается цикл while, который выполняется до тех пор, пока x меньше xk + dx/2.
Внутри цикла выводится форматированная строка с константой goriz.
Вычисляется значение y по формуле y = exp(ax) + sqrt(sqrt(abs(ax + xixixi))) + ln(xi)/2/x.
Если не выполнено условие (y < min) или (y == min), то min присваивается значение y.
Устанавливается значение переменной fmin в true.
Выводится значение y.
Выводится форматированная строка с константой goriz.
Если выполнено условие (fmin == true), то выводится минимальное значение min.
Если условие не выполнено, то выводится сообщение о том, что минимальное значение не найдено.
Выполняется чтение следующего значения с помощью функции readln.
Цикл while повторяется.
Выполняется чтение следующего значения с помощью функции readln.
Цикл while повторяется.