C++ 计算空间三点球圆心

假设三个点为 A(x1, y1), B(x2, y2), C(x3, y3)，则：

求出 AB 和 AC 的中点 M1 和 M2，分别为((x1+x2)/2, (y1+y2)/2) 和 ((x1+x3)/2, (y1+y3)/2)。
求出 AB 和 AC 的垂直平分线 L1 和 L2，分别为(y1-y2)x + (x2-x1)y + (x1y2-x2y1) = 0 和 (y1-y3)x + (x3-x1)y + (x1y3-x3y1) = 0。
求出 L1 和 L2 的交点，即为圆心 O。O 的坐标为

((L1 的 b * L2 的 c - L2 的 b * L1 的 c) / (L1 的 a * L2 的 b - L2 的 a * L1 的 b), (L2 的 a * L1 的 c - L1 的 a * L2 的 c) / (L1 的 a * L2 的 b - L2 的 a * L1 的 b))，

其中 L1 的 a、b、c 和 L2 的 a、b、c 分别为 L1 和 L2 的一般式中的系数。

求出 OA 的长度 R，即为 O 点到 A 点的距离。R 的值为

sqrt((x1-O 的 x)^2 + (y1-O 的 y)^2)。

以 O 为圆心，R 为半径，画出圆即可。