function main;
echo on
rho_b1=0:5:35;
rho_b2=0:0.1:35;
for i=1:length(rho_b1)
    smld_err_prb(i)=ss_Pe96(rho_b1(i));
end
for i=1:length(rho_b2)
    temp=10^(rho_b2(i)/10);
    if temp>2
        theo_err_rate(i)=1/(exp(1)*temp);
    else
        theo_err_rate(i)=1/2*exp(-temp/2);
    end
end
semilogy(rho_b1,smld_err_prb,'-b*',rho_b2,theo_err_rate,'-r');