%---
%  Q distribution
%
%---

EbNodB = [-2:0.01:10];

EbNo = 10 .^(EbNodB/10);

BER = 0.5*erfc(sqrt(EbNo));

semilogy(EbNodB, BER);
grid
xlabel('Eb/No (dB)')
ylabel('Bit Error Rate')
axis([-2 10 1E-6 1])