Differential equation for a ratio of consecutive Bessel functions - MathOverflow most recent 30 from http://mathoverflow.net 2013-05-25T07:02:43Z http://mathoverflow.net/feeds/question/36350 http://www.creativecommons.org/licenses/by-nc/2.5/rdf http://mathoverflow.net/questions/36350/differential-equation-for-a-ratio-of-consecutive-bessel-functions Differential equation for a ratio of consecutive Bessel functions J. M. 2010-08-22T02:46:03Z 2011-12-02T02:24:51Z <p>My attempts to search via Google seem to be failing, so I thought of asking here.</p> <p>All the derivatives of the function</p> <p><code>$r_n(z):=\frac{J_n(z)}{J_{n-1}(z)}$</code></p> <p>where <code>$J_n(z)$</code> is the Bessel function of the first kind are expressible in terms of <code>$r_n(z)$</code>, for instance <code>$\frac{\mathrm{d}}{\mathrm{d}z}r_n(z)=r_n(z)^2-\frac{2n-1}{z}r_n(z)+1$</code> . I've been trying to derive a (linear?) differential equation (hopefully just second-order) that might be satisfied by <code>$r_n(z)$</code>, but my manipulative ability does not seem to be up to snuff.</p> <p>Probably my problem can be resolved in two ways:</p> <ol> <li>Are there any papers where generating functions/differential equations of ratios of Bessel functions have been studied? ; or</li> <li>How can I derive a differential equation for <code>$r_n(z)$</code> with the knowledge that all higher derivatives are expressible in terms of <code>$r_n(z)$</code>?</li> </ol> <p>(On the other hand, the <em>difference</em> equation for <code>$r_n(z)$</code> (and thus its continued fraction representation) is easily derived, so no problem for me there.)</p> <p>I will be interested in any input. Thanks!</p> <p><strong>EDIT:</strong></p> <p>Per Pietro's request, I now tip my hand and reveal my reason for interest: I saw <a href="http://dx.doi.org/10.1016/S0010-4655%2898%2900193-3" rel="nofollow">this paper</a> many years ago on a neat method for computing the first few roots of the Bessel function of the first kind. Some time later, I came across J.F. Traub's "Iterative Methods for the Solution of Nonlinear Equations", where he shows the <a href="http://books.google.com/books?id=se3YdgFgz4YC&amp;pg=PA78" rel="nofollow">construction of iteration functions involving derivatives</a> and can be constructed to have quadratic, cubic... convergence. (Newton's method is but the first member of this family). I also came across this <a href="http://www.ams.org/journals/mcom/1958-12-061/S0025-5718-1958-0099752-6/S0025-5718-1958-0099752-6.pdf" rel="nofollow">short note</a> by D.J. Hofsommer on how one might profitably exploit the methods derived by Traub if the function of interest satisfies a simple differential equation (Essentially, one just constructs the Newton correction <code>$u=\frac{f(z)}{f^{\prime}(z)}$</code>, and the high-order iteration functions are merely a series in powers of <code>$u$</code>). That got me wondering on how one might recursively generate iteration functions with increasing order of convergence for the case of finding the roots of the Bessel function. (On another note, I was able to successfully use the ideas of Traub and Hofsommer for the generation of Gaussian quadrature rules, e.g. Legendre, Lobatto, Radau, and was hoping things might be just as successful for Bessel function root-finding).</p> http://mathoverflow.net/questions/36350/differential-equation-for-a-ratio-of-consecutive-bessel-functions/36372#36372 Answer by Gerald Edgar for Differential equation for a ratio of consecutive Bessel functions Gerald Edgar 2010-08-22T12:37:11Z 2010-08-22T12:37:11Z <p>As noted by Pietro, the nonlinear differential equation you gave is a Riccati equation. Next, there is a standard method (due to Ince?) to convert a Riccati equation into a second-order linear equation.</p> <p><a href="http://en.wikipedia.org/wiki/Riccati_equation" rel="nofollow">http://en.wikipedia.org/wiki/Riccati_equation</a></p> <p>Try it, and see if that is what you want.</p> http://mathoverflow.net/questions/36350/differential-equation-for-a-ratio-of-consecutive-bessel-functions/36374#36374 Answer by Jacques Carette for Differential equation for a ratio of consecutive Bessel functions Jacques Carette 2010-08-22T13:32:42Z 2010-08-22T13:32:42Z <p>No, you will not find a linear ordinary differential equation (with polynomial coefficients) for $r_n$. This is because $J_{n-1}$ has infinitely many zeroes which are not cancelled out by the zeroes of $J_n$, so that $r_n$ has infinitely many poles. Holonomic functions, aka solutions of LODEs with polynomial coefficients, can only have a finite number of singularities.</p> <p>The transformation between LODEs and Riccati equations involves transformations between an equation for a quantity $y(z)$ [your $r_n$] and an LODE for $-u'/u$. So what you'll get is indeed an LODE for a single Bessel $J$ function!</p> http://mathoverflow.net/questions/36350/differential-equation-for-a-ratio-of-consecutive-bessel-functions/82425#82425 Answer by Chad for Differential equation for a ratio of consecutive Bessel functions Chad 2011-12-02T02:24:51Z 2011-12-02T02:24:51Z <p>I do not know how helpful this would be to you however it was very helpful to understand the physics and numerics of Bessel. if you are studying elastic wave propagation. The solution of the differential equations of potential wave is the cylindrical Bessel:</p> <p>$r^2 \frac{d^2 R}{dr^2} + r \frac{dR}{dr} + (r^2 - \alpha^2)R = f(r)$</p> <p>for an arbitrary real integer number α (the '''order''' of the Bessel function). In solving problems in cylindrical coordinate systems, Bessel functions are of integer order (α = ''n''). Since this is a second-order differential equation, there must be two [[linearly independent]] solutions. My solutions use Bessel J(n,.) and Hankel H(n,.) (as previously mentioned)</p> <p>The potential is assumed for each media to be:</p> <p>$\phi=\left(a_{1}J_{n}(K*r)+a_{2}*H_{n}(K*r)\right)*e^{in\theta} ,$</p> <p>$\psi(t) = \left(a_{3}J_{n}(k*r)+a_{4}*H_{n}(k*r)\right)*e^{in\theta} \,$</p> <p>However for numerical stability: They (ref. 1) normalize the potential for each layer and for each nth iteration The potential will have Hankel function equal to 1 at the inner radius, while Bessel J will be multiplied by Hankel at outer radius. </p> <p>$\phi=\left(a_{1}J_{n}(K*r)*H_{n}(K*r_{out})+a_{2}*\frac{H_{n}(K*r)}{H_{n}(K*r_{in})}\right)*e^{in\theta},$</p> <p>$\psi(t)=\left(a_{3}J_{n}(k*r)*H_{n}(k*r_{out})+a_{4}*\frac{H_{n}(k*r)}{H_{n}(K*r_{in})}\right)*e^{in\theta},$</p> <p>For more details: <a href="http://mathoverflow.net/questions/82303/bessel-functions-in-wave-propagation-and-scattering" rel="nofollow">http://mathoverflow.net/questions/82303/bessel-functions-in-wave-propagation-and-scattering</a> Reference: David C. Ricks and Henrik Schmidt, "A numerically stable global matrix method for cylindrically layered shells excited by ring forces" 1994</p>