close all
clear all

x= -3:0.002:3;
y = normpdf(x);
figure(1) %----------------------------------------------------------------
plot(x,y, 'k-','linewidth',3)
plot([-3 3],[0 0],'k-')
axis([-3 3 -0.2 0.45])
xx = -2:0.002:2;
yy = normpdf(xx);
hold on
fill([-2 xx 2], [0 yy 0], 'y')
plot(x,y, 'k-','linewidth',3)
set(gca,'XTick',[],'YTick',[])
%-----------------marker dots-----------------------------------------
plot([0 0],[0 0.4],'k:','linewidth',2)
plot(0,0,'-mo', 'LineWidth',2,'MarkerEdgeColor','k',...
                'MarkerFaceColor',[.49 1 .63],'MarkerSize',10)
 plot(-2,0,'-mo', 'LineWidth',2,'MarkerEdgeColor','k',...
                'MarkerFaceColor',[.49 1 .63],'MarkerSize',10)  
 plot(2,0,'-mo', 'LineWidth',2,'MarkerEdgeColor','k',...
                'MarkerFaceColor',[.49 1 .63],'MarkerSize',10)       
%--------------------------------------------------------------------------
% text(0.05,0.03,'\fontsize{25} \mu')
% text(-1.05,0.02,'\fontsize{25} a')
% text(1.95,0.02,'\fontsize{25} b')

text('Interpreter','latex',...
	'String','$$\mu$$',...
	'Position',[0.15 .03],...
	'FontSize',25)
text('Interpreter','latex',...
	'String','$$\mu - z_{1-\alpha/2} \frac{\sigma}{\sqrt{n}}$$',...
	'Position',[-2.9   -.07],...
	'FontSize',23)

  xw = 0.36:0.02:0.65; yw = normpdf(xw); %whiten!
  plot(xw, yw, 'w-', 'linewidth',4)

text('Interpreter','latex',...
	'String','$$\mu + z_{1-\alpha/2} \frac{\sigma}{\sqrt{n}}$$',...
	'Position',[0.7  -.07],...
	'FontSize',23)

text('Interpreter','latex',...
	'String','$$\bar X \sim {\cal N}\left(\mu,\frac{\sigma^2}{n}\right)$$',...
	'Position',[0.2 .34],...
	'FontSize',20)

text('Interpreter','latex',...
	'String','$$1 - \alpha$$',...
	'Position',[-0.6  0.2],...
	'FontSize',30)


%print -depsc 'C:\Springer\Estim\Estimeps\ci01.eps'

