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

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