function x = Gaussian_Elimination(A,b)
s = length(A);
cc=1;
for j = 1:(s-1)
    for i = s:-1:j+1
        m = A(i,j)/A(j,j);
        m_arr(cc)="m"+i+j+"="+m;
        cc=cc+1;
        A(i,:) = A(i,:) - m*A(j,:);
        b(i) = b(i) - m*b(j);
        disp([A b]);
    end
end
disp(m_arr);
x = zeros(s,1);
x(s) = b(s)/A(s,s);               
for i = s-1:-1:1                    
    sum = 0;
    for j = s:-1:i+1                
        sum = sum + A(i,j)*x(j);    
    end 
    x(i) = (b(i)- sum)/A(i,i);
end