Найти координаты точек пересечения прямой.

procedure check(x, y : real);
   begin
      write('Точка (', x:0:4, ', ', y:0:4, ') в ');
      if (x > 0) and (y > 0) then writeln('I четверти.');
      if (x < 0) and (y > 0) then writeln('II четверти.');
      if (x < 0) and (y < 0) then writeln('III четверти.');
      if (x > 0) and (y < 0) then writeln('IV четверти.')
   end;
 
var k, b, r, d : real; {входные параметры и дискриминант}
   x1, y1, x2, y2 : real; {координаты первой и второй точки}
begin
   write('k = '); readln(k);
   write('b = '); readln(b);
   write('r = '); readln(r);
   d := sqr(2 * k * b) - 4 * (1 + sqr(k)) * (sqr(b) - sqr(r));
   if d < 0 then writeln('Прямая не пересекает окружность.')
   else if d = 0 then begin
      x1 := (-k * b) / (1 + sqr(k));
      y1 := k * x1 + b;
      writeln('Прямая касается окружности.');
      check(x1, y1)
   end else begin
      x1 := (-2 * k * b - sqrt(d)) / (2 * (1 + sqr(k)));
      x2 := (-2 * k * b + sqrt(d)) / (2 * (1 + sqr(k)));
      y1 := k * x1 + b;
      y2 := k * x2 + b;
      writeln('Прямая пересекает окружность в двух точках:')
      check(x1, y1);
      check(x2, y2)
   end
end.