* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
628
Решение задач линейной алгебры и оптимизации
Прямое решение этим методом выполняется одной из двух команд, отличаю щихся формой записи:
> x := LinearSolve(A, B, method='LU'); x := LinearSolve(
, method='LU');
Проверим решение данной системы уравнений:
> A.x-B;
В данном случае решение точно (в пределах точности вычислений по умол чанию). Можно также выполнить решение, проведя отдельно LU декомпозицию, что делает наглядным алгоритм решения и операции подстановки:
> P,L,U:=LUDecomposition(A);
> V2:=Transpose(P).B;
> V3:=ForwardSubstitute(L,V2);