3604.sagews

p = var('p')
q = var('q')

Hb(p)= -p.log(2)*p-(1-p).log(2)*(1-p) # Binary Entropy

show (plot(-p.log(2)*p-(1-p).log(2)*(1-p), xmin=0, xmax=1, gridlines=true))

Ht(p,q)= -(p.log(2)*p)-(q.log(2)*q)-(1-p-q).log(2)*(1-p-q) # Tertiary Entropy

show(plot3d(Ht(p,q),(p,0,1),(q,0,1)))