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

a := sqr(b1) + sqr(a1);
b := 2 * b1 * c1 - 2 * a1 * x * b1 - 2 * sqr(a1) * y;
c := sqr(c1) - 2 * a1 * x * c1 + sqr(a1) * x + sqr(a1) * sqr(y) - sqr(a1) * sqr(r);

var
  x, y, r, a, b, c, a1, b1, c1, d, yo, xo, yo2, xo2: double;
  
begin
  readln(x, y, r, a1, b1, c1);
  a := sqr(b1) + sqr(a1);
  b := 2 * b1 * c1 - 2 * a1 * x * b1 - 2 * sqr(a1) * y;
  c := sqr(c1) - 2 * a1 * x * c1 + sqr(a1) * x + sqr(a1) * sqr(y) - sqr(a1) * sqr(r);
  d := sqr(b1) - 4 * a1 * c1;
  if d < 0 then
    writeln(0)
  else
    if d = 0 then
    begin
      yo := -b / (2 * a);
      xo := -((b1 * yo + c1) / a1);
      writeln(1);
      writeln(xo, ' ', yo);
    end
    else
    begin
      yo := (-b + sqrt(d)) / (2 * a);
      xo := -((b1 * yo + c1) / a1);
      yo2 := (-b - sqrt(d)) / (2 * a);
      xo2 := -((b1 * yo2 + c1) / a1);
      writeln(2);
      writeln(xo, ' ', yo);
      writeln(xo2, ' ', yo2);
    end;
end.

var x0,y0,r:real;//координаты центра окружности, радиус
    a,b,c:real; //коэффициенты прямой ax+by+c=0
    x1,y1,x2,y2:real;//координаты точек пересечения
    k,s:real;//коэффициенты прямой y=kx+s
    x,a1,b1,c1,d:real; //вспомогательные переменные
begin
clrscr;
repeat
write('Радиус окружности r>0 r=');
readln(r);
until r>0;
writeln('Координаты центра окружности');
readln(x0,y0);
repeat
writeln('Коэффициенты a,b,c a и b одновремено не равны 0');
readln(a,b,c);
until(a<>0)or(b<>0);
if b=0 then //прямая вертикальна
 begin
  x:=-c/a;
  if (x0-r>x)or(x0+r<x) then write('Окружность и прямая не пересекаются')
  else if (x0-r=x)or(x0+r=x) then  write('Прямая касается окружности в точке (',x:0:2,';',y0:0:2,')')
  else
   begin
    y1:=y0+sqrt(r*r-sqr(x0-x));
    y2:=y0-sqrt(r*r-sqr(x0-x));
    write('Прямая пересекает окружность в 2-х точках (',x:0:2,';',y1:0:2,')  (',x:0:2,';',y2:0:2,')');
   end;
 end
else //другие случаи
 begin
  k:=-a/b;//преобразуем уравнение ax+by+c=0 к виду y=kx+s
  s:=-c/b;
  //и решим систему
  //(x - x0)^2 + (y - y0)^2 = r^2 (окружность)
  //y=kx+s {прямая}
  //после приведения найдем коэффициенты уравнения a1x^2+b1x+c1=0
  a1:=1+k*k;
  b1:=-2*x0+2*k*s-2*k*y0;
  c1:=s*s+x0*x0-2*s*y0+y0*y0-r*r;
  d:=b1*b1-4*a1*c1;
  if d<0 then write('Прямая и окружность не пересекаются')
  else if d=0 then
   begin
    x1:=-b1/(2*a1);
    y1:=k*x1+s;
    write('Прямая касается окружности в точке (',x1:0:2,';',y1:0:2,')')
   end
  else
   begin
    x1:=(-b1+sqrt(d))/(2*a1);
    x2:=(-b1-sqrt(d))/(2*a1);
    y1:=k*x1+s;
    y2:=k*x2+s;
    writeln('Прямая пересекает окружность в 2-х точках:');
    write('(',x1:0:2,';',y1:0:2,')  (',x2:0:2,';',y2:0:2,')');
   end;
 end
end.