Place the 'omega' parameter to successive_over_relaxation() first -- a more natural...
[octave.git] / iterative / gauss_seidel.m
index a844142a73ae8be1a67223d0c3b68e4a49c583e2..62826fb215ada0c4f8a43935b6a44451d1bcb883 100644 (file)
@@ -13,9 +13,9 @@ function [x, iterations, residual_norms] = ...
 
   if (nargout > 2)
     [x, iterations, residual_norms] = ...
-      successive_over_relaxation(A, b, omega, x0, tolerance, max_iterations);
+      successive_over_relaxation(omega, A, b, x0, tolerance, max_iterations);
   else
     [x, iterations] = ...
-      successive_over_relaxation(A, b, omega, x0, tolerance, max_iterations);
+      successive_over_relaxation(omega, A, b, x0, tolerance, max_iterations);
   end
 end