% A =[6 -3 1; 4 -2 1; 1 3 6];
% b=[15 16 1]'; 
% x0=[0 0 0]'; 
% omega=1.536;
% tolerantion=0.01;
% max_iter=200;
% sor(A, b, x0, omega, tolerantion, max_iter);

sor([6 -3 1; 4 -2 1; 1 3 6],[15 16 1]',[0 0 0]',  1.536  ,0.01,200);
