[BACK]
 Computer Modeling in Engineering & Sciences

DOI: 10.32604/cmes.2022.017669

ARTICLE

Lacunary Generating Functions of Hybrid Type Polynomials in Viewpoint of Symbolic Approach

1Mathematics Section, Women' s College, Aligarh Muslim University, Aligarh, 202002, India
2Department of Mathematics, Aligarh Muslim University, Aligarh, 202002, India
3Department of Economics, Faculty of Economics, Administrative of Social Sciences, Hasan Kalyoncu University, Gaziantep, TR-27410, Turkey
*Corresponding Author: Serkan Araci. Email: serkan.araci@hku.edu.tr
Received: 29 May 2021; Accepted: 19 August 2021

Abstract: In this paper, we introduce mon-symbolic method to obtain the generating functions of the hybrid class of Hermite-associated Laguerre and its associated polynomials. We obtain the series definitions of these hybrid special polynomials. Also, we derive the double lacunary generating functions of the Hermite-Laguerre polynomials and the Hermite-Laguerre-Wright polynomials. Further, we find multiplicative and derivative operators for the Hermite-Laguerre-Wright polynomials which helps to find the symbolic differential equation of the Hermite-Laguerre-Wright polynomials. Some concluding remarks are also given.

Keywords: Hermite-laguerre polynomials; laguerre-wright polynomials; hermite-laguerre-wright polynomials; hermite-mittag-leffler functions

1  Introduction

Recently, it has been realized that the symbolic method of operational as well as umbral nature provides powerful tool for the study of special functions [1,2]. Babusci et al. [3] developed the formalism of the symbolic method of operational nature for obtaining the generating functions of the Laguerre polynomials. By using symbolic method, the properties of the Laguerre polynomials could have accordingly been reduced to those of a Newton binomial containing an operator treated, in all the manipulations, as an ordinary algebraic quantity. Then Dattoli et al. [4] exploited symbolic method of umbral nature for obtaining the generating functions of the Hermite polynomials.

The Hermite and Laguerre polynomials, being orthogonal polynomial sequences, arise in various fields of engineering, mathematics and physics. The Hermite polynomials have applications in signal processing, probability, combinatorics, numerical analysis, quantum mechanics, random matrix theory, etc. These polynomials are also an example of Appell sequence, which obeys the umbral calculus. The Laguerre polynomials appear in quantum mechanics of the Morse potential, solution of the Schro¨dinger equation for a one-electron atom, oscillator system in phase space and 3D isotropic harmonic oscillator (see [58]), etc. The study of Laguerre polynomials becomes more convenient when these are expressed in terms of Laguerre-Wright polynomials. The symbolic method has emerged as a powerful tool to solve the problems involving certain special functions related to different branches of engineering and sciences for example fluid mechanics, physics, engineering and dynamical processes, etc. Most of the special functions do not show certain characteristics by classical methods and hence the importance of symbolic and umbral techniques are realized. These methods are widely used to establish numerous properties such as exponential and ordinary lacunary generating functions, symbolic definitions, symbolic differential equations, symbolic multiplicative and derivative operators etc. of known and new special functions and to evaluate various types of integrals involving these functions. A remarkable application of symbolic method is the convolution of two special functions to introduce and study new special functions. Motivated by the importance and applications of the Hermite and Laguerre-Wright polynomials and the usefulness of symbolic method, in this paper, we introduce and study the Hermite-Laguerre-Wright polynomials using symbolic method.

It has been shown that the concept related to monomiality techniques of the classical and generalized polynomials can be exploited to derive certain properties of families of polynomials, including Hermite and Laguerre polynomials used in pure and applied mathematics see for example [9,10].

The concept and formalism behind the monomiality principle and the symbolic methods can be exploited to introduce certain generalized as well as hybrid special polynomials and functions and to simplify the derivation of the properties of known and newly introduced special functions. The key element of this work is the introduction of the mon-symbolic method, which is the combination of monomiality techniques and symbolic method.

We recall that the ordinary Laguerre polynomials Ln(x) are defined by means of the following series definition [11]:

Ln(x)=n!r=0n(x)rr!2(nr)!.(1)

Babusci et al. have proved that the symbolic method of defining special functions can be exploited for obtaining several properties of certain special functions, which cannot be easily established by other well known methods [3]. So this method filled a huge vacuum in the study of special functions.

Babusci defined a symbolic shift operator c^z as [3]:

c^zα:f(z)f(z+α),(2)

which satisfies the property c^zαc^zβ=c^zβc^zα=c^zα+β. In particular c^zβc^zβ=c^zβc^zβ=1.

Clearly, for f(z)=1Γ(1+z),

c^zα1Γ(1+z)|z=0=1Γ(1+α).(3)

In view of Eq. (1), the symbolic definition of Ln(x) is given as [3]:

Ln(x)=(1xc^z)n1Γ(1+z)|z=0.(4)

Interchanging the variables xy and yx in the definition of Laguerre-Wright polynomials [3] the symbolic definition of the 2-variable Laguerre-Wright polynomials (2VLWP) Λn(α,β)(x,y) are as follows [3]:

Λn(α,β)(x,y)=(yc^zβx)n1Γ(1+z)|z=α,(5)

which on simplifying, gives the following series definition for the 2VLWP Λn(α,β)(x,y) [3]:

Λn(α,β)(x,y)=n!r=0n(x)rynrr!(nr)!Γ(βr+α+1).(6)

The 2-variable associated Laguerre polynomials (2VALP) Ln(α)(x,y) are specified by means of the following series definition [12]:

Ln(α)(x,y)=r=0nΓ(1+α+n)(x)rynrr!(nr)!Γ(1+α+r),

which in view of Eq. (6), gives [3]

Ln(α)(x,y)=Γ(1+α+n)n!Λn(α)(x,y),(7)

where

Λn(α)(x,y):=Λn(α,1)(x,y).(8)

In view of Eqs. (5), (7) and (8), symbolic definition of the 2VALP Ln(α)(x,y) is as follows [3]:

Ln(α)(x,y)=Γ(1+α+n)n!(yc^zx)n1Γ(1+z)|z=α.(9)

For y =1, Eq. (5) gives the symbolic definition of the Laguerre-Wright polynomials (LWP) Λn(α,β)(x) in terms of operator c^z is as follows:

Λn(α,β)(x)=(1c^zβx)n1Γ(1+z)|z=α,(10)

simplifying Eq. (10), we get the following series definition of LWP Λn(α,β)(x):

Λn(α,β)(x)=n!r=0n(x)rr!(nr)!Γ(βr+α+1).(11)

Next, we recall that the associated Laguerre polynomials (ALP) Ln(α)(x) are defined by the following series definition [11]:

Ln(α)(x)=r=0nΓ(1+α+n)(x)rr!(nr)!Γ(1+α+r).(12)

In view of Eqs. (11) and (12), we have:

Ln(α)(x)=Γ(1+α+n)n!Λn(α)(x),(13)

where

Λn(α)(x):=Λn(α,1)(x),(14)

which in view of Eq. (10), gives:

Λn(α)(x)=(1c^zx)n1Γ(1+z)|z=α.(15)

In view of Eqs. (13) and (15), the symbolic definition of the ALP Ln(α)(x) is given as:

Ln(α)(x)=Γ(1+α+n)n!(1c^zx)n1Γ(1+z)|z=α.(16)

The symbolic definition of the Bessel-Wright function Wβ,α(x) is given as [3]:

W(β,α)(x)=ec^zβx1Γ(1+z)|z=α1,(17)

which on simplifying, gives the following series definition of Wβ,α(x) [13,14]:

W(β,α)(x)=r=0xrr!Γ(βr+α).(18)

The symbolic definition of the Mittag-Leffler function Eβ,~α(x) is given as [3]:

Eβ,α(x)=[11xc^zβ]1Γ(1+z)|z=α1,(19)

which on simplifying, gives the following series definition of Eβ,~α(x) [13,14]:

Eβ,α(x)=r=0xrΓ(βr+α).(20)

The study of Hermite polynomials help to solve the classical boundary-value problems in the parabolic regions, through the use of parabolic coordinates and in quantum mechanics as well as in other areas of sciences. We recall that the 2-variable Hermite Kampe de Feriet polynomials (2VHKdFP) Hn(x, y) are defined by means of the following generating function and series definition [15]:

n=0Hn(x,y)tnn!=exp(xt+yt2)(21)

and

Hn(x,y)=n!k=0[n2]xn2kykk!(n2k)!,(22)

respectively.

Next, we recall that the nth-order Bessel function Jn(x) is defined by the following series definition [16]:

Jn(x)=r=0(1)r(x2)n+2rr!Γ(1+n+r).(23)

In view of Eq. (23) for n = 0, the 0th-order Bessel function Jn(x) is defined as [17]:

J0(2x)=r=0(1)rxr(r!)2.(24)

According to the monomi ality principle proposed by Steffensen [18] and developed by Dattoli [10], a polynomial set {pn(x)}n=0 is called quasi-monomial if there exist two operators multiplicative operator M^ and derivative operator P^, respectively, such that [10]:

M^{pn(x)}=pn+1(x)(25)

and

P^{pn(x)}=npn1(x).(26)

The multiplicative operator M^ and derivative operator P^ satisfy the following commutation relation:

[P^,M^]=P^M^M^P^=1^,

thus, the operatorM^and P^ display a weyl group structure [10]. Several characteristics of polynomial pn(x) can be obtained by using the operators M^ and P^.

Some characteristics are as follows:

(i)   If M^ and P^ have differential realizations, then the polynomial pn(x) satisfy the following differential equation:

M^P^{pn(x)}=npn(x).(27)

(ii)   Assuming here and in the following p0(x) = 1, then pn(x) can be explicitly constructed as:

pn(x)=M^n{1}(28)

which gives the series definition for pn(x).

(iii)   Consequently the generating function of pn(x) can be obtained as:

exp(tM^){1}=n=0pn(x)tnn!(|t|<).

By induction, Eq. (25) gives:

M^rpn(x)=pn+r(x).(29)

We recall that the 2VHKdFP Hn(x, y) is quasi-monomial with respect to the following multiplicative and derivative operators [10]:

M^H=x+2yDx(30)

and

P^H=Dx,(31)

respectively.

In view of Eq. (22), it can be easily verified that [10]:

Hn(ax,a2y)=anHn(x,y).(32)

Dattoli et al. [4] used the transition of 2VHKdFP Hn(x, y) from monomiality to umbral interpretation. The newton binomial realizes the umbral image of 2VHKdFP Hn(x, y). Dattoli et al. [4] have proved that the operational method become a fairly powerful tool once complimented with a notation of umbral nature. The umbral approach to the 2VHKdFP Hn(x, y) is particularly useful for a straightforward derivation of the relevant properties.

Dattoli redefined the 2VHKdFP Hn(x, y) using umbral approach of symbolic method as [4]:

Hn(x,y)=(x+h^y)nϕ0,(33)

where h^y denotes umbra, which acts on the vacuum ϕ0 in the following manner [4]:

h^yrϕ0=yr2r!Γ(r2+1)|cos(rπ2)|,(34)

which for r =0, gives ϕ0 =1. Thus Eq. (33) can be rewritten as:

Hn(x,y)=(x+h^y)n{1}.(35)

In view of Eqs. (28) and (35), we get the following symbolic multiplicative operator of 2VHKdFP:

M^H=(x+h^y).(36)

The use of umbral formalism looks much promising to develop a new technique to study the theory of special polynomials and special functions as well. Hybrid special functions as well as polynomials and their applications has been recognized by Dattoli and his co-workers [2,9,12].

We recall that the 2-variable Hermite-Laguerre polynomials (2VHLP) HLn(x, y) are defined by means of the following series definition [9]:

HLn(x,y)=n!r=0n(1)rHr(x,y)(r!)2(nr)!(37)

and the Hermite-Bessel-Wright function (HBWF) HWβ,α(x, y) is defined by means of the following series definition [3]:

HW(β,α)(x,y)=r=0Hr(x,y)r!Γ(βr+α).(38)

Motivated by the work of Babusci and his co-authors on lacunary generating functions for Laguerre polynomials [3] and application of the Laguerre polynomials [12,19,20], in this paper, we introduce certain generating functions for the 2-variable Hermite-Laguerre polynomials and some new families of polynomials. In Section 2, we introduce the Hermite-associated Laguerre polynomials and obtain ordinary generating function via mon-symbolic approach. Further, we define some new families of polynomials such as Hermite-Laguerre-Wright polynomials and find their exponential and ordinary generating functions. In Section 3, we derive the double and the triple lacunary 2-variable Hermite-Laguerre polynomials and find double lacunary generating functions for the 2-variable Hermite-Laguerre polynomials and the Hermite-Laguerre-Wright polynomials. In Section 4, we proposed an idea to find multiplicative and derivative operators for Hermite-Laguerre-Wright polynomials which helps to find symbolic differential equation for the same polynomials.

2  Generating Functions of the Hermite-Associated Laguerre and Hermite-Laguerre-Wright Polynomials

The most interesting example to understand the flexibility and usefulness of the symbolic method is derivation of generating functions of special polynomials and special functions as well. In this paper, we consider two types of generating functions, namely exponential and ordinary.

In this section, we introduce the Hermite-associated Laguerre polynomials and the Hermite-Laguerre-Wright polynomials by using the mon-symbolic method. Also, we obtain certain generating functions of these polynomials.

In view of Eqs. (15) and (36), we introduce the 1-parameter Hermite-Laguerre-Wright polynomials (1PHLWP) HΛn(α)(x,y) by means of the following symbolic definition:

HΛn(α)(x,y):=Λn(α)(x+h^y)=(1c^zH)n1Γ(1+z)|z=α,(39)

where

c^zH:=c^z(x+h^y).(40)

Since, in view of Eqs. (35) and (40), we have:

c^zHα1Γ(1+z)|z=0=(c^z(x+h^y))α1Γ(1+z)|z=0=Hα(x,y)Γ(α+1)(41)

and

c^zαc^zHβ1Γ(1+z)|z=0=c^zα+βHβ(x,y)1Γ(1+z)|z=0.

Binomially expanding the right hand side of Eq. (39) and then using Eqs. (35) and (40), we have:

HΛn(α)(x,y)=r=0n(nr)(1)rc^zrHr(x,y)1Γ(1+z)|z=α.

Simplifying, we get the following series definition of the 1PHLWP HΛn(α)(x,y):

HΛn(α)(x,y)=n!r=0n(1)rHr(x,y)r!(nr)!Γ(r+α+1).(42)

Now, in view of Eq. (39), we define the Hermite-associated Laguerre polynomials (HALP) HLn(α)(x,y) by replacing x with (x+h^y) in Eq. (13) as:

HLn(α)(x,y)=Γ(1+α+n)n!HΛn(α)(x,y).(43)

Next, we define the Hermite-Laguerre-Wright polynomials (HLWP) HΛn(α,β)(x,y) by replacing x with (x+h^y) in Eq. (10) as:

HΛn(α,β)(x,y)=(1c^zβ(x+h^y))n1Γ(1+z)|z=α,(44)

which on using Eq. (35), gives the following series definition of HΛn(α,β)(x,y):

HΛn(α,β)(x,y)=n!r=0n(1)rHr(x,y)r!(nr)!Γ(βr+α+1).(45)

In view of Eqs. (42) and (45), it is clear that:

HΛn(α,1)(x,y)=HΛn(α)(x,y).(46)

In the recent years, Dattoli used the monomiality principle to find the generating functions for some special polynomials and special functions as well [910,12]. Also, Babusci et al. [3] used the concept of symbolic method to find the generating function of the Laguerre Polynomials. In this paper, we combine the symbolic method with the monomiality technique to obtain the generating functions of the HALP HLn(α)(x,y) and the HLWP HΛn(α,β)(x,y).

Now, we establish following result for the ordinary generating function of the HALP HLn(α)(x,y):

Theorem 2.1 The ordinary generating function for HALP HLn(α)(x,y) is given by:

n=0,HLn(α)(x,y)tn=1(1t)(α+1)ext1t+yt2(1t)2.(47)

Proof. From Eq. (43), we have:

n=0HLn(α)(x,y)tn=n=0Γ(1+α+n)n!HΛn(α)(x,y)tn,(48)

which on using Eqs. (44) and (46) in the right hand side, it gives:

n=0HLn(α)(x,y)tn=n=0Γ(1+α+n)n!(1c^z(x+h^y))n1Γ(1+z)|z=αtn.(49)

Using the following relation between Gamma function and pochhammer symbol:

(α)n=Γ(α+n)Γ(α)(50)

in Eq. (49), we find:

n=0HLn(α)(x,y)tn=Γ(1+α)n=0(α+1)nn!(ttc^z(x+h^y))n1Γ(1+z)|z=α,

which on using the following series expansion:

n=0(α)ntnn!=1(1t)α,(|t|<1)(51)

gives:

n=0HLn(α)(x,y)tn=Γ(1+α)1(1(ttc^z(x+h^y)))α+11Γ(1+z)|z=α.

Again, using Eqs. (50) and (51), we have:

n=0HLn(α)(x,y)tn=Γ(1+α)1(1t)α+1n=0Γ(α+1+n)Γ(α+1)(t1t)nc^zn(x+h^y)n1Γ(1+z)|z=α1n!.

Using Eq. (35) in the right hand side of the above equation, we have:

n=0HLn(α)(x,y)tn=1(1t)α+1n=0(t1t)nHn(x,y)n!,

which on using Eqs. (21) and (32) gives assertion (47).

Next, we establish following result for the exponential generating function of the HLWP HΛn(α,β)(x,y):

Theorem 2.2 The exponential generating function for HLWP HΛn(α,β)(x,y) is given by:

n=0tnn!HΛn(α,β)(x,y)=etHW(β,α+1)(tx,t2y).(52)

Proof. Using Eq. (44), we have:

n=0tnn!HΛn(α,β)(x,y)=n=0tnn!(1c^zβ(x+h^y))n1Γ(1+z)|z=α=e(1c^zβ(x+h^y))t1Γ(1+z)|z=α.(53)

From [t,c^zβ(x+h^y)t]=0 and Weyl decoupling identity given by [21]:

eA^+B^=eA^eB^ek2,k=[A^,B^](kC),(54)

We find:

n=0tnn!HΛn(α,β)(x,y)=etec^zβt(x+h^y)1Γ(1+z)|z=α.(55)

Now, expanding the second exponential in the right hand side of the above equation, Eq. (55) yields:

n=0tnn!HΛn(α,β)(x,y)=etr=0c^zβr(t)r(x+h^y)rr!1Γ(1+z)|z=α,

which on simplifying, gives:

n=0tnn!HΛn(α,β)(x,y)=etr=0(t)rHr(x,y)r!Γ(βr+α+1).(56)

Using Eq. (32) in Eq. (56), we have:

n=0tnn!HΛn(α,β)(x,y)=etr=0Hr(xt,yt2)r!Γ(βr+α+1),(57)

which in view of Eq. (38), yields assertion (52).

Since, in view of Eqs. (43) and (46), HΛn(α,1)(x,y)n! := HLn(α)(x,y)Γ(1+α+n).

Therefore, from Theorem 2.2, we get the following result.

Corollary 2.1 The ordinary generating function for the polynomial HLn(α)(x,y)Γ(1+α+n) is given by:

n=0tnHLn(α)(x,y)Γ(1+α+n)=etHW(1,α+1)(tx,t2y).(58)

Now, we proceed to obtain the ordinary generating function for the HLWP HΛn(α,β)(x,y). For obtaining the ordinary generating function for the HLWP HΛn(α,β)(x,y), we define the Hermite-Mittag-Leffler function (HMLF) HEβ, α(x, y) by replacing x with (x+h^y) in Eq. (19) as:

HEβ,α(x,y)=[11c^zβ(x+h^y)]1Γ(1+z)|z=α1,(59)

which on simplifying and then using Eq. (35), gives the following series definition for the HEβ, α(x, y):

HEβ,α(x,y)=r=0Hr(x,y)Γ(βr+α).(60)

Now, we establish following result for the ordinary generating function of the HLWP HΛn(α,β)(x,y):

Theorem 2.3 The ordinary generating function for HLWP HΛn(α,β)(x,y) is given by:

n=0tnHΛn(α,β)(x,y)=1(1t)HEβ,α+1(xt1t,yt2(1t)2).(61)

Proof. From Eq. (44), we have:

n=0tnHΛn(α,β)(x,y)=n=0tn(1c^zβ(x+h^y))n1Γ(1+z)|z=α,(62)

which on using Eq. (51) for α = 1, gives:

n=0tnHΛn(α,β)(x,y)=11t+c^zβt(x+h^y)1Γ(1+z)|z=α,=1(1t)[1+c^zβt(x+h^y)1t]1Γ(1+z)|z=α.(63)

In view of Eq. (51) for α = 1, we find:

n=0tnHΛn(α,β)(x,y)=1(1t)r=0c^zβr(t1t)r(x+h^y)r1Γ(1+z)|z=α.

Using Eqs. (32) and (35) in the above equation, we have:

n=0tnHΛn(α,β)(x,y)=1(1t)r=0Hr(xt1t,yt2(1t)2)Γ(βr+α+1),

which in view of Eq. (60), gives assertion (61).

Now, we list the following examples of the special polynomials introduced or discussed in this paper.

In the next section, we obtain the lacunary generating functions of the Hermite-Laguerre polynomials HLn(x, y) and the HLWP HΛn(α,β)(x,y).

3  Lacunary Generating Functions

Babusci et al. redefined the Laguerre polynomials of degrees 2n and 3n (nN) by means of the series definitions and obtained the generating functions for these polynomials by using the symbolic method [3]. The generating functions of the Laguerre polynomials of degrees 2n and 3n (nN) are named double and triple lacunary generating functions, respectively. Recently, Dattoli et al. [4] obtained the double lacunary generating function of the 2VHKdFP Hn(x, y).

In this section, we define the 2-variable Hermite-Laguerre polynomials HL2n(x, y) and HL3n(x, y) by using the symbolic method. Also, we obtain double lacunary generating functions for the 2VHLP and the HLWP.

In the case when α = 0 in Eq. (43), the associated Laguerre polynomials Ln(α)(x) reduce to the Laguerre polynomials Ln(x) as follows:

HLn(x,y)=HΛn(0)(x,y),

which in view of Eqs. (39) and (40), gives the following symbolic definition of the 2VHLP HLn(x,y):

HLn(x,y)=(1c^zH)n1Γ(1+z)|z=0.(64)

Now, we give the following theorem for the 2-variable Hermite-Laguerre polynomials HL2n(x, y) and HL3n(x, y).

Theorem 3.1 The series definition for HL2n(x, y) and HL3n(x, y) are given by:

HL2n(x,y)=r=0n(nr)(1)rHχn(r)(x,y)(65)

and

HL3n(x,y)=k,r=0n(nr)(nk)(1)k+rHχn(k+r)(x,y).(66)

Proof. In view of Eq. (64), we have:

HL2n(x,y)=(1c^zH)2n1Γ(1+z)|z=0,(67)

which on simplifying, gives:

HL2n(x,y)=(1c^zH)n(1c^zH)n1Γ(1+z)|z=0=r=0n(nr)(1)r(x+h^y)r(1c^zH)n1Γ(1+z)|z=r.

Using Eq. (39) in the right hand side of the above equation, we find:

HL2n(x,y)=r=0n(nr)(1)r(x+h^y)rHΛn(r)(x,y).(68)

Therefore, in view of Eqs. (35) and (42), we have:

(x+h^y)rHΛn(r)(x,y)=n!k=0n(1)kHr+k(x,y)k!(nk)!Γ(k+r+1).

If we denote xrΛn(r)(x) by the polynomial χn(r)(x), then by replacing x with (x+h^y), we find:

Hχn(r)(x,y)=(x+h^y)rHΛn(r)(x,y),(69)

which on using Eq. (39), gives the following symbolic definition of the Hχn(r)(x,y):

Hχn(r)(x,y)=(x+h^y)r(1c^zH)n1Γ(1+z)|z=r.(70)

Thus, using Eq. (69) in Eq. (68), we get the assertion (65).

By the same way, we can prove it for HL3n(x, y), given by (66).

Now, we proceed to obtain the double lacunary generating function of the 2VHLP HLn(x, y). For this, we denote the product of two 2VHKdFP Hn(x, y) and Hn(z, w), by Hn2(x,y;z,w) which are called Bi-Hermite polynomials.

Thus, in view of Eq. (33), we have:

Hn2(x,y;z,w):=Hn(x,y)Hn(z,w)Hn2(x,y;z,w)=[(x+h^y)(z+h^w)]nϕ0,yϕ0,w.(71)

Now, we define the Bi-Hermite-Bessel function H2J0(x, y; z, w) as:

H2J0(x,y;z,w)=J0(2(x+h^y)(z+h^w))ϕ0,yϕ0,w,

which on using Eq. (24), gives:

H2J0(x,y;z,w)=r=0(1)r[(x+h^y)(z+h^w)]r(r!)2ϕ0,yϕ0,w.

In view of Eq. (33), we get the following series definition of H2J0(x, y; z, w):

H2J0(x,y;z,w)=r=0(1)rHr(x,y)Hr(z,w)(r!)2.(72)

Now, we are in a position to state the following theorem:

Theorem 3.2 The double lacunary generating function for 2VHLP HLn(x, y) is given by:

n=0tnn!HL2n(x,y)=etH2J0(x,y;2t,t),(73)

where H2J0(x, y; 2t, t) denotes the Bi-Hermite-Bessel function.

Proof. In view of Eq. (67), we have:

n=0tnn!HL2n(x,y)=n=0tnn!(1c^z(x+h^y))2n1Γ(1+z)|z=0,

which becomes:

n=0tnn!HL2n(x,y)=et(1c^z(x+h^y))21Γ(1+z)|z=0.(74)

Simplifying the above equation, we find:

n=0tnn!HL2n(x,y)=ete2c^z(x+h^y)t+c^z2(x+h^y)2t1Γ(1+z)|z=0.

Using Eq. (21) in the above equation, we obtain:

n=0tnn!HL2n(x,y)=etr=0c^zrr!Hr(2(x+h^y)t,(x+h^y)2t)1Γ(1+z)|z=0,

which on using Eqs. (32) and (35), it gives:

n=0tnn!HL2n(x,y)=etr=0(1)rHr(x,y)Hr(2t,t)r!2.(75)

In view of Eqs. (72) and (75), we get assertion (73).

Similarly, to obtain the double lacunary generating function for the HLWP HΛn(α,β)(x,y), we define the Bi-Hermite-Wright function (BHWF) H2Wβ,α(x, y; z, w) as:

H2W(β,α)(x,y;z,w)=W(β,α)((x+h^y)(z+h^w))ϕ0,yϕ0,w.

Using Eq. (18) in the right hand side of the above equation, we have:

H2W(β,α)(x,y;z,w)=r=0[(x+h^y)(z+h^w)]rr!Γ(βr+α)ϕ0,yϕ0,w,

which on using Eq. (33), gives the following series definition of H2W(β,α)(x,y;z,w):

H2W(β,α)(x,y;z,w)=r=0Hr(x,y)Hr(z,w)r!Γ(βr+α).(76)

We, now state the following theorem:

Theorem 3.3 The double lacunary generating function for HLWP HΛn(α,β)(x,y) is given by

n=0tnn!HΛ2n(α,β)(x,y)=etH2W(β,α+1)(x,y;2t,t).(77)

Proof. Using Eq. (44), we have

n=0tnn!HΛ2n(α,β)(x,y)=n=0tnn!(1c^zβ(x+h^y))2n1Γ(1+z)|z=α,

which can be written as:

n=0tnn!HΛ2n(α,β)(x,y)=et(1c^zβ(x+h^y))21Γ(1+z)|z=α=ete((2t)c^zβ(x+h^y)+tc^z2β(x+h^y)2)1Γ(1+z)|z=α,

which on using Eq. (21), gives:

n=0tnn!HΛ2n(α,β)(x,y)=etr=0(c^zβ)rr!Hr(2t(x+h^y),t(x+h^y)2)1Γ(1+z)|z=α.(78)

In view of Eqs. (32) and (35), we have:

n=0tnn!HΛ2n(α,β)(x,y)=etr=0Hr(x,y)Hr(2t,t)r!Γ(βr+α+1).(79)

Using Eq. (76) in the right hand side of (79), we get assertion (77).

In the next section, we introduce the symbolic multiplicative and derivative operators and establish the symbolic differential equation for HLWP HΛn(α,β)(x,y) by using the monomiality method.

4  Monomiality Principle and Symbolic Differential Equation

In this section, we develop the theory of symbolic multiplicative and derivative operators for the HLWP HΛn(α,β)(x,y), which helps to find the symbolic differential equation for the HLWP HΛn(α,β)(x,y).

We obtain the following symbolic-differential recurrence relation for HLWP HΛn(α,β)(x,y):

Theorem 4.1 The HLWP HΛn(α,β)(x,y) satisfies the following recurrence relation:

HΛn1(α,β)(x,y)xc^zHβΛn1(α,β)(x,y)2yDxc^zHβΛn1(α,β)(x,y)=HΛn(α,β)(x,y).(80)

Proof. Using Eqs. (30) and (36) in Eq. (44), we have:

HΛn1(α,β)(x,y)=(1c^zβ(x+2yDx))n11Γ(1+z)|z=α.(81)

Operating (1c^zβ(x+2yDx)) on both sides of Eq. (81), and then using Eqs. (36) and (44), we obtain:

(1c^zβ(x+2yDx))HΛn1(α,β)(x,y)=HΛn(α,β)(x,y),(82)

which on simplifying, gives assertion (80).

Now, we obtain the symbolic multiplicative and derivative operators for the HLWP HΛn(α,β)(x,y).

Theorem 4.2 The HLWP HΛn(α,β)(x,y) is quasi-monomial with respect to the following symbolic multiplicative and derivative operators:

MΛ=(1c^zβ(x+2yDx))(83a)

and

PΛ=c^zβDx,(83b)

respectively.

Proof. Using Eqs. (30) and (36) in Eq. (44), we have:

HΛn(α,β)(x,y)=(1c^zβ(x+2yDx))n1Γ(1+z)|z=α,(84)

Now, operating (1c^zβ(x+2yDx)) in the both sides of Eq. (84), we have:

(1c^zβ(x+2yDx))HΛn(α,β)(x,y)=(1c^zβ(x+2yDx))n+11Γ(1+z)|z=α,(85)

which on using (84), gives:

(1c^zβ(x+2yDx))HΛn(α,β)(x,y)=HΛn+1(α,β)(x,y).(86)

In view of Eqs. (25) and (86), we obtain assertion (83a).

Now, differentiating Eq. (45) partially with respect to x, we have:

DxHΛn(α,β)(x,y)=n!r=1n(1)rHr1(x,y)(r1)!(nr)!Γ(βr+α+1).(87)

Replacing r by r + 1 and operating c^zβ on both sides of the above equation, we have:

c^zβDxHΛn(α,β)(x,y)=n(n1)!r=0n1(1)rHr(x,y)r!(n1r)!Γ(βr+α+1).(88)

Using Eq. (45), we get:

c^zβDxHΛn(α,β)(x,y)=nHΛn1(α,β)(x,y),(89)

which in view of Eq. (26), we get assertion (83b).

Further, we establish the following theorem for the HLWP HΛn(α,β)(x,y):

Theorem 4.3 The HLWP HΛn(α,β)(x,y) satisfies the following symbolic differential equation:

2yDxH2Λn(α,β)(x,y)+(xc^zβ)DxHΛn(α,β)(x,y)nHΛn(α,β)(x,y)=0.(90)

Proof. Using multiplicative and derivative operators given by Eqs. (83) and (84) in Eq. (27), we find:

(1c^zβ(x+2yDx))c^zβDx{HΛn(α,β)(x,y)}=nHΛn(α,β)(x,y),(91)

which on simplification, gives assertion (90).

5  Concluding Remarks

It has been realized that the symbolic method serves as a useful tool to introduce several special polynomials and their lacunary forms. In this paper, we used the symbolic method to find the generating functions of different polynomials. In this section, we define the 2-parameter, 2-variable Laguerre polynomials (2P2VLP) by using symbolic approach.

We introduce the 2-parameter, 2-variable Laguerre polynomials (2P2VLP) Ln(α,β)(x,y) by means of the symbolic definition.

In view of Eqs. (13) and (14), we define the 2P2VLP Ln(α,β)(x,y) as:

Ln(α,β)(x,y)=Γ(1+α+n)n!Λn(α,β)(x,y),(92)

which on using Eq. (5), gives the following symbolic definition of the 2P2VLP Ln(α,β)(x,y):

Ln(α,β)(x,y)=Γ(1+α+n)n!(yc^zβx)n1Γ(1+z)|z=α.(93)

Using Eq. (11) in Eq. (93), we get the following series definition of the 2P2VLP Ln(α,β)(x,y):

Ln(α,β)(x,y)=Γ(1+α+n)r=0n(x)rr!(nr)!Γ(βr+1+α).(94)

Finding generating functions of this polynomial is an open problem for further research.

The results established in this paper can also be obtained by using the umbra of 2VHKDFP Hn(x, y) but the method will not be purely symbolic. We conclude this paper with the fact that the symbolic method makes study of special functions easier than classical techniques.

We are now in a position to conclude our paper by investigating the special cases for our main equations which was given below as Table 1.

Acknowledgement: The authors wish to express their appreciation to the reviewers for their helpful suggestions which greatly improved the presentation of this paper.

Funding Statement: The authors received no specific funding for this study.

Conflicts of Interest: The authors declare that there are no conflicts of interest regarding the publication of this paper.

## References

1.  Babusci, D., Dattoli, G., Duchamp, G. H. E., Górska, K., Penson, K. A. (2012). Definite integrals and operational methods. Applied Mathematics and Computation, 219(6), 3017–3021. DOI 10.1016/j.amc.2012.09.029. [Google Scholar] [CrossRef]

2.  Dattoli, G. (2000). Advance special functions and integration methods. Proceedings of the Workshop, Melfi (PZItaly. [Google Scholar]

3.  Babusci, D., Dattoli, G., Górska, K., Penson, K. (2017). Lacunary generating functions for the laguerre polynomials. Séminaire Lotharingien de Combinatoire, 76, B76b. [Google Scholar]

4.  Dattoli, G., Germano, B., Martinelli, M., Ricci, P. (2015). Lacunary generating functions of hermite polynomials and symbolic methods. Journal of Mathematics, 4(1), 16–23. [Google Scholar]

5.  Arfken, G. W. H. (2000). mathematical methods for physicists. San Diego: Academic Press. [Google Scholar]

6.  Bhrawy, A. H., Alghamdi, M. A. (2013). The operational matrix of caputo fractional derivatives of modified generalized laguerre polynomials and its applications. Advances in Difference Equations, 2013(1), 1–19. DOI 10.1186/1687-1847-2013-307. [Google Scholar] [CrossRef]

7.  Karaseva, I. A. (2011). Fast calculation of signal delay in rc-circuits based on laguerre functions. Russian Journal of Numerical Analysis and Mathematical Modelling, 26(3), 295–301. DOI 10.1515/rjnamm.2011.016. [Google Scholar] [CrossRef]

8.  Karaseva, I. A. (2020). Laguerre-type exponentials, laguerre derivatives and applications. A survey. Mathematics, 8(11), 2054. DOI 10.3390/math8112054. [Google Scholar] [CrossRef]

9.  Dattoli, G., Lorenzutta, S., Cesarano, C. (2001). Generalized polynomials and new families of generating functions. Annali Dellâ Universita di Ferrara, 47(1), 57–61. [Google Scholar]

10. Dattoli, G. (2000). Hermite-bessel and laguerre-bessel functions: A by-product of the monomiality principle. Advanced Special Functions and Applications, 1171(1), 147–164. [Google Scholar]

11. Andrews, L. C. (1985). Special functions for applied mathematics and engineering. New York: MacMillan. [Google Scholar]

12. Dattoli, G. (2000). Generalized polynomials, operational identities and their applications. higher transcendental functions and their applications. Journal of Computational and Applied Mathematics, 118(1--2), 111–123. DOI 10.1016/S0377-0427(00)00283-1. [Google Scholar] [CrossRef]

13. Babusci, D., Dattoli, G., Górska, K., Penson, K. (2013). Symbolic methods for the evaluation of sum rules of bessel function. Journal of Mathematical Physics, 54(7), 73501. DOI 10.1063/1.4812325. [Google Scholar] [CrossRef]

14. Podlubny, I. (1999). Fractional differential equations. San Diego: Academic press. [Google Scholar]

15. Appell, P., de Fériet, J. K. (1926). Functions hypergéométriques et hypersphériques. Polynomes d'Hermite. Paris: Gauthier-Villars. [Google Scholar]

16. Rainville, E. D. (1960). Special functions, vol. 5. New York: The Macmillan Company. [Google Scholar]

17. Babusci, D., Dattoli, G., Licciardi, S., Sabia, E. (2019). Mathematical methods for physicists. Singapore: World Scientific. [Google Scholar]

18. Steffensen, J. F. (1941). The poweroid an extension of the mathematical notion of power. Acta Mathematica, 73, 333–366. DOI 10.1007/BF02392231. [Google Scholar] [CrossRef]

19. Dattoli, G., Torre, A. (1998). Operatorial methods and two variable laguerre polynomials. Atti della Accademia delle Scienze di Torino. Classe di Scienze Fisiche, Matematiche e Naturali, 132, 408. [Google Scholar]

20. Tricomi, F. G. (1959). Funzioni speciali. Gheroni. [Google Scholar]

21. Dattoli, G., Ottaviani, P., Torre, A., Vázquez, L. (1997). Evolution operator equations: Integration with algebraic and finitedifference methods. applications to physical problems in classical and quantum mechanics and quantum field theory. La Rivista del Nuovo Cimento (1978–1999), 20(2), 3. DOI 10.1007/BF02907529. [Google Scholar] [CrossRef]