Arxiv:1911.10319V2 [Math.CA] 29 Dec 2019 Roswl Egvni H Etscin.Atraiey Ecnuse Can We Alternatively, Sections

arXiv:1911.10319v2 [math.CA] 29 Dec 2019 roswl egvni h etscin.Atraiey ecnuse can we Alternatively, sections. next the in given be will proofs a rl xrse n[3 et 2]. ce the Sect. accepted of [13, tacitly moment in was second-order expressed function the firmly non-elementary that was a opinion is The operators ” Zeller p.2]: [14, in said Here rpitsbitdt ... to submitted Preprint ucin nsre rfiiesm feeetr yegoercfun hypergeometric elementary the of Keywords: with sums concerned finite are or series applications in exp Other are functions t operators polylogarithms. these Zeller and of Meyer-König moments and tions order of a higher moments they The that order functions. prove second and the functions sequently, hypergeometric some consider We Abstract prtr a eepesdas expressed be hypergeometr can elementary operators of an sums operators finite type or series Zeller Meyer-König in and of functions Heun moments the to tions Introduction 1. 00MSC: 2010 functions Heun polylogarithms; operators, lmnayhpremti ucin,Hu ucin,and functions, Heun functions, hypergeometric Elementary 1 a hsi upiig eas ti lmnayt rv that prove to elementary is it because surprising, is This lsdfr of form closed A n[] ...Akmd rvdta h eododrmmn fteM the of moment order second the that proved Alkemade J.A.H. [4], In hypergeometric elementary of families some to devoted is paper This orsodn author corresponding uinBaaUiest fSbu eateto Mathematic of Department Sibiu, of University Blaga Lucian e ahmtc,Sr eoadmlin.2 ljNpc,Roma Cluj-Napoca, 28 nr. Memorandumului Str. Mathematics, b r ehia nvriyo ljNpc,Fclyo Automatio of Faculty Cluj-Napoca, of University Technical ( t := ) yegoercfntos lmnayfntos Meyer-König functions, elementary functions, hypergeometric 30,3C0 33,41A36 33E30, 33C90, 33C05, t r , t ∈ [0 M M , n n 1], O501 ii,Rmna -al [email protected] e-mail: Romania, Sibiu, RO-550012 ( ( e e 2 2 r ; ; x x ≥ 2 = ) ,soigta ti neeetr ucin a ie n[4.I is It [14]. in given was function, elementary an is it that showing ), F ,and 0, 1 (1 x n-ai Acu Ana-Maria , 2 ;3; 2; fMZoperators MKZ of + 2 F x x 1 (1 n ( = ) ,b a, − 1 + x − ; ) c 2 2 ; x x 2 1a − F eoe h yegoercfunction. hypergeometric the denotes ) 2 1 onRasa Ioan , ( (1 x , log(1 + 2; n nomtc,Sr r .Rtu No.5-7, Ratiu, I. Dr. Str. Informatics, and s n optrSine eatetof Department Science, Computer and n n 2; + i,emi:[email protected] e-mail: nia, 2 F − b 1 esdi em feeetr func- elementary of terms in ressed x (1 x ) x , )) , ctions. ic 90”Tesm opinion same The 1960.” since 2; cfunctions. ic eeeetr ucin.Con- functions. elementary re p prtr r elementary are operators ype , oteepnino certain of expansion the to d ∈ n 2; + xaso fcranHeun certain of expansion [0 ertdMeyer-König and lebrated , ucin,wt applica- with functions, 1) n elrtype Zeller and x yrKoi n Zeller eyer-König and . seeetr:such elementary: is ) aur ,2020 1, January moments (1.2) (1.1) and [12, 2.8(24)]:

m+c−a−b (c)m(1 − x) a+b−c (m) 2F1(a,b; c + m; x)= (1 − x) 2F1(a,b; c; x) , (1.3) (c − a)m(c − b)m where (r)m := r(r + 1) ... (r + m − 1), m ≥ 1, and (r)0 := 1. They lead to

n−1 (3)n−1(1 − x) (n−1) 2F1(1, 2; n + 2; x)= (2F1(1, 2; 3; x)) (2)n−1(1)n−1 n +1 (n−1) = − (1 − x)n−1 x−1 + x−2 log(1 − x) , (1.4) (n − 1)! which is obviously an elementary function. Thus Mn(e2; x) given by (1.1) is elementary. On the other hand, [8, (32)] provides a closed form of the (elementary) hypergeometric function 2F1(1,m; m +2k + 1; x), where m ≥ 1 and k ≥ 0 are integers. In Section 2 we consider the more general function 2F1(m,n; p; x) with m, p positive integers, p ≥ m + 1, and n ∈ R. We present a closed form of this function, showing that it is elementary. A general formula for 2F1(a,b; a + b + l; x) can be found in [12, 2.3(2)], but for arbitrary parameters it is not always an elementary function. Generalizing [8, (32)], we give in Section 3 two closed forms for the elementary function 2F1(1,m; m + l + 1; x) with m ≥ 1, l ≥ 0 integers. Three closed forms for 2F1(1, 2; n + 2; x) are presented in Section 4. Using a result of [2] we give in Section 5 another proof that the second MKZ moment is elementary. The same result of [2] is instrumental in Section 6, where we show that the higher order MKZ moments can be expressed as finite sums of functions, but containing polylogarithms. In the final section we consider two modifications of the MKZ operators and show that their second order moments are elementary functions; here we use Theorem 2.1 and the fact that 2F1(1, 3; n + 3; x) is elementary. Let us mention that the moments of the MKZ operators are involved in studying the asymptotic behavior of the iterates of these operators. For a qualitative result, see [13]. Generalized qualitative and quantitative results can be found in [23, Sect. 2], [7] and [6, Sect.1.3]. Section 8 is devoted to some applications involving Heun equations. The power-series solutions of such an equation are governed by three-term recurrence relations between the successive coefficients, instead of two-term recurrence relations which appear in the hypergeometric case. There are few Heun equations for which the exact solutions are known in terms of simpler functions. This is why expansions of the solutions of Heun equations using other functions instead of powers were considered by many authors; see [16], [17], [19], [24], and the references therein. In particular, expansions in terms of hypergeometric functions and conditions for deriving finite-sum solutions can be found in [16]. Using results from [16], we describe a family of Heun equations for which the solutions are expanded in series of elementary hypergeometric functions. Moreover, we present conditions under which such a series terminates and so the solutions are elementary functions. Other results and applications will be given in a forthcoming paper.

2 2. The function 2F1(m, n; p; x)

Theorem 2.1. Let m, p be positive integers, p ≥ m +1, and n ∈ R. Then

m−1 (m) x1−p m − 1 F (m,n; p; x)= p−m (−1)k 2 1 (p − m − 1)! k k X=0 p−m−1 j p − m − 1 j 1 · (−1)jxp−m−1−j (−1)i si+k−nds. (2.1) j i j=0 i=0 1−x X X Z

Consequently, 2F1(m,n; p; x) is an elementary function. Proof. According to the deﬁnition,

∞ (m)j (n)j j 2F1(m,n; p; x)= x . (p)j j! j=0 X This implies

∞ p−1 (p−m) (m)j n + j − 1 j+p−1 (p−m) x 2F1(m,n; p; x) = x (p)j j j=0 X ∞ n + j − 1 xm−1 = (m) xm−1 xj = (m) . (2.2) p−m j p−m (1 − x)n j=0 X p−1 For the function x 2F1(m,n; p; x), the Taylor polynomial of degree p − 2 ≥ p − 1 − m at 0 vanishes, and so

t 1 1 (p−m) F (m,n; p; x)= (x − t)p−m−1 tp−1 F (m,n; p; t) dt 2 1 xp−1 (p − m − 1)! 2 1 Z0 1 1 t = (x − t)p−m−1(m) tm−1(1 − t)−nd. xp−1 (p − m − 1)! p−m Z0 Therefore,

m−1 1 (m) m − 1 x F (m,n; p; x)= p−m (−1)k (x − t)p−m−1(1 − t)k−ndt. 2 1 xp−1 (p − m − 1)! k k 0 X=0 Z Now use the binomial formula for (x − t)p−m−1 and set s := 1 − t in order to verify that the last integral equals

p−m−1 j p − m − 1 j 1 (−1)j xp−m−1−j (−1)i si+k−nds. j i j=0 i=0 1−x X X Z This concludes the proof of (2.1). Remark 2.1. The proof of Theorem 2.1 is truly elementary: it uses only the deﬁnition of the hypergeometric function and Taylor’s formula. Another proof has been recently given in [3].

3 For m = 1, Theorem 2.1 yields Corollary 2.1. Let p ≥ 2 be an integer, and n ∈ R. Then p−2 p − 1 p − 2 1 F (1,n; p; x)= (x − 1)p−2−i si−nds. (2.3) 2 1 xp−1 i i=0 1−x X Z Proof. Setting m = 1 in (2.1), we get p−2 j p − 1 p − 2 j 1 F (1,n; p; x)= (−1)j xp−2−j (−1)i si−nds 2 1 xp−1 j i j=0 i=0 1−x X X Z p−2 p−2 p − 1 p − 2 p − 2 − i 1 = (−1)i (−1)jxp−2−j si−nds xp−1 i j − i i=0 j=i 1−x X X Z p−2 p − 1 p − 2 1 = (x − 1)p−2−i si−nds. xp−1 i i=0 1−x X Z

Remark 2.2. Let λ, µ, and ν be nonnegative integers. Theorem 2.1 may be considered also in analogy to the following result (see [12, p. 69], [11, (23)]):

2F1(ν +1,ν + µ + 1; ν + µ + λ + 2; x) (−1)µ+1(ν + µ + λ + 1)! dν+µ dλ log(1 − x) = (1 − x)µ+λ . (2.4) λ!ν!(ν + µ)!(µ + λ)! dxν+µ dxλ x Setting m = ν +1, n = ν + µ +1, p = ν + µ + λ +2, (2.4) leads to

2F1(m,n; p; x) p − 1 p − 2 p − 2 dn−1 dp−n−1 log(1 − x) = (−1)n−m+1 (1 − x)p−m−1 , (2.5) (p − 2)! n − 1 m − 1 dxn−1 dxp−n−1 x for all integers m,n,p with p ≥ n +1 ≥ m +1 ≥ 2.

3. The function 2F1(1,m; m + l + 1; x) Theorem 3.1. Let m ≥ 1 and l ≥ 0 be integers. Then (m) (−1)l+1 F (1,m; m + l + 1; x)= l+1 (1 − x)l log(1 − x) 2 1 xm+l l! m+l−1 m + l m + l − 1 (−1)m+l−1−i + (1 − x)m+l−1−i − (1 − x)l , xm+l i i − m +1 i ,i6 m− =0X= 1 (3.1) and also (m) (−1)l+1 F (1,m; m + l + 1; x)= l+1 (1 − x)l log(1 − x) 2 1 xm+l l! l+1 l j−1 i m−2 l+i+1 (m)l+1 (−1) j (−1) l x + m+l x − . (3.2) x  l! j − i i (i + 1)l  j=1 i=0 i=0 +1  X X X    4 Proof. (3.1) is a consequence of (2.3). To prove (3.2), let us replace in (2.2) m by 1, n by m, an p by m + l + 1; thus we have

m−1 (l+1) x 1 xn+l F (1,m; m + l + 1; x) = (m) = (m) − 1 − x −···− xm−2 . 2 1 l+1 1 − x l+1 1 − x It follows that

l m−2 l+1+i m+l l+1 (1 − x) x x 2F1(1,m; m + l + 1; x) = (m)l+1 (−1) log(1 − x) − − Ql(x) , l! (i + 1)l ( i=0 +1 ) X where Ql is a polynomial of degree less or equal to l. m+l The right hand side should be divisible by x , hence Ql(x) is the Taylor polynomial of degree (1 − x)l l, at 0, for the function (−1)l+1 log(1 − x). It is not diﬃcult to obtain l!

l j−1 (−1)l (−1)i l Q (x)= xj . l l! j − i i j=1 i=0 X X Putting all things together we get (3.2). Remark 3.1. Equations (3.1) and (3.2) extend Eq. (32) in [8].

4. The function 2F1(1, 2; n + 2; x)

For this function, which is involved in the representation (1.1), we have three explicit expressions: Theorem 4.1. Let n ≥ 1 be an integer.Then

n− n +1 1 n − j F (1, 2; n + 2; x) = (−1)n (1 − x)n−1 [x + n log(1 − x)] − (−1)j xj (1 − x)n−j−1 , 2 1 xn+1  j  j=1  X  (4.1)   n +1 F (1, 2; n + 2; x) = (−1)n (1 − x)n−1 [x + n log(1 − x)] 2 1 xn+1 n (−1)i n + (1 − x)n−i − (1 − x)n−1 , (4.2) i − 1 i i=2 ) X n−1 j−1 n(n + 1) (−1)i n − 1 n +1 F (1, 2; n + 2; x) = (−1)n (1 − x)n−1 log(1 − x)+ xj − . 2 1 xn+1  j − i i  x j=1 i=0  X X  (4.3)   Proof. (4.1) can be obtained if we continue the calculation in (1.4). (4.2) is a consequence of (3.1), and (4.3) is obtained from (3.2).

Remark 4.1. Other expressions for 2F1(1,m; m + l + 1; x) and 2F1(1, 2; n + 2; x) can be obtained from (2.5).

5 5. The second MKZ moment is elementary: another proof

From [2, (6) and (7)] we know that

r r (−n)j ∞ tj−1e−nt M e (x)=1+(1 − x)n+1 dt. n r j (j − 1)! (1 − xe−t)n+1 j=1 0 X Z On the other hand,

∞ ∞ ∞ −n − 1 tj−1e−nt(1 − xe−t)−n−1dt = tj−1e−nt (−1)kxke−ktdt k Z0 Z0 k=0 ∞ X ∞ n + k ∞ n + k Γ(j) = xk tj−1e−(n+k)tdt = xk , k k (n + k)j k 0 k X=0 Z X=0 so that r r M e (x)=1+(1 − x)n+1 (−n)jf (x), (5.1) n r j n,j j=1 X where ∞ n + k 1 f (x) := xk, n ≥ 1, j ≥ 0. (5.2) n,j k (n + k)j k X=0 Let us remark that 1 f (0) = . (5.3) n,j nj Moreover, ∞ n + k 1 (xnf (x))′ = xn−1+k = xn−1f , n,j k (n + k)j−1 n,j−1 k X=0 and so 1 x f (x)= tn−1f (t)dt, j ≥ 1. (5.4) n,j xn n,j−1 Z0 1 We have f (x)= , and (5.4) yields n,0 (1 − x)n+1 1 f (x)= . (5.5) n,1 n(1 − x)n

From (5.4) and (5.5) it is easy to obtain

n−1 (−1)n−1 n − 1 (−1)i 1 − (1 − x)i f (x)= − log(1 − x) . (5.6) n,2 nxn i i (1 − x)i ( i=1 ) X

Now (5.1), (5.5) and (5.6) show that Mne2(x) is an elementary function.

6 6. Higher order moments of the MKZ operators

We recall the deﬁnition (see [20]) of the dilogarithm

∞ x log(1 − t) xk Li (x) := − dt = , |x|≤ 1, (6.1) 2 t k2 0 k Z X=1 and the polylogarithm of order n,

∞ x Li (t) xk Li (x) := n−1 dt = , n ≥ 3. (6.2) n t kn 0 k Z X=1 Using (5.4) and (5.6) we get

n−1 (−1)n−1 x n − 1 (−1)i 1 − (1 − t)i log(1 − t) f (x)= − dt n,3 nxn i i t(1 − t)i t 0 ( i=1 ) Z X n−1 i−2 (−1)n−1 x n − 1 (−1)i 1 log(1 − t) = (1 − t)l−i + − dt, nxn i i 1 − t t 0 ( i=1 l ! ) Z X X=0 so that ﬁnally,

n−1 i−2 (−1)n−1 n − 1 (−1)i 1 − (1 − x)i−l−1 f (x)= −log(1−x) + Li (x) . (6.3) n,3 nxn i i (i−l−1)(1−x)i−l−1 2 ( i=1 " l # ) X X=0 n Theorem 6.1. Let n ≥ 2. For each j ≥ 2, the functions x fn,j(x) are linear combinations, with constant coeﬃcients, of n functions as follows: 1) For 2 ≤ j ≤ n:

1 − (1 − x)i , i =1,...,n − j + 1; log(1 − x); Li (x),...,Li (x). (1 − x)i 2 j−1

2) For j = n +1: log(1 − x); Li2(x),...,Lin(x).

3) For j ≥ n +2: Lij−n(x),...,Lij−1(x).

Proof. We use induction on j. The ﬁrst step, involving j = 2, j = 3, is accomplished by using (5.6) and (6.3). Then the induction is continued with the recurrence relation (5.4) and (6.1), (6.2). Remark 6.1. Using (5.1) and Theorem 6.1 we get information about the structure of the moments Mner. In particular, we see once more that Mne2 is elementary. The higher order moments contain polylogarithms.

7 7. Other operators

n+k n+1 k Let mn,k(x) := k (1 − x) x , x ∈ [0, 1], and for f ∈ L1[0, 1], ∞ (n + k + 1)(n + k + 2) 1 L (f; x) := m (t)f(t)dt m (x). n n +1 n,k n,k k 0 X=0 Z The operators Ln were investigated in [10]; they coincide with the operators Mn,2 from [15]; details can be found in [9]. It was proved in [10] (see also [9, p. 123]) that 2x(1 − x)2 L (e ; x)= x2 + F (1, 3; n + 3; x). (7.1) n 2 n +2 2 1

Closed forms of the elementary functions 2F1(1, 3; n+3; x) can be obtained from (3.1) and (3.2) taking m = 3, l = n − 1. In particular, we see that Ln(e2; x) is an elementary function. It would be interesting to study from this point of view the higher order moments of Ln and, more generally, of other modifications of the MKZ operators presented in [9]. Remark 7.1. In [5, Chapter 7] the generalized Meyer-König and Zeller operators are defined as follows: ∞ n + r + k − 1 k + β M α,βf (x) := (1 − x)n+r xkf , n,r k n + k + α k=0 X x ∈ [0, 1), α,β ∈ R, α ≥ β ≥ 0, r ∈ Z,n,n + r ∈ N, and it is proved that β n + r − β M α,βe (x)= F (1, α − r; n +1+ α; x)+ x F (1, α − r +1,n +2+ α; x). n,r 1 n + α 2 1 n +1+ α 2 1 (7.2)

α,β If α ≥ 0 is an integer, we can apply Theorem 2.1 to conclude that Mn,r e1 is an elementary function. Moreover, suppose that α ≥ 0 is an integer and r = α +1. Then (7.2) yields β β M α,β e (x)= + 1 − x. (7.3) n,α+1 1 n + α n + α j Let gj(x) := (1 − x) , j = 0, 1,.... Using a trick invented by Professor Ulrich Abel (see [1]), we can write ∞ (n + α − β)j (n + k + α)! xk M α,β g (x)= (1 − x)n+α+1 n,α+1 j (n + α)! k! (n + k + α)j k X=0 ∞ (n+α) (n + α − β)j xn+k+α = (1 − x)n+α+1 (n + α)! (n + k + α)j k ! X=0 (n + α − β)j = (1 − x)n+α+1 (Li (x))(n+α) . (n + α)! j

8 m m Combined with e = (−1)j g , this leads to m j j j=0 X m (1 − x)n+α+1 m M α,β e (x)= (−1)j (n + α − β)j (Li (x))(n+α) . n,α+1 m (n + α)! j j j=0 X α,β Now it is easy to recover (7.3) and to see that Mn,α+1e2 is an elementary function.

8. Elementary Heun functions

The general Heun equation in canonical form is (see [18], [21], [22]): d2u γ δ ε du αβx − q + + + + u =0, (8.1) dx2 x x − 1 x − a dx x(x − 1)(x − a) where the parameters satisfy α + β +1= γ + δ + ε. Under the hypotheses that αβ 6= 0 and γ + ε is not zero or a negative integer, an expansion of the solution u(x) as a series of hypergeometric functions is investigated in [16]. Let m and p be positive integers, p ≥ m + 1, and n ∈ R, n 6= 0. In (8.1) choose

α = m, β = n, γ = p +1 − m − n, δ = m + n − p +1, 1 1 ε = m + n − 1, a = , q = (mn − (m + n − p)(m + n − 1)) . 2 2 Theorem 8.1. With the above parameters, the solution u(x) of (8.1) is given by m + n − 1 p − m p − n ∞ 2 2 2 u(x)= k k k F (m,n; p +2k; x), (8.2) p p +1 2 1 k X=0 k! 2 k 2 k and 2F1(m,n; p +2k; x) are elementary functions. Moreover, if there exists a positive integer r such that m + n = 3 − 2r or n − p = 2r − 2, then the series (8.2) reduces to a ﬁnite sum and so u(x) is an elementary function. 1 Proof. It is easy to check that αβ 6= 0, γ + ε = p ≥ 2, a = , γ + δ = 2, q = aαβ + a(1 − δ)ε, and 2 so the conditions [16, (43)] are satisﬁed. According to [16, (45)], the solution of the Heun equation (8.1) is given by

∞ ε γ+ε−α γ+ε−β 2 k 2 k 2 u(x)= k F (α, β; γ + ε +2k; x), γ+ε 1+γ+ε 2 1 k k! 2 k 2 k X=0 which leads immediately to (8.2). Using Theorem 2.1 we see that 2F1(m,n; p +2k; x) is an elementary function for each k ≥ 0. Let m + n =3 − 2r, or n − p =2r − 2, with r an integer, r ≥ 1. Then m + n − 1 p − n =0, or =0, 2 2 r r

9 and so the series in (8.2) reduces to the sum for k ∈{0, 1,...,r − 1}; this means that in this case u(x) is an elementary function for which explicit expressions can be obtained using results from Sections 2-4. Remark 8.1. Usually, a Heun function is normalized by u(0) = 1. For the equation (8.1), [16, (46)] yields γ + ε − α γ + ε − β ε γ + ε 1+ γ + ε u(0) = F , , ; , ;1 , 3 2 2 2 2 2 2 where 3F2 is Clausen’s generalized hypergeometric function. With the parameters used in Theorem 8.1 we get p − m p − n m + n − 1 p 1+ p u(0) = F , , ; , ;1 . 3 2 2 2 2 2 2 Thus u(x)/u(0) satisﬁes the normalization condition.

Acknowledgement. The work of the ﬁrst author was ﬁnanced by Lucian Blaga University of Sibiu & Hasso Plattner Foundation research grants LBUS-IRG-2019-05.

References

[1] U. Abel, Elementary function representations for the moments of the Meyer-König and Zeller operators, arXiv: 1911.02563v1 [math.GM] [2] U. Abel, The moments for the MeyerKönig and Zeller operators, J. Approx. Theory 82(3), 352361 (1995). [3] U. Abel, Unpublished manuscript, 2019. [4] J.A.H. Alkemade, The second moment for the MeyerKönig and Zeller operators, J. Approx. Theory 40(3), 261273 (1984). [5] J.A.H. Alkemade, On a differential equation and the Laplace- Stieltjes transform in the theory of linear positive operators, Thesis, Delft University Press, 1984. [6] F. Altomare, M. Cappelletti Montano, V. Leonessa, I. Raa, Markov Operators, Positive Semi- groups and Approximation Processes, 61, De Gruyter Studies in Mathematics, 2014. [7] F. Altomare, I. Rasa, Lipschitz contractions, unique ergodicity and asymptotics of Markov semigroups, Bollettino U. M. I. (9), 5(2012), 1-17. [8] A. B˘arar, G. Mocanu, I. Ra¸sa, Heun functions related to entropies, RACSAM (2019), 113: 819-830. [9] J. Bustamante, M.A. Jiménez Pozo, Meyer-König and Zeller operators and some of their modifications, Jaen J. Approx. 5(2), 2013, 101-178. [10] W. Chen, The second moment for the Meyer-König and Zeller-type operators, J. Xiamen Univ. Natur. Sci. 29(2), 1990, 119-123 (in Chinese). [11] J. Detrich, R.W. Conn, Finite sum evaluation of the Gauss hypergeometric function in an important special case, Math. Comput. 33, 788791 (1979).

10 [12] A. Erdelyi, Higher Transcendental Functions, Vol. I, McGraw-Hill Book Company, New York, 1953. [13] I. Gavrea, M. Ivan, On the iterates of positive linear operators preserving the affine functions, J. Math. Anal. Appl. 372 (2010) 366-368. [14] I. Gavrea, M. Ivan, An elementary function representation of the second-order moment of the Meyer-König and Zeller operators, Mediterr. J. Math. (2018), 15:20. [15] H. Heilmann, Commutativity of Durrmeyer-type modifications of Meyer-König and Zeller and Baskakov operators, in: B. Bojanov (ed.), Constructive Theory of Functions, Proc. Internat. Conf., Varna, Bulgaria, 2002, pp. 295-301, DARBA, Sofia. [16] T.A. Ishkhanyan, T.A. Shahverdyan, A.M. Ishkhanyan, Expansions of the solutions of the general Heun equation governed by two-term recurrence relations for coefficients, Advances in High Energy Physics 2018, Article ID 4263678, (2018). [17] A. Ishkhanyan, K.A. Suominen, New solutions of Heuns general equation, J. Phys. A Math. Gen. 36, L81L85 (2003). [18] G. Kristensson, Second Order Differential Equations. Special Functions and Their Classifica- tion, Springer, 2010. [19] C. Leroy, A.M. Ishkhanyan, Expansions of the solutions of the confluent Heun equation in terms of the incomplete Beta and the Appell generalized hypergeometric functions, Integral Transform. Spec. Funct. 26(6), 451459 (2015). [20] L. Lewin, Polylogarithms and associated functions, North-Holland, Amsterdam, 1981. [21] R.S. Maier, The 192 solutions of the Heun equation, Math. Comp. 76, 811-843 (2007). [22] R.S. Maier, On reducing the Heun equation to the hypergeometric equation, J. Differ. Equ. 213, 171203 (2005).

[23] I. Rasa, C0-Semigroups and iterates of positive linear operators: asymptotic behaviour, Ren- diconti del Circolo Matematico di Palermo, Ser. II, Suppl. 82(2010),123-142. [24] V.A. Shahnazaryan, T.A. Ishkhanyan, T.A. Shahverdyan, A.M. Ishkhanyan, New relations for the derivative of the conﬂuent Heun function, Armen. J. Phys. 5, 146156 (2012).