plot C log (b/a) + (N-C) log ((1-b)/(1-a)) where a = 0.1, N = 10, C = 4, b between 0 and 1