Вычисление таблицы значений функции u(x, y) при изменении значений аргументов

program Project1;
var
  upr: real; {Промежуточная переменная для хранения значения вычислений }
  xn,yn: real; {Переменные для хранения значения левой границы вычислений}
  xk,yk: real; {Переменные для хранения значения правой границы вычислений}
  x,y: real;  {Переменные для хранения значения аргумента}
  hx,hy: real; {Переменные для хранения значения шага вычислений}
  u: real; {переменная для хранения значения результата вычислений}
begin
  {Заданные значения аргументов и границ}
  xn:= 2;
  xk:= 2.5;
  hx:= 0.5;
  yn:= -0.5;
  yk:= 0.3;
  hy:= 0.4;
  {-----------------------------------------------}
  x:= xn;        {Присвоение аргументу X значения левой границы вычислений }
  while x <= xk + hx/2 do    {цикл с условием}
    begin
      y:= yn;        {Присвоение аргументу Y значения левой границы вычислений }
      while y <= yk + hy/2 do    {цикл с условием}
        begin
          if (x*x*y <= -1) then
            begin
              u:= sin(x + y) + x;    {Вычисление первой формулы}
              writeln('x = ',x:7:4, ';  y = ', y:7:4, ';  u(x,y) = ',u:7:4);
            end
          else
            if (x*x*y >= 1) then
              begin
                u:= sin(y);    {Вычисление второй формулы}
                writeln('x = ',x:7:4, ';  y = ', y:7:4, ';  u(x,y) = ',u:7:4);
              end
            else
              begin
                upr:= (x*x*x*x*x*x*x + y*y*y);
                if upr > 0 then    {Проверка аргументов на область определения}
                  begin
                    u:= 4*x*x + ln(upr);  {Вычисление третьей формулы}
                    writeln('x = ',x:7:4, ';  y = ', y:7:4, ';  u(x,y) = ',u:7:4);
                  end
                else
                  writeln('x = ',x:7:4, ';  y = ', y:7:4, ';  функция не определена');
              end;
            y:= y + hy;         {Приращение аргумента Y}
        end;
      x:= x + hx;         {Приращение аргумента X}
    end;
  readln;
end.