function [x,pdfx,cdfx]=pdf_cdf_gen(x0, M)
up=max(x0);
dn=min(x0);
N=max(size(x0));
delta=(up-dn)/M;
pdfx=zeros(M+1,1);

for n=1:N
    d0=floor((x0(n)-dn)/delta)+1;
    pdfx(d0)=pdfx(d0)+1;
end
x=dn:delta:up;
pdfx=pdfx/(N*delta);

for m=1:M+1
    cdfx(m)=sum(delta*pdfx(1:m));
end