%Metoda Gauss-Jordan
function X=Gauss_Jordan(A)
[m n]=size(A);

X= eye(n);

%triunghiularizare
for p=1:n-1
    for i=p+1:n
        t=A(i,p)/A(p,p);
        A(i,p)=0;
        A(i,p+1:n)=A(i,p+1:n)-t*A(p,p+1:n);
    end
    X(i,1:n)=X(i,1:n)-t*X(p,1:n);
end
%rezolvare sisteme triunghiulare
for i=n:-1:1
    for  k=1:n
            suma=A(i,i+1:n)*X(i+1:n,k);
            X(i,k)=(X(i,k)-suma)/A(i,i);
    end
end