Решение систем линейных алгебраических уравнений. Нормы векторов и матриц. Метод Гаусса (схема единственного деления). LU-разложение матрицы.

% Найти LU-разложение матрицы

% Введём матрицу
A = [3,4,-9,5;-15,-12,50,-16;-27,-36,73,8;9,12,-10,-16];

% Выполним LU-разложение
[L,U, P] = lu(A);
% L - нижнетреугольнаЯ, U - верхнетреугольнаЯ, P - матрица перестановок. P * A = L * U
P * A - L * U

>>
ans =
1.0e-014 *