function Gauss_elim_method(A,b)
n = max(size(A)); 
for j=2:n
   for i=j:n
       m = A(i,j-1) / A(j-1,j-1);
       A(i,:) = A(i,:) - A(j-1,:) * m;
       b(i) = b(i) - m * b(j-1);
    end
end
disp('Upper-triangular form of the system is:');
[A b]
x = zeros(n,1);      
x(n) = b(n) / A(n,n);    
for j=n-1:-1:1
   x(j) = (b(j) - A(j,j+1:n) * x(j+1:n)) / A(j,j);
end
disp('The solution is vector:');
x