%hartsfield 
%
clear all
close all


 disp('Hartsfield Dirichlet Example')
  lw = 3; 
  set(0, 'DefaultAxesFontSize', 16);
  fs = 15;
  msize = 10;
  

% need hartsfield.mat
% This mat file has two columns first 1-40000 40001-80000
% and the second lambda(1:40000) P[10](1:40000).
load hartsfield
p10=a(1:40000,:);
lambda=a(40001:80000,:);
figure(1)
hist(lambda(:,2), 100)
print -deps 'C:\NPBook\NPBayes\hartsfield1.eps'
figure(2)
hist(p10(:,2), 100)
print -deps 'C:\NPBook\NPBayes\hartsfield2.eps'