test.m

error = [0.238095238095238	0.342465753424658	0.342465753424658	0.280898876404494	0.238095238095238	0.206611570247934	0.206611570247934	0.163398692810458	0.100401606425703	0.0798722044728435	0.0511247443762781	0.0291715285880980	0.00581801256690714	0.00211846453690365	0.000329185594838370	4.80332313107500e-05	4.01244371143339e-06	1.79584251682300e-07	0	0	0];
EbN0dB = -20:0;
semilogy(EbN0dB, errors, 'Color', 'Red'); 
hold on
ylim([10^(-7) ,1])
semilogy(-20:20, qfunc(sqrt(2*db2pow(-20:20))), 'Color', 'Blue')


hold off

xlabel("Eb/N0");
ylabel("BER");

title("Bit Error Rate x Eb/N0");
legend({'BER', 'Theoretical error'})