Метод Гаусса-Жордана

Хммм... Насколько помню - это метод решения СЛАУ основанный на последовательном получении нулей в столбцах (кроме одной строки).
Например есть матрица, описывающая СЛАУ из 3х уравнений:
a11 a12 a13 b1
a21 a22 a23 b2
a31 a32 a33 b3
Задача получить матрицу вида
d11 d12 d13 c1
0 d22 d23 c2
0 d32 d33 c3
Это реализуется следующим образом: умножаем 2ю строку исходной матрицы на (-a11/a21) и складываем с первой. Таким образом d22=(-a11/a21)*a22+a12 и d23==(-a11/a21)*a23+a13. Аналогично поступаем с третьей строкой, умножая на (-a11/a31) и складываю с первой.
Следующим этапом получаем 0 в элементе d32. Для этого первую строку оставляем в покое, третью умножаем на (-d22/d32) и складываем со второй.
В результате получим треугольную матрицу по которой и находим решение СЛАУ.

Т.е. надо получить нули под главной диагональю?

Ага. Метод работает. Проверил на системке. Но это точно метод Гаусса-Жордана? Где-то слышал, что надо матрицу приводить к единичной, и тогда в расширенной матрице 4-ый столбец(система из 3 уравнений) будет решением СЛАУ.

Так и есть

[QUOTE=Ginza9]Где-то слышал, что надо матрицу приводить к единичной, и тогда в расширенной матрице 4-ый столбец(система из 3 уравнений) будет решением СЛАУ.[/QUOTE]

Что теперь стоит привести её к единичной? Аналогично тому, как матрица приводилась к верхне-триугольной, приводишь её к диагональной, а потом к единичной -> ответ готов :)

Это точно метод Жордана-Гаусса :)
К единичной матрицу приводить не обязательно, ибо не всегда удастся.