function [out]=Puncturing(datain,length,rate);
temp=reshape(datain,3,(length/3));
if isequal(rate==[0 0 0 0 0])
     R=1/3;
elseif isequal(rate,[0 0 0 0 1])
     R=1/2;
elseif isequal(rate,[0 0 0 1 0])
     R=1/3;
elseif isequal(rate,[0 0 0 1 1])
     R=1/2;
elseif isequal(rate,[0 0 1 0 0])
     R=5/8;
elseif isequal(rate[0 0 1 0 1])
     R=1/2;
elseif isequal(rate==[0 0 1 1 0])
     R=5/8;
elseif isequal(rate==[0 0 1 1 1])
     R=3/4;
end


if(R==1/3)
    dout=datain;
elseif(R==1/2)
    dou=[temp(1,:);temp(3,:)];
    dout=reshape(dou,1,(length*3/2));
elseif(R==5/8)
    for i=1:(length/15)
        dout((i-1)*8+1)=datain((i-1)*15+1);
        dout((i-1)*8+2)=datain((i-1)*15+2);
        dout((i-1)*8+3)=datain((i-1)*15+6);
        dout((i-1)*8+4)=datain((i-1)*15+7);
        dout((i-1)*8+5)=datain((i-1)*15+8);
        dout((i-1)*8+6)=datain((i-1)*15+12);
        dout((i-1)*8+7)=datain((i-1)*15+13);
        dout((i-1)*8+8)=datain((i-1)*15+14);
    end
elseif(R==3/4)
    for i=1:(length/9)
        dout((i-1)*4+1)=datain((i-1)*9+1);
        dout((i-1)*4+2)=datain((i-1)*9+2);
        dout((i-1)*4+3)=datain((i-1)*9+6);
        dout((i-1)*4+4)=datain((i-1)*9+9);
    end
end
out=dout;