Gauss_Jordan 소거법

function x=gauss_jordan2(A)

% A: 계수행렬

% I: 단위행렬

% x: inv(A)의 계산값

[m1,k1]=size(A);% [m2,k2]=size(b); x=zeros(m1,k2);

I=eye(m1);

% 상삼각 행력 형성

for i=1:m1-1

    m=-A(i+1:m1,i)/A(i,i);   

    A(i+1:m1,:)=A(i+1:m1,:)+m*A(i,:);

    I(i+1:m1,:)=I(i+1:m1,:)+m*I(i,:);

end

%A

%I

% 하삼각 행렬 형성

for i=1:m1-1

    for j=i+1:k1

    m=-A(i,j)/A(j,j);   

    A(i,:)=A(i,:)+m*A(j,:);

    I(i,:)=I(i,:)+m*I(j,:);

    end

    %A

    %I

end

%A

%I

for i=1:m1

    x(i,:)=I(i,:)/A(i,i);

    A(i,i)=A(i,i)/A(i,i);

end