clear;
clc;
rho_in_dB=42;
BPH=2;
number_of_states=16;
tim=125;
pb=0;
p=zeros(1,tim);
for i=1:tim
    %p(i)=log_map_simu4(rho_in_dB,BPH,number_of_states);
    %p(i)=va_demod(rho_in_dB,BPH,number_of_states);
    %p(i)=log_map_rayleigh(rho_in_dB,BPH,number_of_states);
    %p(i)=va_demod_rayleigh(rho_in_dB,BPH,number_of_states);
    p(i)=symbysym(rho_in_dB,BPH,number_of_states);
    pb=pb+p(i);
    i
end
pb=pb/tim;
sprintf('pb=%f',pb)