function [D] = qhls(A )
A=input('请输入一个矩阵')
[r,c]=size(A);n=length(A);q=zeros(1,n);D=0;E=1;
if r~=c
disp('请注意：因为行数不等于列数，无法计算矩阵行列式，请重新输入！')
return
end
s=1;
for p=1:n-1
k=p;
y=A(k,p);
while ((y==0)&&(k<=n))
k=k+1;
y=A(k,p);
end
if y==0
D=0
return;
else
if (k~=p)
q=A(p,:);
A(p,:)=A(k,:);A(k,:)=q;
s=-s;
end
end
for u=p+1:n
m=A(u,p)/A(p,p);
A(u,p:n)=A(u,p:n)-m*A(p,p:n);
end
end
for e=1:n
E=E*(A(e,e) );
end
D=s*E;
fprintf(' 矩阵的行列式为：end
D=%3.2f\\n',D);
因篇幅问题不能全部显示，请点此查看更多更全内容