function out_i=interleaving(in,N_cbps,N_bpsc)                                                                                                                                                                                                                                                                    
 s=max(N_bpsc/2,1);
 for k=1:1:N_cbps
     out_i((N_cbps/16)*mod(k-1,16)+floor((k-1)/16)+1)=in(k);
 end
 for i=1:1:N_cbps
     out_j(s*floor((i-1)/s)+(mod(i-1+N_cbps-floor(16*(i-1)/N_cbps),s))+1)=out_i(i);
 end