Вот, реализовал метод сопряженных направлений решения СЛАУ, описанный в работах по приведенной ссылке.
Вот что получилось:
Имеем СЛАУ с несимметричной неположительно определенной матрицей
Задается итерационный процесс:
нулевой шаг:

,

,

.

-тый шаг:
Я не очень разобрал, этот метод является собственной разработкой автора или это стандартная его реализация? Дело в том, что он приводится наряду со стандартными методами Гаусса и Халецкого, но с другой стороны выбор чисел

и

, если я правильно понял, предложен автором.