arxiv: v9 [math-ph] 6 Nov 2014

Transcription

1 Generalization of the three-term recurrence formula and its applications Yoon Seok Choun Baruch College, The City University of New York, Natural Science Department, A506, 7 Lexington Avenue, New York, NY 000 arxiv: v9 math-ph] 6 Nov 204 Abstract The history of linear differential equations is over 350 years By using Frobenius method and putting the power series expansion into linear differential equations, the recursive relation of coefficients starts to appear There can be between two and infinity number of coefficients in the recurrence relation in the power series expansion During this period mathematicians developed analytic solutions of only two term recursion relation in closed forms Currently the analytic solution of three term recurrence relation is unknown In this paper I will generalize the three term recurrence relation in the linear differential equation This paper is 2nd out of 0 in series Special functions and three term recurrence formula 3TRF The next paper in series deals with the power series expansion in closed forms of Heun function 3] The rest of the papers in the series show how to solve mathematical equations having three term recursion relations and go on producing the exact solutions of some of the well known special functions including: Mathieu, Heun, Biconfluent Heun and Lame equations See section 7 for all the papers and short descriptions in the series Keywords: Three-term recurrence relation; Linear ordinary differential equation; Polynomial; Infinite series 2000 MSC: 33E30, 34A30, 39A05, 39A0 Introduction Mathieu functions6], is an example of three term recurrence relation appears in physical problems involving elliptical shapes7] or periodic potentials, were introduced by Mathieu 8682] when he investigated the vibrating elliptical drumhead Mathieu function has been described using numerical approximations Whittaker 943], Frenkel and Portugal 200] Whittaker tried to obtain the analytic solution of Mathieu equation using Floquet s theorem and he reached the conclusion using three term recursive relation He did not represent the solution of Mathieu equation in closed forms because of its three term Correspondence to: Baruch College, The City University of New York, Natural Science Department, A506, 7 Lexington Avenue, New York, NY 000 address: YoonChoun@baruchcunyedu; ychoun@gccunyedu; ychoun@gmailcom Yoon Seok Choun Preprint submitted to Elsevier November 7, 204

2 recursive coefficients He argued that While the general character of the solution from the function theory point of view is thus known, its Mathieu] actual analytical determination presents great difficulties The chief impediment is that the constant µ the Mathieu exponent cannot readily be found in terms of a and q 3]: a and q see in Ref3] are corresponding toλand 8q see in Ref5] In my opinion, Mathieu functions are difficult to represent in analytic closed forms and in its integral forms because of the three recursive coefficients Heun function, is an example of three term recurrence relation, generalizes all well-known special functions such as: Spheroidal Wave, Lame, Mathieu, and hypergeometric 2 F, F and 0F functions The Heun equation has four kinds of confluent forms such as Confluent Heun, Doubly-Confluent Heun, Biconfluent Heun and Triconfluent Heun equations Mathieu equation is the special case of Confluent Heun equation According to Hortacsu, Heun differential equation8, 9] starts to appear in astronomy or in general relativity0] inevitably Heun equation is a general equation of all well-known special functions; Mathieu, Lame and Coulomb spheroidal equations The power series expansion on Heun function can not be described as two term recursion relation any more The coefficients in a power series expansions of Heun equation have a recursive relation between a 3-term 5] For the past 350 years, we only have constructed the power series expansion in closed forms using the two term recurrence relation in linear differential equation However, since 930, we do not have simple problems to solve in theoretical particle physics and scientists and mathematicians doing research on this field have to tackle more difficult problems, either with more difficult metrics or in higher dimensions Most of the difficult problems must include three term or more Hortacsu 20 5] In this paper I will construct three term recurrence relation in the form of power series expansion for the cases of polynomial and infinite series In the following papers see section 7 using the method described in this paper I will show how to obtain exact solution in the power series expansion, 2 in the integral formalism and 3 in the generating function of any linear ordinary differential equations having three term recurrence relation Using the method described in this paper one might be able to obtain exact solutions of higher term recurrence formula; four, five,, m th, term recurrence formulas If this is possible, then one would be able to generalize all homogeneous and inhomogeneous linear differential equations Most problems in nature turns out to be nonlinear For simplification purposes we usually try to linearize these nonlinear system The systems are linearized using certain methods of simplification resulting in better approximation of physical systems All physics theories E & M, Newtonian mechanics, quantum mechanic, QCD, supersymmetric field theories, string theories, general relativity, etc, generally involve solutions of linear differential equations Unfortunately in some instances there are no analytic solutions to these physical systems Hopefully the methods and procedures described in these papers see section 7 will help obtain better analytic solutions to linearized nonlinear systems 2 Lame equation, Frobenius method and three term recurrence relation There are no exact solutions in closed forms for second order or higher order ordinary linear differential equations consisting of three term recurrence relation: some of the examples are Lame function, the generalized Lame function, Heun s equation, GCH function4], Mathieu function, etc Two term recurrence relation are solvable Some examples are: Legendre function, hypergeometric function, Kummer function, Bessel function, etc Let s think about an example 2

3 which has no analytic solution in closed forms in detail d 2 y dt 2+ 2 t a + t b + t c dy dt αα+t+β y=0 2 tt at bt c 2 is Lame differential equation Ifαis not positive integer, the solution of it is called as the generalized Lame function Replace t by x+a in 2 By using the function yx as Frobenious series in it, yx= c n t a n+λ = c n x n+λ 22 λ is an indicial root Plug 22 in 2 And its recurrence relation is where, c n+ = A n c n + B n c n ; n 23 K n = A n + B n K n ; n 24a K n = c n+ c n K n = c n c n 24b A n = αα+a+β} 22 2a b cn+λ a ba cn+λ+n+λ+ 2 24c B n = α+2n+λ }α 2n+λ+2} 2 2 a ba cn+λ+n+λ+ 2 24d All other differential equations having no analytic solution in closed forms can be described as in 23 If I get a formula of 23 type, I can apply it to all other functions having no analytic solution in closed forms to diverse areas such as Lame function, generalized Lame function, Mathieu function, Heun function, GCH function4], etc 3 Infinite series Assume that c = A 0 c is a necessary boundary condition The three term recurrence relation in all differential equations having no analytic solution in closed forms follow 3 a i n=a i B n = where a i is positive integer including is also a necessary condition Every differential equations having no analytic solution in closed forms also take satisfied with 32 3

4 My definition of B i, j,k,l refer to B i B j B k B l Also, A i, j,k,l refer to A i A j A k A l For,, 2, 3,, 23 gives c 0 c = A 0 c 0 c 2 = B + A 0, c 0 c 3 = A 0 B 2 + A 2 B + A 0,,2 c 0 c 4 = B,3 + A 0, B 3 + A 0,3 B 2 + A 2,3 B + A 0,,2,3 c 0 c 5 = A 0 B 2,4 + A 2 B,4 + A 4 B,3 + A 0,,2 B 4 + A 0,,4 B 3 + A 0,3,4 B 2 + A 2,3,4 B + A 0,,2,3,4 c 0 c 6 = B,3,5 + A 0, B 3,5 + A 0,3 B 2,5 + A 0,5 B 2,4 + A 2,3 B,5 + A 2,5 B,4 + A 4,5 B,3 +A 0,,2,3 B 5 + A 0,,2,5 B 4 + A 0,,4,5 B 3 + A 0,3,4,5 B 2 + A 2,3,4,5 B + A 0,,2,3,4,5 c 0 c 7 = A 0 B 2,4,6 + A 2 B,4,6 + A 4 B,3,6 + A 6 B,3,5 + A 0,,2 B 4,6 + A 0,,4 B 3,6 + A 0,,6 B 3,5 + A 0,3,4 B 2,6 +A 0,3,6 B 2,5 + A 0,5,6 B 2,4 + A 2,3,4 B,6 + A 2,3,6 B,5 + A 2,5,6 B,4 + A 4,5,6 B,3 +A 0,,2,3,4 B 6 + A 2,3,4,5,6 B + A 0,3,4,5,6 B 2 + A 0,,4,5,6 B 3 + A 0,,2,5,6 B 4 + A 0,,2,3,6 B 5 + A 0,,2,3,4,5,6 c 0 c 8 = B,3,5,7 + A 0, B 3,5,7 + A 0,3 B 2,5,7 + A 0,5 B 2,4,7 + A 0,7 B 2,4,6 + A 2,3 B,5,7 + A 2,5 B,4,7 + A 2,7 B,4,6 +A 4,5 B,3,7 + A 4,7 B,3,6 + A 6,7 B,3,5 + A 0,,2,3 B 5,7 + A 0,,2,5 B 4,7 + A 0,,2,7 B 4,6 + A 0,,4,5 B 3,7 +A 0,,4,7 B 3,6 + A 0,,6,7 B 3,5 + A 0,3,4,5 B 2,7 + A 0,3,4,7 B 2,6 + A 0,3,6,7 B 2,5 + A 0,5,6,7 B 2,4 + A 2,3,4,5 B,7 +A 2,3,4,7 B,6 + A 2,3,6,7 B,5 + A 2,5,6,7 B,4 + A 4,5,6,7 B,3 + A 0,,2,3,4,5 B 7 + A 0,,2,3,4,7 B 6 +A 0,,2,3,6,7 B 5 + A 0,,2,5,6,7 B 4 + A 0,,4,5,6,7 B 3 + A 0,3,4,5,6,7 B 2 + A 2,3,4,5,6,7 B + A 0,,2,3,4,5,6,7 c 0 33 In 33 the number of individual sequence c n follows Fibonacci sequence:,,2,3,5,8,3,2,34,55, The sequence c n consists of combinations A n and B n in 33 First observe the term inside parentheses of sequence c n which does not include any A n s: c n with even index c 0, c 2, c 4, a Zero term of A n s c 0 c 2 = B c 0 c 4 = B,3 c 0 c 6 = B,3,5 c 0 34 c 8 = B,3,5,7 c 0 When a function yx, analytic at x=0, is expanded in a power series, we write yx= c n x n+λ = y m x=y 0 x+y x+y 2 x+ 35 m=0 4

5 where y m x= c m l x l+λ 36 l=0 λ is the indicial root y m x is sub-power series that have sequence c n including m term of A n s in 33 For example y 0 x has sequences c n including zero term of A n s in 33, y x has sequences c n including one term of A n s in 33, y 2 x has sequences c n including two term of A n s in 33, etc Substitute 34 in 36 putting m=0 y 0 x=c 0 B 2i0 + x2n+λ 37 Observe the terms inside parentheses of sequence c n which include one term of A n s in 33: c n with odd index c, c 3, c 5, b One term of A n s c = A 0 c 0 c 3 = A 0 c 5 = A 0 c 7 =,6 A 0,6,8 c 9 = A 0 c = } +A2 B B c0 2 } +A2 B B 2 + A 4 B B,3 c0 2,4,6 + A2 B 2 +A2,6,8 B 2,6,8 +A 8 B B,3,5,7 2,4,6,8 }c 0 A 0,6,8,0 +A 8,6,8,0 B 2,4,6,8 +A2,6,8,0 B + A 4,6 B 2,4 B,3 + A 6,6 B 2,4,6 B,3,5 } c0 B + A 4,6,8 B 2,4,6,8,0 B B 2 + A 4 B 2,4,6,8,0 B,3,5,7 + A 0 B B,3,5,7,9 2,4,6,8,0 }c 0 B,3 + A 6,6,8 B 2,4,6 B,3,5,6,8,0 B,3 + A 6 B B,3,5 2,4, is simply c 2n+ = c 0 n i 0 A 2i 0 Substitute 39 in 36 putting m= n i 0 y x=c 0 A 2i 0 B 2i + B 2i2 +2 i =0 i 2 =i 0 B 2i + B 2i2 +2 i =0 i 2 =i 0 39 x2n++λ 30 Observe the terms inside parentheses of sequence c n which include two terms of A n s in 33: c n with even index c 2, c 4, c 6, 5

6 c Two terms of A n s c 2 = A 0, c 0 c 4 = A B,3 B,3 } 0 A B +A3 B,3 + B,3 } } A2 A3 B 2 B B,3 c 0 c 6 = A 0 A B,3,5 B B,3,5 B,3,5 } +A3 B +A5,3 B,3,5 c 8 = +A B,3,5 B,3,5 } B,3,5 } } 2 A3 B 2 B B,3 + A 5 B 2 B B,3,5 + A4 A5 B 2,4 B B,3,3,5 c 0 A 0 A B,3,5,7 B,3,5,7 B,3,5,7,6 B,3,5,7 } B +A3 B +A5,3 B +A7,3,5 B,3,5,7 +A B,3,5,7 B,3,5,7,6 B,3,5,7 } 2 A3 B 2 B B,3 + A 5 B 2 B B,3,5 + A 7 B 2 B B,3,5,7 B,3,5,7,6 B,3,5,7 } +A 4 A5 B 2,4 B B,3,3,5 + A 7 B 2,4 B B,3,3,5,7 +A 6 A7,6 B 2,4,6 B,3,5,7 B,3,5,7 B,3,5 } } c 0 3 is simply n n c 2n+2 = c 0 A i 0 2i 0 A 2i + B 2i2 + i =i 0 i 2 =0 n i 0 n = c 0 A 2i 0 B 2i + A 2i 2 + i =0 i 2 =i 0 i B 2i3 +2 B 2i4 +3 i 3 =i 0 i 4 =i i 2 B 2i3 +2 B 2i4 +3 i 3 =i 0 i 4 =i Substitute 32 in 36 putting m=2 n i 0 y 2 x=c 0 A 2i 0 B 2i + i =0 i 2 =i 0 n A 2i 2 + i 2 B 2i3 +2 B 2i4 +3 i 3 =i 0 i 4 =i 2 x2n+2+λ 33 Observe the terms inside parentheses of sequence c n which include three terms of A n s in 33: c n with odd index c 3, c 5, c 7, 6

7 d Three terms of A n s c 3 = A 0,,2 c 0 c 5 = A 0 A A2 B4 + A4 B3 B4 ] B 4 c 7 = c 9 = +A 3 A 4 B 2 B3 B4 } B 3 B 4 + A2 A3 A4 B B3 B4 ]} } B 2 B 3 B 4 c 0 A 0 A A2 B4,6 + A4 B3 B4,6 B 4 + A6 B3,5 B4,6 +A 3 A4 B 2 B3 B4,6 B 3 B 4 + A6 B 2 B3,5 B4,6 ] B 3 B 4,6 + A5 A6 B 2,4 B3,5 B4,6 ]} B 3,5 B 4,6 +A 2 A3 A4 B B 2 B3 B 3 B4,6 B 4 + A6 B B 2 B3,5 B 3 B4,6 B 4,6 ] B3,5 B4,6 ]} B3,5 B4,6 ]} } +A 5 A6 B B 2 B 3,5 B 4,6 + A4 A5 A6 B,3 B 2,4 B 3,5 B 4,6 c 0 A 0 A A2 + A4 B3 B 4 + A6 +A 3 A4 B 2 B3 +A 5 A6 B 2,4 B 3 B 4 + A6 B 2 B3,5 B 3,5 +A 2 A3 A4 B B 2 B3 B 3 B 4 B3,5 B 4,6 + A8 B 2,4 B 4,6 ] B 3 B 4,6 + A8 B 2 B 3 B3,5,7 ] B 3,5 B 4,6,8 + A7 A8 B 2,4,6 + A6 B B 2 B3,5 B 3 B 4,6 + A8 B B 2 B3,5,7 B 3 +A 5 A6 B B 2 B3,5 B 3,5 B 4,6 + A8 B B 2 B3,5,7 B 3,5 B 4,6,8 ] + A7 A8 B,6 B 2 +A 4 A5 A6 B,3 B 2,4 B3,5 B 3,5 B 4,6 + A8 B,3 B 2,4 B3,5,7 B 3,5 B 4,6,8 ] B3,5 B 4,6 + A8 B3,5,7 ] B 4,6,8 B3,5,7 ]} B 3,5,7 B 4,6,8 ] B 4,6,8 B3,5,7 ]} B 3,5,7 B 4,6,8 +A 7 A8 B,3,6 B 2,4 B3,5,7 B 3,5,7 B 4,6,8 ]} + A6 A7 A8 B,3,5,6 B 2,4,6 B3,5,7 B 3,5,7 B 4,6,8 ]} } c 0 B3,5,7 ] B 4,6,8 34 is simply n n c 2n+3 = c 0 A n 2i 0 A i 0 2i + A 2i 2 +2 i =i 0 i 2 =i i 3 =0 n i 0 n = c 0 A i 2 2i 0 B 2i + A 2i 2 + i =0 i 2 =i 0 i i 2 B 2i3 + B 2i4 +2 B 2i5 +3 B 2i6 +4 i 4 =i 0 i 5 =i i 6 =i 2 B 2i3 +2 i 3 =i 0 i 4 =i 2 n A 2i 4 +2 i 4 B 2i5 +3 B 2i6 +4 i 5 =i 2 i 6 =i

8 Substitute 35 in 36 putting m=3 n i 0 y 3 x = c 0 A 2i 0 i =0 n i 4 A 2i 4 +2 i 4 =i 2 B 2i + i 2 =i 0 B 2i5 +3 B 2i6 +4 i 5 =i 2 i 6 =i 4 n i 2 A 2i 2 + B 2i3 +2 i 3 =i 0 x2n+3+λ 36 By repeating this process for all higher terms of A s, we obtain every y m x terms where m 4 Substitute 37, 30, 33, 36 and including all y m x terms where m 4 into 35 Theorem The general expression of yx for infinite series is yx = y 0 x+y x+y 2 x+y 3 x+ n i 0 = c 0 B 2i0 + x2n+λ + A 2i 0 i =0 n i 0 N n + N=2 A 2i 0 B 2i + i =0 k= x2n+n+λ i 2N =i 2N B 2i2N +N+ B 2i + B 2i2 +2 i 2 =i 0 i 2k x2n++λ A 2i2k +k B 2i2k+ +k+ i 2k =i 2k i 2k+ =i 2k 37 4 Polynomial which makes B n term terminated Now let s investigate the polynomial case of yx Assume that B n is terminated at certain value of n Then each y i x where i=0,, 2, will be polynomial Examples of these are Heun s equation, GCH function4], Lame function, etc First B 2k+ is terminated at certain value of k I choose eigenvalueβ 0 which B 2k+ is terminated whereβ 0 = 0,, 2, B 2k+2 is terminated at certain value of k I choose eigenvalueβ which B 2k+2 is terminated whereβ = 0,, 2, Also B 2k+3 is terminated at certain value of k I choose eigenvalueβ 2 which B 2k+3 is terminated whereβ 2 = 0,, 2, By repeating this process I obtain B 2βi +i+= 0 where i,β i = 0,, 2, 4 In general, the two term recurrence relation for polynomial has only one eigenvalue: for example, the Laguerre function, confluent hypergeometric function, Legendre function, etc But three term recurrence relation has infinite eigenvalues which areβ i, where i,β i = 0,, 2, First observe the term in sequence c n which does not include any A n s in 34: c n with even index c 0, c 2, c 4, a Asβ 0 =0, then B =0 in 34 c 0 42 b Asβ 0 =, then B 3 =0 in 34 c 0 43 c 2 = B c 0 8

9 c Asβ 0 =2, then B 5 =0 in 34 c 0 Substitute 42,43 and 44 in 36 putting m=0 β 0 y 0 x=c 0 c 2 = B c 0 c 4 = B,3 c 0 44 B 2i + i =0 x2n+λ 45 Observe the terms inside curly brackets of sequence c n which include one term of A n s in 38: c n with odd index c, c 3, c 5, a Asβ 0 =0, then B =0 in 38 c = A 0 c 0 c 3 = A 0 B 2 c 0 c 5 = A 0 B 2,4 c 0 c 7 = A 0 B 2,4,6 c 0 46 c 9 = A 0 B 2,4,6,8 c 0 As i= in 4, Substitute 46 in 36 putting m= by using 47 β y 0 x=c 0A 0 B 2β +2= 0 whereβ = 0,, 2, 47 B 2i +2 i =0 x2n++λ 48 In 48 y 0 x is sub-power series, having sequences c n including one term of A n s in 33 as β 0 =0, for the polynomial case in which makes B n term terminated b Asβ 0 =, then B 3 =0 in 38 c = A 0 c 0 c 3 =A 0 B 2 +A 2 B }c 0 c 5 =A 0 B 2,4 +A 2 B 4 B }c 0 c 7 =A 0 B 2,4,6 +A 2 B 4,6 B }c 0 49 c 9 =A 0 B 2,4,6,8 +A 2 B 4,6,8 B }c 0 The first term in curly brackets of sequence c n in 49 is same as 46 Then, its solution is equal to 48 Substitute 47 into the second term in curly brackets of sequence c n in 49 Its power 9

10 series expansion including the first and second terms in curly brackets of sequence c n in 49, analytic at x = 0, is β y x=c 0 A β 0 B 2i +2 + A 2B B 2i +2 x2n++λ 40 i =0 In 40 y x is sub-power series, having sequences c n including one term of A n s in 33, as β 0 = for the polynomial case in which makes B n term terminated c Asβ 0 =2, then B 5 =0 in 38 c = A 0 c 0 c 3 =A 0 B 2 +A 2 B }c 0 n= i = c 5 =A 0 B 2,4 + A 2 B 4 B + A 4 B,3 }c 0 c 7 =A 0 B 2,4,6 +A 2 B 4,6 B + A 4 B 6 B,3 }c 0 c 9 =A 0 B 2,4,6,8 + A 2 B 4,6,8 B + A 4 B 6,8 B,3 }c 0 4 c =A 0 B 2,4,6,8 +A 2 B 4,6,8 B + A 4 B 6,8,0 B,3 }c 0 The first and second term in curly brackets of sequence c n in 4 is same as 49 Then its power series expansion is same as 40 Substitute 47 into the third term in curly brackets of sequence c n in 4 Its power series expansion including the first, second and third terms in curly brackets of sequence c n in 4, analytic at x=0, is β y 2 x = c 0 A β 0 B 2i +2 + A 2B B 2i +2 i =0 n= i = β + A 4 B,3 B 2i +2 x2n++λ 42 n=2 i =2 In 42 y 2 x is sub-power series, having sequence c n including one term of A n s in 33 as β 0 =2, for the polynomial case in which makes B n term terminated By repeating this process for allβ 0 = 3, 4, 5,, I obtain every y j x terms where j 3 According to 48, 40, 42 and every y j x where j 3, the general expression of y x for allβ 0 is β 0 i 0 β y x=c 0 A i 2 2i 0 B 2i + x2i 2++λ whereβ 0 β 43 i =0 i 2 =i 0 i 3 =i 0 B 2i3 +2 Observe the terms of sequence c n which include two terms of A n s in 3: c n with even index c 2, c 4, c 6, 0

11 a Asβ 0 =0, then B =0 in 3 c 2 = A 0, c 0 c 4 = A 0 A B 3 +A 3 B 2 }c 0 c 6 = A 0 A B 3,5 +A 3 B 2 B 5 +A 5 B 2,4 }c 0 c 8 = A 0 A B 3,5,7 +A 3 B 2 B 5,7 +A 5 B 2,4 B 7 +A 7 B 2,4,6 }c 0 c 0 = A 0 A B 3,5,7,9 + A 3 B 2 B 5,7,9 + A 5 B 2,4 B 7,9 +A 7 B 2,4,6 B A 9 B 2,4,6,8 }c 0 i Asβ =0, then B 2 =0 in 44 c 2 = A 0, c 0 c 4 = A 0 A B 3 c 0 c 6 = A 0 A B 3,5 c 0 c 8 = A 0 A B 3,5,7 c 0 45 c 0 = A 0 A B 3,5,7,9 c 0 As i=2 in 4, Substitute 45 in 36 putting m=2 by using 46 β 2 y 0,0 2 x=c 0A 0 A B 2β2 +3= 0 whereβ 2 = 0,, 2, 46 B 2i +3 i =0 x2n+2+λ 47 In 47 y 0,0 2 x is sub-power series, having sequences c n including two term of A n s in 33 as β 0 =0 andβ =0, for the polynomial case in which makes B n term terminated ii Asβ =, then B 4 =0 in 44 c 2 = A 0, c 0 c 4 = A 0 A B 3 + A 3 B 2 }c 0 c 6 = A 0 A B 3,5 +A 3 B 2 B 5 }c 0 c 8 = A 0 A B 3,5,7 + A 3 B 2 B 5,7 }c 0 48 c 0 = A 0 A B 3,5,7,9 +A 3 B 2 B 5,7,9 }c 0 The first term in curly brackets of sequence c n in 48 is same as 45 Then its power series expansion is equal to 47 Substitute 46 into the second term in curly brackets of sequence

12 c n in 48 Its power series expansion including the first and second terms in curly brackets of sequence c n, analytic at x=0, is y 0, 2 x=c 0A 0 A β 2 B 2i +3 i =0 β 2 + A 3B 2 n= B 2i +3 i = x2n+2+λ 49 In 49 y 0, 2 x is sub-power series, having sequences c n including two term of A n s in 33 as β 0 =0 andβ =, for the polynomial case in which makes B n term terminated iii Asβ =2, then B 6 =0 in 44 c 2 = A 0, c 0 c 4 = A 0 A B 3 + A 3 B 2 }c 0 c 6 = A 0 A B 3,5 +A 3 B 2 B 5 +A 5 B 2,4 }c 0 c 8 = A 0 A B 3,5,7 + A 3 B 2 B 5,7 +A 5 B 2,4 B 7 }c c 0 = A 0 A B 3,5,7,9 +A 3 B 2 B 5,7,9 +A 5 B 2,4 B 7,9 }c 0 The first and second term in curly brackets of sequence c n in 420 is same as 48 Then its power series expansion is same as 49 Substitute 46 into the third term in curly brackets of sequence c n in 420 Its power series expansion including the first, second and third terms in curly brackets of sequence c n in 420, analytic at x=0, is y 0,2 2 x = c 0A 0 A β 2 i =0 β 2 + A 5 B 2,4 n=2 B 2i +3 B 2i +3 i =2 β 2 + A 3B 2 B 2i +3 n= i = x2n+2+λ 42 In 42 y 0,2 2 x is sub-power series, having sequences c n including two term of A n s in 33 as β 0 =0 andβ =2, for the polynomial case in which makes B n term terminated By repeating this process for allβ = 3, 4, 5,, we obtain every y 0, j 2 x terms where j 3 According to 47, 49, 42 and every y 0, j 2 x where j 3, the general expression of y0 2 x for the case ofβ 0= 0 replacing the index n by i 0 is β y 0 2 x=c i 0 β 2 0A 0 A i 2 2i 0 + B 2i +2 x2i 2+2+λ 422 i =0 i 2 =i 0 i 3 =i 0 B 2i3 +3 In 422 y 0 2 x is sub-power series, having sequences c n including two term of A n s in 33 as β 0 =0, for the polynomial case in which makes B n term terminated 2

13 b Asβ 0 =, then B 3 =0 in 3 c 2 = A 0, c 0 c 4 = A 0 A B 3 + A 3 B 2 ] + A 2 A3 B ]} c0 c 6 = A 0 A B 3,5 +A 3 B 2 B 5 +A 5 B 2,4 ] +A 2 A3 B 5 B + A 5 B 4 B ]} c0 c 8 = A 0 A B 3,5,7 + A 3 B 2 B 5,7 +A 5 B 2,4 B 7 + A 7 B 2,4,6 ] +A 2 A3 B 5,7 B + A 5 B 4 B 7 B + A 7 B 4,6 B ]} c0 423 c 0 = A 0 A B 3,5,7,9 +A 3 B 2 B 5,7,9 +A 5 B 2,4 B 7,9 +A 7 B 2,4,6 B 9 +A 9 B 2,4,6,8 ] +A 2 A3 B 5,7,9 B + A 5 B 4 B 7,9 B + A 7 B 4,6 B 9 B + A 9 B 4,6,8 B ]} c0 The first square brackets including A 0 inside curly brackets in sequence c n in 423 is same as 44 Then its power series expansion is same as 422 Observe the second square brackets including A 2 inside curly brackets in sequence c n in 423 i Asβ =, then B 4 =0 in the second square brackets including A 2 inside curly brackets in sequence c n in 423 c 4 = A 2 B A3 } c 0 c 6 = A 2 B A3 B 5 } c0 c 8 = A 2 B A3 B 5,7 } c0 424 c 0 = A 2 B A3 B 5,7,9 } c0 Its power series expansion of 424 by using 46, analytic at x=0, is β 2 y, 2 x=c 0A 2 B A 3 B 2i +3 x2n+2+λ 425 n= In 425 y, 2 x is sub-power series, for the second square brackets inside curly brackets in sequence c n including two term of A n s in 423 asβ 0 = andβ =, for the polynomial case in which makes B n term terminated ii Asβ =2, then B 6 =0 in second square brackets inside curly brackets in sequence c n in- i = 3

14 cluding A 2 in 423 c 4 = A 2 B A3 } c 0 c 6 = A 2 B A3 B 5 + A 5 B 4 } c 0 c 8 = A 2 B A3 B 5,7 + A 5 B 4 B 7 } c0 426 c 0 = A 2 B A3 B 5,7,9 + A 5 B 4 B 7,9 } c0 The first term in curly brackets of sequence c n in 426 is same as 424 Then its power series expansion is same as 425 Substitute 46 into the second term in curly brackets of sequence c n in 426 Its power series expansion including the first and second terms in curly brackets of sequence c n in 426, analytic at x=0, is y,2 2 x=c 0A 2 B A 3 β 2 n= B 2i +3 i = β 2 + A 5B 4 n=2 B 2i +3 i =2 x2n+2+λ 427 In 427 y,2 2 x is sub-power series, for the second square brackets inside curly brackets in sequence c n including two term of A n s in 423 asβ 0 = andβ =2, for the polynomial case in which makes B n term terminated By using similar process as I did before, the solution forβ =3 with B 8 =0 for the second square brackets inside curly brackets in sequence c n including A 2 in 423 is β 2 y,3 2 x = c 0A 2 B A β 2 3 B 2i +3 + A 5B 4 B 2i +3 n= i = n=2 i =2 β 2 + A 7 B 4,6 B 2i +3 x2n+2+λ 428 n=3 i =3 By repeating this process for allβ = 4, 5, 6,, we obtain every y, j 2 x terms where j 4 for the second square brackets inside curly brackets in sequence c n including two term of A n s in 423 According to 422, 425, 427, 428 and every y, j 2 x where j 4, the general expression of y 2 x for allβ 0= replacing the index n by i 0 is β y 2 x = c 0 A i 0 β 2 0 A i 2 2i 0 + B 2i +2 B 2i i =0 i 2 =i 0 i 3 =i 0 β i 0 β 2 + A 2 B A i 2 2i 0 + B 2i +2 x2i 2+2+λ i 0 = i = i 2 =i 0 i 3 =i 0 B 2i3 +3 4

15 c Asβ 0 =2, then B 5 =0 in 3 c 2 = A 0, c 0 c 4 = A 0 A B 3 + A 3 B 2 ] + A 2 B A3 ]} c 0 c 6 = A 0 A B 3,5 + A 3 B 2 B 5 + A 5 B 2,4 ] + A 2 B A3 B 5 + A 5 B 4 ] +A 4 B,3 A5 ]} c 0 c 8 = A 0 A B 3,5,7 + A 3 B 2 B 5,7 +A 5 B 2,4 B 7 + A 7 B 2,4,6 ] +A 2 B A3 B 5,7 + A 5 B 4 B 7 + A 7 B 4,6 ] + A 4 B,3 A5 B 7 + A 7 B 6 ]} c 0 c 0 = A 0 A B 3,5,7,9 + A 3 B 2 B 5,7,9 + A 5 B 2,4 B 7,9 + A 7 B 2,4,6 B 9 + A 9 B 2,4,6,8 ] +A 2 B A3 B 5,7,9 + A 5 B 4 B 7,9 + A 7 B 4,6 B 9 + A 9 B 4,6,8 ] 430 +A 4 B,3 A5 B 7,9 + A 7 B 6 B 9 + A 9 B 6,8 ]} c 0 By repeating similar process from the above, the general expression of y 2 2 x for allβ 0= 2 in 430 is β y 2 2 x = c 0 A i 0 β 2 0 A i 2 2i 0 + B 2i +2 B 2i3 +3 i =0 i 2 =i 0 i 3 =i 0 β i 0 β 2 +A 2 B A i 2 2i 0 + B 2i +2 B 2i i 0 = i = i 2 =i 0 i 3 =i 0 β i 0 β 2 + A 4 B,3 A i 2 2i 0 + B 2i +2 x2i 2+2+λ i 0 =2 i =2 i 2 =i 0 i 3 =i 0 B 2i3 +3 Again by repeating this process for allβ 0 = 3, 4, 5,, I obtain every y j 2 x terms where j 3 Then I have general expression y 2 x for allβ 0 of two term of A n s according 422, 429, 43 and y j 2 x terms where j 3 β 0 i 0 β y 2 x = c 0 A i 2 β 2 2i 0 B 2i + A i 4 2i 2 + B 2i3 +2 B 2i5 +3 x2i 4+2+λ i =0 i 2 =i 0 i 3 =i 0 i 4 =i 2 i 5 =i 2 whereβ 0 β β By using similar process for the previous cases of zero, one and two term of A n s, the function 5

16 y 3 x for the case of three term of A i s is β 0 i 0 y 3 x = c 0 A 2i 0 i =0 β 2 A 2i 4 +2 i 4 =i 2 β B 2i + i 2 =i 0 i 4 β 3 B 2i5 +3 i 5 =i 2 i 6 =i 4 i 2 A 2i 2 + B 2i i 3 =i 0 i 6 x2i 6+3+λ i 7 =i 4 B 2i7 +4 whereβ 0 β β 2 β 3 By repeating this process for all higher terms of A n s, I obtain every y m x terms where m>3 Substitute 45, 43, 432, 433 and including all y m x terms where m>3 into 35 Theorem 2 The general expression of yx for the polynomial case which makes B n term terminated in the three term recurrence relation is yx = y 0 x+y x+y 2 x+y 3 x+ β 0 i 0 β 0 i 0 β i 2 = c 0 B 2i + x2i 0+λ + A 2i 0 B 2i + B 2i3 +2 x2i 2++λ i =0 i =0 i 2 =i 0 i 3 =i 0 β 0 i 0 N β k i 2k + N=2 A 2i 0 B 2i + A 2i2k +k B 2i2k+ +k+ i =0 k= i 2k =i 2k i 2k+ =i 2k β N i 2N x2i 2N+N+λ 434 i 2N =i 2N i 2N+ =i 2N B 2i2N+ +N+ On the above,β i β j only if i j where i, j,β i,β j N 0 Theorem 3 For infinite series, replacingβ 0,β,β k andβ N by in 434 yx = y 0 x+y x+y 2 x+y 3 x+ i 0 = c 0 B 2i + x2i 0+λ i 0 i 2 + A 2i 0 B 2i + i =0 i =0 i 2 =i 0 i 0 N i 2k + N=2 A 2i 0 B 2i + A 2i2k +k i =0 k= i 2k =i 2k i 2N x2i 2N+N+λ i 2N =i 2N i 2N+ =i 2N B 2i2N+ +N+ i 3 =i 0 B 2i3 +2 i 2k+ =i 2k B 2i2k+ +k+ x2i 2++λ is exactly equivalent to is the another general expression of yx for the infinite series 5 Recurrence relation and Fibonacci numbers of higher order 5 Three term recurrence relation and Fibonacci sequence The Fibonacci sequence is:,, 2, 3, 5, 8, 3, 2, 34, 55, 89, 44, 5 6

17 The numbers in the sequence follows the recursive relation with seed values c n+ = c n + c n : n 52 c 0 = c = 53 The power series of the generating function of the Fibonacci sequence is c n x n x = x x 2 54 As we know, the three term recurrence relation is with seed values c n+ = A n c n + B n c n ; n 55 c = A 0 c is the expansion of 55 If coefficients A n = B n = in 55 and c 0 = c = in 56, then c n follows Fibonacci sequence Lucas series is another sequence generated by the three term recurrence relation with constant coefficients A n = B n = in 55 and c 0 = 2, c = in 56 You can think of these two sequences having constant coefficients as the most basic three term recurrence relation In contrast, Heun and Mathieu equations coefficients A n and B n are defined to be non-constant second order polynomial in denominator and numerator: see 3, 4a-4c in Ref3] and 4, 5a-5c in Ref5] The generating function of the Fibonacci sequence corresponds to infinite series of three term recurrence relation If all coefficients A n = B n = in 37 and 435, then it is equivalent to Two term recurrence relation There is an algebraic number sequence in which is I call 57 the identity sequence, and it s recurrence relation is with seed values,,,,,, 57 c n+ = c n : n 0 58 The power series of the generating function of the identity sequence is c 0 = 59 c n x n = x 50 If B n =0 in 55, then three term recurrence relation becomes a two term recurrence relation c n+ = A n c n ; n 0 5 7

18 Some examples are the Legendre function, Kummer function, hypergeometric function, Bessel function, etc The number of each sequence c n in 5 is,,,,,, is equivalent to 57 As I put A n = in 5, it becomes 58 You can think of this sequence having constant coefficients as the most basic two term recurrence relation The power series expansion of 5 for the infinite series is yx= c n x n+λ = c 0 A i xn+λ 53 And polynomial case of 5 is α 0 α 0 yx= c n x n+λ = c 0 A i xn+λ 54 The generating function of the identity sequence corresponds to infinite series of two term recurrence relation If all coefficients A n = in 53, then it definitely equivalent to Four term recurrence relation and Tribonacci sequence Now, let s think about four term recurrence relation in ordinary differential equation The four term recurrence formula is i=0 i=0 c n+ = A n c n + B n c n + C n c n 2 ; n 2 55 with seed values c = A 0 c 0 c 2 = A 0 A + B c 0 56 And the number of each of sequence c n in 55 is the following way: 57 is Tribonacci number, and it s recurrence relation is,, 2, 4, 7, 3, 24, 44, 57 c n+ = c n + c n + c n 2 ; n 2 58 with seed values c 0 = 0 c = c 2 = 59 If A n = B n = C n = in 55, it s exactly equivalent to Tribonacci recurrence relation Four term recurrence relation non-constant coefficients A n, B n and C n is the more general form than Tribonacci recurrence relation constant coefficients A n, B n and C n 8

19 54 Five term recurrence relation and Tetranacci sequence And five term recurrence relation in ordinary differential equation is with seed values c n+ = A n c n + B n c n + C n c n 2 + D n c n 3 ; n c = A 0 c 0, c 2 = A 0 A + B c 0, c 3 = A 0 A A 2 + A 0 B 2 + A 2 B + C 2 c 0 52 And the number of each of sequence c n in 520 is the following way: 522 is Tetranacci number, and it s recurrence relation is with seed values,, 2, 4, 8, 5, 29, 56, 08, 208, 522 c n+ = c n + c n + c n 2 + c n 3 ; n c 0 = 0 c = c 2 = c 3 = If A n = B n = C n = D n = in 520, it becomes exactly equivalent to Tetranacci recurrence relation From the above multi-term recurrence relation is more general form than n-nacci recurrence relation 6 Conclusion In this paper I show how to generalize three-term recurrence relation for polynomials and infinite series analytically In the next papers I will work out the analytic solution for the three term recurrence relation such as Heun equation and its confluent form, Lame equation, Mathieu equation I will derive the power series expansion, the integral form and the generating function of Heun, Mathieu, etc functions see section 7 for more details 7 Series Special functions and three-term recurrence formula 3TRF This paper is 2nd out of 0 Approximative solution of the spin free Hamiltonian involving only scalar potential for the q q system ] - In order to solve the spin-free Hamiltonian with light quark masses we are led to develop a totally new kind of special function theory in mathematics that generalize all existing theories of confluent hypergeometric types We call it the Grand Confluent Hypergeometric Function Our new solution produces previously unknown extra hidden quantum numbers relevant for description of supersymmetry and for generating new mass formulas 2 Generalization of the three-term recurrence formula and its applications 2] - Generalize three term recurrence relation in linear differential equation Obtain the exact solution of the three term recurrence relation for polynomials and infinite series 3 The analytic solution for the power series expansion of Heun function 3] - Apply three term recurrence formula to the power series expansion in closed forms of Heun function infinite series and polynomials including all higher terms of A n s 9

20 4 Asymptotic behavior of Heun function and its integral formalism, 4] - Apply three term recurrence formula, derive the integral formalism, and analyze the asymptotic behavior of Heun function including all higher terms of A n s 5 The power series expansion of Mathieu function and its integral formalism, 5] - Apply three term recurrence formula, analyze the power series expansion of Mathieu function and its integral forms 6 Lame equation in the algebraic form 6] - Applying three term recurrence formula, analyze the power series expansion of Lame function in the algebraic form and its integral forms 7 Power series and integral forms of Lame equation in Weierstrass s form and its asymptotic behaviors 7] - Applying three term recurrence formula, derive the power series expansion of Lame function in Weierstrass s form and its integral forms 8 The generating functions of Lame equation in Weierstrass s form 8] - Derive the generating functions of Lame function in Weierstrass s form including all higher terms of A n s Apply integral forms of Lame functions in Weierstrass s form 9 Analytic solution for grand confluent hypergeometric function 9] - Apply three term recurrence formula, and formulate the exact analytic solution of grand confluent hypergeometric function including all higher terms of A n s Replacingµandεω by and q, transforms the grand confluent hypergeometric function into Biconfluent Heun function 0 The integral formalism and the generating function of grand confluent hypergeometric function 20] - Apply three term recurrence formula, and construct an integral formalism and a generating function of grand confluent hypergeometric function including all higher terms of A n s References ] Frenkel, D and Protugal, R, Algebraic Methods to Compute Mathieu Functions, J Phys A: Math Gen 34, 200,354 2] Mathieu, E, Mémoire sur Le Mouvement Vibratoire d une Membrane de forme Elliptique, J Math Pure Appl 3, 868, p ] Whittaker, ET, On the General Solution of Mathieu s Equation, Proc Edinburgh Math Soc 32, 94, p ] Choun, YS, Approximative solution of the spin free Hamiltonian involving only scalar potential for the q q system, arxiv: ] M Hortacsu, Heun Functions and their uses in Physics, arxiv:0047 6] Gutirrez-Vega, J C, Rodrguez-Dagnino,R M, Meneses-Nava, M A and Chvez-Cerda, S, Mathieu functions, a visual approach, Amer J Phys 73, 2003, p ] Troesch, B A and Troesch, H R, Eigenfrequencies of an elliptic membrane, Math Comp 2724, 973, p ] Birkandan, T, Hortacsu, M, Examples of Heun and Mathieu functions as solutions of wave equations in curved spaces, J Phys A: Math Theor 40, 2007, 05 9] Birkandan, T, Hortacsu, M Dirac equation in the background of the Nutku helicoid metric, J Phys A: Math Theor 48, 2007, ] Aliev, AN, Hortacsu, M, Kalayci, J, and Nutku, Y, Gravitational instantons from minimal surfaces, Class Quantum Grav 6, 999, 63 ] Link to arxiv: ] Link to arxiv:

21 3] Link to arxiv: ] Link to arxiv: ] Link to arxiv: ] Link to arxiv: ] Link to arxiv: ] Link to arxiv: ] Link to arxiv: ] Link to arxiv: