% c6exp.m - Chapter 6: exponential distribution illustration
%
help c6exp
disp('This M-file illustrates 3 probability distributions...')

      rand('seed',0); n=20000;% Reset the 'seed' value
      m=1.2; s=1.2;           % mean and std
      x = unirand(0,2.4,n,1); % Uniform in (0,2.4)

     randn('seed',0);         % Reset the 'seed' value
      y = normrand(m,s,n,1);  % Normal N(1.2,1.2)

      rand('seed',0);         % Reset the 'seed' value
      z = exprand(m,n,1);     % Exponential

      subplot(311); hist(x,90)% Plot x
      axis([0 2.4 0 225])     % Cut the right half
      hold on; plot([m m],[0 225],':r'); hold off

      subplot(312); hist(y,20)% Plot y
      hold on; plot([m m],[0 4000],':r'); hold off

      subplot(313); hist(z,40)% Plot z
      axis([0 10 0 6852])     % Cut the right half
      hold on; plot([m m],[0 6852],':r'); hold off

%%% The following lines are for "cumsum" %%%
t=cumsum(z); 
fprintf('For n = %d, t first and last = %10.4f, %10.4f\n',n,t(1),t(n));
format compact; tm=hist(z,40); sum(tm);tm=max(tm)
