Suprarelaxare

function x = suprarelaxare2(A, x0, b, w, maxiter, epsilon)
    [n n] = size(A);
    xc = zeros(n,1);
    xp = x0;
    ok = 1;
    for k = 1 : maxiter
        for i = 1 : n
            ok = 1;
            s1 = 0;
            s2 = 0;
            for j = 1: i-1
                s1 = s1 + A(i,j)*xc(j);
            endfor
            for j = i+1 : n
                s2 = s2 + A(i,j)*xp(j);
            endfor
            xc(i) = w *(b(i) - s1 - s2)/A(i,i) + xp(i)*(1-w);;
            if norm(xc - xp) < epsilon
                ok = 0;
                break;
            endif
            xp = xc;
        endfor
    endfor
    x = xc;
    if ok == 0
        disp('Am ajuns la convergenta inainte de maxiter!');
        disp(poz);
    endif
endfunction