mygaussx.m

来自「线性规划的程序」· M 代码 · 共 10 行

function A=mygaussx(A,r,c)
%以A(r,c)为轴，对列高斯消去
m=min(size(A));                         %取循环次数，行数
A(r,:)=A(r,:)./A(r,c);
for i=1:m
    if i~=r
        A(i,:)=A(i,:)-A(i,c).*A(r,:);
    end
end