DIFFERENTIAL EQUATION AND RECURRENCE RELATIONS OF THE SHEFFER-APPELL POLYNOMIAL SEQUENCE: A MATRIX APPROACH H. M. SRIVASTAVA, SAIMA JABEE AND MOHAMMAD SHADAB∗ Abstract. Motivated by the effective impact of the Pascal functional and the Wron- skian matrices, we investigate several identities and differential equation for the Sheffer- Appell polynomial sequence by using matrix algebra. The matrix approach, which we have used in this article, is convenient to derive the generating functions of the Sheffer- Appell polynomial sequence. By means of examples, we apply and also illustrate our results to an extended class of polynomial sequences. 1. introduction Sequences of polynomials play an important rˆole in many problems of pure and ap- plied mathematical sciences such as those occurring in approximation theory, statistics, combinatorics and analysis (see, for example, [10, 11, 12, 13]). The class of Sheffer se- quences is one of the most important classes of polynomial sequences. A polynomial ∞ sequence {sn(x)}n=0 is called a Sheffer polynomial sequence [3, 4, 12, 15] if and only if its generating function has the following form: ∞ n xH(y) y A(y)e = sn(x) , (1.1) n n! X=0 where A(y)= A0 + A1y + · · · , and arXiv:1903.09620v1 [math.CA] 21 Mar 2019 2 H(y)= H1y + H2y + · · · , with A0 6= 0 and H1 6= 0. Let us recall an alternate definition of the Sheffer sequences in terms of a pair of gen- erating functions (l(y), h(y)) (see, for example, [8]): 2010 Mathematics Subject Classification. Primary 15A15, 15A24, 33C45; Secondary 65Q30. Key words and phrases. Sheffer-Appell polynomial sequences; Pascal functional; Wronskian matrices; Differential equation; Recurrence relations, Generalized Laguerre polynomials; Miller-Lee type Appell polynomials; Orthogonal polynomialds. *Corresponding author. 1 2 H.M.SRIVASTAVA,SAIMAJABEEANDMOHAMMADSHADAB Let h(y) be a delta series and let l(y) be an invertible series, defined as follows: ∞ yn h(y)= hn (h(0) = 0; h(1) 6= 0) (1.2) n n! X=0 and ∞ yn l(y)= ln (l(0) 6= 0). (1.3) n n! X=0 Then there exists a unique sequence of Sheffer polynomials sn(x) satisfying the orthog- onality conditions: k hl(y)h(y) |sn(x)i = n!δn,k (∀ n, k ≧ 0), (1.4) where δn,k is the Kronecker delta. Roman [12, p. 18, Theorem 2.3.4] introduced the exponential generating function of sn(x) as follows: ∞ n 1 −1 y exh (y) s x . −1 = n( ) (1.5) h (y) n n! X=0 The Sheffer sequence for the pair (l(y),y) is called an Appell sequence for l(y). In fact, Roman [12] characterized Appell sequences in several ways: {An(x)}n∈N is an Appell set if either d α (x) = nα (x) n ∈ N dx n n or if there exists an exponential generating function of the form (see also the recent works [9, 16]): ∞ n xy y A(y)e = αn(x) , (1.6) n n! X=0 where N denotes the set of positive integers and 1 A(y)= . l(y) We also note that, for H(y)= y, the generating function (1.1) of the Sheffer polynomials sn(x) reduces to the generating function (1.6) of the Appell polynomials αn(x). The polynomials defined as the discrete convolution of known polynomials are used to (A) investigate new families of special functions. For example, the polynomial hn (x) given by ∞ (A) hn (x)= Akhn−k(x), (1.7) r X=0 n is known as a discrete Appell convolution by setting hn(x)= x in the above equation. SHEFFER-APPELL POLYNOMIAL SEQUENCE: A MATRIX APPROACH 3 In the year 2015, Subuhi et al. [7] introduced the determinantal definition and other properties of the Sheffer-Appell polynomials. The Sheffer-Appell polynomial sequences are combination of the families of the Sheffer and the Appell polynomials sequences. Now, in order to recall the definition of the generalized Pascal functional matrix of an analytic function (see [18]), let ∞ yr F = h(y)= αr αr ∈ C ( r! ) r=0 X C be the set of power series possessing the -algebra. Then the generalized Pascal func- tional matrix [Pn(h(y))], which is a lower triangular matrix of order (n + 1) × (n + 1) for h(y) ∈ F, is defined by i h(i−j)(y) (i ≧ j) j! Pn[h(y)]ij = (1.8) 0 (otherwise), for all i, j = 0, 1, 2, · · · ,n. Hereh(i)(y) is the ith order derivative of h(y). We next recall the nth order Wronskian matrix of several analytic functions h1(y), h2(y), · · · , hm(y) of order (n + 1) × m as follows: h1(y) h2(y) h3(y) · · · hm(y) ′ ′ ′ ′ h1(y) h2(y) h3(y) · · · hm(y) Wn[h1(y), h2(y), · · · , hm(y)] = . (1.9) . .. (n) (n) (n) (n) h1 (y) h2 (y) h3 (y) · · · hm (y) We also record here some properties and relationships between the Wronskian matrices and the generalized Pascal functional matrices as they are the main tool of our work (see, for example, [19, 20]). Property I. For h(y), l(y) ∈ C, Pn[h(y)] and Wn[h(y)] are linear, that is, Pn[uh(y)+ vl(y)] = uPn[h(y)] + vPn[l(y)] and Wn[uh(y)+ vl(y)] = uWn[h(y)] + vWn[l(y)], (1.10) where u, v ∈ C. Property II. For h(y), l(y) ∈ C, Pn[h(y)l(y)] = Pn[h(y)]Pn[l(y)] = Pn[l(y)]Pn[h(y)]. (1.11) Property III. For h(y), l(y) ∈ C, Wn[h(y)l(y)] = Pn[h(y)]Wn[l(y)] = Pn[l(y)]Wn[h(y)]. (1.12) 4 H.M.SRIVASTAVA,SAIMAJABEEANDMOHAMMADSHADAB Property IV. For h(y), l(y) ∈ C, with h(0) = 0 and h′(0) 6= 0, W l h y W , h y , h2 y , , hn y −1W l y , n[ ( ( ))] y=0 = n 1 ( ) ( ) · · · ( ) y=0 Ωn n[ ( )] y=0 (1.13) h i where 0!, 1!, · · · ,n! are the diagonal entries in the diagonal matrix given by Ωn = diag[0!, 1!, 2!, · · · ,n!]. 2. The Sheffer-Appell polynomial sequence and its differential equation He and Ricci ([6]; see also [14]) derived some recurrence relations and differential equation for the Appell polynomial sequence. Further, Youn and Yang ([20]; see also [1]) obtained some identities and differential equation for the Sheffer polynomial sequence by using matrix algebra. Here, in this paper, we study some recursive formulas and differential equation for the Sheffer-Appell polynomial sequence by using matrix algebra. The Sheffer-Appell polynomial sequence, which is denoted by sAn(x), is defined as the discrete Appell convolution of the Sheffer polynomials sn(x). The generating function of the Sheffer-Appell polynomials sAn(x) is given by ∞ n 1 −1 y exh (y) A x , −1 = s n( ) (2.1) l(h (y))l(y) n n! X=0 where h−1(y) is the compositional inverse of h(y), that is, we have (see [19]) h−1(h(y)) = (h(h−1(y)) = y. Thus, if the following generating function in (2.1): 1 −1 exh (y) l(h−1(y))h(y) is analytic, then (by using Taylor’s expansion theorem), we obtain k d 1 xh−1(y) sAk(x)= − e (k ≧ 0). (2.2) dyk l(h 1(y))l(y) y=0 The Sheffer-Appell polynomial sequence sAn(x) in vector form for the pair (l(y), h(y)) is denoted by sA~ n(x) and it is defined by ~ T sAn(x) = [sA0(x), sA1(x), · · · , sAn(x)] , (2.3) which can also be expressed as follows: ~ T 1 xh−1(y) sAn(x) = [sA0(x), sA1(x), · · · , sAn(x)] = Wn − e . (2.4) l(h 1(y))l(y) y=0 SHEFFER-APPELL POLYNOMIAL SEQUENCE: A MATRIX APPROACH 5 Lemma. Let sAn(x) be the Sheffer-Appell polynomial sequence for the pair (l(y), h(y)). Then T −1 Wn[sA0(x), sA1(x), · · · , sAn(x)] Ωn W , h−1 y , h−1 y 2, , h−1 y n = n[1 ( ( )) ( ( )) · · · ( ( )) ] y=0 −1 1 1 xy · Ωn Pn Pn Pn [e ] . (2.5) l(y) l(h(y)) y=0 y=0 y=0 Proof. Let us begin with the equation (2.4), that is, ~ 1 xh−1(y) sAn(x)= Wn − e . (2.6) l(h 1(y))l(y) y=0 Applying Property IV in the equation (2.6), we get ~ −1 −1 2 −1 n −1 1 xy sAn(x)= Wn[1, (h (y)), (h (y)) , · · · , (h (y)) ] Ωn Wn e .(2.7) l(y)l(h(y)) y=0 y=0 In view of the following result: xy 2 n T Wn[e ] = [1, x, x , · · · ,x ] , (2.8) y=0 the (2.7) becomes −1 −1 2 −1 n sA~ n(x)= Wn[1, (h (y)), (h (y)) , · · · , (h (y)) ] y=0 −1 1 1 2 n T · Ωn Pn Pn [1, x, x , · · · ,x ] . (2.9) l(y) l(h(y)) y=0 y=0 Now, by taking the kth order derivative of both sides of the equation (2.9) with respect to x and dividing the resulting equation by k!, we obtain 1 [ A(k)(x), A(k)(x), · · · , A(k)(x)]T k! s 0 s 1 s n −1 −1 2 −1 n = Wn[1, (h (y)), (h (y)) , · · · , (h (y)) ] y=0 −1 1 1 · Ωn Pn Pn l(y) l(h(y)) y=0 y=0 T k + 1 k + 1 n · 0, · · · , 0, 1, x, x2, · · · , xn−k . (2.10) " k ! k ! k! # 6 H.M.SRIVASTAVA,SAIMAJABEEANDMOHAMMADSHADAB Hence, clearly, the right-hand side and left-hand side of the equation (2.10) are the kth columns of −1 −1 2 −1 n −1 Wn[1, (h (y)), (h (y)) , · · · , (h (y)) ] Ωn , y=0 1 1 xy Pn Pn Pn[e ] l(y) l(h(y)) y=0 y=0 y=0 and T −1 Wn[sA0(x), sA1(x), · · · , sAn(x)] Ω , respectively.