load 1.dat;
Pf=X1;
len=length(Pf);
sum=0;
temp=0;
k=6;
%for n=10:1:20
n=8;
Nmax=0;
 for j=1:1:len
        sum=0;
        for l=k:1:n
            temp=nchoosek(n,l)*((1-Pf(j))^l)*(Pf(j)^(n-l));
            sum=sum+temp ;
        end
     Pt(j)=sum; //全部數(shù)據(jù)正確接受的概率  
     Pf_fec(j)=exp(log(1-Pt(j))/(Nmax+1));//frame loss rate
     Pw_fec(j)=1-Pt(j); //packet loss rate    
    end
fid = fopen('2.dat','w');
for o=1:1:len   
  fprintf(fid,'%f\n',Pf_fec(o));
  %fprintf(fid,'%f\n',Pw(o));
end
fclose(fid);