%
% Written by:
% -- 
% John L. Weatherwax                2008-02-20
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

p = [ 0.25, 0.5, 0.75 ]; 

ns  = 10.^[ 2:7 ]; 
for ni=1:length(ns),
  n  = ns(ni); 
  x  = randn(1,n); 
  q1 = sample_quantiles( x, p ); 
  q2 = sample_quantiles_linInterp( x, p ); 
  fprintf('for n=%10d...\n',n); 
  [q1(:).'; q2(:).']
end