function x_0=my_sample(vahy,x_1)
N=size(vahy,2);
F= (vahy*triu(ones(N,N),0));
u=rand(1,N)'; 
index=min(N-sum((ones(1,N)'*F-u*ones(1,N))>0,2)+1,N); 
x_0=x_1(:,index);
end