QR – разложение, Fortran.

1) Написать две программы, которые для N×N матрицы A вычисляют QR разложение, используя элементарные матрицы Хаусхолдера, и находят A-1, используя найденное разложение. Q – ортогональная матрица, R – верхняя треугольная матрица.

2) Оформить отчет по программе.

Содержание отчета.
1. Формулировка задачи (конкретный вариант).
2. Расчётные формулы.
3. Тестовые примеры для демонстрации работы программ.
4. Тексты программ и результаты их работы.
Результаты работы должны содержать:
- печать исходных данных,
- результаты промежуточных вычислений,
- невязки для вычисленных разложений матриц, систем линейных алгебраических уравнений и обратных матриц,
- разность между точным и приближенным результатом.

Задачу требуется выполнить в программе Code::Blocks.

QR – разложение и подробное описание прикрепил в файле.