function f=find(Q,data,adj,open)
m=15-length(open);
adj=adj~=0;
tempopen=open;
s=0;
while m>0
    weijianods=sum(adj(tempopen,:))>0; %%%%%%  high
    for i=1:length(tempopen)
        wejianods(tempopen(i))=0
    end
    n=m;
    for i=1:15
        if weijianods(i)~=0
            weijadress(n)=i;
        end
    end
    high=min(data(weijiadress,1))+1.2;   %%%%%%%   high
    for i=i:length(tempopen)
        if (high-data(tempopen(i),1))>0
            s=s+(high-data(tempopen(i),1))*high;
        end
    end
    if Q-s<0
        f=tempopen;
        return;
    else
        nextdress=weijadress(1);
        for i=1:length(weijadress)
            if data(nextdress,4)
            
    m=m-1;
end