%  Rayleigh distribution example 6.2

sig=3;
r=[0:1500]/100;
p=r.*exp(-r.^2/(2*sig^2))/sig^2.;
figure;
plot(r,p,'k')
axis([0 15 0 .25])
xlabel('Envelope r (mV)')
ylabel('p(r)')
legend('\sigma=.005')
hold on
rr=[r(500) r(501:1501)];
pp=[0 p(501:1500) 0];
fill(rr,pp,[.8 .8 .8])