Метод Гаусса: объяснить что происходит в каждой строчке кода

 /* метод Гаусса для решения системы уравнений, заданной 
 матрицей коэффициентов matr.
 Метод рекурсивно приводит матрицу к треугольному виду,
 параметр line определяет номер уравнения (ряда), «ниже» которого
 необходимо обнулить коэффициенты под номером line */
 
float* gauss(float** matr, int rows, int columns, int line) {
    if (line == rows - 1) {
        float* result = new float[rows];
        
        for (int i = rows - 1; i >= 0; i--) {
            float temp = matr[i][columns - 1];
            for (int j = i + 1; j < rows; j++) {
                temp = temp - matr[i][j] * result[j];
            }
            if (matr[i][i] == 0) {
                result[i] = 0;
            }
            else {
                result[i] = temp / matr[i][i];
            }
        }
        return result;
    }
    cout << " "  << endl;
    float koef = matr[line][line];
    for (int i = line + 1; i < rows; i++) {
        float koef2 = matr[i][line];
        for (int j = 0; j < columns; j++) {
            matr[i][j] = matr[i][j] * koef - koef2 * matr[line][j];
        }
    }
    
    print(matr, rows, columns);
    printf("\n");
    return gauss(matr, rows, columns, line + 1);
}

float *koefs = gauss(matrix, 4, 5, 0); // решение системы уравнений методом Гаусса
 
    for (int i = 0; i < 4; i++) { // вывод коэффициентов на экран
        cout << koefs[i] << endl;