solve d/df p*log(1 + f*b) + q * log(1 - f) = 0 for f