bryn

function[xnext,iteration] = jarcobi(A,b,e,kmax)

 x=zeros(size(b));
 [n,m] = size(A);

 iteration = kmax;

 nondia = A - eye(n).*(A);
 dia = diag(A);

 x = (b-nondia*x)./dia;

 for k =1:kmax
     fprintf('Loop iteration %i\n', k)

     xnext = (b-nondia*x)./dia;

     diff = sum((((xnext-x))./x).^2);
     fprintf('The error is %f\n', diff)

     if diff<e
         iteration = k
         break
     end
     x = xnext;
 end
end