Roots of the equation $I_1(b x) - x I_0(b x) = 0$

91 Views Asked by At

I'm interested in the roots of the equation:

$I_1(bx) - x I_0(bx) = 0$

Where $I_n(x)$ is the modified Bessel function of the first kind and $b$ is real positive constant.

More specifically, I'm interested in the behaviour of the largest non-negative root for $x \geq 0$. I have solved the problem numerically in a crude fashion to see what the function looks like, but as far as I can tell there isn't a closed form.

If the value of the largest root is given by $x(b)$, then based on the numerical solution I think the following is true:

$x(b) < 0 \; $ for $\; b < 2$

$\lim_{b\to\infty} x(b) = 1$

Can anyone give me some pointers to the possible approaches I might take to approximating $x(b)$ for $b>2$ ?

Thanks!

2

There are 2 best solutions below

4
On BEST ANSWER

Notice that the equation also can be rewritten as $$ \frac{I_1(b x)}{I_0(b x)} = x $$ Using the asymptotic series expansion for large $b$ and some fixed $x>0$ we get: $$ 1 - \frac{1}{2 x b} - \frac{3}{8 b^2 x^2} - \frac{1}{8 b^3 x^3} + \mathcal{o}\left(b^{-1}\right) = x $$ resulting in $$ x(b) = 1 - \frac{1}{2 b} - \frac{3}{8 b^2} - \frac{9}{16 b^3} + \mathcal{o}\left(b^{-3}\right) $$ Confirming $\lim_{b \to + \infty} x(b) = 1$.

1
On

I'm a little confused by your notation. Here's what I see:

Assume $b>0$. I can tell you why there is a discontinuity at $b=2$, because

$$I_1(b x) - x I_0(b x) = \left (\frac{b}{2}-1 \right ) x + O(x^3)$$

so when $b=2$, this expression is very close to zero near $x=0$ anyway. Away from $b=2$, I get a plot for the zero like this:

enter image description here

The abcissa is $b$, while the ordinate is the zero for that value of $b$. I do see the limiting behavior you see, i.e., the zero tends to $1$ as $b \rightarrow \infty$.

For $b<2$, $I_1(b x) - x I_0(b x) < 0$ and there are no zeroes.