gauss.m

数值分析

function [x,det,index]=Gauss(A,b)
[n,m]=size(A);nb=length(b);
if n~=m
    error('A的行数不等于列数');
    return;
end
if m~=nb
    error('A的列数不等于b的行数');
    return;
end
index=1;det=1;x=zeros(n,1);
for k=1:n-1
    a_max=0;
    for i=k:n
        if abs(A(i,k))>a_max
            a_max=abs(A(i,k));r=i;
        end
    end
    if a_max<1e-10
        index=0;return;
    end
    if r>k
        for j=k:n
            z=A(k,j);A(k,j)=A(r,j);A(r,j)=z;
        end
        z=b(k);b(k)=b(r);b(r)=z;det=-det;
    end
    for i=k+1:n
        m=A(i,k)/A(k,k);
        for j=k+1:n
            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-10
    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