function x=indexcat(y,val);
%indexcat Extract indices for y being equal to val if val is a scaler
%         or extract indices for y being val(1)<y<=val[2] if val is 2x1
if rows(val)==1 & cols(val)==1;
   x=find(y==val);
else;
   x=find((y>val(1)).*(y<val(2)));
end;
return;