matlab中将任意矩阵转换成上三角矩阵的源码

function B=gauss(A)
B=A;
[m,n]=size(A);
for i=1:min([m,n])
j=find(B(i:end,i),1);
if isempty(j)
break;
else
B([i,i+j-1],:)=B([i+j-1,i],:);
B(i+1:end,i:end)=B(i+1:end,i:end)-ones(m-i,1)*B(i,i:end)/B(i,i).*(B(i+1:end,i)*ones(1,n+1-i));
end
end