rref2.m

来自「有兴趣学系matlab的朋友」· M 代码 · 共 27 行

function y=rref2(A)
[m n]=size(A);
i=1;
while i<=m
    k=1;
    while k<=n
        if any(A(i:m,k)')==0
            k=k+1;
        else
            break;
        end
    end
    if k==n|k==n+1
        break;
    end
    C=exchange(A(i:m,k:n));
    A(i:m,k:n)=C;
    l=1;
    while l<=m
        if A(l,k)~=0 & l~=i
            A(l,k:n)=add(A(i,k:n),A(l,k:n));
        end
       l=l+1;    
    end    
    i=i+1;   
end
y=exchange(A);