将矩阵A高斯顺序消元成上三角矩阵，但是出不来结果，哪条语句出问题了

function A=gauss(A)
[n n]=size(A);
for k=1:n
    for i=k+1:n
        A(k,i)=A(k,i)/A(k,k);
    end
    for j=k+1:n
        for i=k+1:n
            A(j,i)=A(j,i)-A(k,i)*A(j,i);
        end
    end
end