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

% Решить систему Ax=b методом простой итерации
A = [3 -2 0; -2 3 0; 0 0 3];
b = [-21;24;15];

e = eig(A);
% Вычислим итерационный параметр
tau = 2 / (min(e) + max(e));

r = 1;

% Выполняем метод простой итерации с начальным приближением x0
x0 = [0;0;0];
y = x0;
while r > 100 * eps
x = y - tau * (A * y - b);
r = norm(x-y);
y = x;
end
% Проверим полученное значение
A * x - b

>>
ans =
1.0e-013 *