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