4407.sagews

def H(p,q):
    res=-(p*log(p,2) + q*log(q,2) + (1-p-q)*log((1-p-q),2))
    return res

H(0.3,0.3)

1.57095059445467

p, q= var('p q')

W = plot3d(H(p,q),(p,-1,1),(q,-1,1), frame=False)

show(W, figsize=8)