Алгоритм Гаусса: что означает ноль в столбце ?

/*
Решение системы N линейных уравнений с N неизвестными методом Гаусса с выбором главного элемента.
n - кол-во уравнений и неизвестных;
Matrix - полная матрица системы, размерность N строк на N+1 столбцов (последний столбец - свободных членов).
    В ходе решения массив будет изменён !
X - одномерный массив решений (вектор X{x1,,xn}).
Возврашает 1 при неудаче, 0 - решение получено.
*/
int Simq(int n, double **Matrix, double *X)
{
    const double eps = 1e-38;
    double * tmp;
    double max,v,u;
    int i,j,l,k1;
 
    if (n < 1) return 1;
    //ПРЯМОЙ ХОД
    for (i = 0; i < n; i++)//идём по диагонали (i - текущий индекс и строки и столбца)
        {
         //если в текущем столбце есть элемент больше, чем он в текущей строке, то запоминаем номер строки (k1)
        max = fabs(Matrix[i][i]); k1 = i;
        for (l = i+1; l < n; l++)
            if (fabs(Matrix[l][i]) > max) {max = fabs(Matrix[l][i]); k1 = l;}
         //если максимальный элемент не в текущей  строке, то меняем местами строки
        if (k1 != i) //(путём обмена указателей на строки)
            {
            tmp = Matrix[k1];
            Matrix[k1] = Matrix[i];
            Matrix[i] = tmp;
            }
        if (max < eps) return 1; //Что означает ?
        u = Matrix[i][i];
        //первый(идём по диагонали !) элемент в текущей строке теперь максимальный (среди эл. тек. столбца)
        //и после следующей операции должен стать равен 1
         //делим все элементы текущей строки на первый(текущий) элемент.
        for (j = i+1 ; j <= n; j++) 
            Matrix[i][j] /= u;
        for (l = i+1; l < n; l++)//Из следующей строки вычитаем текущую, умноженную на
            {
            v=Matrix[l][i]; //элемент в текущем столбце следующей строки.
            for (j = i+1; j <= n; j++) 
                Matrix[l][j] -= Matrix[i][j]*v;
            }
        }
    //получили верхнюю треугольную матрицу и решение xn
    //прим.: в матрице значения по и ниже главной диагонали не рассчитываются, для уменьшения кол-ва операций,
    //ибо их значения (0 и 1) и так известны.
    //ОБРАТНЫЙ ХОД
    X[n-1] = Matrix[n-1][n];
    for (signed int i = n-2; i >= 0; i--)
        { 
        X[i] = Matrix[i][n];//запоминаем свободный член
        //вычитаем из него уже определённые x домноженные на соответствующий коэфф. в матрице.
        for (j = i+1; j < n; j++) X[i] -= Matrix[i][j]*X[j];
        };
    return 0;
};