clc;clear;
a=[0.25,0.25,0.20,0.15,0.10,0.05];%初始概率，也可為輸入
k=length(a);y=0;
s=zeros (k,1);b=zeros(k,1);p=zeros(8,1);
for m=1:k
   s(m)=y;
   y=y+a(m);                     %求和，
   b(m)=ceil(-log2(a(m)));       %取整，求位數
   z=zeros(b(m),1);
   x=s(m);
   p=b2d10 (x);                  %調用子函數
   
for r=1:b(m)
    z(r)=p(r);
   end
  
  disp('輸出結果為:')            %結果顯示
  disp('初始概率'),  disp(a(m))  
  disp('求和結果'),  disp(s(m))
  disp('編碼位數'),  disp(b(m)) 
  disp('最終編碼'),  disp(z') 

end