function [P_CRC]=mydigui(t,cont,k,i,D,p,P_CRC)
%

    
if i==cont   
for k(i)=k(i-1)+1:n
        if t(k(i))=='1'
            t(k(i))='0';
        else
            t(k(i))='1';
        end
      
            rt=t;
           [q,l,rr,m]=mod2_div(rt,length(rt),p,16);
            RR=bin2dec(rr);
            if RR==0
                P_CRC(cont)=P_CRC(cont)+(pb^cont)*(1-pb)^(n-cont);
            end
end
reture(P_CRC);
elseif i==1
    for k(i)=1:n
        if t(k(i))=='1'
            t(k(i))='0';
        else
            t(k(i))='1';
        end
        P_CRC=mydigui(t,cont,k,i+1,D,p,P_CRC);
    end
else
    for k(i)=k(i-1)+1:n
        if t(k(i))=='1'
            t(k(i))='0';
        else
            t(k(i))='1';
        end
        P_CRC=mydigui(t,cont,k,i+1,D,p,P_CRC);
    
end