arxiv: v1 [math.nt] 4 Nov 2018

Transcription

1 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES REVISITED arxiv:8.0455v math.nt] 4 Nov 08 YAMILET QUINTANA ), WILLIAM RAMÍREZ), AND ALEJANDRO URIELES ) Abstract. This paper is concerned with the generalized Euler polynomial matrix E α) x) and the Euler matrix E. Taing into account some properties of Euler polynomials and numbers, we deduce product formulae for E α) x) and determine the inverse matrix of E. We establish some explicit expressions for the Euler polynomial matrix Ex), which involving the generalized Pascal, Fibonacci and Lucas matrices, respectively. From these formulae we get some new interesting identities involving Fibonacci and Lucas numbers. Also, we provide some factorizations of the Euler polynomial matrix in terms of Stirling matrices, as well as, a connection between the shifted Euler matrices and Vandermonde matrices.. Introduction The classical Euler polynomials E n x) and the generalized Euler polynomials E α) n x) of real or complex) order α, are usually defined as follows see, for details,, 5, 0, ]):.) and ) α e z e xz + n0 E α) n x)zn n!, z < π, α :,.) E n x) : E ) n x), n N 0, where N 0 : N {0} and N {,,3,...}. The numbers E n α) : E n α) 0) are called generalized Euler numbers of order α, n N 0. It is well-nown that the classical Euler numbers are defined by the generating function.3) e z +e z z n ε n n!. The sequence {ε n } n 0 counts the numbers of alternating n-permutations. Let recall us that a permutation σ of a set of n elements or n-permutation), is said alternating if and only if the n differences σi + ) σi) for i,,...,n have alternating signs cf. 5, p. 58]). From.) n0 00 Mathematics Subect Classification. B68, B83, B39, 05A9. Key words and phrases. Euler polynomials, Euler matrix, generalized Euler matrix, generalized Pascal matrix, Fibonacci matrix, Lucas matrix. ) Partially supported by the grants Impacto Caribe IC ) from Universidad del Atlántico, Colombia, and DID-USB S-IC-CB-003-6) from Decanato de Investigación y Desarrollo. Universidad Simón Bolívar, Venezuela. ) Partially supported by the grant Impacto Caribe IC ) from Universidad del Atlántico, Colombia.

2 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES and.3) it is easy to chec that the connection between the classical Euler numbers and the Euler polynomials is given by the formula ).4) ε n n E n, n N 0. So, the numbers E n : E n 0) also are nown in the literature as Euler numbers cf., e.g., 4, 0]). The first six generalized Euler polynomials are E α) 0 x), E α) x) x α, Eα) x) x αx+ αα ), 4 E α) 3 x) x 3 3α x + 3αα ) x 3α α ), 4 8 E α) 4 x) x 4 αx 3 + 3αα ) x α α 3) E α) 5 x) x 5 5α x4 + 5αα ) x 3 5α α 3) 4 α α 3 0α +5α+0). 3 x+ αα3 6α +3α 6), 6 x + 5αα )α 5α ) x 6 Recent and interesting wors dealing with these polynomials, Appell and Apostol type polynomials, their properties and applications in several areas as such as combinatorics, number theory, numerical analysis and partial differential equations, can be found by reviewing the current literature on this subect. For a broad information on old literature and new research trends about these classes of polynomials we strongly recommend to the interested reader see 5, 8 0, 4, 3 6]. From the generating relation.), it is fairly straightforward to deduce the addition formula:.5) E α+β) n x+y) And, it follows also that ) n E α) x)eβ) n y)..6) E n α) x+)+e n α) x) E n α ) x). Since E n 0) x) x n, maing the substitution β 0 into.5) and interchanging x and y, we get ) n.7) E n α) x+y) E α) y)xn. And, as an immediate consequence, we have.8).9) E n x+y) E n x) ) n E y)x n, ) n E x n.

3 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 3 Using.4),.8) and the well-nown relation E n x) ) n E n x), it is possible to deduce the following connection formula between E n and the classical Euler numbers ε n : n n ) n εn, if n is odd,.0) E n 0, if n is even. Inspired by the article 30] in which the authors introduce the generalized Bernoulli matrix and establish some algebraic properties of the Bernoulli polynomial and Bernoulli matrices, in the present article we focus our attention on the algebraic and differential properties of the generalized Euler matrix. It is worthwhile to mention that the authors of 30] to point out that their proposed methodology can be used for obtaining similar properties in the setting of the generalized Euler matrix. However, the authors of 30] do not actually write out any proof about this statement. The outline of the paper is as follows. Section has an auxiliary character and provides some bacground as well as some results which will be used throughout the paper. Maing use of the some identities above, we introduce the generalized Euler matrix in Section 3. Then, we study some interesting particular cases of this matrix, namely, the Euler polynomial matrix, the Euler matrix and the specialized Euler matrix. The main results of this section are Theorems 3., 3., 3.3 and 3.4, because these theorems contain the information concerning the product formula for the Euler matrix, an explicit expression for the inverse matrix of the specialized Euler matrix, the factorization of the Euler matrix via the generalized Pascal matrix of first ind, and a useful factorization for the inverse matrix of a particular horizontal sliding of the Euler polynomial matrix, respectively. Also, some consequences of these results are showed see for instance, Corollaries 3., 3., 3.3 and 3.4). Section 4 shows several factorizations of the generalized Euler matrix in terms the Fibonacci and Lucas matrices, respectively cf. Theorems 4. and 4.). Also, some new identities involving Fibonacci and Lucas numbers are given in this section. Finally, in Section 5 we provide some factorizations of the Euler polynomial matrix in terms of Stirling matrices, and the shifted Euler matrices and their connection with Vandermonde matrices are given.. Bacground and previous results Throughout this paper, all matrices are in M n+ R), the set of all n+)-square matrices over the real field. Also, for i, any nonnegative integers we adopt the following convention ) i 0, whenever > i. In this section we recall the definitions of the generalized Pascal matrix, the Fibonacci matrix and the Lucas matrix, as well as, some properties of these matrices. Definition.. Let x be any nonzero real number. The generalized Pascal matrix of first ind Px] is an n+) n+) matrix whose entries are given by see, 8]): i ) x i, i,.) p i, x) In, 8, 9] some properties of the generalized Pascal matrix of first ind are showed, for example, its matrix factorization by special summation matrices, its associated differential equation and its

4 4 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES bivariate extensions. The following proposition summarizes some algebraic and differential properties of Px]. Proposition.. Let Px] be the generalized Pascal matrix of first ind and order n+. Then the following statements hold. a) Special value. If the convention 0 0 is adopted, then it is possible to define.) P0] : I n+ diag,,...,), where I n+ denotes the identity matrix of order n+. b) Px] is an invertible matrix and its inverse is given by.3) P x] : Px]) P x]. c), Theorem ] Addition theorem of the argument. For x,y R we have.4) Px+y] Px]Py]. d), Theorem 5] Differential relation Appell type polynomial entries). Px] satisfies the following differential equation.5) D x Px] LPx] Px]L, where D x Px] is the matrix resulting from taing the derivative with respect to x of each entry of Px] and the entries of the n+) n+) matrix L are given by p i, 0), i, l i, 0, otherwise, +, i +, e) 8, Theorem ]) The matrix Px] can be factorized as follows..6) Px] G n x]g n x] G x], where G x] is the n+) n+) summation matrix given by ] In 0,,...,n, 0 S G x] x] S n x], n, being S x] the +) +) matrix whose entries S x;i,) are given by x i, i, S x;i,) 0 i, ). 0, > i, Another necessary structured matrices in what follows, are the Fibonacci and Lucas matrices. Below, we recall the definitions of each one of them.

5 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 5 Definition.. Let {F n } n be the Fibonacci sequence, i.e., F n F n +F n for n with initial conditions F 0 0 and F. The Fibonacci matrix F is an n+) n+) matrix whose entries are given by ]: F i +, i + 0,.7) f i, 0, i + < 0. Let F be the inverse of F and denote by f i, the entries of F. In ] the authors obtained the following explicit expression for F..8) fi,, i,, i +, +, Definition.3. Let {L n } n be the Lucas sequence, i.e., L n+ L n+ + L n for n with initial conditions L and L 3. The Lucas matrix L is an n+) n+) matrix whose entries are given by 3]: L i +, i 0,.9) l i, Let L bethe inverse of L and denote by l i, the entries of L. In 3, Theorem.] the authors obtained the following explicit expression for L..0) li,, i, 3, i +, 5 ) i i, i +, For x any nonzero real number, the following relation between the matrices Px] and L was stated and proved in 3, Theorem 3.]..) Px] L Gx] H x]l,

6 6 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES where the entries of the n+) n+) matrices Gx] and H x] are given by ) ) i i x g i, x) x x i+ 3x i +5 ) i+ i m i,+ ) ], ) i i h i, x) x x i+ 3x )+ ) i +5xi n i+,+3 respectively, with i ) ) x, i, m i, x) : 0, i <, )], x i ) i ) x, i, and n i, x) : 0, i <. 3. The generalized Euler matrix Definition 3.. The generalized n+) n+) Euler matrix E α) x) is defined by 3.) E α) i, x) i α) ) E i x), i, While, Ex) : E ) x) and E : E0) are called the Euler polynomial matrix and the Euler matrix, respectively. In the particular case x, we call the matrix E : E ) specialized Euler matrix. It is clear that.6) yields the following matrix identity: 3.3) E α) x+)+e α) x) E α ) x). Since E 0) x) Px], replacing α by in 3.3) we have 3.4) Ex+)+Ex) Px]. Then, putting x 0 in 3.4) and taing into account.), we get E)+E I n+. Analogously, E + E ) P ]. From.4) it follows that the entries of the specialized Euler matrix E are given by i ) i ε i, i, 3.5) e i,

7 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 7 From.0) it follows that the entries of the Euler matrix E are given by i ) Ei, i > and i odd, 3.6) E i,, i, The next result is an immediate consequence of Definition 3. and the addition formula.5). Theorem 3.. The generalized Euler matrix E α) x) satisfies the following product formula. 3.7) E α+β) x+y) E α) x)e β) y) E β) x)e α) y) E α) y)e β) x). Proof. We proceed as in the proof of 30, Theorem.], maing the corresponding modifications. Let A α,β) i, x,y) be the i,)-th entry of the matrix product E α) x)e β) y), then by the addition formula.5) we have ) ) A α,β) i i, x,y) E α) i x) E β) y) i i ) ) i E α) i x) E β) y) i ) i ) ) i E α) i i x)eβ) y) i ) i E α) i x)eβ) y) ) i E α) i x+y), which implies the first equality of 3.7). The second and third equalities of 3.7) can be derived in a similar way. Corollary 3.. Let x,...,x ) R. For α real or complex parameters, the Euler matrices E α ) x) satisfies the following product formula,,...,. 3.8) E α +α + +α ) x +x + +x ) E α ) x )E α ) x ) E α ) x ). Proof. The application of induction on gives the desired result. If we tae x x x x and α α α α, then we obtain the following simple formula for the powers of the generalized Euler matrix, and consequently, for the powers of the Euler polynomial and Euler matrices. Corollary 3.. The generalized Euler matrix E α) x) satisfies the following identity. 3.9) E x)) α) E α) x).

8 8 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES In particular, 3.30) Ex)) E ) x), E E ). Remar 3. Note that Theorem 3. and Corollaries 3. and 3. are respectively, the analogous of Theorem. and Corollaries. and.3 of 30] in the setting of Euler matrices. Let D be the n+) n+) matrix whose entries are defined by + ) i ) i ) i, i, 3.3) d i, Theorem 3.. The inverse matrix of the specialized Euler matrix E is given by Furthermore, E D. E ) )] D. Proof. Taing into account.4) and 3.5), it is possible to deduce + ) ) ) ) n n + ) ) ) n E n ε n δ n,0, where δ n,0 is the Kronecer delta cf., e.g., 8, pp ]). So, we obtain that the i,)-th entry of the matrix product DE may be written as i ) i + ) i ) ) ) i E i i ) i i i ) i ) i + ) i ) ) E ) i + ) i ) ) E ) i i δ i,0, and consequently, DE I n+. Similar arguments allow to show that ED I n+, and hence E D. Finally, from the identity E D and 3.30) we see that E )] ) E ) E ) D. This last chain of equalities finishes the proof.

9 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 9 It is worthwhile to mention that the calculation of E strongly depends on the use of inverse relations derived from exponential generating functions cf. 8, Chap. 3, Sec. 3.4]). This tool can be applied in order to determine E, but it not wors for determining of E. This fact and.0) suggest that methodology proposed in 30] does not suffice to finding an explicit formula for E. The next result establishes the relation between the generalized Euler matrix and the generalized Pascal matrix of first ind. Theorem 3.3. The generalized Euler matrix E α) x) satisfies the following relation. 3.3) E α) x+y) E α) x)py] Px]E α) y) E α) y)px]. In particular, 3.33) Ex+y) Px]Ey) Py]Ex), 3.34) Ex) Px]E, 3.35) E x+ ) Px]E, and 3.36) E P Proof. The substitution β 0 into 3.7) yields Since E 0) x) Px], we obtain ] E. E α) x+y) E α) x)e 0) y) E 0) x)e α) y) E α) y)e 0) x). E α) x+y) Px]E α) y). A similar argument allows to show that E α) x+y) E α) x)py] and E α) x+y) E α) y)px]. Next, the substitution α into 3.3) yields 3.33). From the substitutions y 0 and y into 3.33), we obtain the relations 3.34) and 3.35), respectively. Finally, the substitution x into 3.35) completes the proof. Remar 3. Note that the relation 3.33) is the analogous of 30, Eq. 3)] in the context of Euler polynomial matrices and, the counterpart of 3.34) is 30, Eq. 4)]. However, the relation 3.35) is slightly different from 30, Eq. 4)], since it involves an Euler polynomial matrix with shifted argument and the specialized Euler matrix. More precisely, the relation 30, Eq. 4)] reads as Bx) Px]B, consequently, this relation expresses to the Bernoulli polynomial matrix Bx) in terms of the matrix product between the generalized Pascal matrix of first ind Px] and the Bernoulli matrix B. While, on the left hand side of 3.35) appears an Euler polynomial matrix with shifted argument, and the matrix product on the right hand side of 3.35) contains to the specialized Euler matrix E. The following example shows the validity of Theorem 3.3.

10 0 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES Example 3.. Let us consider n 3. It follows from the definition 3. that E and E x+ ) x 0 0 x x x x 3x 3 4 3x On the other hand, from 3.35) and a simple computation we have E x+ ) x x x 0 x 3 3x x 0 0 x 3x x x }{{}}{{} x 3x 3 4 3x Px] E The next theorem follows by a simple computation. Theorem 3.4. The inverse of the Euler polynomial matrix E x+ ) can be expressed as 3.37) E x+ E )] P x] DP x]. In particular, 3.38) E DP ]. Proof. Using.3), 3.35) and Theorem 3. the relation 3.37) is deduced. The substitution x into 3.37) yields 3.38). Example 3.. Let us consider n 3. From the definition 3. and a standard computation we obtain E x+ )] x 0 0 x 4 x 0 x x 3x 3 4 3x On the other hand, from 3.37) we have E x+ )] x 0 0 x + 4 x 0. x x 3x x x x x 0 x 0 0 x x 3 3x 4 x 0. 3x x }{{}}{{} x 3x x D P x]

11 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES Hence, when x, we get E }{{}}{{} D P ] At this point, an apart mention deserves the recent wor ] since it states an explicit formula to the inverse matrix of the q-pascal matrix plus one in terms of the q-analogue of the Euler matrix E. As a consequence of the relations.6),.), and Theorems 3.3 and 3.4, we obtain the following corollaries. Corollary 3.3. The Euler polynomial matrix E x+ ) and its inverse can be factorized by summation matrices as follows. E x+ ) G n x]g n x] G x]e, E x+ DG n x]g n x] G x]. )] In particular, E G n ]G n ] G ] E, E DG n ] ] ] G n G. Corollary 3.4. For x any nonzero real number, the Euler polynomial matrix E x+ ) and its inverse can be factorized, respectively, in terms of the Lucas matrix L and its inverse as follows. E x+ ) L Gx]E H x]le, E x+ DGx]) )] L DL H x]). In particular, E L G ] E H ] LE, E D G ]) L DL H ]).

12 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES We end this section showing others identities, which can be easily deduced from the content of this paper. So, we will omit the details of their proofs. D x Ex+y) LPx]Ey), D x E D x Ex) LPx]E, x+ ) LPx]E, D x E x+ )] DLP x]. 4. Generalized Euler polynomial matrices via Fibonacci and Lucas matrices For 0 i, n and α a real or complex number, let M α) x) be the n+) n+) matrix whose entries are given by cf. 30, Eq. 8)]): ) ) ) 4.39) m α) i i, x) E α) i i x) E α) i i x) E α) i x). We denote Mx) M ) x) and M M0). Similarly, let N α) x) be the n+) n+) matrix whose entries are given by cf. 30, Eq. 3)]): ) ) ) 4.40) ñ α) i i, x) E α) i i x) E α) i i + x) E α) i + x). We denote N x) N ) x) and N N 0). From the definitions of M α) x) and N α) x), we see that m α) 0,0 m α) 0, m α),0 m α), m α) i,0 x) mα), x) ñα) 0,0 x) ñα) 0, x) 0,, x) ñα), x) Eα) 0 x), x) ñα),0 x) Eα) x) E α) 0 x) x α, x) ñα), x) 0,, x) ñα) i,0 x) Eα) i x) E α) i x) Eα) i x), i.

13 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 3 For 0 i, n and α a real or complex number, let L α) x) be the n+) n+) matrix whose entries are given by 4.4) ˆlα,) i, x) ) ) i E α) i i x) 3 E α) We denote L x) L ) x) and L L 0). Similarly, let L α) 4.4) ˆlα,) i, x) i x)+5 i ) i i x) be the n+) n+) matrix whose entries are given by ) ) i E α) i i ) i x) 3 E α) i i + x)+5 ) We denote L x) L ) x) and L L 0). From the definitions of L α) x) and L α) x), we see that ˆlα,) α,) i,i x) ˆl i,i x), i 0, ˆlα,) α,) 0, x) ˆl 0, x) 0,, + ˆlα,) α,),0 x) ˆl,0 x) E α) x) 3E α) 0 x) x α 3, ˆlα,) α,), x) ˆl, x) 0,, ˆlα,) i,0 x) E α) i x) 3E α) i x)+5 i )i i E α) x), i, ˆlα,) i,0 x) E α) i x) 3iE α) i x)+5 i ) ) i α) E i x), i, ) E α) x). ˆlα,) i, x) ie α) 7i ) i x) E α) i x)+5 i )i i E α) x), i 3. E α) i x). Thefollowingresultsshowseveralfactorizations of E α) x)intermsoffibonacciandlucasmatrices, respectively. Theorem 4.. The generalized Euler polynomial matrix E α) x) can be factorized in terms of the Fibonacci matrix F as follows. 4.43) E α) x) FM α) x), or, 4.44) E α) x) N α) x)f. In particular, 4.45) FMx) Ex) N x)f, 4.46) FM E N F,

14 4 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES and 4.47) FM ) E N ) F. Proof. Since the relation 4.43) is equivalent to F E α) x) M α) x), it is possible to follow the proof given in 30, Theorem 4.], maing the corresponding modifications, for obtaining 4.43). The relation 4.44) can be obtained using a similar procedure. The relations 4.45), 4.46) and 4.47) are straightforward consequences of 4.43) and 4.44). Also, the relations 4.43) and 4.44) allow us to deduce the following identity: M α) x) F N α) x)f. As a consequence of Theorems 3.4 and 4., we can derive simple factorizations for the inverses of the polynomial matrices M x+ ) and N x+ ) : Corollary 4.. The inverses of the polynomial matrices M x+ ) and N x+ ) can be factorized as follows. 4.48) M x+ DP x]f, )] 4.49) N x+ )] FDP x]. In particular, 4.50) M DP 4.5) M ] F, and N FDP )] DF, and N ], )] FD. An analogous reasoning as used in the proof of Theorem 4. allows us to prove the results below. Theorem 4.. The generalized Euler polynomial matrix E α) x) can be factorized in terms of the Lucas matrix L as follows. 4.5) E α) x) L L α) x), or, 4.53) E α) x) L α) x)l. In particular, 4.54) L L x) Ex) L x)l, 4.55) L L E L L, and 4.56) L L ) x) E L ) x)l.

15 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 5 Also, the relations 4.5) and 4.53) allow us to deduce the following identity: L α) x) L L α) x)l. Corollary 4.. The inverses of the polynomial matrices L x+ as follows. 4.57) L x+ )] DP x]l, ) and L x+ ) can be factorized 4.58) L x+ )] L DP x]. In particular, 4.59) L DP 4.60) ] L, and L L DP L )] DL, and ], L )] L D. Remar 4. It is worthwhile to mention that if we consider a C, b C\{0} and s 0,, then Theorems 4. and 4., as well as, their corollaries have corresponding analogous forms for generalized Fibonacci matrices of type s, F a,b,s), and for generalized Fibonacci matrices U a,b,0) with second order recurrent sequence U n a,b) subordinated to certain constraints. The reader may consult 7] in order to complete the details of this assertion. We finish this section with some new identities involving the Fibonacci numbers, the Lucas numbers and the generalized Euler polynomials and numbers. Theorem 4.3. For 0 r n and α any real or complex number, we have ) ] n E α) n r r x) F r +)α+ n r+ + r +)x F n r + r+ ){ E α) r r r x) F n r+ + n x α ) ] F n r n ){ n + E α) n n x) E α) + n E α) r r x)+ Eα) r ]}F x) n + n x)+ + E α) n x) ]}F r+. Proof. We proceed as in the proof of 30, Theorem 4.], maing the corresponding modifications. From 4.39), it is clear that m α) r +)α+ r,r x), mα) r+,r x) r +)x,

16 6 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES and, for : ){ ]} m α),r x) E α) r r r x) E α) r r x)+ Eα) r x). Next, it follows from 4.43) that ) n E α) n r x) Eα) n,r r x) r F n + m α),r x) F n r+ +F n r m α) r+,r x)+ F n r+ + r +)x + r+ r+ r +)α+ ){ E α) r r r x) F n + m α),r x) ] F n r E α) r r x)+ Eα) r ]}F x) n +. This chain of equalities completes the first part of the proof. The second one is obtained in a similar way, taing into account the following identities: ñ α) n,n x), ñα) n,n x x) n α ), and, for 0 n : ){ ñ α) n n, x) E α) n n x) E α) + n ) n n r n x)+ + Corollary 4.3. For 0 r n and α any real number, we have ) E α) r +)x α)+ n r x) )r F n r+ + ) r+ + r+ ) r n ) n +r n + ] F n r E α) n x) ]}. ){ E α) r r x)+ E α) r r x) Eα) r ]}F x) n + ) r F n r+ + ) r+ n x α ) ] F n r ){ E α) n n x)+ E α) + n + Proof. Replacing x by α x in 4.6) and applying the formula E α) n x) )n E α) n α x) x)+ n E α) n x) ]}F r+.

17 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 7 to the resulting identity, we obtain the first identity of Corollary 4.3. An analogous reasoning yields the second identity. Analogous reasonings to those used in the proofs of Theorem 4.3 and Corollary 4.3 allow us to prove the following results. Theorem 4.4. For any real or complex number α, we have the following identities E n α) x) L n+ + x α ) 3 L n + E α) x) 3Eα) )L x) n + 4.6) whenever n. 4.6) whenever n ) s s L n + E s α) x), s0 ne α) n x) L n +x α 3)L n ) s s sl n + E α) s x), 3 s Corollary 4.4. The following identities hold. 4.63) ) n E α) n x) L n+ whenever n. 4.64) whenever n x α +3 ) L n + E α) x) 3 )Eα) )L x) n + ) E α) x)+3eα) )L x) n + ) s s L n + E s α) x), s0 ) n ne α) n x) L n +α x 3)L n + ) E α) x)+3 )Eα) )L x) n ) s sl n + E α) s x), 3 s By 4.6), 4.6), 4.63) and 4.64) we obtain the following interesting identities involving Lucas and Euler numbers. For n : E n L n+ 7 ) L n E 3E +5 ) s s E s )L n +, s0

18 8 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES ) n E n L n+ 5 L n + For n 3: ) E +3E )L n + +5 ne n L n 4L n ) ) s s L n + E s. s0 E 3 )E )L n ) n ne n L n α x 3)L n ) s s se s L n +, 3 s ) E +3 )E )L n ) s sl n + E s. 3 s Another similar combinatorial identities may be obtained using the results of 3]. We leave to the interested reader the formulation of them. 5. Euler matrices and their relation with Stirling and Vandermonde matrices Let sn,) and Sn,) be the Stirling numbers of the first and second ind, which are respectively defined by the generating functions 0, Chapther, Section.6]: 5.65) sn,)z zz ) z n+), log+z))! z n n sn,) zn, z <, n! Sn,)zz ) z +), e z )! n Sn,) zn n!. In combinatorics, it is well-nown that the value sn, ) represents the number of permutations of n elements with disoint cycles. While, the Stirling numbers of the second ind Sn,) give the number of partitions of n obects into non-empty subsets. Another way to compute these numbers is by means of the formula see 6, Eq. 5.)] or 8, p. 6]: Sn,)! ) ) l l n, n. l l0

19 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 9 A recent connection between the Stirling numbers of the second ind and the Euler polynomials is given by the formula see 7, Theorem 3., Eq. 3.3)]): ) ] n + n 5.66) E n x) ) n ) l l )! l Sn +,l) x. l Proceeding as in the proofof 7, Theorem 3.], onecan findasimilar relation to theprevious one but connecting Stirling numbers of the first ind and a particular class of generalized Euler polynomials. Theorem 5.. Let us assume that α m N. Then, the connection between the Stirling numbers of the first ind and the generalized Euler polynomial E n m) x) is given by the formula: 5.67) E n m) x) ) n n n sn,) m) x), Proof. By Leibniz s theorem for differentiation we have r ) m ] r ) ) r m ] ) z r e z e xz r + e z + z r exz ) r ) ) r m ] ) e z x r e xz + ) m r ) r e z e x+)z m) x r + e z +), 0 where in the last expression m) denotes the falling factorial with opposite argument m. Combining this with the r-th differentiation on both sides of the generating function in.) reveals that nr E n m) x) z n r n r)! e z + ) m e x+)z r Further taing z 0 and employing 5.65) give r r E r m) x x) ) m) r r r r ) r r 0 ) r m) x r e z +). r ) r m) r x) sr,) m) x). Finally, changing r by n the proof of the formula 5.66) is complete. Definition 5.. For the Stirling numbers si,) and Si,) of the first ind and of the second ind respectively, define S and S to be the n+) n+) matrices by { { si,), i, Si,), i, 5.68) S i, and S 0, otherwise, i,

20 0 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES The matrices S and S are called Stirling matrix of the first ind and of the second ind, respectively see 3]). In order to obtain factorizations for Euler matrices via Stirling matrices we will need the following matrices: Let S n be the factorial Stirling matrix, i.e., the n n matrix whose i,)-th entry is given by S i,,n :!S i,, i and otherwise 0. For m N, let S m) be the n+) n+) matrix whose i,)-entries are defined by 5.69) S m) i, i ) i si,) m), i, Let C and D be the n+) n+) matrices whose i,)-entries are defined by i ) ) i i ) Si,,i, i, 5.70) Ci, i ) ) i i ) Si,+,i, i, 5.7) Di, The next theorem shows the corresponding factorizations of the generalized Euler matrix E m), m N, in terms of the Stirling matrices, when the expressions 5.66) and 5.67) are incorporated. Theorem 5.. For m N, the generalized Euler matrix E m) x) can be factorized as follows. 5.7) E m) x) S m) Px]. In the case of the Stirling matrix of the second ind, we have 5.73) Ex) C + D)Px]. Furthermore, 5.74) S ) C + D. Proof. For m N and i, let A m) i, x) be the i,)-th entry of the matrix product Sm) Px], then i i ) ) ] A m) i, x) S m) i i, p i,x) si,r) m) r x r0 i ) i ] i i ) i si,r) m) r x) i ) i i i ) i r0 r0 si,r) m) r ]x) ) i E m) i x),

21 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES The last equality is an immediate consequence of 5.67), and 5.7) follows from the previous chain of equalities. In order to prove 5.73), we proceed in a similar way to the previous one, taning into account 5.66), 5.70), 5.7), and maing the corresponding modifications. Finally, the substitution m into 5.7) yields 5.74). Definition 5.. The n+) n+) shifted Euler polynomial matrix Ẽx) are given by 5.75) Ẽ i, x) E i +x), 0 i, n. Let us consider the Vandermonde matrix: x +x +x n+x V x) : x +x) +x) n+x) x n +x) n +x) n n+x) n In 4, Theorem.], the following factorization for the Vandermonde matrix V x) was stated. 5.76) V x) ] S n ) n+ x)p T : ] S n ) n+ x)p]) T, where n+ x)p]) T represents the LU-factorization of a lower triangular matrix whose i,)-th entry is x i ), if i and otherwise 0. The relation between the shifted Euler polynomial matrix Ẽx) and the matrices V x) and S n is contained in the following result. Theorem 5.3. The shifted Euler polynomial matrix Ẽx) can be factorized in terms of the Vandermonde matrix V x) and consequently, in terms of the factorial Stirling matrix S n as follows: 5.77) Ẽx) E V x), 5.78) Ẽx) E] S n ) n+ x)p T. Proof. Let Ẽi,x) be the i,)-th entry of the shifted Euler polynomial matrix Ẽx). Then, using.9) we get i ) i i Ẽ i, x) E i +x) E i +x) E i, V, x). Hence, 5.77) follows from this chain of equalities. The relation 5.78) is a straightforward consequence of 5.76). Remar 5. Note that the relations 5.77) and 5.78) are the analogous of 30, Eqs. 37), 38)], respectively, in the context of Euler polynomial matrices. Finally, in the present paper, all matrix identities have been expressed using finite matrices. Since such matrix identities involve lower triangular matrices, they have a resemblance for infinite matrices. We state this property briefly as follows.

22 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES Let E α) ) x), E x+, E x), E, E, D, Ẽ x), P x], F, F, L, G x], H x], M α) ) x), M x+ α) ), N x), N x+, V and S m), be the infinite cases of the matrices E α) x), E x+ ), Ex), E, E, D, Ẽx), Px], F, F, L, Gx], H x], M α) x), M x+ ), N α) x), N x+ ), V and S m) respectively. Then the following identities hold. E α ) x) E α) x+)+e α) x), E α) x+y) E α) x)p y] P x]e α) y) E α) y)p x], E x+y) P x]e y) P y]e x), E x) P x]e, E x+ ) P x]e, E x+ )] D P x], E x+ ) L G x]e H x]l E, E α) x) F M α) α) x) N x)f, M α) x) F N α) x)f, M x+ )] D P x]f, N x+ )] F D P x], E m) x) S P x], Ẽ x) E V x). References ] T. Araawa, T. Ibuiyama, M. Kaneo, Bernoulli Numbers and Zeta Functions, Springer, New Yor 04). ] G.S. Call, D.J. Velleman, Pascal s matrices, Amer. Math. Monthly 00, ). 3] G.-S. Cheon, J.-S. Kim, Stirling matrix via Pascal matrix, Linear Algebra Appl. 39, ). 4] G.-S. Cheon, J.-S. Kim, Factorial Stirling matrix and related combinatorial sequences, Linear Algebra Appl. 357, ). 5] L. Comtet, Advanced Combinatorics. The Art of Finite and Infinite Expansions, D. Reidel Publishing Co., Dordrecht 974). 6] T. Ernst, A comprehensive treatment of q-calculus, Birhäuser, Boston 0). 7] B.-N. Guoa, F. Qi, Explicit formulae for computing Euler polynomials in terms of Stirling numbers of the second ind, J. Comput. Appl. Math. 7, ). 8] Y. He, S. Araci, H.M. Srivastava, Some new formulas for the products of the Apostol type polynomials, Adv. Differ. Equ. 06Article ID 87), 8 06).

23 EULER MATRICES AND THEIR ALGEBRAIC PROPERTIES 3 9] Y. He, S. Araci, H.M. Srivastava, M. Acigöz, Some new identities for the Apostol-Bernoulli polynomials and the Apostol-Genocchi polynomials, Appl. Math. Comput. 6, ). 0] P. Hernández-Llanos, Y. Quintana, A. Urieles, About extensions of generalized Apostol-type polynomials, Results Math. 68, ). ] G.I. Infante, J.L. Ramírez, A. Şahin, Some results on q-analogue of the Bernoulli, Euler and Fibonacci matrices, Math. Rep. Bucur.) 969) No. 4, ). ] G.-Y. Lee, J.-S. Kim, S.-G. Lee, Factorizations and eigenvalues of Fibonacci and symmetric Fibonacci matrices, Fibonacci Quart. 403), 03 00). 3] G.-Y. Lee, J.-S. Kim, S.-H. Cho, Some combinatorial identities via Fibonacci numbers, Discrete Appl. Math. 303), ). 4] Q.-M. Luo, H.M. Srivastava, Some generalizations of the Apostol-Bernoulli and Apostol-Euler polynomials, J. Math. Anal. Appl. 308), ). 5] N.E. Nørlund, Vorlesungen über Differenzenrechnung, Springer-Verlag, Berlin 94) reprinted 954), in German). 6] Á. Pintér, H.M. Srivastava, Addition theorems for the Appell polynomials and the associated classes of polynomial expansions, Aequationes Math. 85, ). 7] Y. Quintana, W. Ramírez, A. Urieles, On an operational matrix method based on generalized Bernoulli polynomials of level m, Calcolo 553) 9 pages, 08). 8] J. Riordan, Combinatorial Identities, Wiley, New Yor, 968). 9] H.M. Srivastava, M.A. Boutiche, M. Rahmani, A class of Frobenius-type Eulerian polynomials, Rocy Mountain J. Math. 483), ). 0] H.M. Srivastava, J. Choi, Zeta and q-zeta Functions and Associated Series and Integrals, Elsevier, London 0). ] H.M. Srivastava, I. Kucuoğlu, Y. Simse, Partial differential equations for a new family of numbers and polynomials unifying the Apostol-type numbers and the Apostol-type polynomials, J. Number Theory 8, ). ] H.M. Srivastava, H.L. Manocha, A Treatise on Generating Functions, Ellis Horwood Ltd., West Sussex 984). 3] H.M. Srivastava, M. Mased-Jamei, M. Reza Beyi, A parametric type of the Apostol-Bernoulli, Apostol-Euler and Apostol-Genocchi polynomials, Appl. Math. Inform. Sci. 5), , 08). 4] H.M. Srivastava, M.A. Özarslan, C. Kaanuglu, Some generalized Lagrange-based Apostol- Bernoulli, Apostol-Euler and Apostol-Genocchi polynomials, Russian J. Math. Phys. 0, ). 5] H.M. Srivastava, M.A. Özarslan, B. Yilmaz, Some families of differential equations associated with the Hermite-based Appell polynomials and other classes of Hermite-based polynomials, Filomat 84), ). 6] H.M. Srivastava, Á. Pintér, Remars on some relationships between the Bernoulli and Euler polynomials, Appl. Math. Lett. 7, ). 7] P. Stanimirović, J. Niolov, I. Stanimirović, A generalization of Fibonacci and Lucas matrices, Discrete Appl. Math. 564), ). 8] Z.Z. Zhang, The linear algebra of generalized Pascal matrix, Linear Algebra Appl. 50, ).

24 4 Y. QUINTANA, W. RAMÍREZ, AND A. URIELES 9] Z.Z. Zhang, M.X. Liu, An extension of generalized Pascal matrix and its algebraic properties, Linear Algebra Appl. 7, ). 30] Z. Zhang, J. Wang, Bernoulli matrix and its algebraic properties, Discrete Appl. Math. 54, ). 3] Z. Zhang, Y. Zhang, The Lucas matrix and some combinatorial identities, Indian J. Pure Appl. Math. 385), ). Departamento de Matemáticas Puras y Aplicadas, Edificio Matemáticas y Sistemas MYS), Apartado Postal: 89000, Caracas 080 A, Universidad Simón Bolívar, Venezuela address: yquintana@usb.ve Departamento de Ciencias Básicas, Universidad de la Costa - CUC, Barranquilla, Colombia. address: wramirez4@cuc.edu.co Programa de Matemáticas, Universidad del Atlántico, Km 7 Vía Pto. Colombia, Barranquilla, Colombia. address: aleandrourieles@mail.uniatlantico.edu.co