arxiv:math/0603516v1 [math.ca] 1 Mar 006 THE FOURTH-ORDER TYPE LINEAR ORDINARY DIFFERENTIAL EQUATIONS W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ Abstract. This note reports on the recent advancements in the search for explicit representation, in classical special functions, of the solutions of the fourth-order linear ordinary differential equations named Bessel-type, Jacobi-type, Laguerre-type, Legendre-type. Contents 1. Introduction 1. Bessel-type differential equation 3. Laguerre-type differential equation 3 4. Legendre-type differential equation 5 5. Jacobi-type differential equation 6 6. Maple files 10 7. Acknowledgements 10 References 10 1. Introduction The four fourth-order type linear ordinary differential equations are named after the four classical second-order equations of Bessel, Jacobi, Laguerre, Legendre. The three fourth-order equations Jacobi-type, Laguerre-type, Legendre-type, and their associated orthogonal polynomials, were first defined by H.L. Krall in 1940, see [10] and [11], and later studied in detail by A.M. Krall in 1981, see [9], and by Koornwinder in 1984, see [1]. The structured definition of the general-even order Bessel-type special functions is dependent upon the Jacobi and Laguerre classical orthogonal polynomials, and the Jacobi-type and Laguerre-type orthogonal polynomials; see [8]. The properties of the fourth-order Bessel-type functions have been studied in the papers [3], [5] and [6]. This short paper reports on the recent progress that has been made in the search for representation of solutions of these four type differential equations in terms of classical special functions. Date: 1 March 006. 000 Mathematics Subject Classification. Primary: 33C15, 33D15, 33F10; Secondary: 33C05, 34B05. Key words and phrases. Ordinary linear differential equations, special functions, Bessel-type, Jacobi-type, Laguerre-type, Legendre-type equations. 1
W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ For the Bessel-type equation two independent solution were obtained, see [8] and [3], and then a complete set of four independent solutions in the papers [14] and [4]; these solutions are dependent upon the classical Bessel functions J r,y r,i r,k r for r = 0,1. Subsequently, using the techniques developed in [13], solutions of the three Jacobi-type, Laguerre-type, Legendre-type equations have been obtained in terms of the solutions of the corresponding classical differential equations. In all these four cases the verification of these representations as solutions of the type differential equations, is a formal process using the computational methods in the computer program Maple. Copies of the corresponding Maple.mws files, for the Jacobi-type, Laguerretype, Legendre-type equations, are available; see Section 6 below. There are four subsequent sections devoted to reporting on the form of these solutions, for each of the type fourth-order differential equations. The four type differential equations are all written in the Lagrange symmetric formally self-adjoint) form with spectral parameter Λ or ; if further properties of the solutions of these equations, involving the theory of special functions, are to be studied, then these spectral notations may have to be changed.. Bessel-type differential equation The fourth-order Bessel-type differential equation takes the form.1) xy x)) 9x 1 +8M 1 x)y x)) = Λxyx) for all x 0, ) where M 0, ) is a positive parameter and Λ C, the complex field, is a spectral parameter. The differential equation.1) is derived in the paper [8, Section 1, 1.10a)]. The manyspectralpropertiesofthisequationintheweighted HilbertfunctionspaceL 0, );x) are studied in the papers [3], [5] and [6]; for a collected account see the survey paper [4]. Our knowledge of the special function solutions of the Bessel-type differential equation.1) is now more complete than at the time the paper [8] was written. However, the results in [8, Section 1, 1.8a)], with α = 0, show that the function defined by.) J 0,M x) := [1+M/) ]J 0 x) M/) x) 1 J 1 x) for all x 0, ), is a solution of the differential equation.1), for all C, and hence for all Λ C, and all M > 0. Here: i) the parameter M > 0 ii) the parameter C iii) the spectral parameter Λ and the parameter M, in the equation.1), and the parameters M and, in the definition.), are connected by the relationship.3) Λ Λ,M) = +8M 1 ) for all C and all M > 0 iv) J 0 and J 1 are the classical Bessel functions of the first kind), see [15, Chapter III]. Similar arguments to the methods given in [8] show that the function defined by.4) Y 0,M x) := [1+M/) ]Y 0 x) M/) x) 1 Y 1 x) for all x 0, ),
FOURTH-ORDER TYPE EQUATIONS 3 is also an independent solution of the differential equation.1), for all C, and hence for all Λ C and all M > 0; here, again, Y 0 and Y 1 are classical Bessel functions of the second kind), see [15, Chapter III]. These earlier studies of the fourth-order differential equation.1) failed to find any explicit form of two linearly independent solutions, additional to the solutions J 0,M and Y 0,M. However, results of van Hoeij, see [13] and [14], using the computer algebra program Maple have yielded the required two additional solutions, here given the notations of I 0,M and K 0,M, with explicit representation in terms of the classical modified Bessel functions I 0,K 0 and I 1,K 1. These two additional solutions are defined as follows, where as far as possible we have followed the notation used for the solutions J 0,M and Y 0,M, i) given C, with arg) [0,π), M 0, ) and using the principal value of, define.5) c c,m) := +8M 1 and d d,m) := 1+M/).6) ii) define the solution, for all x 0, ), I 0,M x) := di 0 cx)+ 1 cmx 1 I 1 cx) := [1+M/) ]I 0 x ) +8M 1 + M +M/) )x 1 I 1 x ).7) +8M 1 iii) define the solution, for all x 0, ),.8) K 0,M x) := dk 0 cx)+ 1 cmx 1 K 1 cx).9) := [1+M/) ]K 0 x ) +8M 1 + M +M/) )x 1 K 1 x ) +8M 1. Remark.1. We have 1) The four linearly independent solutions J 0,M,Y 0,M,I 0,M,K 0,M provide a basis for all solutions of the original differential equation.1), subject to the Λ, ) connection given in.3). ) These four solutions are real-valued on their domain 0, ) for all R. 3) The domain 0, ) of the solutions J 0,M and I 0,M can be extended to the closed half-line [0, ) with the properties J 0,M 0) = I 0,M 0) = 1 for all R and all M 0, ). 3. Laguerre-type differential equation The Laguerre-type differential equation was discovered by H.L. Krall, see [10] and [11], and subsequently studied by other authors, see [9], [1] and [7]. The differential equation may be written in two forms; here the parameter A 0, ) and the spectral parameter C:
4 W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ 3.1) 3.) 1) The Frobenius form: x y 4) x)+ x +4x)y 3) x)+x A+6)x)y x) +A+)x A)y x) = yx) for all x 0, ). ) The Lagrange symmetric form formally self-adjoint form): x exp x)y x)) A+)x+)exp x)y x)) = exp x)yx) for all x 0, ). The spectral properties of the equation 3.) are considered in the weighted Hilbert function space L 0, );exp )). The van Hoeij method of searching for factors and solutions of this differential equation depends on the use of the computer algebra program Maple; see [13] and [14]. Following the solutions of the fourth-order Bessel-type obtained in [8, Section ], and then in [4] on using the methods of [14], it has proved possible to express all solutions of the equations 3.1) and 3.) in terms of the confluent hypergeometric functions. The special functions involved are the Whittaker functions M κ,µ z) and W κ,µ z) as defined in the compendium [1, Chapter 13, Sections 13.1.31 to 13.1.34]. Define the functions Φ and Ψ, for the range of variables given in [1, Chapter 13]: 3.3) Φκ,µ;z) := M κ,µ z) 3.4) Ψκ,µ;z) := W κ,µ z). Given the parameters A 0, ) and C define 3.5) Γ,A) := 4A +4A+1+4. We designate the four linearly independent solutions of 3.1), equivalently 3.), again the parameters A 0, ) and C, by L r,a;x) for all x 0, ) and r = 1,,3,4. These solutions are defined as follows, where d/dx denotes differentiation with respect to the independent variable x, L 1 3.6) 3.7) 3.8) L L 3 L 1,A;x) := 1 + 1Γ,A)) x 1/ exp 1x)Φ A 1Γ,A),0;x) d x 1/ exp 1 dx x)φ A 1Γ,A),0;x)). L,A;x) := 1 1Γ,A)) x 1/ exp 1x)Φ A+ 1Γ,A),0;x) d x 1/ exp 1 dx x)φ A+ 1 Γ,A),0;x)). L 3,A;x) := 1 + 1 Γ,A)) x 1/ exp 1 x)ψ A 1 Γ,A),0;x) d x 1/ exp 1 dx x)ψ A 1Γ,A),0;x)).
FOURTH-ORDER TYPE EQUATIONS 5 L 4 3.9) L 4,A;x) := 1 1Γ,A)) x 1/ exp 1x)Ψ A+ 1Γ,A),0;x) d x 1/ exp 1 dx x)ψ A+ 1Γ,A),0;x)). 4. Legendre-type differential equation The Legendre-type differential equation was discovered by H.L. Krall, see [10] and [11], and subsequently studied by other authors, see [9], [1] and [7]. The differential equation may be written in two forms; here the parameter A 0, ) and the spectral parameter C: 1) The Frobenius form: 4.1) x 1) y 4) x)+8xx 1)y 3) x)+4a+1)x 1)y x) +8Axy x) = yx) for all x 1,+1). ) The Lagrange symmetric form formally self-adjoint form): 4.) 1 x ) y x)) 8+4A1 x ))y x)) = yx) for all x 1,+1). The spectral properties of the equation 4.) are considered in the Hilbert function space L 1,+1). The van Hoeij method of searching for factors and solutions of this differential equation depends on the use of the computer algebra program Maple; see [13] and [14]. Following the solutions of the fourth-order Bessel-type obtained in [8, Section ], and then in [4] on using the methods of [14], it has proved possible to express all solutions of the equations 4.1) and 4.) in terms of the confluent hypergeometric functions. The special functions involved are the Legendre functions P ν µ z) and Qµ ν z) as defined in the compendium [1, Chapter 8, Section 8.1]. Define the functions Pν,x) and Qν,x), for the range of variables given in [1, Chapter 8], by: 4.3) Pν,x) := P 0 ν x) 4.4) Qν,x) := Q 0 νz). Given the parameters A 0, ) and C define, noting the use of the ± symbol, 4.5) Γ ±,A) := ± 4A 4A+1+ and Ω ±,A) := 5 8A+4Γ ±,A) We designate the four linearly independent solutions of 4.1), equivalently 4.), again the parameters A 0, ) and C, by 4.6) Le r,a;x) for all x 1,+1) and r = 1,,3,4.
6 W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ These solutions are defined as follows, where d/dx denotes differentiation with respect to the independent variable x, 4.7) Le 1,A;x) := 1 1+Ω+,A))xP 1 9 8A+4Γ+,A)+4Ω +,A) 1,x ) + [ +3 4A+4A )+Ω +,A) 3Γ +,A)+Ω +,A)Γ +,A)++4A+4A )x ] +4A+4A ) 1 d dx P1 9 8A+4Γ+,A)+4Ω +,A) 1,x). 4.8) 1 1+Ω,A))xP 1 Le,A;x) := ) 9 8A+4Γ,A)+4Ω,A) 1,x + [ +3 4A+4A )+Ω,A) 3Γ,A)+Ω,A)Γ,A)++4A+4A )x ] +4A+4A ) 1 d dx P1 9 8A+4Γ,A)+4Ω,A) 1,x). 4.9) 1 1+Ω+,A))xQ 1 Le 3,A;x) := ) 9 8A+4Γ+,A)+4Ω +,A) 1,x + [ +3 4A+4A )+Ω +,A) 3Γ +,A)+Ω +,A)Γ +,A)++4A+4A )x ] +4A+4A ) 1 d dx Q1 9 8A+4Γ+,A)+4Ω +,A) 1,x). 4.10) 1 1+Ω,A))xQ 1 Le 4,A;x) := ) 9 8A+4Γ,A)+4Ω,A) 1,x + [ +3 4A+4A )+Ω,A) 3Γ,A)+Ω,A)Γ,A)++4A+4A )x ] +4A+4A ) 1 d dx Q1 9 8A+4Γ,A)+4Ω,A) 1,x). 1) It has been verified, using Maple 9, that Le 1,Le,Le 3,Le 4 are solutions of the differential equations 4.1) and 4.) for the domain x 1,+1), and the parameter values A 0, ) and C; this verification follows from the original results of van Hoeij. ) The factor +4A+4A ) 1 which appears in the derivative term of all four of the solutions Le r, for r = 1,,3,4, presents an apparent singularity in these solutions when = 4A+4A ); however Maple 9 determines that the other factors [ +3 4A+4A )+Ω ±,A) 3Γ ±,A)+Ω ±,A)Γ ±,A)++4A+4A )x ] both have a zero at this value = 4A+4A ). 5. Jacobi-type differential equation This note concerns information about the van Hoeij solutions to the fourth-order Jacobitype linear ordinary differential equation.
FOURTH-ORDER TYPE EQUATIONS 7 The Jacobi-type differential equation was discovered by H.L. Krall, see [10] and [11], and subsequently studied by other authors, see [9], [1] and [7]. The differential equation may be written in two forms; here the parameters A 0, ) and α 1, ), and the spectral parameter C: 1) The Frobenius form: 1 x ) y 4) x) 1 x )α+4)x+α)y 3) x) +1+x)4A 5.1) α +α +9α+14)x+ 4A α +α 3α 10))y x) +4Aα α +8A α +α +6α+4)x+4Aα α +α +6α+4)y x) = yx) for all x 1,+1). ) The Lagrange symmetric form formally self-adjoint form): J[y]x) := 1 x) α+ 1+x) y x)) 5.) 1 x) α+1 4A α +α+)x+4a α +α+6)y x)) = 1 x) α yx) for all x 1,+1). The spectral properties of the equation 5.) are considered in the weighted Hilbert function space L 1,+1);1 x) α ). The van Hoeij method of searching for factors and solutions of this differential equation depends on the use of the computer algebra program Maple; see [13] and [14]. Following the solutions of the fourth-order Bessel-type obtained in [8, Section ], and then in [4] on using the methods of [14], it has proved possible to express all solutions of the equations 5.1) and 5.) in terms of hypergeometric functions. The special functions involved are the JacobiP functions; these functions are defined in terms of the hypergeometric functions p F q ; for certain values of the parameters involved in this definition, they reduce to the classical Jacobi polynomials P n α,β). The notation JacobiP is not to be found in [1, Chapter ] but is given in the program Maple9; here we shorten the notation JacobiP and define the function JP as follows, in terms of the classical Γ and F 1 functions: 5.3) JacobiPν,α,β;z) JPν,α,β;z) Γν +α+1) := Γν +1)Γα+1) F 1 ν,ν +α+β +1;α+1;1 z)/). Here, the independent variables ν,α,β,z are in the complex field C, but for the results in this note we make the restrictions ν C, α 1, ), β = 0, z R. When ν = n N 0 then we have the connection, see [1, Section.] and [, Lecture,.)], P n α,β) x) = JPn,α,β;x) for all α,β 1, ) and all x 1,+1). The Jacobi-type orthogonal polynomials were discovered by H.L. Krall, see [10] and [11], and later studied by A.M. Krall, T.H. Koornwinder and L.L. Littlejohn. These polynomials are denoted by Pn,A α x) for α 1, ),A 0, ),n N 0 and x 1,+1); each Pn,A α is real-valued on 1,+1) and of degree n. These polynomials have the properties:
8 W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ i) For all n N 0 the polynomial P α n,a is a solution of the differential equation, see 5.) 5.4) J[Pn,A α ]x) = nα,a)1 x) α Pn,A α x) for all x 1,+1) where the eigenvalue n α,a) is determined by 5.5) n α,a) = nn+α+1)n +α+1)n+4a α +α) for all n N 0. 5.6) ii) The collection {Pn,A α ) : n N 0} is a complete, orthogonal set in the Lebesgue- Stieltjes function space L 1,+1);µ Jac ), where the Borel measure µ Jac is determined by the monotonic non-decreasing function ˆµ Jac, for all α 1, ) and all A 0, ), ˆµ Jac x) := 1/ for all x, 1) := A α+1) [α+1 1 x) α+1 ] for all x [ 1,+1] := Aα+1 α+1) for all x +1, ). The orthogonality of the set {Pn,A α : n N 0} then appears as Pm,A α x)pα n,a x) dµ Jacx) 5.7) [ 1,+1] = 1 Pα m,a 1)Pα n,a 1)+ 1 +1 A Pm,A α x)pα n,a x)1 x)α dx 1 = 0 for all m,n N 0 with m n. The van Hoeij solutions for the general differential equation 5.) begin with solutions to the set of special equations 5.4) when the general spectral parameter takes the value of one of the eigenvalues { n α,a) : n N 0 }. Given n N 0 define, for all x 1,+1), with JP given by 5.3): 5.8) S 1,n x) := nα+a α +n+n )JPn,α,0;x) +1 x) d dx JPn,α,0;x). Then the van Hoeij methods determine that S 1,n is a solution of 5.1) for this value of n, i.e. for = n α,a). Remark 5.1. From the definition 5.8) it would appear the this solution S 1,n is linearly dependent upon the Jacobi-type polynomial P α n,a. To proceed further, given C, let {ρ r ) C : r = 1,,3,4} denote the roots of the quartic polynomial 5.9) ρρ+α+1)ρ +α+1)ρ+4a α +α) = 0.
FOURTH-ORDER TYPE EQUATIONS 9 Maple gives explicit representations for these roots as follows: ρ 1 α,a,) = 1 α 1+ ) α +1 8A α +ξα,a,) ρ α,a,) = 1 α 1 ) α +1 8A α +ξα,a,) 5.10) ρ 3 α,a,) = 1 α 1+ ) α +1 8A α ξα,a,) ρ 4 α,a,) = 1 α 1 ) α +1 8A α ξα,a,), where ξα,a,) := α +8αA α +16A α +4. The van Hoeij solutions to the Jacobi-type equation, see 5.10), 5.11) J[y]x) = 1 x) α yx) for all x 1,+1) now develop as follows. With C, replace n in the solution S 1,n with the first root ρ 1 ) taken from 5.10); this yields a solution J 1 in terms of the function JP, for all α 1, ) and all A 0, ), 5.1) J 1 α,a,;x) := 1 [α ξα,a,)]jpρ 1α,A,),α,0;x) 1 x) d dx JPρ 1α,A,),α,0;x) for all x 1,+1). Similarly, for C, there is a solution J derived from the third root ρ 3 ), for all α 1, ) and all A 0, ), J α,a,;x) := 1 5.13) [α+ξα,a,)]jpρ 3α,A,),α,0;x) 1 x) d dx JPρ 3α,A,),α,0;x) for all x 1,+1). Given n N 0 define, for all x 1,+1), there is a second solution S,n, similar to S 1,n of 4.1), defined by 5.14) S,n x) := [n+1)α+a α+1 +n+n )]x 1) α JP 1 n, α,0;x) x 1) 1 α d dx JP 1 n, α,0;x). Then the van Hoeij methods determine that S,n is a solution of 5.1) for this value of n, i.e. for = n α,a). With C, replace n in the solution S,n with the first root ρ 1 ) taken from 5.10); this yields a solution J 3 in terms of the function JP, for all α 1, ) and all A 0, ), 5.15) J 3 a,a,;x) := 1 [α+ξα,a,)] 1+x) α JP 1 ρ 1 α,a,), α,0;x) + 1+x) 1 α d dx JP 1 ρ 1α,A,), α,0;x) for all x 1,+1).
10 W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ Similarly, for C, there is a solution J 4 derived from the third root ρ 3 ), for all α 1, ) and all A 0, ), 5.16) J 4 α,a,;x) := 1 [α ξα,a,)] 1+x) α JP 1 ρ 3 α,a,), α,0;x) + 1+x) 1 α d dx JP 1 ρ 3α,A,), α,0;x) for all x 1,+1). 6. Maple files The formal validity of these solution results were originally obtained by van Hoeij, following the methods given in [13] and [14], using the computer program Maple. These results have been subsequently confirmed by David Smith in the Maple program at the University of Birmingham. The Maple.mws files for the following three cases are named as follows: Jacobi-type-solutions.mws Laguerre-type-solutions.mws Legendre-type-solutions.mws Copies of these files may be obtained on application to David Smith at smithd@for.mat.bham.ac.uk 7. Acknowledgements The authors indicate here that the results reported on in this note follow from the combined efforts of the many collaborators whose names are given in the references below. The author Norrie Everitt is particularly grateful for the advice and help received from Dr David Smith, School of Mathematics at the University of Birmingham, on the application and use of the computer program Maple. The author Mark van Hoeij is supported by National Science Foundation grant 0511544. References [1] M. Abramowitz and I.A. Stegun. Handbook of mathematical functions. Dover Publications, New York: 1965.) [] R. Askey. Orthogonal Polynomials and Special Functions. SIAM, Philadelphia, USA: 1975.) [3] J. Das, D.B. Hinton, L.L. Littlejohn and C. Markett. The fourth-order Bessel-type differential equation. Applicable Analysis. 83 004), 35-36. [4] W.N. Everitt. Fourth-order Bessel-type special functions: a survey. Submitted for proceedings of the International Conference on Difference Equations, Special Functions and Applications; Technical University Munich, Germany: 5 to 30 July 005.) [5] W.N. Everitt, H. Kalf, L.L. Littlejohn and C. Markett. Additional properties of the fourth-order Besseltype differential equation. Math. Nachr. 78 005), 1538-1549. [6] W.N. Everitt, H. Kalf, L.L. Littlejohn and C. Markett. The fourth-order Bessel equation; eigenpackets and a generalised Hankel transform. Submitted to Integral Transforms and Special Transforms.) [7] W.N. Everitt, A.M. Krall, L.L. Littlejohn and V.P. Onyango-Otieno. Differential operators and the Laguerre type polynomials. SIAM J. Math. Anal. 3 199), 7 736. [8] W.N. Everitt and C. Markett. On a generalization of Bessel functions satisfying higher-order differential equations. Jour. Computational Appl. Math. 54 1994), 35-349.
FOURTH-ORDER TYPE EQUATIONS 11 [9] A.M. Krall. Orthogonal polynomials satisfying fourth order differential equations. Proc. Roy. Soc. Edinburgh A) 87 1981), 71-88. [10] H.L. Krall. Certain differential equations for Tchebycheff polynomials. Duke Math. J. 4 1938) 705-718. [11] H.L. Krall. On orthogonal polynomials satisfying a certain fourth order differential equation. The Pennsylvania State College Studies: No. 6. The Pennsylvania State College, State College, PA.: 1940.) [1] T.H.Koornwinder.Orthogonalpolynomialswithweightfunction1 x) α 1+x) β +Mδx+1)+Nδx 1). Canad. Math. Bull. 7 ) 1984), 05-14. [13] M. van Hoeij. Formal solutions and factorization of differential operators with power series coefficients. J. Symbolic Comput. 4 1997), 1-30. [14] M. van Hoeij. Personal contribution. International Conference on Difference Equations, Special Functions and Applications; Technical University Munich, Germany: 5 to 30 July 005.) [15] G.N. Watson. A treatise on the theory of Bessel functions. Cambridge University Press; second edition: 1950.) W.N. Everitt, School of Mathematics and Statistics, University of Birmingham, Edgbaston, Birmingham B15 TT, England, UK E-mail address: w.n.everitt@bham.ac.uk D. J. Smith, School of Mathematics and Statistics, University of Birmingham, Edgbaston, Birmingham B15 TT, England, UK E-mail address: smithd@for.mat.bham.ac.uk M. van Hoeij, Department of Mathematics, Florida State University, 08 James J. Love Building, Tallahassee, FL 3306-4510, USA E-mail address: hoeij@math.fsu.edu