A = [4,1,0;1,4,1;0,1,4];
b = [0.5;0.5;0.5];
gauss_elimination(A, b)
