gauss.m

MEXP100-matlab源码

function u=gauss(a,n)
for k=1:n-1
    for i=k+1:n
        for j=k+1:n+1
            if abs(a(k,k))>1e-6
                a(i,j)=a(i,j)-a(i,k)/a(k,k)*a(k,j);
            else
                msgbox('Gauss消去失败','计算结果');
                pause;
                exit;
            end
        end
    end
end
x(n)=a(n,n+1)/a(n,n);
for i=n-1:-1:1
    s=0;
    for j=i+1:n
        s=s+a(i,j)*x(j);
    end
    x(i)=(a(i,n+1)-s)/a(i,i);
end
u=x;