Next Article in Journal
Maximal Product and Symmetric Difference of Complex Fuzzy Graph with Application
Next Article in Special Issue
Some New Results for the Kampé de Fériet Function with an Application
Previous Article in Journal
Evaluation of the Heterogeneity of Wheat Kernels as a Traditional Model Object in Connection with the Asymmetry of Development
Previous Article in Special Issue
Certain Generalizations of Quadratic Transformations of Hypergeometric and Generalized Hypergeometric Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lagrange-Based Hypergeometric Bernoulli Polynomials

1
Department of Mathematics, College of Science, University of Ha’il, Ha’il 2440, Saudi Arabia
2
Departamento de Matemáticas, Universidad Carlos III de Madrid, Avenida de la Universidad 30, Leganés, 28911 Madrid, Spain
3
Department of Mathematics and Natural Sciences, Prince Mohammad Bin Fahd University, P.O. Box 1664, Al Khobar 31952, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(6), 1125; https://doi.org/10.3390/sym14061125
Submission received: 15 April 2022 / Revised: 22 May 2022 / Accepted: 24 May 2022 / Published: 30 May 2022
(This article belongs to the Special Issue Recent Advances in Special Functions and Their Applications)

Abstract

:
Special polynomials play an important role in several subjects of mathematics, engineering, and theoretical physics. Many problems arising in mathematics, engineering, and mathematical physics are framed in terms of differential equations. In this paper, we introduce the family of the Lagrange-based hypergeometric Bernoulli polynomials via the generating function method. We state some algebraic and differential properties for this family of extensions of the Lagrange-based Bernoulli polynomials, as well as a matrix-inversion formula involving these polynomials. Moreover, a generating relation involving the Stirling numbers of the second kind was derived. In fact, future investigations in this subject could be addressed for the potential applications of these polynomials in the aforementioned disciplines.

1. Introduction

As it is well-known, diverse differential equations can only be treated by utilizing families of special polynomials that provide novel viewpoints of mathematical analysis. Moreover, these special polynomials yield the derivation of other useful identities in a fairly straightforward manner and allow the consideration of new families of special polynomials. In addition, it is important that any polynomial has explicit formulas, symmetric identities, summation formulas, and relations with other polynomials.
The Lagrange polynomials g n α 1 , , α r x 1 , , x r in the variables x 1 , , x r and complex parameters α j ( j = 1 , , r ) are defined by means of the following generating function:
j = 1 r 1 x j z α j = n = 0 g n α 1 , , α r x 1 , , x r z n ,
where | z | < min { | x 1 | 1 , , | x r | 1 } and 1 α j : = 1 for j = 1 , , r . This class of multivariate polynomials (also known as the class of Chan–Chyan–Srivastava polynomials) was introduced in [1].
It is clear that the Lagrange polynomials g n α 1 , , α r x 1 , , x r provide a natural extension of the class of bivariate Lagrange polynomials:
( 1 x z ) α ( 1 y z ) β = n = 0 g n ( α , β ) ( x , y ) z n , | z | < min | x | 1 , | y | 1 ,
where α and β are complex numbers.
As is well known, these bivariate polynomials appear in some statistics problems (cf., e.g., [2] (p. 267), and (Chs. 1,7,8)) in [3] and can be expressed as follows:
g n ( α , β ) ( x , y ) = ( y x ) n P n ( α n , β n ) x + y x y , x y , x n P n ( α + β 1 , β n ) 1 , x = y ,
where P n ( α , β ) ( x ) is nth classical Jacobi polynomial given by (cf., (Equation (4.3.2)) in [4]).
P n ( α , β ) ( x ) : = k = 0 n n + α k n + β n k x + 1 2 k x 1 2 n k .
The seminal idea underlaying recent studies about special polynomials related to Lagrange polynomials (1) has been to make appropriate modifications for the generating functions associated with these polynomials by mixing generating functions that follow directly from a multiparameter and multivariate extension of Carlitz theorem (cf., (Ch. 7, Sec. 7.6)) in [5] and obtaining similar algebraic and/or differential properties for them (see, for instance, [6,7,8,9]).
Following the same methodology, one can consider for a fixed natural number m the hypergeometric Bernoulli polynomials (also called generalized Bernoulli polynomials of level m) defined by means of the following generating function [5,8,10,11,12,13]:
z m e x z e z l = 0 m 1 z l l ! = n = 0 B n [ m 1 ] ( x ) z n n ! , | z | < 2 π ,
and define the Lagrange-based hypergeometric Bernoulli polynomials in variables x , x 1 , , x r , and complex parameters α j ( j = 1 , , r ) as follows:
j = 1 r 1 x j z α j z m e x z e z l = 0 m 1 z l l ! = n = 0 B n [ m 1 , α 1 , , α r ] x | x 1 , , x r z n ,
where | z | < min { 2 π , | x 1 | 1 , , | x r | 1 } and 1 α j : = 1 for j = 1 , , r .
It is clear that this new class of special polynomials generalizes to the families of Lagrange-based Bernoulli polynomials (cf., (Equation (7))) in [9] and the hypergeometric Bernoulli polynomials and, hence, to the classical Bernoulli polynomials. Furthermore, if T n , λ , k α 1 , , α r ; α x 1 , x r ; x , y and T n , β , k α 1 , , α r x | x 1 , , x r ; a , b denote, respectively, the Lagrange-based Apostol type Hermite (cf., (Equation (2.1))) in [14] and the Lagrange-based unified Apostol-type polynomials (cf., (Equation (2.1)), in [15] given by the following:
j = 1 r 1 x j z α j 2 k z λ e z + ( 1 ) k + 1 α e x t + y t 2 = n = 0 T n , λ , k α 1 , , α r ; α x 1 , x r ; x , y z n , a n d
2 1 k z k β b e z a b e x z j = 1 r 1 x j z j α j = n = 0 T n , β , k α 1 , , α r x | x 1 , x r ; a , b z n .
then it is not difficult to see from (4) that the following is the case.
B n [ 0 , α 1 , , α r ] x | x 1 , , x r = T n , 1 , 0 α 1 , , α r ; 1 x 1 x r ; x , 0 ,
and
B n [ 0 , α 1 , α 2 ] x | x 1 , x 2 = T n , β , 1 α 1 , α 2 x | x 1 , x 2 ; a , 0 .
Inspired by the recent articles [13,14,15,16,17,20] in which the authors introduce the ( p , q ) -Hermite and ( p , q ) -Bernstein polynomials, the generalized Lagrange-based Apostol-type polynomials, the generalized Lagrange-based Apostol type Hermite polynomials, and Laguerre-based Hermite-Bernoulli polynomials associated with bilateral series and studied several analytic/numerical aspects of generalized Bernoulli polynomials of level m and the generalized mixed type Bernoulli–Gegenbauer polynomials, respectively, in the present article, we focus our attention on some algebraic and differential properties of polynomials B n [ m 1 , α 1 , , α r ] x | x 1 , , x r and its corresponding matrix-inversion formula.
Moreover, it is worthy to mention that the use of the Cauchy product of a power series is the technique behind these formulations. This approach is not a novelty; however, it has been useful for generating new families of special polynomials (satisfying or not Appell-type conditions), even those explored very recently. In this regard, we refer the interested reader to [17,18] and the references cited therein for a detailed exposition about some very recent trends in this broad field.
The paper is organized as follows. Section 2 contains the basic background about the Lagrange polynomials and the hypergeometric Bernoulli polynomials and some other auxiliary results that will be used throughout the paper. In Section 3, we prove some relevant algebraic and differential properties of the Lagrange-based hypergeometric Bernoulli polynomials (4) (Theorem 1), as well as their relation with the Stirling numbers of second kind (Theorem 2). Finally, we derive matrix-representation formulas for these polynomials (Theorems 3 and 4)

2. Background and Previous Results

Throughout this paper, let N , N 0 , R + , and C denote, respectively, the set of all natural numbers, the set of all nonnegative integers, the set of all positive real numbers, and the set of all complex numbers, and P n denotes the linear space of polynomials with real coefficients and a degree less than or equal to n. For λ C and k N 0 , we use notations ( λ ) k and λ ( n ) for the rising and falling factorials, respectively:
( λ ) 0 = 1 , ( λ ) k = λ ( λ + 1 ) ( λ + k 1 ) ,
and the following.
λ ( 0 ) = 1 , λ ( k ) = λ ( λ 1 ) ( λ k + 1 ) .
Moreover, as usual, the numbers given by the following:
B n [ m 1 ] : = B n [ m 1 ] ( 0 ) , f o r   a l l n 0 .
denote the hypergeometric Bernoulli numbers (or generalized Bernoulli numbers of level m N ). It is clear that if m = 1 in (3) then we recover the definition of the classical Bernoulli polynomials B n ( x ) , and classical Bernoulli numbers, respectively, i.e., B n ( x ) = B n [ 0 ] ( x ) , and B n = B n [ 0 ] , respectively, for all n 0 .
It is worth noticing here that there exist many families of special polynomials (both univariate and multivariate) generalizing the classical Bernoulli polynomials: for instance, those one drawing on the formalism and techniques of exponential operators or unified versions ((usually including Apostol-type generalizations and their further reductions), we refer the interested reader to [19,21,22] for more details).
Clearly, (1) yields the following explicit representation (cf., ( [1], Equation (6))).
g n α 1 , , α r x 1 , , x r = k 1 + + k r = n α 1 + k 1 1 k 1 α r + k r 1 k r x 1 k 1 x r k r .
Hence, using α j + k j 1 k j = ( α j ) k j k j ! , for j = 1 , , r , we obtain the following equivalent identity.
g n α 1 , , α r x 1 , , x r = k 1 + + k r = n ( α 1 ) k 1 ( α r ) k r x 1 k 1 k 1 ! x r k r k r ! ,
g n ( 0 , , 0 ) x 1 , , x r = k 1 + + k r = n x 1 k 1 k 1 ! x r k r k r ! .
Moreover, from (4), it is clear that the following is the case.
B n [ m 1 , α 1 , , α r ] x | x 1 , , x r = k = 0 n B k [ m 1 ] ( x ) k ! g n k ( α 1 , , α r ) ( x 1 , , x r ) = k = 0 n g k ( α 1 , , α r ) ( x 1 , , x r ) B n k [ m 1 ] ( x ) ( n k ) ! .
Thus, if we take α 1 = = α r = 0 and combine (4), (3), and (6), we have the following:
B n [ m 1 , 0 , , 0 ] x | x 1 , , x r = B n [ m 1 ] ( x ) n ! = k = 0 n g k 0 , , 0 x 1 , , x r B n k [ m 1 ] ( x ) ( n k ) ! = k = 0 n k 1 + + k r = n x 1 k 1 k 1 ! x r k r k r ! B n k [ m 1 ] ( x ) ( n k ) ! ,
for any n 0 .
Moreover, from (5) and (7), we obtain that if x = 0 , then (4) induces the following multivariate polynomials.
B n [ m 1 , α 1 , , α r ] x 1 , , x r : = B n [ m 1 , α 1 , , α r ] 0 | x 1 , , x r , n 0 .
Thus, the substitution x 1 = = x r = y , with y 0 into (8) yields the following univariate polynomials.
B n [ m 1 , α 1 , , α r ] y : = B n [ m 1 , α 1 , , α r ] y , , y , n 0 .
In this case, j = 1 r ( 1 x j z ) α j = j = 1 r ( 1 y z ) α j = ( 1 y z ) ( α 1 + + α r ) , and from (4), we can deduce the following (see also [1] (Equation (36))).
g n ( α 1 , , α r ) ( y , , y ) = ( α 1 + + α r ) n n ! y n , α 1 + + α r 0 , n 0 , 0 , α 1 + + α r = 0 , n 1 , 1 , α 1 + + α r = 0 , n = 0 .
Notice the following.
g n ( α 1 , , α r ) ( 0 , , 0 ) = lim y 0 g n ( α 1 , , α r ) ( y , , y ) = lim y 0 ( α 1 + + α r ) n n ! y n = 1 , i f n = 0 , 0 , i f n 1 .
Consequently, (9) takes the following form:
B n [ m 1 , α 1 , , α r ] y = k = 0 n ( α 1 + + α r ) k k ! B n k [ m 1 ] ( n k ) ! y k ,
whenever n 0 and α 1 + + α r 0 .
On the other hand, for x 1 = = x r = y , with y 0 , the polynomials B n [ m 1 , α 1 , , α r ] y can be described by means of the following generating function:
1 y z ( α 1 + + α r ) z m e z l = 0 m 1 z l l ! = n = 0 B n [ m 1 , α 1 , , α r ] ( y ) z n ,
where | z | < min { 2 π , | y | 1 } and 1 α j : = 1 for j = 1 , , r . Now, if we assume that α 1 + + α r = 0 , then (12) becomes the following:
z m e z l = 0 m 1 z l l ! = n = 0 B n [ m 1 , α 1 , , α r ] ( y ) z n ,
and from this last equality, it is easily deducible that the following is the case.
B n [ m 1 , α 1 , , α r ] ( y ) = B n [ m 1 ] n ! , n 0 , y 0 .
Then, the following is the case.
B n [ m 1 , α 1 , , α r ] ( 0 ) = lim y 0 B n [ m 1 , α 1 , , α r ] ( y ) = B n [ m 1 ] n ! , n 0 .
Finally, from (11), (13), and (14), we conclude that for each n 0 , we have the following.
B n [ m 1 , α 1 , , α r ] y = k = 0 n ( α 1 + + α r ) k k ! B n k [ m 1 ] ( n k ) ! y k .
It is clear that for n N , the univariate polynomials in (9) are different from the hypergeometric Bernoulli polynomials (3) (cf., e.g., [7,8,11,12,13]). Furthermore, it is not difficult to see from summation Formula (11) that it is possible to find explicit expressions for polynomials B n [ m 1 , α 1 , , α r ] ( y ) when α 1 + + α r 0 . Indeed, the first ones are given by the following.
B 0 [ m 1 , α 1 , , α r ] ( y ) = m ! , B 1 [ m 1 , α 1 , , α r ] ( y ) = m ! ( α 1 + + α r ) y 1 m + 1 = B 1 [ m 1 ] ( α 1 + + α r ) y , B 2 [ m 1 , α 1 , , α r ] ( y ) = m ! ( α 1 + + α r ) ( 1 + α 1 + + α r ) 2 y 2 ( α 1 + + α r ) m + 1 y + 2 ( m + 1 ) 2 ( m + 2 ) .

3. The Polynomials B n [ m 1 , α 1 , , α r ] x | x 1 , , X R and Their Properties

Now, we can proceed to investigate some relevant properties of the Lagrange-based hypergeometric Bernoulli polynomials.
Theorem 1.
For a fixed m N , let B n [ m 1 , α 1 , , α r ] x | x 1 , , x r be the Lagrange-based hypergeometric Bernoulli polynomials in the variables x , x 1 , , x r , and parameters α j C ( j = 1 , , r ). Then, the following statements hold:
(a) 
Summation formulas. For every n 0 , we have the following.
B n [ m 1 , α 1 , , α r ] x + y | x 1 , , x r = k = 0 n B n k [ m 1 , α 1 , , α r ] x | x 1 , , x r y k k ! = k = 0 n B k [ m 1 , α 1 , , α r ] y | x 1 , , x r x n k ( n k ) ! .
In particular, the following obtains.
B n [ m 1 , α 1 , , α r ] x | x 1 , , x r = k = 0 n B k [ m 1 , α 1 , , α r ] x 1 , , x r x n k ( n k ) ! .
(b) 
Differential relations (Appell-type polynomial sequences). For n , j 0 with
0 j n and any nonzero x 1 , , x r , we have the following.
j x j B n [ m 1 , α 1 , , α r ] x | x 1 , , x r = B n j [ m 1 , α 1 , , α r ] x | x 1 , , x r .
(c) 
Representation formulas. If at least an α j is nonzero, j = 1 , r , then the Lagrange polynomials g n α 1 , , α r x 1 , , x r can be expressed in terms of multivariate polynomials (8) as follows.
g n α 1 , , α r x 1 , , x r = k = 0 n B n k [ m 1 , α 1 , , α r ] x 1 , , x r ( n k ) ! ( m + k ) ! .
(d) 
Inversion formula. If α 1 + + α r 0 , then the following is the case.
( α 1 + + α r ) n n ! y n = k = 0 n B n k [ m 1 , α 1 , , α r ] ( y ) ( n k ) ! ( m + k ) ! , n 0 .
(e) 
Integral formulas. For any nonzero x 1 , , x r , we have the following.
y 0 y 1 B n [ m 1 , α 1 , , α r ] x | x 1 , , x r d x = B n + 1 [ m 1 , α 1 , , α r ] y 1 | x 1 , , x r B n + 1 [ m 1 , α 1 , , α r ] y 0 | x 1 , , x r = k = 0 n + 1 B k [ m 1 ] ( y 1 ) B k [ m 1 ] ( y 0 ) k ! g n + 1 k ( α 1 , , α r ) ( x 1 , , x r ) .
In particular, we have the following.
B n [ m 1 , α 1 , , α r ] x | x 1 , , x r = 0 x B n 1 [ m 1 , α 1 , , α r ] t | x 1 , , x r d t + B n [ m 1 , α 1 , , α r ] x 1 , , x r .
Proof. 
Since (a), (b), and (e) are straightforward consequences of (4) and a suitable use of the Fundamental Theorem of Calculus, respectively, we shall omit their proof. Thus, we focus our efforts on the proof of (c) and (d).
Assume that at least an α j is nonzero, j = 1 , r . By (4), (8), and direct calculations, we have the following.
z m j = 1 r 1 x j z α j = e z l = 0 m 1 z l l ! n = 0 B n k [ m 1 , α 1 , , α r ] x 1 , , x r z n , = n = 0 z n + m ( n + m ) ! n = 0 B n k [ m 1 , α 1 , , α r ] x 1 , , x r z n .
Or equivalently, we also have the following.
z m j = 1 r 1 x j z α j = n = 0 k = 0 n B n k [ m 1 , α 1 , , α r ] x 1 , , x r ( n k ) ! ( m + k ) ! z n + m .
Now, from (1), we have the following.
z m j = 1 r 1 x j z α j = n = 0 g n α 1 , , α r x 1 , , x r z n + m .
Hence, comparing the coefficients of z n + m on the right hand side of (19) and (20), we obtain (17).
Finally, assume that α 1 + + α r 0 and take x 1 = = x r = y . Then, the substitution of (10) into (17) and the use of (9) yield (18). □
The combination of (2) and (17) provides the following connection formula between Lagrange-based hypergeometric Bernoulli polynomials and Jacobi polynomials:
k = 0 n B n k [ m 1 , α , β ] x , y ( n k ) ! ( m + k ) ! = ( y x ) n P n ( α n , β n ) x + y x y , x y , x n P n ( α + β 1 , β n ) 1 , x = y ,
where P n ( α , β ) ( x ) is nth classical Jacobi polynomial.
Moreover, notice that the inversion formula (18) immediately implies the following.
Proposition 1.
If α 1 + + α r 0 , then for a fixed m N and each n 0 , the set B 0 [ m 1 , α 1 , , α r ] ( y ) , B 1 [ m 1 , α 1 , , α r ] ( y ) , , B n [ m 1 , α 1 , , α r ] ( y ) is a basis for P n .
With respect to the study of zeros of the nth Lagrange-based hypergeometric Bernoulli polynomial B n [ m 1 , α 1 , , α r ] x | x 1 , , x r when x 1 , , x r are fixed, relatively little is known. For instance, it is possible to use the Hurwitz theorem (see (Chapter I, p. 22)) in [4] for obtaining the fact that the complex zeros of B n [ m 1 , α 1 , , α r ] x | x 1 , , x r must move further away from the origin as n proceeds to infinity, because the functions to which they converge only have real zeros. In Figure 1, the plots for the zeros of B 11 [ m 1 , α 1 , , α r ] ( y ) and B 50 [ m 1 , α 1 , , α r ] ( y ) are shown for prescribed values of m and α 1 , , α r .
There is another relation of Lagrange-based hypergeometric Bernoulli polynomials with Lagrange polynomials and hypergeometric Bernoulli numbers in terms of Stirling numbers of the second kind S ( n , k ) , for which its generating function is given by the following.
e z 1 k k ! = n = k S ( n , k ) z n n ! .
Theorem 2.
For a fixed m N and n 0 , we have the following.
B n [ m 1 , α 1 , , α r ] x | x 1 , , x r = j = 0 n s = 0 n j k = 0 s g j α 1 , , α r x 1 , , x r S ( n , k ) x ( k ) s ! B n s j [ m 1 ] ( n s j ) ! .
Proof. 
Using (1), (4), and the Abel binomial theorem, we obtain the following.
n = 0 B n [ m 1 , α 1 , , α r ] x | x 1 , , x r z n = j = 1 r 1 x j z α j n = 0 B n [ m 1 ] z n n ! e z 1 + 1 x = j = 1 r 1 x j z α j n = 0 B n [ m 1 ] z n n ! n = 0 x n n ! n = k S ( n , k ) z n n ! = j = 1 r 1 x j z α j n = 0 B n [ m 1 ] z n n ! n = 0 k = 0 n S ( n , k ) x ( k ) z n n ! = n = 0 j = 0 n s = 0 n j k = 0 s g j α 1 , , α r x 1 , , x r S ( n , k ) x ( k ) s ! B n s j [ m 1 ] ( n s j ) !
Therefore, comparing the coefficients on both sides, we obtain (21). □
When α 1 = = α r = 0 , expression (21) reduces to a relation of hypergeometric Bernoulli polynomials with their numbers in terms of Stirling numbers of the second kind (see, e.g., (Proposition 5)) in [10].
The results in [13,17] allow us to obtain a matrix form of B s [ m 1 , α 1 , , α r ] x | x 1 , , x r , s = 0 , 1 , , n , as follows.
Part (a) of Theorem 1 yields the following:
B s [ m 1 , α 1 , , α r ] x | x 1 , , x r = G r [ m 1 , α 1 , , α r ] ( x 1 , , x r ) B [ m 1 ] ( x ) ,
where the following is the case:
G s [ m 1 , α 1 , , α r ] ( x 1 , , x r ) = g s [ m 1 , α 1 , , α r ] ( x 1 , , x r ) g 0 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) 0 ,
and the null entries of matrix G s [ m 1 , α 1 , , α r ] ( x 1 , , x r ) appear ( n r ) -times, and the matrix B [ m 1 ] ( x ) is given by B [ m 1 ] ( x ) = B 0 [ m 1 ] ( x ) B 1 [ m 1 ] ( x ) B r [ m 1 ] ( x ) r ! B n [ m 1 ] ( x ) n ! T .
Then, by (22), the matrix of the following:
B [ m 1 , α 1 , , α r ] x | x 1 , , x r = B 0 [ m 1 , α 1 , , α r ] x | x 1 , , x r B n [ m 1 , α 1 , , α r ] x | x 1 , , x r T ,
can be expressed as follows:
B [ m 1 , α 1 , , α r ] x | x 1 , , x r = G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) B [ m 1 ] ( x ) ,
where G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) is the following ( n + 1 ) × ( n + 1 ) matrix.
G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) = g 0 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) 0 0 g 1 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) g 0 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) 0 g 2 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) g 1 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) 0 g n [ m 1 , α 1 , , α r ] ( x 1 , , x r ) g n 1 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) g 0 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) .
The following theorem summarizes the ideas described above.
Theorem 3.
For a fixed m N , let B n [ m 1 , α 1 , , α r ] x | x 1 , , x r n 0 be the sequence of Lagrange-based hypergeometric Bernoulli polynomials in variables x , x 1 , , x r and parameters α 1 , , α r C . Then, matrix B [ m 1 , α 1 , , α r ] x | x 1 , , x r has the following matrix form.
B [ m 1 , α 1 , , α r ] x | x 1 , , x r = G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) B [ m 1 ] ( x ) .
Remark 1.
Note that according to (22), the rows of matrix G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) are precisely the matrices G s [ m 1 , α 1 , , α r ] ( x 1 , , x r ) for s = 0 , , n . Furthermore, matrix G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) is an ( n + 1 ) × ( n + 1 ) lower triangular matrix for each x 1 , , x r R \ { 0 } such that the following is the case.
det G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) = g 0 [ m 1 , α 1 , , α r ] ( x 1 , , x r ) n + 1 = 1 .
Therefore, G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) is an invertible matrix for each x 1 , , x r R \ { 0 } .
Remark 2.
Using (Equation (8)) in [13], we can deduce the following:
B [ m 1 ] ( x ) = M [ m 1 ] T ( x ) = B 0 [ m 1 ] n ! 0 0 0 0 B 1 [ m 1 ] ( n 1 ) ! B 0 [ m 1 ] n ! 0 0 0 B 2 [ m 1 ] 2 ! ( n 2 ) ! B 1 [ m 1 ] ( n 1 ) ! B 0 [ m 1 ] n ! 0 0 B 3 [ m 1 ] 3 ! ( n 3 ) ! B 2 [ m 1 ] 2 ! ( n 2 ) ! B 1 [ m 1 ] ( n 1 ) ! B 0 [ m 1 ] n ! 0 B n [ m 1 ] n ! B n 1 [ m 1 ] ( n 1 ) ! B n 2 [ m 1 ] 2 ! ( n 2 ) ! B n 3 [ m 1 ] 3 ! ( n 3 ) ! B 0 [ m 1 ] n ! T ( x ) ,
where T ( x ) = 1 x x r x n T . Again, matrix M [ m 1 ] is an ( n + 1 ) × ( n + 1 ) has a lower triangular matrix satisfying the following.
det M [ m 1 ] = B 0 [ m 1 ] n ! n + 1 = 1 ( n ! ) n + 1 .
Hence, M [ m 1 ] is an invertible matrix.
Finally, on the account of Theorem 3 and Remark 2, we can deduce the following result.
Theorem 4.
For fixed m N and x 1 , , x r R \ { 0 } , let B n [ m 1 , α 1 , , α r ] x | x 1 , , x r n 0 be the sequence of Lagrange-based hypergeometric Bernoulli polynomials in the variable x, parameters α 1 , , α r C . Then, the following matrix-inversion formula holds:
T ( x ) = N [ m 1 ] Q [ m 1 ] B [ m 1 , α 1 , , α r ] x | x 1 , , x r ,
where Q [ m 1 ] and N [ m 1 ] are the inverse matrices of G [ m 1 , α 1 , , α r ] ( x 1 , , x r ) and M [ m 1 ] , respectively.

4. Concluding Remarks

The main goal of our research has been to introduce Lagrange-based hypergeometric Bernoulli polynomials and to investigate some algebraic and analytic properties of these polynomials. We derived summation formulas, differential relations, and integral formulas for them. In addition, an interesting matrix-inversion formula (cf. Theorems 3 and 4) and a generating relation involving the Stirling numbers of the second kind have been derived for these polynomials (see Theorem 2).
In our study, we have obtained some formulas for some classical special numbers and recovered some well-known identities in the literature (Theorems 1 and 2). We have used the techniques of the theory of generating functions, mainly variants of Bernoulli generating functions into our investigation of some new identities for some special numbers. However, there is a different approach for the study of generating functions based on the use of the theory of zeta functions, which provides a new description of special polynomials and special numbers in terms of special values of certain zeta functions such as the Riemann zeta function (cf. [23,24,25] where this approach is adopted). To the best of our knowledge, a unified approach to the study of special numbers remains a work in progress (cf., e.g., [21] or more recently, [22], and references thereof).
We would also like to mention that if we consider the following broad class of generating functions:
G ( x 1 , x 2 , , x m , z ) : = θ ( z ) exp k = 1 m x k z k ,
where θ ( z ) is a function having an explicit Laurent expansion near z = 0 , then (4) becomes a special case of (25).
Hence, some generalizations of this research might be considered probably in the future. Finally, the results of this article might potentially be used in mathematics, in mathematical physics, and/or in engineering.

Author Contributions

Conceptualization, Y.Q. and W.A.K.; methodology, S.A., Y.Q., A.I. and W.A.K.; validation, Y.Q. and W.A.K.; formal analysis, S.A., Y.Q., A.I. and W.A.K.; investigation, S.A., Y.Q., A.I. and W.A.K.; writing—original draft preparation, Y.Q., A.I. and W.A.K.; writing—review and editing, S.A., Y.Q., A.I. and W.A.K.; visualization, S.A., Y.Q., A.I. and W.A.K.; supervision, Y.Q., A.I. and W.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical and privacy issues.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chan, W.-C.C.; Chyan, C.-J.; Srivastava, H.M. The Lagrange polynomials in several variables. Integral Transform. Spec. Funct. 2001, 12, 139–148. [Google Scholar] [CrossRef]
  2. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F. Higher Transcendental Functions; McGraw Hill: New York, NY, USA, 1953; Volume 3, p. 267. [Google Scholar]
  3. Mathai, A.M.; Haubold, H.J. Special Functions for Applied Scientists; Springer Science Business Media: New York, NY, USA, 2008; pp. 1–309. [Google Scholar]
  4. Szegö, G. Orthogonal Polynomials; American Math. Soc: Providence, RI, USA, 1939; pp. 68–69. [Google Scholar]
  5. Srivastava, H.M.; Manocha, H.L. A Treatise on Generating Functions; Ellis Horwood Ltd.: West Sussex, UK, 1984; pp. 378–384, 441–442. [Google Scholar]
  6. Bretti, G.; Ricci, P.E. Multidimensional extensions of the Bernoulli and Appell polynomials. Taiwanese J. Math. 2004, 8, 415–428. [Google Scholar] [CrossRef]
  7. Bretti, G.; Natalini, P.; Ricci, P.E. Generalizations of the Bernoulli and Appell polynomials. Abstr. Appl. Anal. 2004, 7, 613–623. [Google Scholar] [CrossRef] [Green Version]
  8. Natalini, P.; Bernardini, A. A generalization of the Bernoulli polynomials. J. Appl. Math. 2003, 2003, 155–163. [Google Scholar] [CrossRef] [Green Version]
  9. Srivastava, H.M.; Özarslan, M.A.; Kaanoğlu, C. Some generalized Lagrange-based Apostol-Bernoulli, Apostol-Euler and Apostol-Genocchi polynomials. Russian J. Math. Phys. 2013, 10, 110–120. [Google Scholar] [CrossRef]
  10. Chakraborty, K.; Komatsu, T. Generalized hypergeometric Bernoulli numbers. Rev. R. Acad. Cienc. Exactas Fís. Nat. Ser. A Math. RACSAM 2021, 115, 1–14. [Google Scholar] [CrossRef]
  11. Hassen, A.; Nguyen, H.D. Hypergeometric Bernoulli polynomials and Appell sequences. Int. J. Number Theory 2008, 4, 767–774. [Google Scholar] [CrossRef] [Green Version]
  12. Howard, F.T. Some sequences of rational numbers related to the exponential function. Duke Math. J. 1967, 34, 701–716. [Google Scholar] [CrossRef]
  13. Quintana, Y.; Ramírez, W.; Urieles, A. On an operational matrix method based on generalized Bernoulli polynomials of level m. Calcolo 2018, 55, 1–29. [Google Scholar] [CrossRef]
  14. Khan, W.A.; Abouzaid, M.S.; Abusufian, A.H.; Nisar, K.S. Some new classes of generalized Lagrange-based Apostol type Hermite polynomials. J. Inequal. Spec. Funct. 2019, 10, 1–11. [Google Scholar]
  15. Khan, W.A. On generalized Lagrange-based Apostol-type and related polynomials. Kragujevac J. Math. 2022, 46, 865–882. [Google Scholar]
  16. Duran, U.; Acikgoz, M.; Esi, A.; Araci, S. The Lagrange polynomials in several variables. Appl. Math. Inf. Sci. 2018, 12, 227–231. [Google Scholar] [CrossRef]
  17. Quintana, Y. Generalized mixed type Bernoulli-Gegenbauer polynomial. Kragujevac J. Math. 2020, 47, 245–257. [Google Scholar]
  18. Cesarano, C.; Ramírez, W.; Khan, S. A new class of degenerate Apostol-type Hermite polynomials and applications. Dolomites Res. Notes Approx. 2022, 15, 1–10. [Google Scholar]
  19. Cesarano, C.; Parmentier, A. A note on Hermite-Bernoulli Polynomials. In Nonlocal and Fractional Operators; Beghin, L., Garrappa, R., Mainardi, F., Eds.; Springer Nature: Cham, Switzerland, 2021; Volume 26, pp. 101–109. [Google Scholar]
  20. Khan, W.A.; Araci, S.; Acikgoz, M.; Esi, A. Laguerre-based Hermite-Bernoulli polynomials associated with bilateral series. Tbilisi Math. J. 2018, 11, 111–121. [Google Scholar]
  21. Hernández-Llanos, P.; Quintana, Y.; Urieles, A. About extensions of generalized Apostol-type polynomials. Results Math. 2015, 68, 203–225. [Google Scholar] [CrossRef]
  22. Navas, L.; Ruiz, F.J.; Varona, J.L.; Bernardini, A. Existence and reduction of generalized Apostol-Bernoulli, Apostol-Euler and Apostol-Genocchi polynomials. Arch. Math. 2019, 55, 157–165. [Google Scholar]
  23. Bayad, A. Special values of Lerch zeta function and their Fourier expansions. Adv. Stud. Contemp. Math. 2011, 21, 1–4. [Google Scholar]
  24. Bayad, A.; Hajli, M. On the multidimensional zeta functions associated with theta functions, and the multidimensional Appell polynomials. Math. Methods Appl. Sci. 2020, 43, 2679–2694. [Google Scholar] [CrossRef]
  25. Hajli, M. On a formula for the regularized determinant of zeta functions with application to some Dirichlet series. Q. J. Math. 2020, 71, 843–865. [Google Scholar] [CrossRef]
Figure 1. Plot for the zeros of B 11 [ m 1 , α 1 , α 2 , α 3 , α 4 ] ( y ) (a) and plot for the zeros of B 50 [ m 1 , α 1 , , α r ] ( y ) (b), when m = 6 and α 1 = 1 / 4 , α 2 = 1 , α 3 = 2 , α 4 = 3 .
Figure 1. Plot for the zeros of B 11 [ m 1 , α 1 , α 2 , α 3 , α 4 ] ( y ) (a) and plot for the zeros of B 50 [ m 1 , α 1 , , α r ] ( y ) (b), when m = 6 and α 1 = 1 / 4 , α 2 = 1 , α 3 = 2 , α 4 = 3 .
Symmetry 14 01125 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Albosaily, S.; Quintana, Y.; Iqbal, A.; Khan, W.A. Lagrange-Based Hypergeometric Bernoulli Polynomials. Symmetry 2022, 14, 1125. https://doi.org/10.3390/sym14061125

AMA Style

Albosaily S, Quintana Y, Iqbal A, Khan WA. Lagrange-Based Hypergeometric Bernoulli Polynomials. Symmetry. 2022; 14(6):1125. https://doi.org/10.3390/sym14061125

Chicago/Turabian Style

Albosaily, Sahar, Yamilet Quintana, Azhar Iqbal, and Waseem A. Khan. 2022. "Lagrange-Based Hypergeometric Bernoulli Polynomials" Symmetry 14, no. 6: 1125. https://doi.org/10.3390/sym14061125

APA Style

Albosaily, S., Quintana, Y., Iqbal, A., & Khan, W. A. (2022). Lagrange-Based Hypergeometric Bernoulli Polynomials. Symmetry, 14(6), 1125. https://doi.org/10.3390/sym14061125

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop