給定三角形的三點 求外心

Point heart(int a,int b,int c)
{
    double a1=p[b].x-p[a].x;
    double b1=p[b].y-p[a].y;
    double c1=(p[b].x*p[b].x+p[b].y*p[b].y-p[a].x*p[a].x-p[a].y*p[a].y)/2;
 
    double a2=p[c].x-p[b].x;
    double b2=p[c].y-p[b].y;
    double c2=(p[c].x*p[c].x+p[c].y*p[c].y-p[b].x*p[b].x-p[b].y*p[b].y)/2;
    double  x=(b2*c1-b1*c2)/(a1*b2-a2*b1);
    double y=(c1*a2-a1*c2)/(b1*a2-a1*b2);
    return Point(x,y);


}