function bpskerr_rate=trng_gauss
close all;
clear all;
clc
echo on
SNRindb1=0:2:12;
SNRindb2=0:2:12;
for i=1:length(SNRindb1)
    smld_err_prb(i)=cm_bpske(SNRindb1(i));   %simulated error rate
     echo off;
end
echo on
for i=1:length(SNRindb2)
    SNR=exp(SNRindb2(i)*log(10)/10);         %signal to noise ratio
    theo_err_prb(i)=erfc(sqrt(SNR/2))/2;       %theoretical symbol error rate
    echo off;
end
echo on

%plotting commands follows
semilogy(SNRindb1,smld_err_prb,'*');
xlabel('SNR(dB)---------->');
ylabel('BER-------------->');
hold on
semilogy(SNRindb2,theo_err_prb);
xlabel('SNR(dB)---------->');
ylabel('BER-------------->');
legend(' simulated result',' theoretical result');
hold off