三点求圆心模板

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