Определение попадания точки с заданными координатами X и Y в закрашенную область

CONST R=100; AA=170; BB=80;
VAR p1p2,p2p3,p3p1,p4p5,p5p6,p6p7,p7p4:BOOLEAN;
   X,Y,X0ok,Y0ok,X0el,Y0el:REAL;
BEGIN
   {Здесь должны быть вычислены все необходимые
   параметры используемых далее выражений:
   X0ok,Y0ok,X0el,Y0el и т.д.}
p1p2:=((X-X0el)*(X-X0el))/(AA*AA)+((Y-Y0el)*(Y-Y0el))/(BB*BB) > 1;
p2p3:=X>0;
p3p1:=(X0ok)*(X0ok)+(Y-Y0ok)*(Y-Y0ok) < R*R;
p4p5:=((X-X0el)*(X-X0el))/(AA*AA)+((Y-Y0el) *(Y-Y0el))/(BB*BB) < 1;
p5p6:=X>0;
p6p7:=(X-X0ok)*(X-X0ok)+(Y-Y0ok)*(Y-Y0ok) > R*R;
p7p4:=Y<0;
IF (p1p2 AND p2p3 AND p3p1) OR (p4p5 AND p5p6 AND p6p7 AND p7p4)
THEN
WRITELN ('Точка принадлежит эакрашенной области')
ELSE
WRITELN ('Точка не принадлежит закрашенной области');
END.

const
  R = 24; //радиус окружности
  a = 35; //полуоси эллипса
  b = 17;
  X3 = 25; //координаты соответствующих точек
  X4 = -28;
  Y5 = -31;
  Y6 = 6;
  eps = 1E-14; //точность для метода дихотомии 10^(-14)
var
  X0C, Y0C, X0E, Y0E, X1, Y1, Y1E, X2, Y2, X, Y, XL, XR: real;
  in_c, in_e, out_e, u_l, left, right, in_a, in_b: boolean;
begin
  //центр окружности
  X0C := X3 - R;
  Y0C := sqrt(sqr(R) - sqr(X0C)) + Y5;
  //центр эллипса
  X0E := X4 + a;
  Y0E := Y6 - b * sqrt(1 - sqr(X0E / a));
  //пересечения прямой:
  //пересечение с осью абсцисс
  X2 := sqrt(sqr(R) - sqr(Y0C)) + X0C;
  Y2 := 0;
  //пересечение с осью ординат
  XL := X0C - R; //начальная левая граница интервала поиска корня
  XR := 0; //начальная правая граница интервала поиска корня
  repeat //цикл определения координаты точки пресечения окружности и эллипса
    X1 := (XL + XR) / 2; //центр интервала
    Y1 := Y0C - sqrt(sqr(R) - sqr(X1 - X0C)); //ордината точки на окружности
    Y1E := Y0E - b * sqrt(1 - sqr((X1 - X0E) / a)); //ордината точки на эллипсе
    if Y1 < Y1E then XR := X1 else XL := X1; //определение новых границ интервала
  until (abs(Y1 - Y1E) < eps) and (XR - XL < eps); //выход из цикла при достижении заданной точности
  //ввод координат точки M
  writeln('Введите координаты точки:');
  write('x = ');
  readln(X);
  write('y = ');
  readln(Y);
  //вычисление логических переменных
  in_c := sqr(X - X0C) + sqr(Y - Y0C) <= R * R; //внутри окружности
  in_e := sqr(b * (X - X0E)) + sqr(a * (Y - Y0E)) <= sqr(a * b); //внутри эллипса
  out_e := sqr(b * (X - X0E)) + sqr(a * (Y - Y0E)) >= sqr(a * b); //снаружи эллипса
  u_l := (Y - Y1) / (Y2 - Y1) <= (X - X1) / (X2 - X1); //под прямой
  left := X <= 0; //слева от оси ординат
  right := X >= 0; //справа от оси ординат
  in_a := in_c and out_e and left; //в области A
  in_b := in_c and in_e and u_l and right; //в области B
  //печать результатов
  if in_a then writeln('Точка принадлежит области A')
  else writeln('Точка не принадлежит области A');
  if in_b then writeln('Точка принадлежит области B')
  else writeln('Точка не принадлежит области B');
  readln
end.