Решение систем алгебраических уравнение итерационными методами.

Решить систему [Graphics:1.gif] методом простой итерации, где  [Graphics:2.gif]  и  [Graphics:3.gif].

Введем матрицу и вектор

[Graphics:4.gif]

[Graphics:5.gif]

[Graphics:6.gif]

[Graphics:7.gif]

Вычислим итерационный параметр

[Graphics:8.gif]

[Graphics:9.gif]

Нам нужно будет вычислять нормы векторов и матриц, поэтому загрузим пакет [Graphics:10.gif]

[Graphics:11.gif]

Проверим условие сходимости метода

[Graphics:12.gif]

[Graphics:13.gif]

Запишем процедуру нахождения решения системы методом простых итераций с точностью [Graphics:14.gif] и начальным приближением [Graphics:15.gif]

[Graphics:16.gif]

Найдем решение при начальном приближении [Graphics:17.gif]

[Graphics:18.gif]

[Graphics:19.gif]

Ту же самую процедуру можно записать более компактным способом

[Graphics:20.gif]

Найдем решение при начальном приближении [Graphics:21.gif]

[Graphics:22.gif]

[Graphics:23.gif]

Можно найти решение с помощью встроенной функции [Graphics:24.gif]

[Graphics:25.gif]

[Graphics:26.gif]

[Graphics:27.gif]