%%credset.eps

clear all
close all

set(0, 'DefaultAxesFontSize', 16);
fs = 16;
lw = 2.5;  
msize = 6;
zz = 102.8 + 1.96 * sqrt(48) %end point of CS
lev = normpdf(zz, 102.8, sqrt(48)) %level of cutting line
% lev =0.0084
aa = linspace(70, 130, 300);
plot(aa,normpdf(aa, 102.8, sqrt(48)),'k-','linewidth',3)
hold on
plot(aa, lev * ones(size(aa)), 'b--','linewidth',lw)
%--------------------------
plot([98 - norminv(0.975)*sqrt(98) 98 + norminv(0.975)*sqrt(98)],[0,0],'r-',...
    'LineWidth',40)
   %78.5973  117.4027
plot([102.8 - norminv(0.975) * sqrt(48),102.8 + norminv(0.975) * sqrt(48)],[0,0], 'b-','LineWidth',30)


  print -depsc 'C:\Springer\Bayes\Bayeseps\credset.eps'

