% LMFig070301_4th.m
xmax=6.5;
r=0:0.01:xmax;
y=fpdf(r,3,5);
r95=finv(0.95,3,5);
plot(r,y,'--r','LineWidth',3)
axis([0 xmax 0 0.71])
ax1=gca;
xlabel('r','FontSize',20),
ylabel('Probability Density','FontSize',20);
set(ax1,'xtick',[0:5 r95 6],'FontSize',18);
set(ax1,'ytick',0:0.2:0.7,'FontSize',14)
hold on;
xf=r95:0.01:xmax;
yf=fpdf(xf,3,5);
fill([r95 xf xmax]',[0 yf 0]',[.8 .8 1])
text(r95+0.05,0.05,'Area=0.05','FontSize',18);
title('Figure 7.3.1','FontSize',22)
figure(gcf);pause
hold off;