* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
Работа с пакетом LinearAlgebra и алгоритмами NAG
> x:=BackwardSubstitute(U,V3);
629
> A.x-B;
8.3.5. Решение системы линейных уравнений методом QR декомпозиции
Выполним теперь решение для тех же исходных данных методом QR декомпози ции, обозначив метод в функции LinearSolve:
> x := LinearSolve(A, B, method='QR');
> A.x-B;
Другой, более явный, но и более громоздкий метод решения представлен ниже:
> Q,R := QRDecomposition(A);
> V2:=Transpose(Q).B;