invers matriks program matlab

clc;

clear;

A= input ('matriks awal =');
disp('matriks awal =');
disp(A);
[n m]=size (A);
I= eye(n);
disp('proses invers')
for k=1:n
    dummy=A(k, k);
    for j=1:n
        A(k,j)=A(k, j)/dummy;
        I(k,j)=I(k, j)/dummy;
    end
    for i=1:n
        if i~=k
            faktor=A(i, k);
            for j=1:n
                A(i,j)=A(i,j)+(-faktor*A(k,j));
                I(i,j)=I(i,j)+(-faktor*I(k,j));
            end
        end
        P=[A I];
        disp(P)
    end
end