p1=[0.23003,0.14465,0.06858,0.02236,0.00449,0.000745,0.000088];
p2=[0.24111,0.154130,0.068910,0.02254,0.0053,0.000875,0.000141];
r=[0,1,2,3,4,5,6];
r1=[0:0.1:6];
for i=1:length(r1)
    p3(i)=log_map_bound(2,16,r1(i));
end
r2=[5:11];
for i=1:length(r2)
    p4(i)=awgn_theo(64,2,r2(i));
end