plot p * log(p/q) + (1-p) * log((1-p)/(1-q)) for 0<p<1 and 0<q<1