用C#写一个矩阵求逆算法

以下是一个简单的矩阵求逆算法的C#实现：

public static double[,] Inverse(double[,] matrix)
{
    int n = matrix.GetLength(0);
    double[,] inverse = new double[n, n];
    double[,] identity = Identity(n);

    // 高斯-约旦消元法求解
    for (int i = 0; i < n; i++)
    {
        double pivot = matrix[i, i];
        for (int j = 0; j < n; j++)
        {
            matrix[i, j] /= pivot;
            identity[i, j] /= pivot;
        }
        for (int j = 0; j < n; j++)
        {
            if (j != i)
            {
                double factor = matrix[j, i];
                for (int k = 0; k < n; k++)
                {
                    matrix[j, k] -= factor * matrix[i, k];
                    identity[j, k] -= factor * identity[i, k];
                }
            }
        }
    }

    return identity;
}

private static double[,] Identity(int n)
{
    double[,] identity = new double[n, n];
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            identity[i, j] = (i == j) ? 1.0 : 0.0;
        }
    }
    return identity;
}

该算法使用高斯-约旦消元法求解逆矩阵。在求解过程中，我们同时对单位矩阵进行相同的变换，最终得到的单位矩阵就是原始矩阵的逆矩阵。

你可以调用`Inverse`方法并传入一个3*3矩阵，该方法会返回该矩阵的逆矩阵。注意，该方法只适用于可逆的3*3矩阵。