Вычислить значение выражения

program laba4;
uses crt, Math;
var x,y:float;
begin
clrscr();
writeln('Vedite x');
readln(x);
if x>0 then
begin
y := arctan(sqrt(1-sqr(power(x,arctan(((x+1)/sin(x)/cos(x))/sqrt(1-sqr((x+1)/sin(x)/cos(x)))))))/power(x,arctan(((x+1)/sin(x)/cos(x))/sqrt(1-sqr((x+1)/sin(x)/cos(x))))));
writeln('result= ',y:5:6);
readln;
end
else
writeln('ne verno odz');
end.

uses crt,math;
//введем дополнительно переменную а, угол в интервале (pi/2;3pi/2)
function odz(x,a:real):boolean;
begin
//поскольку x должен быть больше 0 (при х=0 tg=0, деление на 0, а при x<0 нельзя возводить х в 
//вещественную степень) берем интервал угла первый положительный пригодный для тангенса
//если угол в пределах (pi/2;3pi/2)
//и тангенс не равен 0
//и значение синуса и косинуса=[-1;1]
// то Х в ОДЗ
odz:=(a>pi/2)and(a<3*pi/2)and(tan(a)<>0)and(abs((x+1)/tan(a))<=1)
      and(abs(power(x,arcsin((x+1)/tan(a))))<=1)
end;
function y(x,a:real):real;
begin
y:=arccos(power(x,arcsin((x+1)/tan(a))))
end;
var x,a:real;
begin
clrscr;
write('x=');
readln(x);
if x=-1 then write('y=0')//частный случай х=-1 y=0
else
 begin
  a:=x;//получим угол в нужном интервале если это возможно
  while a>3*pi/2 do a:=a-pi;
  if not odz(x,a) then write('х не принадлежит ОДЗ')
  else write('y=',y(x,a):0:3);
 end;
readln
end.