求解直线和圆的交点

C#求结直线ax+by+c=0与圆的交点程序

 //直线和圆的交点
  public static void InterC2L(double centerR, double centerC, double R, double a, double b, double c, out List<double> Rows, out List<double> Columns)
  {
       Rows = new List<double>(); //交点1
       Columns = new List<double>();//交点2
	   //点到直线的距离，参考我的其他博文
       double dis = DistanceP2L(centerR, centerC, a, b, c); 
       if (dis <= R)//有交点
       {
           if (a == 0) 
            {
               Columns.Add(-c / b);
               Columns.Add(-c / b);
               double sub = R * R - (Columns[0] - centerC) * (Columns[0] - centerC);
               double kf = Math.Sqrt(sub);
               Rows.Add(-kf + centerR);
               Rows.Add(kf + centerR);
             }
             else if (b == 0)
             {
                Rows.Add(-c / a);
                Rows.Add(-c / a);
                double sub = R * R - (Rows[0] - centerR) * (Rows[0] - centerR);
                double kf = Math.Sqrt(sub);
                Columns.Add(centerC - kf);
                Columns.Add(kf + centerC);
             }
             else
             {
                 double k = -a / b;
                 double _b = -c / b;
                 double equ_a, equ_b, equ_c;
                 equ_a = 1 + k * k;
                 equ_b = -2 * centerR + 2 * k * (_b - centerC);
                 equ_c = centerR * centerR + (_b - centerC) * (_b - centerC) - R * R;
                 
				//求解一元二次方程，参考我的其他博文
                 GenEqu2(equ_a, equ_b, equ_c, out double res1, out double res2);
                 Rows.Add(res1);
                 Rows.Add(res2);

                 double sub = R * R - (Rows[0] - centerR) * (Rows[0] - centerR);
                 double kf = Math.Sqrt(sub);
                 double col1 = -kf + centerC;
                 double col2 = kf + centerC;
				//本应该是直接==，但是程序结算结果又一定的偏差，故用绝对值代替
				//如果有更好的方法或者程序有问题，欢迎留言
                 if (Math.Abs(col1 - (k * Rows[0] + _b)) <= 0.001)
                    Columns.Add(col1);
                 else
                    Columns.Add(col2);

                 sub = R * R - (Rows[1] - centerR) * (Rows[1] - centerR);
                 kf = Math.Sqrt(sub);
                 col1 = -kf + centerC;
                 col2 = kf + centerC;
                 if (Math.Abs(col1 - (k * Rows[1] + _b)) <= 0.001)
                     Columns.Add(col1);
                 else
                     Columns.Add(col2);
            }

       }

 }