%% ================================================================
% John Kerl
% kerl at math dot arizona dot edu
% 2007-09-06
%% ================================================================

hold on;

%% ----------------------------------------------------------------

x = [0.0 : 0.5 : 4.0];
%y = [639, 234, 77, 31, 14, 5]/1000;
%y = [610, 252, 74, 33, 18, 9]/1000;
y = [656, 215, 97, 20, 8, 2, 1, 0, 1]/1000;

plot(x,y,'*');

lambda = 2;
x = [0.0 : 0.5 : 4.0];
z = exp(-lambda * x) - exp(-lambda * (x+0.5));
plot(x,z, 'o');

set(gcf, 'PaperPositionMode', 'Manual');
set(gcf, 'PaperPosition', [0 0 5.0 1.0])
axis([-0.2,4.7,0,1]);
print -depsc exphisto
