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