function j = qrselect(t,s)
%
%  The function QRSELECT returns the index to the element in 
%  the vector S closest to the scalar T.
%
   [a,j] = min(abs(t-s));
end

