Метод гаусса для решений СЛАУ

[Вечно грустная...]

Помогите создать метод который пребразует в нули все элементы ниже главной диагонали(прямой ход)?

Вот нашла один метод:

Помогите пожалуйста разобрать код, на бумажке решила этим методом но код разобрать не получается=)

Не могу понять все эти бесконечные циклы=)

void PryamoiHod(int n, double **a, double *b)
{
        double v;
        for(int k = 0,i,j,im; k < n - 1; k++)
        {
                im = k;
                for(i = k + 1; i < n; i++)
                {
                        if(fabs(a[im][k]) < fabs(a[i][k]))
                        {
                                im = i;
                        }
                }
                if(im != k)
                {
                        for(j = 0; j < n; j++)
                        {
                                v                = a[im][j];
                                a[im][j] = a[k][j];
                                a[k][j]  = v;
                        }
                        v     = b[im];
                        b[im] = b[k];
                        b[k]  = v;
                }
                for(i = k + 1; i < n; i++)
                {
                        v               = 1.0*a[i][k]/a[k][k];
                        a[i][k] = 0;
                        b[i]    = b[i] - v*b[k];
                                                if(v != 0)
                        for(j = k + 1; j < n; j++)
                        {
                                a[i][j] = a[i][j] - v*a[k][j];
                        }
                }
        }
}

[gamer-25]

Вот по другому

#include <iostream.h>
#include <windows.h>


const int ArrSIZE = 5;
void PryamoiHod(int matrix[ArrSIZE][ArrSIZE]);

int main()
{
   int Matrix[ArrSIZE][ArrSIZE] =             {{5,3,7,8,1},
                                               {2,7,9,0,3},
                                               {4,2,8,7,1},
                                               {0,5,7,6,4},
                                               {2,4,9,8,5} };

   PryamoiHod(Matrix);


   cout << endl << endl;
   system("pause");
   return 0;
}

void PryamoiHod(int matrix[ArrSIZE][ArrSIZE])
{
   int row, column;
   
   for (row = 0; row < ArrSIZE; row++)
      for (column = 0; column < ArrSIZE; column++)
         if (row > column)
            matrix[row][column] = 0;

   // Печать матрицы
   for (row = 0; row < ArrSIZE; row++)
   {
      for (column = 0; column < ArrSIZE; column++)
         cout << matrix[row][column] << "  ";
      cout << endl;
   }

}