三点求圆心模板

 void Cir(double ax,double ay,double bx,double by,double cx,double cy,double &x,double &y)
 {
      double a1=atan2(by-ay,bx-ax)+acos(-1.)/2;
      double a2=atan2(cy-by,cx-bx)+acos(-1.)/2;
      ax=(ax+bx)/2,ay=(ay+by)/2;
      bx=(cx+bx)/2,by=(cy+by)/2;
      double r=(sin(a2)*(ax-bx)+cos(a2)*(by-ay))/(sin(a1)*cos(a2)-sin(a2)*cos(a1));
      x=ax+r*cos(a1),y=ay+r*sin(a1);
 }
 // 在这模板下面 同一直线上的点所得得圆心坐标值将超过 0x3fffffff 可定义无穷大的数来排除