end
[p1,i1]=min(p)
index=[(1:i1-1),(i1+1:n)] p=p(index) n=n-1
[p2,i2]=min(p)
index2=[(1:i2-1),(i2+1:n)] p=p(index2);
i2=index(i2)
index=index(index2) p(n)=p1+p2 c=huffman(p)
c{n+1}=strcat(c{n},'1') c{n}=strcat(c{n},'0') index=[index,i1,i2] c(index)=c