plot 1 + (1-epsilon) log2(1-epsilon) + epsilon log2(epsilon) for 0 < epsilon < 1/2