jor.m

来自「various matlab programs to slove various」· M 代码 · 共 19 行

function [x, iter]= jor ( a, b, x0, nmax, toll, omega)
[n,n]=size(a);
iter = 0; r = b - a * x0; r0 = norm(r); err = norm (r); x = x0;
while err > toll & iter < nmax
 iter = iter + 1;
 for i=1:n
   s = 0;
   for j = 1:i-1   
      s = s + a (i,j) * x (j);
   end
   for j = i+1:n
      s = s + a (i,j) * x (j);
   end
   x (i) = omega * ( b(i) - s) / a(i,i) + (1 - omega) * x(i);
 end
 r = b - a * x;
 err = norm (r) / r0;
end
return