Алгоритмы циклической структуры

program ufunc;
 
const xs = 1; xe = 3; hx = 1.3;
      ys = 2; ye = 4; hy = 1.3;
      zs = 0.5; ze = 1.5; hz = 0.8;
 
var u, x, y, z, umax, xmax, ymax, zmax: real;
    xn, yn, zn, xc, yc, zc: integer;
 
begin
  xc := trunc((xe - xs) / hx); yc := trunc((ye - ys) / hy); zc := trunc((ze - zs) / hz);
  writeln('Tabulation of U=-sqrt(y^2+(4x^2)/3)+(5cos^3(xy))/(xz)');
  writeln('   z    y    x         U');
  for zn := 0 to zc do
    begin
      z := zs + zn * hz;
      for yn := 0 to yc do
        begin
          y := ys + yn * hy;
          for xn := 0 to xc do
            begin
              x := xs + xn * hx;
              write('  ', z:3:1, ', ', y:3:1, ', ', x:3:1);
              if x * z = 0
                then writeln('<indeterminacy>')
                else
                  begin
                    u := 5 * cos(x * y) * cos(x * y) * cos(x * y) / x / z - sqrt(y * y + 4 * x * x / 3);
                    writeln('    ', u:10:5);
                    if (zn + yn + xn = 0) or (u > umax)
                      then
                        begin
                          xmax := x; ymax := y; zmax := z; umax := u
                        end
                 end
            end
        end
    end;
  writeln('Maximum of U is');
  writeln('U(', zmax:3:1, ', ', ymax:3:1, ', ', xmax:3:1, ') = ', umax:10:5);
  readln
end.