function [W,IX]=SimilarSearch(D,y,sigma,k)
T=zeros(1,size(D,2));
mu= mean(y)./mean(D);
for i=1:size(D,2)
    x = y-mu(i)*D(:,i);
    a=x'*x;
    b=mean(x);
    c=var(x)-sigma^2;
    T(1,i)=a+abs(b)+abs(c);
end
size(T)
[W,In]= sort(T)
W = W(1:k)
IX=In(1:k);
%G=D(:,IX);