function Y=quantile(Q,X)
% quantile - compute quantiles of X
%
% quantile :: 
%    [[L]->0--1] ~'the L quantiles to compute',
%    [[N,M]] ~'M columns of data'
% -> [[L,M]] ~'the L quantiles for each of M columns'.

j=max(1,round(Q*size(X,1)));
S=sort(X);
Y=S(j,:);


