%  Log-Normal distribution for figure 6.31


sig=.4;
p=.01:.01:8;
p0=2;
pdf=exp(-(log10(p/p0)).^2/(2*sig^2))./(p*sig*sqrt(2*pi));
figure;
plot(p,pdf,'k')
xlabel('Power, p')
ylabel('pdf  p(p)')
set(gca,'xtick',[0 1 2 3 4 5 6 7 8])