Метод простых итераций работает только для 2 уравнений, а для 3 и более зависает

public double[] метод_простых_итераций()
        {
            int N = matrix.GetLength(0);
            double[] TempX = new double[N];
            double[,] A = new double[N, N];
            for (int i = 0; i < N; i++)
            {
                for (int j = 0; j < N; j++)
                {
                    A[i, j] = matrix[i, j];
                }
            }
            double[] F = new double[N];
            double[] X = new double[N];
            for (int i = 0; i < N; i++)
            {
                F[i] = matrix[i, N];
                X[i] = F[i] / A[i, i];
            }
            double norm;
            do
            {
                for (int i = 0; i < N; i++)
                {
                    TempX[i] = F[i];
                    for (int g = 0; g < N; g++)
                    {
                        if (i != g)
                            TempX[i] -= A[i, g] * X[g];
                    }
                    TempX[i] /= A[i, i];
                }
                norm = Math.Abs(X[0] - TempX[0]);
                for (int h = 0; h < N; h++)
                {
                    if (Math.Abs(X[h] - TempX[h]) > norm)
                        norm = Math.Abs(X[h] - TempX[h]);
                    X[h] = TempX[h];
                }
            } while (norm > eps);
            return X;
        }