%本程序的功能是形成因子表
A=input('請輸入矩陣：A=');
[n,m]=size(A);
for i=1:n
    A(i,i)=1./A(i,i);
    for j=i+1:n
        A(i,j)=A(i,j)*A(i,i);
    end
    for k=i+1:n
        for j=i+1:n
            A(k,j)=A(k,j)-A(k,i)*A(i,j);
        end
    end
end
disp('矩陣A的因子表為：A=');
disp(A)