Обратная матрица методом Гаусса, исправить код

        public static double[,] gauss(object matr1)
        {
            double[,] mas = new double[3, 6];
            double[,] matr = ((double[,])matr1);
            double[,] gau = new double[3, 3];
            int ige=0;
            double k;
            for (int i = 0; i <= 2; i++)
            {
                for (int j = 0; j <= 2; j++)
                {
                    mas[i, j] = matr[i, j];
                }
            }
 
            mas[0, 3] = 1; mas[0, 4] = 0; mas[0, 5] = 0;
            mas[1, 3] = 0; mas[1, 4] = 1; mas[1, 5] = 0;
            mas[2, 3] = 0; mas[2, 4] = 0; mas[2, 5] = 1;
            for (int n = 0; n <= 1; n++)
            {
                double ge = mas[n, n];
                for (int i = 1 + n - 1; i <= 2; i++)
                {
                    if (Math.Abs(mas[i, n]) >= ge)
                    {
                        ge = Math.Abs(mas[i, n]);
                        ige = i;
                    }
                }
                for (int j = 0; j <= 5; j++)
                {
                    ge = mas[n, j];
                    mas[n, j] = mas[ige, j];
                    mas[ige, j] = ge;
                }
                for (int i = 2 + n - 1; i <= 2; i++)
                {
                    k = mas[i, n] / mas[n, n];
                    if (k != 0)
                    {
                        for (int j = 0; j <= 5; j++)
                        {
                            mas[i, j] = mas[n, j] - mas[i, j] / k;
                        }
 
                    }
                }
            }
            for (int n = 0; n <= 1; n++)
            {
                for (int i = 2; i >= 1+n; i--)
                {
                    k = mas[i, i] / mas[i - n, i];
                    for (int j = 0; j <= 5; j++)
                    {
                        mas[i - n, j] = mas[i - n, j] - mas[i, j] / k;
                    }
                }
            }
            for (int i = 0; i <= 2; i++)
            {
                for (int j = 0; j <= 2; j++)
                {
                    mas[i, j + 3] = mas[i, j + 3] / mas[i, i];
                }
            }
            for (int i = 0; i <= 2; i++)
            {
                for (int j = 0; j <= 2; j++)
                {
                    gau[i, j] = mas[i, j + 3];
                }
            }
            return gau;
        }

for (int n = 0; n <= 1; n++)
            {
                for (int i = 2; i >= 1+n; i--)
                {
                    k = mas[i, i] / mas[i - n, i];
                    for (int j = 0; j <= 5; j++)
                    {
                        mas[i - n, j] = mas[i - n, j] - mas[i, j] / k;
                    }
                }
            }