mu=3;sigma=0.5;					%正態分布參數設定
x=mu+sigma*[-3:-1,1:3];yf=normcdf(x,mu,sigma);
P=[yf(4)-yf(3),yf(5)-yf(2),yf(6)-yf(1)];%計算 
xd=1:0.1:5;yd=normpdf(xd,mu,sigma);		%計算概率密度函數，供圖示。
%為各區域填色而進行的計算
for k=1:3
	xx{k}=x(4-k):sigma/10:x(3+k);		%用元胞數組存放采樣數不同的數據
	yy{k}=normpdf(xx{k},mu,sigma);	%用元胞數組存放采樣數不同的數據
end  
subplot(1,3,1),plot(xd,yd,'b');hold on
fill([x(3),xx{1},x(4)],[0,yy{1},0],'g')
text(mu-0.5*sigma,0.3,num2str(P(1))),hold off
subplot(1,3,2),plot(xd,yd,'b');hold on
fill([x(2),xx{2},x(5)],[0,yy{2},0],'g')
text(mu-0.5*sigma,0.3,num2str(P(2))),hold off
subplot(1,3,3),plot(xd,yd,'b');hold on
fill([x(1),xx{3},x(6)],[0,yy{3},0],'g')
text(mu-0.5*sigma,0.3,num2str(P(3))),hold off