Составить программу для вычисления значения функции

program tab;
uses
    crt;
var
   a: real;
   b: real;
   h: real;
   y: real;
   procedure swap(var a: real; var b: real);
   var
      buf: real;
   begin
        buf := a;
        a := b;
        b := buf;
   end;
   function ctg(a: real): real;
   begin
        ctg := cos(a) / sin(a);
   end;
begin
     clrscr;
     write('Enter a: ');
     readln(a);
     write('Enter b: ');
     readln(b);
     write('Enter step: ');
     readln(h);
     if (a > b) then
        swap(a, b);
     writeln('+------------+------------+');
     writeln('|     x      |     y      |');
     writeln('+------------+------------+');
     while (a <= b) do
     begin
          y := 1 / 2 * ctg(a) / 4 + 4;
          writeln('|', a:11:4, ' |', y:11:4, ' |');
          a := a + h;
     end;
     writeln('+------------+------------+');
     readkey;
end.