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

已解决，谢谢。