Получить таблицу значений функции

Program tablica_1;
uses Wincrt;
  Var
    Xn,Xk,Xr,dX,x,y : Real;
    i,N: Integer;
Begin
  WriteLn('Лабораторная работа N 4');
  WriteLn('Программирование циклических алгоритмических структур');
  WriteLn('Цикл с параметром. Задание 1');
    Xn:=10; 
Xk:=100;
 dX:=10; 
x:=Xn;
  N:=trunc((Xk-Xn)/dX)+1;
  For i:=1 to N do
   Begin
     Xr:=x*Pi/180;
     x1:=SQRT(sin(Pi-x)/x);
     y:=arctang(x1/SQRT(1-x1*x1))*arctang(x1/SQRT(1-x1*x1))*arctang(x1/SQRT(1-x1*x1));
 WriteLn( 'при x= ' , X:2:0, 'функция y= ', Y:7:4);
   x:=x+dX:
   End;
End.

Program tablica_1;
uses crt;
  Var
    Xn,Xk,Xr,dX,x,y,x1 : Real;
    i,N: Integer;
Begin
  WriteLn('Лабораторная работа N 4');
  WriteLn('Программирование циклических алгоритмических структур');
  WriteLn('Цикл с параметром. Задание 1');
  Xn:=10;
  Xk:=100;
 dX:=10;
 x:=Xn;
  N:=trunc((Xk-Xn)/dX)+1;
  For i:=1 to N do
   Begin
     Xr:=x*Pi/180;
     if sin(Pi-x)/x<0 then WriteLn( 'при x= ' , X:2:0, 'функция y= не определена')
     else
      begin
       x1:=SQRT(sin(Pi-x)/x);
       y:=arctan(x1/SQRT(1-x1*x1))*arctan(x1/SQRT(1-x1*x1))*arctan(x1/SQRT(1-x1*x1));
       WriteLn( 'при x= ' , X:2:0, 'функция y= ', Y:7:4);
      end;
     x:=x+dx;
   End;
End.

 if (sin(Pi-x)/x<0)or(sin(Pi-x)/x>=1) then WriteLn( 'при x= ' , X:2:0, 'функция y= не определена')