MATLAB实验一 解线性方程组的直接法(4)

A(i,j)=A(i,j)-m*A(k,j); end

b(i)=b(i)-m*b(k); end

det=det*A(k,k); end

det=det*A(n,n); % 回代过程

if abs(A(n,n))<1e-20 index=0; return; end

for k=n:-1:1 for j=k+1:n

b(k)=b(k)-A(k,j)*x(j); end

x(k)=b(k)/A(k,k); end