%*************** theoretical figure **********************
ebn0=1:0.1:10;
ber=zeros(1,91);
for ii=1:91
    ber(ii)=0.5*erfc(sqrt(10.^(ebn0(ii)/10)));
end
hold on
plot(ebn0,ber,'-k') ;