Разработать алгоритм и программу вычисления координат критических точек

var
  x0o, x0e, y0o, y0e: real;
  x1, x2, y1, y2: real;
 
const
  r = 20;
  a = 25;
  b = 35;
  dx = 0.001;
  dy = 0.001;
procedure A_1;
var
  x, y: real;
begin
  x := 0;
  y := -1.0983*x+21.8298175565;
    while ((sqr(x-x0e))/sqr(a)+(sqr(y-y0e))/sqr(b)<1) do    
  begin
    x := x - dx;
    y := -1.0983*x+21.8298175565;
  end;
  x1 := x;
  y1 := y;
end;
 
procedure B_1;
var
  x, y: real;
begin
  x := 0;
  y := (-sqrt((sqr(a) * sqr(b) - sqr(b) * sqr(x - x0e)) / sqr(a)) + y0e);
  while ((x-x0o)*(x-x0o)+(y-y0o)*(y-y0o) > r*r)  do
  begin
    x := x + dx;
    y := (-sqrt((sqr(a) * sqr(b) - sqr(b) * sqr(x - x0e)) / sqr(a)) + y0e)
  end;
  x2 := x;
  y2 := y;
end;
 
 
begin
  x0o:=8;
  y0o:=7.373;
  x0e:=4.938;
  y0e:=3;
  A_1;
  B_1;
  writeln('-' * 58);
  writeln('|', 'Наименование', '|', 'Аналитические расчёты', '|', 'Вычислительные методы', '|');
  writeln('|', '    точка   ', '|', '     X    |', '     Y    |', '     X    |', '     Y    |');
  writeln('-' * 58);
  writeln('|', '      ', 'A1', '    ', '|  ', '-5.130', '  |  ', '27.464', '  |  ', x1:6:3, '  | ', y1:6:3, '   |');
  writeln('-' * 58);
  writeln('|', '      ', 'B1', '    ', '|  ', '17.623', '  | ', '-14.613', '  |  ', x2:6:3, '  | ', y2:6:3, '  |');
  writeln('-' * 58);
end.