function [pz,py]=ecdfplot(x,iplot)
%
% computes ECDF of a set of k samples in columns of matrix X 
% and delivers both x values evaluated and ecdf values at these 
% Plots on same frame if iplot=1


[n,r]=size(x); d=1/(n+1); py=[]; pz=[];
lc='brgmcyk';

clf;
for i=1:r
    z=sort(x(:,i)); 
    z=reshape(repmat(z,1,2)',1,2*n)'; 
    y=d*reshape(repmat((0:n)',1,2)',1,2*(n+1))'; y([1 end])=[]; 
	py=[py y]; pz=[pz z]; 
    if (iplot) 
 	plot(z,y,'Color',lc(i)); 
	axis([min(z) max(z) 0 1])
        hold on
    end	
end
hold off
py=unique(py,'rows'); py(1,:)=[]; pz=unique(pz,'rows'); 
%