Пересечения круга и отрезка

bool Circle::cCircle(Line l)
{
    double x1=l.getX();
    double y1=l.getY();
    double x2=l.getX2();
    double y2=l.getY2();
    double x4=((x2-x1)*(y2-y1)*(cy-y1)+x1*(y2-y1)*(y2-y1)+cx*(x2-x1)*(x2-x1))/((y2-y1)*(y2-y1)+(x2-x1)*(x2-x1));
    double y4=(y2-y1)*(x4-x1)/(x2-x1)+y1;
    double s1=dist(cx,cy,x1,y1);
    double s2=dist(cx,cy,x2,y2);
    double s3=dist(cx,cy,x4,y4);
    if ((s1<r) && (s2<r)) {return false;}
    if (y1<y2){
        if (cy+r<=y1 && cy-r>=y2){return false;}
    }
    if (y1>y2){
        if (cy+r<=y2 && cy-r>=y1){return false;}
    }
    if (cx+r<=x1) {return false;}
    if (cx-r>=x2) {return false;}
    if (y1==y2 && x1==x2) {
        if(s1<r){return true;}
    }
    if (r>s1){return true;}
    if (r>s2){return true;}
    if (r>s3) {return true;}
    else {return false;}
}

#include <iostream>
#include <cmath>
using namespace std;
 
bool commonSectionCircle    (   double x1, double y1, double x2, double y2, 
                                double xC, double yC, double R)
{   x1 -= xC;
    y1 -= yC;
    x2 -= xC;
    y2 -= yC;
 
  double dx = x2 - x1;
  double dy = y2 - y1;
 
  //составляем коэффициенты квадратного уравнения на пересечение прямой и окружности.
  //если на отрезке [0..1] есть отрицательные значения, значит отрезок пересекает окружность
  double a = dx*dx + dy*dy;
  double b = 2.*(x1*dx + y1*dy);
  double c = x1*x1 + y1*y1 - R*R;
 
  //а теперь проверяем, есть ли на отрезке [0..1] решения
  if (-b < 0)
    return (c < 0);
  if (-b < (2.*a))
    return ((4.*a*c - b*b) < 0);
 
  return (a+b+c < 0);
}
 
 
int main()
{   double x1,y1,x2,y2, xC,yC,R;
 
    cout<<"x1 y1 x2 y2:\t";
    cin>>x1>>y1>>x2>>y2;
    cout<<"xC yC R:\t";
    cin>>xC>>yC>>R;
 
    cout<<commonSectionCircle(x1,y1,x2,y2,xC,yC,R)<<"\n";
 
    cin.sync();cin.get();
    return 0;
}