Решение систем линейных алгебраических уравнений прямыми методами.

Найти LU-разложение матрицы  [Graphics:1.gif] .

Введем матрицу

[Graphics:2.gif]

[Graphics:3.gif]

С помощью встроенной функции [Graphics:4.gif] найдем LU-разложение. Функция [Graphics:5.gif] выдыет список из трех элементов, первым из которых является совмещенные нижняя диагональная и верхняя диагональная матрицы.

[Graphics:6.gif]

[Graphics:7.gif]

Чтобы выделить верхнюю и нижнюю диагональные матрицы, введем следующие функции

[Graphics:8.gif]

[Graphics:9.gif]

Это нижняя диагональная матрица

[Graphics:10.gif]

[Graphics:11.gif]

Это верхняя диагональная матрица

[Graphics:12.gif]

[Graphics:13.gif]

[Graphics:14.gif]