clear ,clc;
n=4;m=100;
x=rand(m,1);e=0.1;ee=0.01;
t=max(x);
d=t.*ones(m,1);
r=find((d-x)<e);
while size(r,1)<=n-1
    e=e+ee;
    r=find((d-x)<e);
end
s=size(r,1)-n+1;
if s>1
    e=e/10;
   ee=ee/10;
end
rr=find((d-x)<e);
while size(rr,1)<=n-1
    e=e+ee;
    rr=find((d-x)<e);
end
if size(rr,1)<size(r,1)
    x(rr)=1;
else x(r)=1;
end
re=find(x<1);
x(re)=0;
size(r,1),size(rr,1)