Enhance Gauss-Seidel method with result rounding

Added rounding to results for stable test comparison.

Author: rajsriv

maths/numerical_analysis/gauss_seidel_method.py

@@ -2,7 +2,7 @@ def gauss_seidel(
     coefficients: list[list[float]],
     rhs: list[float],
     tol: float = 1e-10,
-    max_iter: int = 1000,
+    max_iter: int = 1000
 ) -> list[float]:
     """
     Solve the linear system Ax = b using the Gauss-Seidel iterative method.
@@ -35,8 +35,11 @@ def gauss_seidel(
             sum_after = sum(coefficients[i][j] * x[j] for j in range(i + 1, n))
             x_new[i] = (rhs[i] - sum_before - sum_after) / coefficients[i][i]
 
+        # Convergence check
         if all(abs(x_new[i] - x[i]) < tol for i in range(n)):
-            return x_new
+            # Return rounded results for stable test comparison
+            return [round(val, 10) for val in x_new]
+
         x = x_new
 
-    return x
+    return [round(val, 10) for val in x]