load 1.dat;
Pf=X1;
len=length(Pf);
sum=0;
temp=0;
k=6;
Tt=0.01;
%for n=10:1:20
n=8;
Nmax=2;
 for j=1:1:len
        sum=0;
        for l=k:1:n
            temp=nchoosek(n,l)*((1-Pf(j))^l)*(Pf(j)^(n-l));%fec 
            sum=sum+temp ;
        end
     Pt(j)=sum;%全部正接受的概率   
     Pf_fec(j)=exp(log(1-Pt(j))/(Nmax+1));%frame loss rate for fec
     Pw_fec(j)=1-Pt(j); %packet loss rate    
     tRTT(j)=Tt*(1-Pf_fec(j)^(Nmax+1))/(1-Pf_fec(j));%delay for fec
     divide=(tRTT(j)*sqrt(2*Pw_fec(j)/3)+4*tRTT(j)*(3*sqrt(3*Pw_fec(j)/8))*Pw_fec(j)*(1+32*Pw_fec(j)*Pw_fec(j)));
     R=1.0/divide;
     Rcur(j)=R/2;
    end
fid = fopen('2.dat','w');
for o=1:1:len   
  %fprintf(fid,'%f\n',Rcur(o));
  fprintf(fid,'%f\n',tRTT(o));
end
fclose(fid);