SOR Method by matlab

norm 부분 미완성 =_=;;
 

A=[4,6,8,2;6,10,16,6;8,16,36,20;2,6,20,15];

b=[48;98;228;134];

n=length(A);

k=1;

R=zeros(4);

w=1;

err=100;

while(err > 10^-4)

    T=0;

    for i=1:1:n

        T=0;

        for j=1:1:i-1

            T=T+A(i,j)*R(j,k+1);

        end

        for j=i+1:1:n

            T=T+A(i,j)*R(j,k);

        end

        R(i,k+1)=(1-w)*R(i,k)+(w/A(i,i))*(b(i)-T);

    end    

    

%     sig1=0;sig2=0;

%     for v1=1:1:4

%         sig1 = sig1 + (abs(R(v1,k+1)-R(v1,k)))^2;

%         sig2 = sig2 + (abs(R(v1,k)))^2;

%     end

%     

%     err = (sig1^1/2)/(sig2^1/2);

    err = norm(R(:,k+1)-R(:,k),2)/norm(R(:,k),2);

    k=k+1;

end

        err

        k