给出了球体计算器的方程，给定球体上的四个点
中心位于点 \( (h,k,l) \) 且半径为 \( r \) 的球体方程由下式给出
\[ (x-h)^2 + (y - k)^2 + (z - l)^2 = r^2 \]
展开上面的方程和分组项，我们得到
\[ x^2 + y^2 + z^2 - 2xh - 2 y k - 2 z l + h^2 + k^2 + l^2 = r^2 \]
还是
\[ x^2 + y^2 + z^2 = 2 x h + 2 y k + 2 z l + r^2 - ( h^2 + k^2 + l^2) \]
设 \( A = 2 h \) 、 \( B = 2 k \) 、 \( C = 2 l \) 和 \( D = r^2 - (h^2 + k^2 + l^2) \) 并将上面重写如下
\[ A x + B y + C z + D = x^2 + y^2 + z^2 \]
对于坐标为 \( (x_1 , y_1 , z_1 ) \) 、 \( (x_2 , y_2 , z_2 ) \) 、 \( (x_3 , y_3 , z_3 ) \) 和 \( (x_4 , y_4 , z_4 ) 的 4 个点 ) \) 要在同一球面上，以下具有未知数 \( A, B , C \) 和 \( D \) 的四个方程组必须有解：
\[ \begin{align*}
A x_1 + B y_1 + C z_1 + D & = x_1^2 + y_1^2 + z_1^2 \\
A x_2 + B y_2 + C z_2 + D & = x_2^2 + y_2^2 + z_2^2 \\
A x_3 + B y_3 + C z_3 + D & = x_3^2 + y_3^2 + z_3^2 \\
A x_4 + B y_4 + C z_4 + D & = x_4^2 + y_4^2 + z_4^2
\end{align*} \]
使用克莱默规则，我们求解 \( A, B , C \) 和 \ (D\)如下：
\[ M = \begin{vmatrix}
&x_1&&y_1&&z_1&&1&\\
\\
&x_2&&y_2&&z_2&&1&\\
\\
&x_3&&y_3&&z_3&&1&\\
\\
&x_4&&y_4&&z_4&&1&\\
\end{vmatrix} \]
以及
\[ A = \dfrac{\begin{vmatrix}
&x_1^2 + y_1^2 + z_1^2&&y_1&&z_1&&1&\\
\\
&x_2^2 + y_2^2 + z_2^2 &&y_2&&z_2&&1&\\
\\
&x_3^2 + y_3^2 + z_3^2&&y_3&&z_3&&1&\\
\\
&x_4^2 + y_4^2 + z_4^2&&y_4&&z_4&&1&\\
\end{vmatrix}}{M} \]
\[ B = \dfrac{\begin{vmatrix}
&x_1&&x_1^2 + y_1^2 + z_1^2&&z_1&&1&\\
\\
&x_2&&x_2^2 + y_2^2 + z_2^2&&z_2&&1&\\
\\
&x_3&&x_3^2 + y_3^2 + z_3^2&&z_3&&1&\\
\\
&x_4&&x_4^2 + y_4^2 + z_4^2&&z_4&&1&\\
\end{vmatrix}}{M} \]
\[ C = \dfrac{\begin{vmatrix}
&x_1&&y_1&&x_1^2 + y_1^2 + z_1^2 &&1&\\
\\
&x_2&&y_2&&x_2^2 + y_2^2 + z_2^2&&1&\\
\\
&x_3&&y_3&&x_3^2 + y_3^2 + z_3^2&&1&\\
\\
&x_4&&y_4&&x_4^2 + y_4^2 + z_4^2&&1&\\
\end{vmatrix}}{M} \]
\[ D = \dfrac{\begin{vmatrix}
&x_1&&y_1&&z_1&&x_1^2 + y_1^2 + z_1^2&\\
\\
&x_2&&y_2&&z_2&&x_2^2 + y_2^2 + z_2^2&\\
\\
&x_3&&y_3&&z_3&&x_3^2 + y_3^2 + z_3^2&\\
\\
&x_4&&y_4&&z_4&&x_4^2 + y_4^2 + z_4^2&\\
\end{vmatrix}}{M} \]
输入四个 3D 点的坐标：