# On the Regular Integral Solutions of a Generalized Bessel Differential Equation.

1. IntroductionThe Bessel differential equation was first considered in connexion with the oscillations of a heavy chain [1] and vibrations of a circular membrane [2] and has had since [3] a vast number of applications supported by an extensive theory [4]. A substantial number of applications arise from the separation of variables in the Laplace operator in cylindrical and spherical coordinates that leads, respectively, to the cylindrical and spherical Bessel functions. The cases of specific interest include two types of cylindrical waves: (i) sound waves as compressible perturbations of a uniform flow; (ii) vortical waves as incompressible perturbations of a uniform flow with superimposed rigid body rotation. Whereas (i) and (ii) separately lead to the original Bessel differential equation, their coupling leads to a generalization. Thus, the consideration of coupled acoustic-vortical waves as rotational compressible perturbations of a uniform mean flow with rigid body rotation leads to the generalized Bessel equation that differs from the original in having an extra term involving a second parameter, namely, the degree [mu], in addition to the order v.

The generalized Bessel differential equation of order v and degree [mu] may have other applications and deserves separate study as it leads to generalizations of the Bessel and Neumann functions. The generalized Bessel differential equation may also be obtained, aside from any physical or engineering motivations, by a purely mathematical argument, starting from the original Bessel differential equation and replacing the coefficients of the dependent variable and its derivative by polynomials of the independent variable; in this case the origin remains a regular singularity of the differential equation and the only other singularity is the point-at-infinity. Thus, solutions exist as Frobenius-Fuchs series [5, 6] with recurrence formula for the coefficients reducing to two terms only in the case of the generalized Bessel differential equations.

The generalized Bessel differential equation has singularities only at the origin and infinity. Since the singularity at the origin is regular, the Frobenius-Fuchs method leads to solutions valid for finite values of the variable. The solutions of the generalized Bessel differential equation around the regular singularity at the origin has (i) indices that are exponents of the leading power depending only on the order; (ii) recurrence relation for the coefficients of the power series expansion depending also on the degree. From (i) follows the familiar situation that generalized Bessel functions (Section 2) specify the general integral for noninteger order, and generalized Neumann functions (Section 3) are needed for integer order. From (ii) it follows that the series expansion for the generalized Bessel (Section 2.1) and Neumann (Section 3.1) functions differ from the original series in having finite products multiplying each term; these finite products can be expressed as ratios of Gamma functions, whose arguments become singular for zero degree. The Wronskians are used to select pairs of linearly independent particular integrals that lead to the general integral for noninteger (Section 2.2) and integer (Section 3.2) order.

2. Generalized Bessel Differential Equation in the Complex Plane

The origin is a regular singularity of the generalized Bessel equation with the same indices (Section 2.1) as the original, leading to generalized Bessel functions whose series expansion differs from the original in the coefficients following the leading term. A linear combination of (Section 2.1) generalized Bessel functions of order [+ or -] v and degree [mu] supplies the general integral (Section 2.2) of the generalized Bessel equation if the order v is not an integer. For v = n (i.e. an integer), the generalized Bessel functions of order [+ or -] n are linearly dependent for any degree [mu], their Wronskian is zero, and generalized Neumann functions are needed (Section 3).

2.1. Generalized Bessel Functions of Arbitrary Order and Degree

Definition 1. The generalized Bessel differential equation is a linear second-order ordinary differential equation (1) in the complex plane, whose order v and degree [mu] are complex numbers,

[mu],v,z [member of] C : [z.sup.2]Q" + z(1 - [mu]/2 [z.sup.2])Q' + ([z.sup.2] - [v.sup.2]) Q = 0. (1)

Remark 2. The original Bessel differential equation corresponds to zero degree [mu] = 0. The choice of [mu]/2 rather than [mu] in the coefficient of Q' simplifies subsequent expressions for the generalized Bessel and Neumann functions that are solutions of (1), much as the choice of [v.sup.2] rather than v in the coefficient of Q simplifies the expressions for the original Bessel and Neumann functions.

Theorem 3. A solution of the generalized Bessel differential equation is the generalized Bessel function of order v and degree [mu] specified by the power series with infinite radius of convergence

[mathematical expression not reproducible] (2)

valid for all [mu].

Remark 4. In the case of the original Bessel function the degree is zero, [mu] = 0, and the last factor in (2) is omitted, leading to the usual series expansion [7] that is alternating for real variable z and has fixed sign for imaginary variable. The last factor that distinguishes the generalized Bessel function for real positive v > 0 is (i) positive for [mu] < 0 and does not change the sign of the terms of the series; (ii) for real [mu] > 0 the sign of the terms of the series is fixed for [mu](l + v/2) > 1,andin particular if [mu]v/2 > 1 terms of the series have the same sign. For example, if [mu] and v are real such that [mu]v/2 > 1 then the series for the generalized Bessel function (2) has fixed sign for real z and alternating sign for imaginary z; this is the reverse of the series for the original Bessel function when the last factor in (2) is omitted.

Corollary 5. If the degree is not zero (3a) the last factor in the coefficients (2) may be written in (3b)

[mu] [not equal to] 0: (3a)

[mathematical expression not reproducible]; (3b)

substitution of (3b) in (2) leads to

[mu] [not equal to] 0: [J.sup.[mu].sub.v](z) (4a)

[mathematical expression not reproducible]. (4b)

This is an alternative expression (4b) for the generalized Bessel function (2) when the degree is not zero (4a) and involves powers of ([mu],[z.sup.2])/4 rather than (-[z.sup.2])/4, hence omitting the alternating sign and inserting the factor [mu].

Proof of Theorem 3. The origin z = 0 is a regular singularity, and the only other singularity is at infinity; thus [8] the solution as a Frobenius series [5, 9] by the Fuchs theorem [6,10] converges in the whole finite complex z-plane

[absolute value of (z)] < [infinity] : [Q.sub.[sigma]](z) = [[infinity].summation over (j=0)][a.sub.j]([sigma])[z.sup.j+[sigma]]. (5)

Substituting (5) into (1) and equating the coefficients of equal powers of z lead to a two-term recurrence formula

[[([sigma] + j).sup.2] - [v.sup.2]][a.sub.j] (a)

= -[1 - [mu]/2([sigma] + j -2)][a.sub.j-2]([sigma]). (6)

Setting (7a) leads to (7b); since [a.sub.0] = 0 would lead to [a.sub.j] = 0 in (6) and Q = 0 in (5), a nontrivial solution requires (7c) leading to the indicial equation (7d) with roots (7e):

j = 0: (7a)

([[sigma].sup.2] - [v.sup.2])[a.sub.0] ([sigma]) = 0; (7b)

[a.sub.0] ([sigma]) [not equal to] 0; (7c)

[[sigma].sup.2] = [v.sup.2]; (7d)

[[sigma].sub.[+ or -]] = [+ or -] v; (7e)

thus, the indices (7e) are specified by the order v as in the original Bessel equation, but the recurrence formula for the coefficients (6),

[mathematical expression not reproducible], (8)

shows that they depend also on the degree [mu]. Applying j times (8) leads to the explicit coefficients

[mathematical expression not reproducible] (9)

where r is the Gamma function [11,12]. For the original Bessel functions, [mu] = 0, the last function is unity leading to the usual expression [4], using for consistency (10a):

[a.sub.0](v) = [2.sup.-v]/[GAMMA](1 + v). (10a)

Thus, the generalized Bessel function of the first kind, order v, and degree [mu] is given by (5) and (10b),

[mathematical expression not reproducible], (10b)

leading to (2).

2.2. Wronskian, Linear Independence of Solutions, and General Integral. The generalized Bessel differential equation (1) is satisfied by the generalized Bessel functions (2) of orders [+ or -] v as follows from the indices (7e). If the functions [J.sub.[+ or -]v] are linearly independent, their linear combination specifies the general integral of the generalized Bessel equation (1). Next the Wronskian is calculated for two solutions of the generalized Bessel differential equation (1), in particular for the pair of solutions [J.sub.[+ or -]v]. It follows that the Wronskian is zero for v = n (i.e., an integer), implying that [J.sup.[mu].sub.[+ or -]n](z) are linearly dependent. The linear relation between [J.sup.[mu].sub.-n](z) and [J.sup.[mu].sub.n](z) is obtained for all integer values of n.

Lemma 6. If [Q.sub.1], [Q.sub.2] are any two solutions of the generalized Bessel differential equation,

[z.sup.2][Q".sub.1] + z(1 - [mu]/2[z.sup.2])[Q'.sub.1] + ([z.sup.2] - [v.sup.2])[Q.sub.1] = 0, (11a)

[z.sup.2][Q".sub.2] + z(1 - [mu]/2[z.sup.2])[Q'.sub.2] + ([z.sup.2] - [v.sup.2])[Q.sub.2] = 0, (11b)

their Wronskian is given by

W(z) = [Q.sub.1](z)[Q'.sub.2](z)-[Q'.sub.1](z)[Q.sub.2](z)

= [W.sub.0]([mu],v)/z exp(1/4 [mu][z.sup.2]) (12a)

where [W.sub.0]([mu], v): (i) is the residue of the Wronskian at its simple pole at the origin,

[mathematical expression not reproducible]; (12b)

(ii) does not depend on the variable z but may depend on the order v and degree [mu]; (iii) depends on the choice of particular integrals {[Q.sub.1](z), [Q.sub.2](z)} in (11a) and (11b).

Proof. Multiplying the second equation (11b) by [Q.sub.1] and multiplying the first equation (11a) by [Q.sub.2] and subtracting lead to

([Q.sub.1][Q'.sub.2] - [Q.sub.2][Q'.sub.1]) z (1 - [mu]/2 [z.sup.2]) = -[z.sup.2]([Q.sub.1][Q.sub".2] - [Q.sub.2][Q".sub.1])

= -[z.sup.2]([Q.sub.1][Q'.sub.2] - [Q.sub.2][Q'.sub.1])'. (13)

Thus, the Wronskian (14a) satisfies a first-order differential equation (14b),

W([Q.sub.1].[Q.sub.2]) [equivalent to] ([Q.sub.1][Q'.sub.2] - [Q.sub.2][Q.sub'.1]) [not equal to] 0: (14a)

W'/W = [mu]z/2 - 1/z. (14b)

The solution of (14b) is (12a) where (12b) is an arbitrary constant of integration [W.sub.0]([mu], v) that is independent of z.

Remark 7. The constant factor [W.sub.0]([mu], v) in the Wronskian (12a) depends on the particular choice of linearly independent solutions.

Lemma 8. The generalized Bessel function of order [+ or -] v and degree [mu] has Wronskian

W ([J.sup.[mu].sub.+v](z), [J.sup.[mu].sub.-v](z)) = -2/[pi]z sin (v[pi]) exp(1/4[mu][z.sup.2]). (15)

Proof. In the case of the generalized Bessel functions (2),

[J.sup.[mu].sub.[+ or -]v](z) ~ [(z/2).sup.[+ or -]v/[GAMMA](1 [+ or

-] v) [1 + O([z.sup.2]), (16a)

d/dz [[J.sup.[mu].sub.[+ or -] v](z)] ~ [+ or -] v/2[(z/2).sup.[+ or 1] v-1]/2 [GAMMA](1 [+ or -] v) [1 + O[([z.sup.2])], (16b)

the limit as z [right arrow] 0 specifies the Wronskian

[mathematical expression not reproducible] (17a)

using the symmetry formula [13] for the Gamma function in (17a). Comparison of (17a) with (12a) in the limit z [right arrow] 0, or with (12b), shows that for the choice of particular integrals (16a) of the generalized Bessel differential equation the constant is given by

[W.sub.0]([mu],v) = -2/[pi] sin (v[pi]) (17b)

and thus depends only on the order v but not on the degree [mu]. Substitution of (17b) in (12a) proves (15).

Theorem 9. The general integral of the generalized Bessel equation (1) is a linear combination (18b) with arbitrary constants [C.sub.[+ or -]] of the functions of the first kind (2) for the indices (7e) provided that the indices (18a) are not integers,

v [not equal to] 0, [+ or -] 1, [+ or -] 2, ... = [+ or -] n: (18a)

Q(Z) = [C.sub.+][J.sup.[mu].sub.v](z) + [C.sub.-][J.sup.[mu].sub.-v](z). (18b)

Proof. The general integral (18b) holds provided that [J.sup.[mu].sub.[+ or -] v](z) are linearly independent. This is the case if their Wronskian (15) is not zero. The Wronskian only vanishes for v = n, an integer, so the general integral (18b) holds (18a) only for noninteger values of v.

The Frobenius-Fuchs method [14] suggests that the two solutions may be linearly dependent and the general integral (18b) fails if the difference of indices (7e) is an integer, that is, if [[sigma].sub.+] - [[sigma].sub.-] = 2v corresponds to [15] order v either (i) an integer (cylindrical Bessel functions) or (ii) an integer plus-one-half (spherical Bessel functions). In the case (ii) of order an integer plus-one-half, the two solutions in (18b) are linearly independent because the Wronskian (15) is not zero so the general integral (18b) holds. It remains to consider the case (i) of order an integer in (1) when the functions [J.sup.[mu].sub.[+ or -] n] must be linearly dependent. This is confirmed by the following relation.

Theorem 10. The generalized Bessel functions of integer order and nonzero degree (19a) are related by (19b),

[mu] [not equal to] 0: (19a)

[J.sup.[mu].sub.-n](z) = [[mu].sup.n] [GAMMA](n/2 - 1/[mu])/[GAMMA](-n/2 - 1/[mu])[J.sup.[mu].sub.n](z). (19b)

Corollary 11. Using the recurrence formula for the Gamma function,

[mathematical expression not reproducible], (20)

the relation (19b) between generalized Bessel functions of order [+ or -] n becomes

[mathematical expression not reproducible]. (21)

Remark 12. This result also holds for the original Bessel function [mu] = 0 when the last factor in (21) is unity.

Proofof Theorem 9. If the order is a positive integer, v = n, in the second solution [J.sup.[mu].sub.-n] the Gamma function [GAMMA](1 -n+j) = [infinity] for j = 0,1 , ..., n - 1 in the denominator of (4b) suppresses the first n terms, leading to (22b),

[mathematical expression not reproducible] (22a)

[mathematical expression not reproducible] (22b)

where the substitution (22a) was made.

3. Solutions of the Generalized Bessel Equation for Any Order and Degree

It has been shown that, for any degree [mu] and integer order, v = n, the general integral (18b) fails (18a), and the generalized Bessel function [J.sup.[mu].sub.-n] that is a regular integral of the first kind of the generalized Bessel differential equation (1) must be replaced by a generalized Neumann function [Y.sup.[mu].sub.n] that is a regular integral of the second kind (Section 3.1) and hence is linearly independent, leading by linear combination of [J.sup.[mu].sub.v](z) and [Y.sup.[mu].sub.v](z) to the general integral (Section 3.2).

3.1. Generalized Neumann Function of Arbitrary Degree and Integer Order

Theorem 13. A solution of the generalized Bessel equation (1) with arbitrary degree [mu] and integer order n is the generalized Neumann function

[Y.sup.[mu].sub.n](z) = 2/[pi] log(z/2)[J.sup.[mu].sub.n](z) + [X.sup.[mu].sub.n](z) + [Z.sup.[mu].sub.n](z) (23)

Consisting of the sum of (i) a constant 2/[pi] multiplying the generalized Bessel function (2) multiplied by a logarithmic singularity; (ii) the preliminary function

[mathematical expression not reproducible] (24)

that has a pole of order n for n = 2,3, ...; (iii) the complementary function

[mathematical expression not reproducible] (25)

involving the digamma function y that has no singularities and starts with the power n.

Remark 14. Substitution of (2), (24), and (25) into (23) specifies explicitly the Neumann function of complex degree [mu] and integer order n:

[mathematical expression not reproducible] (26)

For [mu] = 0 the terms in curly brackets reduce to unity, and the original Neumann [16] function is regained.

Corollary 15. In the case of nonzero degree, [mu] [not equal to] 0, there are alternate expressions for the three terms whose sum is (23) the generalized Neumann function: (i) the logarithmic factor multiplies the generalized Bessel function (2) that has the alternate form (4b) for nonzero degree (4a); (ii) the complementary function (25), using (3b) has the alternate form (27b) for nonzero degree (27a),

[mathematical expression not reproducible] (27a)

[mathematical expression not reproducible]; (27b)

(iii) the preliminary function (24) that for nonzero degree (28a) has the alternate form (28b),

[mu] [not equal to] 0, (28a)

[mathematical expression not reproducible]. (28b)

Proof of Theorem 13. The Frobenius-Fuchs method [17] indicates how a linearly independent function of the second kind can be obtained in the case of coincident indices or indices differing by an integer. For v = -n the terms of the series (2) starting with j = n become infinite due to the zero in the denominator. This can be avoided going back to solution (5) for arbitrary index [sigma] in (8)

[mathematical expression not reproducible]. (29)

Multiplying by [sigma] + n cancels, for v = -n, the zero in the denominator 2l + [sigma] + v = 2l + [sigma] - n when l = n, and the limit [sigma] [right arrow] -n can be taken with a finite result. This leads to a constant multiple of the generalized Bessel function. The substitution of (29) with a factor [sigma]+n in the generalized Bessel equation (1) leads to

[mathematical expression not reproducible]; (30)

since the recurrence formula (8) ensures that all terms vanish except the first, as [sigma] [right arrow] [+ or -] n the r.h.s. (right-hand side) of (30) vanishes, leading in both cases to the generalized Bessel function or a constant multiple. The r.h.s. of (30) also vanishes taking the limit [sigma] [right arrow] -n after differentiation with regard to <r, leading to the solution

[mathematical expression not reproducible]; (31)

this solution involves a logarithmic term arising from [partial derivative]([z.sup.[sigma]])/[partial derivative][sigma] = [z.sup.[sigma]] log z in (31) and thus is linearly independent of the generalized Bessel function (2). It is designated generalized Neumann function and the factor 2/[pi] was inserted for consistency with the usual definition [18] of the original Neumann function for [mu] = 0. When substituting the series (31) into (29) the first n terms are separated to specify a preliminary function

[mathematical expression not reproducible]. (32)

If the derivative [partial derivative]/[partial derivative][sigma] is not applied to ([sigma] + n) this factor remains and leads to zero in the limit [sigma] [right arrow] -n; thus the only nonzero terms arise differentiating [partial derivative]([sigma] + n)/[partial derivative][sigma] = 1, that is, suppressing the factor ([sigma] + n) and taking the limit [sigma] [right arrow] -n in the remaining terms leading to

[mathematical expression not reproducible] (33)

Choosing the leading coefficient,

[mathematical expression not reproducible], (34)

and substituting in (33) specifies the preliminary function (24). The first two factors in (34) were chosen in agreement with the preliminary function for the original Neumann function; the factor in curly brackets is relevant for [mu] = 0 and affects only the generalized Neumann function. The second part of the generalized Neumann function (31) involves the remaining terms of the series (29) stating with j = n for which the factor in curved brackets cancels the same factor in the denominator:

[mathematical expression not reproducible]. (35)

The leading coefficient is chosen (36) to suppress the first n factors in the last product in (35),

[mathematical expression not reproducible], (36)

simplifying (35) to (37c),

k = j - n, (37a)

[beta] = l - n: (37b)

[mathematical expression not reproducible] (37c)

where (37a) and (37b) were used. The differentiation with regard to [sigma] of [(z/2).sup.[sigma]] inserts the factor log(z/2) multiplying the generalized Bessel function of the first kind (2),

[pi]/2 [Y.sup.[mu].sub.n](z) - [X.sup.[mu].sub.n](z)] = log(z/2)[J.sup.[mu].sub.n](z) + [pi]/2[Z.sup.[mu].sub.n](z), (38)

and the remaining term is a complementary function. The complementary function corresponds to differentiation with regard to o after the factor [(z/2).sup.[sigma]] in (37a), (37b), and (37c) and replacing the dummy summation indices (k, [beta]) by (j, l) leads to

[mathematical expression not reproducible] (39)

and then (25) follows where [psi] is the digamma function [8, 13]. Thus, the generalized Neumann function of order n and degree [mu] is given by (38) [equivalent to] (23).

3.2. Generalized Integral for Arbitrary Degree and Integer Order. Using the Wronskian it is proven that the generalized Neumann function is linearly independent from the generalized Bessel function, and thus their linear combination specifies the general integral of the generalized Bessel equation (1) in the case missing from (18a) and (18b), namely, integer order, v = n, and arbitrary degree.

Lemma 16. The Wronskian of the generalized Bessel and Neumann functions of complex degree [mu] and integer order n is

W([J.sup.[mu].sub.n](z), [Y.sup.[mu].sub.n](z)) = 2/[pi]z A([mu],n) exp(1/4[mu][z.sup.2]) (40)

where A([mu], n) is given by

[mathematical expression not reproducible]. (41)

Remark 17. Coefficient (41) appears in relation (19b) between the generalized Bessel functions of any degree [mu] and integer order n:

n [member of[ N : [J.sup.[mu].sub.-n](z) = [(-1).sup.n] A([mu],n)[J.sup.[mu].sub.n](z). (42)

Remark 18. In the case of the original Bessel functions of zero degree coefficient (41) is unity (43a) simplifying (42) to the known [4,11] relation (43b) between Bessel coefficients,

A (0, n) = 1: (43a)

[J.sub.-n](z) = [(-1).sup.n][J.sub.n](z). (43b)

Proof of Lemma 16. The generalized Neumann function (26) has leading term as z [right arrow] 0,

[Y.sup.[mu].sub.n](z) = -A([mu],n)(n - 1)!/[pi] [(z/2).sup.-n][1 + O([z.sup.2])], (44a)

d/dz [[Y.sup.[mu].sub.n](z)] = A([mu],n) n!/2[pi][(z/2).sup.-n-1][1 + O([z.sup.2])], (44b)

involving the constant coefficient (41); this leads to the Wronskian (45a) with the generalized Bessel function (16a) and (16b) and also of integer order:

W([J.sup.[mu].sub.n](z), [Y.sup.[mu].sub.n](z)) = 2/[pi]z A([mu],n)[1 + O([z.sup.2])]; (45a)

[W.sub.0]([mu],n) = 2/[pi] A([mu],n), (45b)

and the Wronskian (45a) agrees with (12a) for (45b), specifying the exact Wronskian (40).

Theorem 19. The general integral of the generalized Bessel equation (1) is (46b) with [C.sub.1] and [C.sub.2] arbitrary constants,

n= 0, [+ or -] 1, [+ or -] 2, ... : (46a)

Q(z) = [C.sub.1][J.sup.[mu].sub.n](z) + [C.sub.2][Y.sup.[mu].sub.n] (z), (46b)

for integer order (46a) and arbitrary degree [mu].

Proof. Both the generalized Bessel (2) and Neumann (26) functions satisfy the generalized Bessel equation (1) for arbitrary degree and integer order (46a); their Wronskian is nonzero so they are linearly independent and their linear combination (46b) specifies the general integral.

4. Conclusions

The generalized Bessel differential equation appears for coupled acoustic-vortical wave problems, which would have satisfied the original Bessel differential equation in the decoupled acoustic or vortical case. The generalized Bessel differential equation (1) with order v and degree [mu] reduces to the original Bessel differential equation for zero degree [mu]. The origin is a regular singularity and the other singularity is at infinity, so the Frobenius-Fuchs method specifies power series solutions valid in the finite complex plane. The indices are specified by the order but the recurrence formula for the coefficients depends also on the degree. Thus, the general integral of the generalized Bessel differential equation is a linear combination of generalized Bessel functions [J.sup.[mu].sub.[+ or -] v] if the order is not an integer; the generalized Bessel functions of integer order [J.sup.[mu].sub.[+ or -] n] are linearly dependent, as shown by an explicit relation between them. It follows that the general integral of the generalized Bessel differential equation for integer order and arbitrary degree requires the introduction of a generalized Neumann function. The generalized Bessel and Neumann functions have expressions that reduce to the original Bessel and Neumann functions for zero degree and also alternate expressions valid only for nonzero degree. The Bessel differential equation can be generalized further (see the Appendix) but in that case the coefficients on the power series solutions would no longer satisfy double recurrence formulas; the latter would be replaced by triple or higher order recurrence formulas.

The present paper is, as far as the authors know, the first concerning the generalized Bessel differential equation (1), and its regular integral solutions of the first and second kind in the finite complex plane that specify, respectively, the generalized Bessel and Neumann functions. There is a vast literature on the original Bessel differential equation and its solutions in terms of the original Bessel and Neumann functions: (i) starting with the first developments [1-3, 16]; (ii) detailed in several monographs [4, 7, 9-11, 13-15, 17, 18]; (iii) continuing with research papers up to the present time on various properties [19-26]. Several generalizations of the original Bessel function have been proposed, usually [26, 27] based on extensions of (2) without the last factor; these generalized Bessel functions are not derived from an explicit differential equation and are not associated with Neumann functions. The present approach is different in that it starts from a generalized Bessel differential equation (1) and leads to generalized Bessel (2) and Neumann (26) functions. The vast literature [1-26, 28] on the original Bessel differential equation and related Bessel and Neumann functions includes many properties that could potentially be extended to the generalized Bessel differential equation and generalized Bessel and Neumann functions and may be the subject of future work.

Appendix

Further Generalizations of the Bessel Differential Equation

The original Bessel differential equation

[z.sup.2]Q" + zQ' + ([z.sup.2] - [v.sup.2]) Q = 0. (A.1)

can be generalized replacing the coefficients of Q' and Q by polynomials of the independent variable,

[mathematical expression not reproducible], (A.2)

where ([[alpha].sub.1], ..., [[alpha].sub.N]) and ([[beta].sub.1], [[beta].sub.3], ..., [[beta].sub.M]) are constants. The origin z = 0 remains a regular singularity of the differential equation (A.2) and the only other singularity is the point-at-infinity. Thus, there exists a solution as a Frobenius-Fuchs series (5), with coefficients satisfying the recurrence formula:

[mathematical expression not reproducible] (A.3)

In order to have a two-term recurrence formula, involving only [a.sub.j]([sigma]) and [a.sub.j-2]([sigma]) that can be solved explicitly, all coefficients must vanish (A.4a),

[[alpha].sub.1] = [[alpha].sub.3] = ... = [[alpha].sub.N] = 0 = [[beta].sub.1] = [[beta].sub.3] = ... = [[beta].sub.M], (A.4a)

[[alpha].sub.2] = -[mu]/2, (A.4b)

except (A.4b). Substitution of (A.4a) and (A.4b) into (A.2) leads to the generalized Bessel differential equation (1).

Note that the more general differential equation (A.2) still has the indices (7a)-(7e). Also, the recurrence formula (A.3) with three terms could be solved using continued fractions. For example, relaxing conditions (A.4a) to

[[alpha].sub.3] = ... = [[alpha].sub.N] = 0 = [[beta].sub.3] = ... = [[beta].sub.M], (A.5)

the differential equation (A.2) becomes

[mathematical expression not reproducible] (A.6)

which reduces to the generalized Bessel equation (1) for [alpha] = 0 = [beta]. If [alpha] [not equal to] 0 or [beta] [not equal to] 0 the solution of (A.6) is still a Frobenius-Fuchs series (5) with recurrence formula for the coefficients,

[mathematical expression not reproducible] (A.7)

that can be solved [28] as the finite descending continued fraction,

[mathematical expression not reproducible], (A.8)

with ([a.sub.0], [a.sub.1]) that are arbitrary constants as the starting values to calculate [a.sub.2],[[alpha].sub.3], .... The method of continued fractions would no longer apply if the recurrence formula for the coefficients would have more than 3 terms.

https://doi.org/10.1155/2018/8919516

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by FCT (Foundation for Science and Technology) through IDMEC (Institute of Mechanical Engineering), under LAETA Pest, OE/EME/LA0022.

References

[1] D. Bernoulli, "Theoremata de oscillationibus corporum filo flexili connexorum et catenae verticaliter suspensae," Novi commentarii Academiae Scientiarum Imperialis Petropolitanae, vol. 6, pp. 108-122.

[2] L. Euler, "De motu vibratorio tympanorum," Novi commentarii Academiae Scientiarum Imperialis Petropolitanae, vol. 10, pp. 243-260, 1764.

[3] F. W. Bessel, "Untersuchung des Theils der planetarischen Storungen, welcher aus der Bewegung der Sonne entstcht," Abhandlungen der Koniglichen Akademie der Wissenschaften in Berlin, pp. 1-52, 1824.

[4] G. N. Watson, A Treatise on The Theory of Bessel Functions, Cambridge University Press, Cambridge, UK, 2nd edition, 1944.

[5] G. Frobenius, "Ueber die Integration der linearen Differentialgleichungen durch Reihen," Journal fur die reine und angewandte Mathematik, vol. 76, pp. 214-235, 1873.

[6] L. Fuchs, "Zur Theorie der linearen Differentialgleichungen mit veranderlichen Coefficienten (Ergonzungen zu der im [66.sup.sten] Bande dieses Journals enthaltenen Abhandlung)," Journal fur die reine und angewandte Mathematik, vol. 68, pp. 354-385, 1868.

[7] N. W. McLachlan, Bessel Functions for Engineers, Oxford University Press, Oxford, UK, 1934.

[8] L. M. Campos, Complex Analysis with Applications to Flows and Fields, CRC Press, Boca Raton, FL, USA, 2011.

[9] A. R. Forsyth, A Treatise on Differential Equations, Macmillan, London, UK, 6th edition, 1956.

[10] E. L. Ince, Ordinary Differential Equations, Dover, Mineola, NY, USA, Reprint edition, 2014.

[11] E. T. Whittaker and G. N. Watson, A course of modern analysis: An Introduction to the General Theory of Infinite Processes and of Analytic Functions with an Account of the Principal Transcendental Functions, Cambridge University Press, Cambridge, UK, 4th edition, 1962.

[12] L. M. Campos, Generalized Calculus with Applications to Matter And Forces, CRC Press, Boca Raton, FL, USA, 2014.

[13] E. T. Copson, An introduction to the theory of functions of a complex variable, Oxford University Press, Oxford, UK, 1935.

[14] E. G. C. Poole, Introduction to The Theory of Linear Differential Equations, vol. 21, Oxford University Press, Oxford, UK, 1936.

[15] Z. X. Wang and D. R. Guo, Special Functions, World Scientific, Farrer Road, Singapore, 1989.

[16] C. G. Neumann, Theorie der Bessel'schen Funktionen. Ein Analogon zur Theorie der Kugelfunctionen, BG Teubner Verlag, Leipzig, Germany, 1867.

[17] A. R. Forsyth, Theory of differential equations, vol. 1, Cambridge University Press, Cambridge, UK, 1890-1906.

[18] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover Publications, New York, NY, USA, 9th edition, 1970.

[19] F. Oberhettinger, "On some expansions for Bessel integral functions," Journal of Research of the National Bureau of Standards, vol. 59, pp. 197-201, 1957.

[20] P. R. Parthasarathy and N. Balakrishnan, "A continued fraction approximation of the modified Bessel function I1t," Applied Mathematics Letters, vol. 3, no. 4, pp. 13-15, 1990.

[21] J. Harrison, "Fast and accurate Bessel function computation," in Proceedings of the 2009 IEEE 19th IEEE Symposium on Computer Arithmetic (ARITH), J. D. Bruguera, M. Cornea, D. DasSarma, and J. Harrison, Eds., pp. 104-113, IEEE, Portland, OR, USA, June 2009.

[22] A. Baricz and T. K. Pogany, "Turan determinants of Bessel functions," Forum Mathematicum, vol. 26, no. 1, pp. 295-322, 2014.

[23] A. Baricz and T. K. Pogany, "On a sum of modified Bessel functions," Mediterranean Journal of Mathematics, vol. 11, no. 2, pp. 349-360, 2014.

[24] G. Pang and S. Tang, "Approximate linear relations for Bessel functions," Communications in Mathematical Sciences, vol. 15, no. 7, pp. 1967-1986, 2017.

[25] P. Martin, J. Olivares, and F. Maass, "Analytic approximation for the modified Bessel function," in Journal of Physics: Conference Series, vol. 936, p. 012020, IOP Publishing, Pafos, Cyprus, 1 edition, 2017.

[26] L. Galue, "A generalized Bessel function," Integral Transforms and Special Functions, vol. 14, no. 5, pp. 395-401, 2003.

[27] G. Dattoli, G. Maino, C. Chiccoli, S. Lorenzutta, and A. Torre, "A unified point of view on the theory of generalized Bessel functions," Computers & Mathematics with Applications, vol. 30, no. 7, pp. 113-125, 1995.

[28] L. M. B. C. Campos, Transcendental Representations with Applications to Solids and Fluids, CRC Press, Boca Raton, FL, USA, 2012.

L. M. B. C. Campos, (1) F. Moleiro (iD), (1) M. J. S. Silva, (2) and J. Paquim (2)

(1) Center for Aeronautical and Space Science and Technology (CCTAE), IDMEC, LAETA, Instituto Superior Tecnico, Universidade de Lisboa, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal

(2) Instituto Superior Tecnico, Universidade de Lisboa, Av. Rovisco Pais 1,1049-001 Lisboa, Portugal

Correspondence should be addressed to F. Moleiro; filipa.moleiro@tecnico.ulisboa.pt

Received 15 July 2018; Accepted 28 October 2018; Published 4 November 2018

Academic Editor: Emilio Turco

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Campos, L.M.B.C.; Moleiro, F.; Silva, M.J.S.; Paquim, J. |

Publication: | Advances in Mathematical Physics |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 5552 |

Previous Article: | Adaptive Fuzzy Synchronization of Fractional-Order Chaotic Neural Networks with Backlash-Like Hysteresis. |

Next Article: | Plasma Parameters Effects on Dust Acoustic Solitary Waves in Dusty Plasmas of Four Components. |

Topics: |