plot -xlog((1-x)/x) - (1-x)log((x)/(1-x)), for x> 0