function sub4(A,i,j,s)
x=0;
y=0;
for k=1:j-1
h=i;
f=i;
B=[A(:,1:k-1);A(:,k+1:j-s)]%不包含结果属性
C=[A(:,1:k-1);A(:,k+1:end)]%包含结果属性
%判断矩阵协调行性
for m=1:i-1%第一层循环
for n=m+1:i%第二层循环
if B(m,:)/B(n,:)==1
h=h-1;
x=h;
end
if C(m,:)/C(n,:)==1
f=f-1;
y=f;
end
end %第二层循环结束
end%第一层循环结束
if x==y
A(:,k);
end
end