* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
630
Решение задач линейной алгебры и оптимизации
> x:=BackwardSubstitute(R,V2);
> A.x-B;
Тут, пожалуй, любопытно, что погрешность вычислений оказалась несколько выше, чем при использовании функции LinearSolve. Однако погрешность не вы ходит за рамки допустимой по умолчанию.
8.3.6. Решение системы линейных уравнений методом декомпозиции Холесски
Выполним решение еще и методом декомпозиции Холесски:
> x:=LinearSolve(A, B, method='Cholesky');
Приведем еще один пример решения системы из четырех линейных уравнений с применением метода декомпозиции Холесски:
> M_temp := Matrix(4, (i,j)->i+i*j-7, shape=triangular[lower]); M :=M_temp.Transpose(M_temp); IsMatrixShape(M,symmetric); IsDefinite(M);