function [p]=smldPe(snr_in_dB)
E=1;
SNR=exp(snr_in_dB*log(10)/10);
sgma=E/sqrt(2*SNR);
N=10000;
for i=1:N
temp=rand;
if(temp<0.5)
dsource(i)=0;
else
dsource(i)=1;
end
end
numoferr=0;
for i=1:N
if(dsource(i)==0)
r0=E+gngauss(sgma);
r1=gngauss(sgma);
else
r0=gngauss(sgma);
r1=E+gngauss(sgma);
end
if(r0>r1)
decis=0;
else
decis=1;
end
if(decis~=dsource(i))
numoferr=numoferr+1;
end
end
p=numoferr/N;
end