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
in the variables
and complex parameters
(
) are defined by means of the following generating function:
where
and
for
. 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
provide a natural extension of the class of bivariate Lagrange polynomials:
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:
where
is
nth classical Jacobi polynomial given by (cf., (Equation (4.3.2)) in [
4]).
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]:
and define the Lagrange-based hypergeometric Bernoulli polynomials in variables
, and complex parameters
(
) as follows:
where
and
for
.
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
and
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:
then it is not difficult to see from (
4) that the following is the case.
and
Inspired by the recent articles [
13,
14,
15,
16,
17,
20] in which the authors introduce the
-Hermite and
-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
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
,
,
, and
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
denotes the linear space of polynomials with real coefficients and a degree less than or equal to
n. For
and
, we use notations
and
for the rising and falling factorials, respectively:
and the following.
Moreover, as usual, the numbers given by the following:
denote the hypergeometric Bernoulli numbers (or generalized Bernoulli numbers of level
). It is clear that if
in (
3) then we recover the definition of the classical Bernoulli polynomials
, and classical Bernoulli numbers, respectively, i.e.,
, and
, respectively, for all
.
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))).
Hence, using
, for
, we obtain the following equivalent identity.
Moreover, from (
4), it is clear that the following is the case.
Thus, if we take
and combine (
4), (
3), and (
6), we have the following:
for any
.
Moreover, from (
5) and (
7), we obtain that if
, then (
4) induces the following multivariate polynomials.
Thus, the substitution
, with
into (
8) yields the following univariate polynomials.
In this case,
, and from (
4), we can deduce the following (see also [
1] (Equation (36))).
Consequently, (
9) takes the following form:
whenever
and
.
On the other hand, for
, with
, the polynomials
can be described by means of the following generating function:
where
and
for
. Now, if we assume that
, then (
12) becomes the following:
and from this last equality, it is easily deducible that the following is the case.
Then, the following is the case.
Finally, from (
11), (
13), and (
14), we conclude that for each
, we have the following.
It is clear that for
, 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
when
. Indeed, the first ones are given by the following.
3. The Polynomials and Their Properties
Now, we can proceed to investigate some relevant properties of the Lagrange-based hypergeometric Bernoulli polynomials.
Theorem 1. For a fixed , let be the Lagrange-based hypergeometric Bernoulli polynomials in the variables , and parameters (). Then, the following statements hold:
- (a)
Summation formulas. For every , we have the following. In particular, the following obtains. - (b)
Differential relations (Appell-type polynomial sequences). For with
and any nonzero , we have the following. - (c)
Representation formulas. If at least an is nonzero, , then the Lagrange polynomials can be expressed in terms of multivariate polynomials (
8)
as follows. - (d)
Inversion formula. If , then the following is the case. - (e)
Integral formulas. For any nonzero , we have the following. In particular, we have the following.
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
is nonzero,
. By (
4), (
8), and direct calculations, we have the following.
Or equivalently, we also have the following.
Now, from (
1), we have the following.
Hence, comparing the coefficients of
on the right hand side of (
19) and (
20), we obtain (
17).
Finally, assume that
and take
. 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:
where
is
nth classical Jacobi polynomial.
Moreover, notice that the inversion formula (
18) immediately implies the following.
Proposition 1. If , then for a fixed and each , the set is a basis for .
With respect to the study of zeros of the
nth Lagrange-based hypergeometric Bernoulli polynomial
when
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
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
and
are shown for prescribed values of
m and
.
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
, for which its generating function is given by the following.
Theorem 2. For a fixed and , we have the following. Proof. Using (
1), (
4), and the Abel binomial theorem, we obtain the following.
Therefore, comparing the coefficients on both sides, we obtain (
21). □
When
, 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
,
, as follows.
Part (a) of Theorem 1 yields the following:
where the following is the case:
and the null entries of matrix
appear
-times, and the matrix
is given by
.
Then, by (
22), the matrix of the following:
can be expressed as follows:
where
is the following
matrix.
The following theorem summarizes the ideas described above.
Theorem 3. For a fixed , let be the sequence of Lagrange-based hypergeometric Bernoulli polynomials in variables and parameters . Then, matrix has the following matrix form. Remark 1. Note that according to (22), the rows of matrix are precisely the matrices for . Furthermore, matrix is an lower triangular matrix for each such that the following is the case. Therefore, is an invertible matrix for each .
Remark 2. Using (Equation (8)) in [13], we can deduce the following:where . Again, matrix is an has a lower triangular matrix satisfying the following. Hence, is an invertible matrix.
Finally, on the account of Theorem 3 and Remark 2, we can deduce the following result.
Theorem 4. For fixed and , let be the sequence of Lagrange-based hypergeometric Bernoulli polynomials in the variable x, parameters . Then, the following matrix-inversion formula holds:where and are the inverse matrices of and , 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:
where
is a function having an explicit Laurent expansion near
, 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.