1 function [x, iterations, residual_norms] = ...
2 gauss_seidel(A, b, x0, tolerance, max_iterations)
9 max_iterations = intmax();
15 [x, iterations, residual_norms] = ...
16 successive_over_relaxation(omega, A, b, x0, tolerance, max_iterations);
19 successive_over_relaxation(omega, A, b, x0, tolerance, max_iterations);