fprintf('Потребление продукции в отрасли...\n')
D = [20, 15, 30;  20, 30, 25; 0, 45, 100]

fprintf('Валовый выпуск продукции в отраси...\n')
X = [200; 150; 100]

fprintf('>>> [1]\n')
fprintf('Объем конечной продукции по отраслям...\n')
Y = X - sum(D')'

fprintf('>>> [2]\n')
fprintf('Технологическая матрица...\n')
A = D ./ X

fprintf('>>> [3]\n')
fprintf('Единичная матрица...\n')
E = eye(3)

fprintf('Существует ли матрица полных затрат?\n')
det(E - A)

fprintf('Матрица полных затрат(МПЗ)...\n')
S = inv(E - A)
fprintf('Не все элементы МПЗ положительны -> МПЗ не продуктивна\n')

fprintf('S * (E - A)) == E \n')
S * (E - A)

fprintf('(E - A) * S) == E \n')
(E - A) * S

fprintf('>>> [4]\n')
fprintf('Вектор валовой продукции X = (S * Y)\n')
Y = [110; 200; 175]
X = S * Y