function CumProb(N)
% Example: CumProb(10000)

x=randn(N,1);
z=hist(x,20)/N;
y=cumsum(z);

bar(z)
grid on
title('PDF')

figure
bar(y)
grid on
axis([0 21 0 1.2])
title('CDF')