Monotonicity Properties for Ratios and Products of Modified Bessel Functions and Sharp Trigonometric Bounds

Let Iν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$I_{\nu }(x)$$\end{document} and Kν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_{\nu }(x)$$\end{document} be the first and second kind modified Bessel functions. It is shown that the nullclines of the Riccati equation satisfied by xαΦi,ν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x^{\alpha } \Phi _{i,\nu }(x)$$\end{document}, i=1,2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i=1,2$$\end{document}, with Φ1,ν=Iν-1(x)/Iν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi _{1,\nu }=I_{\nu -1}(x)/I_{\nu }(x)$$\end{document} and Φ2,ν(x)=-Kν-1(x)/Kν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi _{2,\nu }(x)=-K_{\nu -1}(x)/K_{\nu }(x)$$\end{document}, are bounds for xαΦi,ν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x^{\alpha } \Phi _{i,\nu }(x)$$\end{document}, which are solutions with unique monotonicity properties; these bounds hold at least for ±α∉(0,1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm \alpha \notin (0,1)$$\end{document} and ν≥1/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu \ge 1/2$$\end{document}. Properties for the product Pν(x)=Iν(x)Kν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{\nu }(x)=I_{\nu }(x)K_{\nu }(x)$$\end{document} can be obtained as a consequence; for instance, it is shown that Pν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{\nu }(x)$$\end{document} is decreasing if ν≥-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu \ge -1$$\end{document} (extending the known range of this result) and that xPν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$xP_{\nu }(x)$$\end{document} is increasing for ν≥1/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu \ge 1/2$$\end{document}. We also show that the double ratios Wi,ν(x)=Φi,ν+1(x)/Φi,ν(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_{i,\nu }(x)=\Phi _{i,\nu +1}(x)/\Phi _{i,\nu }(x)$$\end{document} are monotonic and that these monotonicity properties are exclusive of the first and second kind modified Bessel functions. Sharp trigonometric bounds can be extracted from the monotonicity of the double ratios. The trigonometric bounds for the ratios and the product are very accurate as x→0+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x\rightarrow 0^+$$\end{document}, x→+∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x\rightarrow +\infty $$\end{document} and ν→+∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu \rightarrow +\infty $$\end{document} in the sense that the first two terms in the power series expansions in these limits are exact.


Introduction
Modified Bessel function, and in particular their ratios, are important special functions appearing in countless applications. Bounds for these ratios are needed in a huge number of different scientific and engineering fields, like finite elasticity [15], telecommunications [1], statistics [8], heat transfer [5], information theory [6] and many others. Not surprisingly, this is an active topic of study; see for instance [2,3,7,9,12,13,16,17].
In this paper we obtain new monotonicity properties and bounds for ratios and products of modified Bessel functions, some of them displaying a remarkable accuracy in all three directions as x → 0 + , x → +∞ and ν → +∞; we also extend previous results, in particular for the product of first and second kind modified Bessel functions.
We analyze the monotonicity of the functions x α Φ i,ν (x), with Φ 1,ν (x) = I ν−1 (x)/I ν (x) and Φ 2,ν (x) = −K ν−1 (x)/K ν (x), by considering the Riccati equation satisfied by these functions. It is shown that the nullclines of the Riccati equation are bounds for x α Φ i,ν (x), at least when ±α / ∈ (0, 1) and ν ≥ 1/2. We show that these monotonicity properties are unique for the first and second kind Bessel functions and no other solution of the Riccati equation is both regular and monotonic when ν ≥ 1. The bounds for the ratios of Bessel functions that can be obtained as a consequence of this analysis are described and then applied to the study of the monotonicity and bounds for the product I ν (x)K ν (x). We prove that I ν (x)K ν (x) is decreasing if ν ≥ −1 (enlarging the range of validity considered so far) while xI ν (x)K ν (x) is increasing for ν ≥ 1/2. Upper and lower bounds for the product are also made available.
In a similar way, the monotonicity properties of the double ratios W i,ν (x) = Φ i,ν+1 (x)/Φ i,ν (x) are established and proved to be unique for the first and second kind modified Bessel functions. New sharp trigonometric bounds for both the first and second kind modified Bessel functions ratios are obtained from this analysis. These bounds, both for the ratios and the products, are shown to be very accurate in the three limits x → 0 + , x → +∞ and ν → +∞, in the sense that at least the first two terms of the power series expansions of the ratios and products, in any of these limits, is given exactly by our new bounds.
The main tool for proving these results is the analysis of the qualitative properties of the first order differential equations satisfied by the ratios and double ratios of Bessel functions. For the case of the single ratios, this analysis is similar to that of [12,13]; we summarize some of these results in Sect. 2, we discuss how the monotonicity properties are unique for the first and second kind functions (and therefore the bounds are sharp only for such functions) and we prove the monotonicity properties for the product P ν (x) and the corresponding bounds. In Sect. 3 we study the monotonicity of the double ratio by considering the first order differential equation satisfied by this ratio. In this analysis, the nullclines of the differential equation satisfied by the double ratio, which are solutions of an algebraic cubic equation, will be shown to provide very sharp bounds for the simple ratios Φ i,ν (x) (similarly as happened in [14] for Parabolic Cylinder functions) and then, as a consequence, for the double ratio W i,ν (x) and the product P ν (x).

Bounds from the Riccati Equation
The starting point in the analysis is the difference-differential system [10, 10.29.2] which is satisfied by I ν (x) and e iπν K ν (x) 1 , which as a consequence also satisfy is also a solution of (1), but it is not independent of I ν (x) for integer ν.
For proving the results in this paper, the only information which will be needed as input is the difference-differential system (1) together with information on the sign of the function ratios and first derivatives as x → 0 + and x → +∞; this information will single out two of the solutions of the system (1), specifically the regular solution at x = 0 (I ν (x) = I ν (x)) and the recessive solution as x → +∞ (I ν (x) = e iπν K ν (x)).
We first briefly review how bounds for the ratios of first and second kind modified Bessel functions can be obtained by analyzing the nullclines of the Riccati equations satisfied by these ratios (Theorem 1), as done in [12,13]. The analysis of the bounds for the first and second kind Bessel functions will be done simultaneously using a same Riccati equation, which differs slightly from the approach in [12,13]. Later in this section we consider the more general case of the general solution of the system (1) and we study the monotonicity properties and bounds for the product P ν (x) = I ν (x)K ν (x).
Starting from the DDE (1) we can obtain the Riccati equation for giving Using {I ν (x), (−1) ν K ν (x)} as a pair of independent solutions of the DDE (1), we can write the solutions Φ ν (x) as , t ∈ (−1, 1], (5) 1 The complex notation is not substantial and we could have also defined a second solution for real ν as (−1) ν Kν (x) Φ t,ν (x) is the general solution of (4). That it is a solution is obvious by construction and that any solution can be written in this form is also clear: for any ν and for each (x, y) there is one and only one solution of (4) such that Φ ν (x) = y and there exists a unique value of t ∈ (−1, 1] such that Φ t,ν (x) = y, and Φ t,ν (x) is precisely this unique solution.
As in [12], we consider a more general Riccati equation by taking We have As we will see, for establishing the bounds on the function ratios it is important that the functionsγ ± a,ν (x) determining the nullclines γ t,a,ν (x) = γ ± a,ν (x) are monotonic. It is easy to prove that: For a ∈ (−1, 0) ∪ (0, 1), let ) has a relative extremum at x e (respectively −x e ), and it is a minimum (respectively maximum) if a < 0 and a maximum (respectively minimum) if a > 0.
We consider next the particular and more important cases t = 0 and t = 1 (that is, γ 0,a,ν (x) = x a I ν−1 (x)/I ν (x) and γ 1,a,ν (x) = −x a K ν−1 (x)/K ν (x)) and prove the monotonicity of both functions; later we consider the general case t ∈ (−1, 1]. Theorem 1. The following monotonicity properties and bounds hold: 1. Properties of γ 0,a,ν (x): (a) If |a| > 1 then for all ν ∈ R, aγ 1,a,ν (x) is increasing and Proof. We don't give a detailed proof for all the cases, but all the results follow from similar arguments. For t = 0 (modified Bessel function of the first kind), the monotonicity property as x → 0 + of γ 0,a,ν (x) (unique for this solutions) is the main input. (7)) and the fact thatγ This would prove the result 1.a for a ≤ −1. Similarly for the rest of 1.
for all a ≥ 0 and ν > 1/2 (which implies that γ 0,a,ν (x) is monotonically decreasing for a ∈ (0, 1) too). In fact, the range of validity as a function of ν increases as a increases from , a ≤ 0, with the range of validity increasing as a decreases. This implies that γ 1,a,ν (x) is also monotonically decreasing for a ∈ (−1, 0). Results Math Now we turn to the general case t ∈ (−1, 1]. We notice that, because as x → +∞ the function I ν (x) is exponentially increasing while K ν (x) is exponentially decreasing we have that, for all real ν and t = 1, On the other hand, if ν ≥ 0 then I ν (x) is regular at the origin, while K ν (0 + ) = +∞ , and therefore, for t = 0 and ν ≥ 1, In other words, the behaviour of the solution which is unique. From this information the next result follows, which will be used to prove that the monotonicity properties for first and second kind Bessel functions are unique (Theorem 2).
Then γ t,a,ν (x) for t ∈ (0, 1) correspond to regular solutions which are inside D, while for t ∈ (−1, 0) they have a vertical asymptote at x * > 0 and their graph is outside D.
Proof. In the first place we notice that the existence and unicity conditions for the solutions of the Riccati equation are fulfilled and that, therefore, given a point (x, y), x > 0, there is only one value of t such that γ t,a,ν (x) = y. Therefore, the integral lines can not cross. Now, taking into account (10) we know the graph of γ t,a,ν (x) approaches the graph of γ 0,a,ν (x) as x → +∞ and, on the other hand, it is easy to check that γ t,a,ν (0 + ) < γ 0,a,ν (0 + ), t = 0 (see (11)) for ν ≥ 1, and use the series given in the Appendix for 0 ≤ ν < 1).
For the case t ∈ (−1, 0), and because I ν (x)/K ν (x) increases monotonically from 0 to +∞ in (0, +∞), there exists a single x ∞ > 0 such that I ν (x ∞ )/K ν (x ∞ ) = − tan(πt/2). Therefore the denominator of (5) is zero at x ∞ , where the function has a vertical asymptote. Because the solution tends and therefore the graph of the solution is outside D.
On the other hand, if t ∈ (0, 1) the denominator of (5) is always positive and γ t,a,ν (x) is continuous and its graph lies below the graph of γ 0,a,ν (x), and therefore is inside the region D.
Theorem 2. If ν ≥ 1, there are no other regular and strictly monotonic solutions of (7) other than γ 0,a,ν (x) and γ 1,a,ν (x). γ 0,a,ν (x) is strictly monotonic except when a ∈ (−1, 0). γ 1,a,ν (x) is strictly monotonic except when a ∈ (0, 1). Proof. As discussed before, the solutions with t ∈ (−1, 0) have a discontinuity and therefore, the only thing left to prove is that the solutions with t ∈ (0, 1) are not monotonic and the particular cases for t = 0 and t = 1 give monotonic solutions if a / ∈ (0, 1) and a / ∈ (−1, 0) respectively. We first observe that Theorem 1 implies that one or two of the nullclines are inside the region D, namely, the graph ofγ (10) and (11) the graph of γ t,a,ν (x), for any t ∈ (0, 1) tends to the upper boundary of D as x → +∞ and to the lower boundary as x → 0 + . Therefore, it crosses the nullcline(s) inside D. More specifically, there is a local maximum if a ≥ 0 becauseγ + a,ν (x) is inside D and a minimum if For a ∈ (0, 1) the solution γ 0,a,ν (x) keeps being monotonic (see Remark 3), but not γ 1,a,ν (x) because the derivative changes sign as can be checked by considering the expansions as x → 0 + and x → +∞ of the Appendix; γ 0,a,ν (x) has a maximum in this case. The rest of solutions, can not be regular and monotonic, by the same arguments as before. The same can be said for a ∈ (−1, 0), changing the roles of γ 0,a,ν (x) and γ 1,a,ν (x); γ 0,a,ν (x) has a minimum in this case.
From the bounds from the Riccati equations and the use of the recurrence relation, most of the know Amos-type inequalities of the form (α + β 2 + x 2 )/x can be established (see [12,13]), with the exception of the Simpson-Spector bound [15], which follows from arguments similar but not identical to the ones considered here for the Riccati equations. We will not be exhaustive in the description of these bounds, and we refer to [7] for a systematic analysis of Amos-type bounds. Here we concentrate on the bounds that can be extracted from the qualitative analysis of first order differential equations. A way to extend the analysis was considered in [12] by iteration of the Riccati equations, and we explore an in Sect. 3 alternative possibility by considering the differential equation satisfied by double ratios, similar to that described in [14] for Parabolic Cylinder functions.
We end this section with an analysis of the monotonicity properties and bounds for the ratio the monotonicity properties discussed so far.

Properties for the Product I ν (x)K ν (x)
We notice that, using the Wronskian relation [10, 10.28 and the recurrence relation (2) we have K ν−1 (z)I ν (z) + K ν (z)I ν−1 (z) = 1/z, and then we obtain the following relation with the product I ν (x)K ν (x): 2 Then, using our previous notation Now we notice that Theorem 1 and Remark 2 estate that both γ 0,−1,ν (x) and −γ 1,−1,ν (x) are increasing functions if ν ≥ −1, and that this proves that I ν (x)K ν (x) is decreasing for ν ≥ −1. This enlarges the range of validity of the result proved in [11], which was later extended to ν ≥ −1/2 in [2]. Here we have just proved this result in a very straightforward way and in the larger range ν ≥ −1. We also prove next that xI ν (x)K ν (x) is increasing for ν ≥ 1/2. We collect both results in a single theorem: Proof. We only need to prove this result for λ = 0, 1; for the rest of values it follows immediately.
Using (13) and the bounds for the ratios of Bessel functions, sharp bounds for the product can be established. We will not be exhaustive in this discussion, as the bounds can be straightforwardly derived. We just give two of these bounds, which are obtained from Theorem 1 (more bounds are available from this same theorem).

Theorem 4. The following two bounds hold:
Proof. For the upper bound use (13) and that Φ 0,ν (x) > λ + 0,ν (x) and Φ 1,ν (x) < −1/λ + 0,ν (x), ν ≥ 1/2. For the lower bound use Φ 0,ν (x) < λ + −1,ν (x), ν ≥ −1, and Φ 1,ν (x) > −1/λ + 1,ν (x), ν ∈ R. These two bounds (as all the bounds that can be extracted from Theorem 1) are sharp as x → +∞. They are not sharp, however, as x → 0 + . This is in contrast with the bounds in [4,Thm. 2], which are sharp as x → 0 + and ν > 0 but not as x → +∞. This is as expected, because we are using bounds for the ratios which are sharp as x → +∞ but not as x → 0 + ; however, upper and lower bounds for the ratios which are sharp in both limits are available (see for instance [12]), and from there it is straightforward to obtain sharp bounds for the product. In particular, the bounds from the iteration of the Riccati equation given in [12] are sharp in both limits. We don't give here such bounds for the product explicitly, which are straightforward applications of previous results, but we will obtain later a new very sharp trigonometric bound which is very accurate in the three limits x → 0 + , x → +∞ and ν → +∞.

Very Sharp Trigonometric Bounds
As in [14], we will study the monotonicity properties of the double ratios and we will establish very sharp trigonometric bounds from these monotonicity properties.
That the double ratios are monotonic both for the first and second kind modified Bessel functions, has been separately shown in two different papers by different methods [16,17]. First, in [16] it was proved that K u (x)K v (x)/ K (u+v)/2 (x) 2 is strictly decreasing for x > 0 and real u, v; integral representation for the product and ratios of modified Bessel functions of the second kind were considered in this analysis. Later, in [17] it was shown that the ratio increasing for x > 0, using the Frobenius series for the Bessel functions. Here, we give a more restricted version of these properties (|u − v| = 2), but we do this in a single analysis for the first and second kind functions, we prove that such monotonicity properties are unique for these two solutions and we obtain bounds for the ratios and products that are sharper and of a different type to those obtained with previous analysis.
Using (2) we have And in terms of which satisfies we have that and differentiating (17) Next we will analyze the qualitative properties of the solutions of the system of equations (18)- (17), and from these we will obtain very sharp trigonometric bounds. Notice that (18) has been obtained by differentiating (17) and using (16) and that, conversely, differentiating (17) and using (18) we obtain that the possible differentiable ψ ν -solutions of the system (18)- (17) are the trivial solution ψ ν (x) = 0 and the solutions of (16), with general solution given by (5). For obvious reasons (the objective is to find properties for modified Bessel functions) we are only considering the latter solutions, in which case the solutions W ν (x) are where α = πt/2, t ∈ (−1, 1]. Remark 5. We notice that, because of (17), W t,ν (x) > − ν 2 x 2 for any real t.
From now on, we will drop the notation W t,ν (x) in favor or W ν (x). We will recover it in Theorem 7.
For analyzing the qualitative properties of the solutions of the system (18)-(17) we need to analyze the nullclines of (18), which determine the monotonicity properties of the solutions.

Properties of the Nullclines
For proving the monotonicity properties and bounds for the double ratio W ν (x), we need to analyze the nullclines in terms of the values of ψ ν (x) which make the right-hand of Eq. (18) zero and then to study the corresponding values of W ν (x) and their monotonicity. We first analyze in Lemma 3 the nullclines in terms of the values of ψ ν (x); after this, the properties for the corresponding values of W ν (x) are analyzed in Lemmas 4 and 5. Once these lemmas are proved, the main results can be estated.

Lemma 3. The cubic equation
has three distinct real roots These solutions can be written where g ν (x) = 3(ν 2 + x 2 ) + 1. Proof ν (x)) and finally a root in (|ν|, +∞) (λ For solving the equation, we transform the cubic f (λ) = 0 to depressed form with the change λ = μ− 1 3 and we have μ 3 −pμ−q = 0 with p = ν 2 +x 2 + 1 3 and q = 2 3 ν 2 − 1 3 x 2 − 2 27 and using the well know trigonometric formula for the solution of a depressed cubic: from where we have that the solutions have the form (21). The solutions are differentiable when the absolute value of the argument of the arccos is smaller than 1, which is equivalent to saying that the discriminant Δ = −(4p 3 + 27q 2 ) is positive, and we have which is positive for x = 0 (see remark 8). Now we expand as x → +∞ for the three values α = 0, ±2π/3 and we get which shows that the ordering λ ν (x) is correct. That the solutions are even functions is immediate given the symmetries of the equation and the monotonicity follows by taking the derivative of (20), from where Consider now x > 0. The numerator does not change sign and neither does the denominator, because if the denominator was zero for some x > 0 then λ ν (x) would not be differentiable for this x, which can not be true. Then, the sign of the denominator for x > 0 is equal to its sign as x → +∞ which, using (22), is positive for λ ν (x) = λ ν (x) > 0 and the other two solutions are negative the monotonicity properties follow. Remark 6. The notation λ (I) ν (x) is used because this solution will be related to a bound for I ν−1 (x)/I ν (x). Similarly, λ  Remark 8. As explained in the previous proof, the solutions λ ν (x) are simple if x = 0. We will be interested in the case x > 0. For x = 0, the solutions are trivially λ ν (x) = −1, ±ν. The identification with the previous notation is given by λ ν (0) = max{−1, −|ν|} (which implies quite curious trigonometric identities) Then, for all x > 0, w Proof. First we observe that the fact that w ν (x). Then, for this value of x we have that both λ ν (x) = λ (I) ν (x) and −λ ν (x) are solutions of (20). Then: and the first two terms suffice to see that w for large x and therefore that this holds for all x > 0.

Remark 9.
That only two monotonic and regular solutions exists is no longer true for smaller values of ν. For instance, it is easy to check that for ν ∈ (0, 1) the solutions satisfying tan α ≤ 2 π sin(πν) are monotonic.
We end this section by obtaining some inequalities for the product of Bessel functions which are a direct consequence of the previous trigonometric bounds, and we propose a conjecture.