rce.m

Matlab的Optimization软件包

function E=rce(A,r,c,p,range)
if nargin==3,p=1;range=size(A,2);end
if nargin==4,if p==1;range=size(A,2);else range=size(A,1);end,end
if p==2
    A=A';
    t=r;r=c;c=t;
end
A(:,c)=A(:,c)/A(r,c);
for i=setdiff(1:range,c)
    A(:,i)=A(:,i)-A(:,c)*A(r,i);
end
if p==2
    A=A';
end
E=A;