function plot_r_vs_db
y=[];
for r=0.001:0.001:0.999
    db = bqpsk_approx(r);
    y = [y db];
end
plot(0.001:0.001:0.999, y);