ps = linspace( 0, 1, 100 ); 

rho1 = 10-9*ps;
rho2 = 5*ps+1; 
rho3 = 13-13*ps;
rho4 = 8-6*ps;
rho5 = 7*ps;
rho6 = 5-2*ps; 
rho7 = 4*ones(1,length(ps));

if( 0 ) 
figure; hold on; 
for ri=1:7,
  plot( ps, eval(['rho',num2str(ri)]) ); 
end
end

RHO = [ rho1; rho2; rho3; rho4; rho5; rho6; rho7 ]; 
figure; plot( ps, RHO.' ); 
xlabel( '\xi' ); ylabel( 'Bayes risk' );
title( 'Bayes risk (for each decsion)' );
saveas( gcf, '../Graphics/chap_8_prob_8.eps', 'epsc' );