% MATLAB script for Illustrative Problem 6, Chapter 9.
echo on
rho_b1=0:5:35;	           		% rho in dB for the simulated error rate 
rho_b2=0:0.1:35; 	   		% rho in dB for theoretical error rate computation
for i=1:length(rho_b1),
  smld_err_prb(i)=ss_pe96(rho_b1(i));	% simulated error rate
end;
for i=1:length(rho_b2),
  temp=10^(rho_b2(i)/10);	
  if (temp>2)
    theo_err_rate(i)=1/(exp(1)*temp);	% theoretical error rate if rho>2
  else
    theo_err_rate(i)=(1/2)*exp(-temp/2);% theoretical error rate if rho<2
  end;
end;
% Plotting command follow