//xishu：3×4陣列//x,y,z方程組的解 private void solutionLinearEqations(double[,] xishu, out double x, out double y, out double z)        {            x = y = z = 0;            double[,] det0 = new double[3, 3];            for (int i = 0; i < 3; i++)                for (int j = 0; j < 3; j++) det0[i, j] = xishu[i, j];            double det00 = determinant(det0);            double[,] detx = new double[3, 3];            detx[0, 0] = xishu[0, 3];            detx[1, 0] = xishu[1, 3];            detx[2, 0] = xishu[2, 3];            for (int i = 0; i < 3; i++)                for (int j = 1; j < 3; j++) detx[i, j] = xishu[i, j];            double detx0 = determinant(detx);            double[,] dety = new double[3, 3];            dety[0, 1] = xishu[0, 3];            dety[1, 1] = xishu[1, 3];            dety[2, 1] = xishu[2, 3];            for (int i = 0; i < 3; i++)                for (int j = 0; j < 3; j++)                {                    if (j != 1)                        dety[i, j] = xishu[i, j];                }            double dety0 = determinant(dety);            double[,] detz = new double[3, 3];            detz[0, 2] = xishu[0, 3];            detz[1, 2] = xishu[1, 3];            detz[2, 2] = xishu[2, 3];            for (int i = 0; i < 3; i++)                for (int j = 0; j < 2; j++) detz[i, j] = xishu[i, j];            double detz0 = determinant(detz);            x = detx0 / det00;            y = dety0 / det00;            z = detz0 / det00;        }