I am wondering if anyone can help me find an analytical solution to the roots of the following function: $$f(b) = c\log \left( \frac{b}{a} \right) + (n-c)\log \left( \frac{1-b}{1-a} \right),$$ $a,b \in (0,1)$ and $n,c \in \mathbb{Z}$ with $c \le n$.
I have done some googling and I found references to using lambert's w function for finding roots in situations where logs are involved; unfortunately, nothing I've read seems like it will solve my problem. I am admittedly completely new to lambert's w function, so it's very possible that I am missing something.
I have proven $f(b)$ is concave and only has two roots (graphical support for this claim). Clearly one root is $b=a$, the other root is my challenge. The context of the actual problem--which I can expound upon if it would be helpful--leads me to believe that this non-trivial root (r) should be a function of $t = \frac{c}{n}$, i.e., $r = h(t)$.
I am open to any ideas or suggestions, even if they do not exactly solve the problem, because they may provide analytical approximations or bounds on $h$ and/or $h^{-1}$, which is my ultimate goal.
Update
@Claude's suggestion to taylor expand moves us forward because we can get a 2-degree polynomial approximation: $$f(b) = f(b_{*}) + \frac{f''(b_{*})(b-b_{*})^2}{2}+R_2(b),$$ where $R_2(b)$ is the remainder/error. We can then use the quadratic equation to find the roots of this 2-degree polynomial, but can we get a bound on $R_2(b)$.
Why Bound The Remainder (Context of My Problem)
So $f(b)$ is a concave function, with 2 roots: trivially when $b = a$ and another non-trivial root that I will call $b_{\max}$. The function $f$ is maximized at $b_{*} = \frac{c}{n}$, i.e., $b_{*} = \arg\max_b f(b)$. I originally wanted to find $b_{\max}$ because I believe there is some other function $h(b_{*}) = b_{\max}$. I believe from this function I should be able to derive an inequality in the vein of $b_{*} \le \frac{b_{\max}}{q}$, where $q$ is a real number or possibly a function of other variables. Intuitively, $b_{\max}$ cannot be "too close" to $b_{*}$.
By doing the taylor expansion approximation as suggested, we will get an approximation of $f$ (called $f_{\text{aprx}}$). From $f_{\text{aprx}}$, we can get an approximation for $b_{\max}$, (called $b_{\text{aprx}}$). However, we need to bound how far away $b_{\text{aprx}}$ is from $b_{\max}$. This is because once we compute our approximation of $h$ (called $h_{\text{aprx}}$)--from deriving the relationship between $b_{*}$ and $b_{\text{aprx}}$ in $f_{\text{aprx}}$--we then compute $b_{*} \le \frac{b_{\text{aprx}}}{q_\text{aprx}}$. However, I need to get a bound on how far away $q_{\text{aprx}}$ if from $q$, which means I need a bound on how far away $b_{\text{aprx}}$ is from $b_{\max}$.