My attempts to search via Google seem to be failing, so I thought of asking here.


All the derivatives of the function


$r\_n(z):=\frac{J\_n(z)}{J\_{n-1}(z)}$


where $J\_n(z)$ is the Bessel function of the first kind are expressible in terms of $r\_n(z)$, for instance $\frac{\mathrm{d}}{\mathrm{d}z}r\_n(z)=r\_n(z)^2-\frac{2n-1}{z}r\_n(z)+1$ . I've been trying to derive a (linear?) differential equation (hopefully just second-order) that might be satisfied by $r\_n(z)$, but my manipulative ability does not seem to be up to snuff.


Probably my problem can be resolved in two ways:


1. Are there any papers where generating functions/differential equations of ratios of Bessel functions have been studied? ; or
2. How can I derive a differential equation for $r\_n(z)$ with the knowledge that all higher derivatives are expressible in terms of $r\_n(z)$?


(On the other hand, the *difference* equation for $r\_n(z)$ (and thus its continued fraction representation) is easily derived, so no problem for me there.)


I will be interested in any input. Thanks!


**EDIT:**


Per Pietro's request, I now tip my hand and reveal my reason for interest: I saw [this paper](http://dx.doi.org/10.1016/S0010-4655%2898%2900193-3) 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 [construction of iteration functions involving derivatives](http://books.google.com/books?id=se3YdgFgz4YC&pg=PA78) 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 [short note](http://www.ams.org/journals/mcom/1958-12-061/S0025-5718-1958-0099752-6/S0025-5718-1958-0099752-6.pdf) 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 $u=\frac{f(z)}{f^{\prime}(z)}$, and the high-order iteration functions are merely a series in powers of $u$). 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).



Differential equation for a ratio of consecutive Bessel functions
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.


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!


