In the world of numerical linear algebra, iterative methods are essential for solving large, sparse systems of linear equations, ( Ax = b ). Among the most famous classical iterative techniques are the Jacobi, Gauss-Seidel, and Successive Over-Relaxation (SOR) methods.

However, you may have encountered a variant called the method. While it sounds more advanced, the "conversion" from MSOR to SOR is not a transformation of results but rather a conceptual and algorithmic simplification.

[ x_i^(k+1) = (1 - \omega) x_i^(k) + \frac\omegaa_ii \left( b_i - \sum_j < i a_ij x_j^(k+1) - \sum_j > i a_ij x_j^(k) \right) ]

Set all ( \omega_i ) in your code to this single ( \omega ). The algorithm becomes: [ x_i^(k+1) = (1 - \omega) x_i^(k) + \frac\omegaa_ii \left( b_i - \sum_j < i a_ij x_j^(k+1) - \sum_j > i a_ij x_j^(k) \right) ]