Табулирование функции 2-х переменных

program one ;
uses CRT;
var
    xn,xk,yn,yk,dx,dy,a,x,y,f : real ;
 
begin
  CLRSCR;
      writeln('Vvedite a');
      readln(a);
      writeln('Vvedite Xn i Xk');
      readln(xn,xk);
      writeln('Vvedite Yn i Yk');
      readln(yn,yk);
      writeln('VVedite dX i dY');
      readln(dx,dy);
x := xn;
y := yn;
   while x < xk do
begin
x:=x+dx;
end;
writeln;
  while y < yk do
begin
x:=xn;
writeln('y=', y:1:1);
  while x<xk do
  begin
f:=(Ln(sin(abs(x*y-a))/Cos(abs(x*y-a)))/ln(10));
write('x= ',x:1:1,'; f= ', f:2:3,';');
          x := x +dx;
        end;
        writeln;
      y:= y +dy;
    end;
  readln;
end.

program one ;
uses crt;
var xn,xk,yn,yk,dx,dy,a,x,y,f:real;
begin
clrscr;
write('Vvedite a=');
readln(a);
writeln('Vvedite Xn i Xk');
readln(xn,xk);
writeln('Vvedite Yn i Yk');
readln(yn,yk);
writeln('VVedite dX i dY');
readln(dx,dy);
x := xn;
while x <= xk do
 begin
  y := yn;
  while y<=yk do
   begin
    if (cos(x*y-a)=0)or(sin(abs(x*y-a))/cos(abs(x*y-a))<=0)then
    writeln('x=',x:6:2,'  y=',y:6:2,' решения нет')
    else
     begin
      f:=ln(sin(abs(x*y-a))/cos(abs(x*y-a)))/ln(10);
      writeln('x=',x:6:2,'  y=',y:6:2,' f=',f:8:3);
     end;
    y:=y+dy;
   end;
  x:=x+dx;
 end;
readln
end.