solve -p log(16) - p log p + (1+p) log(1+p) = 0 for p