All published worksheets from http://sagenb.org
p=[1/8,1/4,1/2,1/23,1/23,7/184];H=0 for i in srange(len(p)): H=H+(-p[i]*log(p[i],2)) print H.n(digits=5)