I'm trying to find the roots $\omega_n$ of, \begin{equation} \tan(\omega) = \frac{2b\omega}{b^2\omega^2 -1} \end{equation}
I know it must be done numerically, so I would use the bisection method. In order to use that I must know the interval to search for each root for any $b \in (0,1)$. The problem is that there is different behaviour of the singularities at $b>1/2$. Does anyone know a way to solve the roots of this?