Решение СЛАУ методом Халецкого

Для удобства рассуждений систему запишем в матричном виде: Ax = b, где A = (a_ij )- квадратная матрица порядка n,

- векторы столбцы. Представим матрицу A в виде произведения нижней треуголной матрицы B=(b_ij ) и верхней треугольной матрицы C=(c_ij ) с единичной диагональю, т.е. A = BC. Тогда искомый вектор x может быть вычислен из цепи уравнений By = b, Cx = y. Так как матрицы В и C - треугольные, то системы легко решаются.

Процедура не проверяет, является ли матрица системы невырожденной. В случае вырожденной матрицы алгоритм неприменим.

Если нашли ошибку в алгоритме - сообщите!

Блоксхемы:

К сожалению, для данного алгоритма нет исходного кода. Это не ошибка и сообщать мне об этом не надо. Возможно, для алгоритма ещё не успели создать исходный код, или же при переносе алгоритма из старой версии библиотеки возникли проблемы с исходником и его написание пришлось отложить до лучших времен. Попробуйте воспользоваться прилагаемыми файлами и блок-схемой или поискать на сайте аналогичный алгоритм, но с исходником.

Бочканов Сергей, Быстрицкий Владимир
Copyright © 1999-2004
При поддержке проекта MANUAL.RU