Differential Equations Associated with Two Variable Degenerate Hermite Polynomials

: In this paper, we introduce the two variable degenerate Hermite polynomials and obtain some new symmetric identities for two variable degenerate Hermite polynomials. In order to give explicit identities for two variable degenerate Hermite polynomials, differential equations arising from the generating functions of degenerate Hermite polynomials are studied. Finally, we investigate the structure and symmetry of the zeros of the two variable degenerate Hermite equations.


Introduction
The classical Hermite numbers H n and Hermite polynomials H n (x) are usually defined by the generating functions: Clearly, H n = H n (0). It can be seen that these numbers and polynomials play an important role in various areas of mathematics and physics, including numerical theory, combinations, special functions, and differential equations. Many interesting properties about them have been explored (see [1][2][3][4][5][6]). For example, in mathematics and physics, the Hermite polynomials are a classical orthogonal polynomial sequence. In probability, they appears as the Edgeworth series; in combinatorics, they arise in the umbral calculus as an example of an Appell sequence; in numerical analysis, they play a role in Gaussian quadrature; and in physics, they give rise to the eigenstates of the quantum harmonic oscillator. The polynomials H n (x) satisfy the Hermite differential equation: dH(x) dx + 2nH(x) = 0, n = 0, 1, 2, . . . .
The special polynomials of two variables provided a new means of analysis for the solution of large classes of partial differential equations often encountered in physical problems. Most of the special polynomials of mathematical physics and their generalization have been proposed by physical problems. As another application of the Hermite differential equation for H n (x, y), we recall that the two variable Hermite polynomials H n (x, y) defined by the generating function (see [2]): are the solution of heat equation: ∂ ∂y H n (x, y) = ∂ 2 ∂x 2 H n (x, y), H n (x, 0) = x n . (2) Observe that H n (2x, −1) = H n (x).
Motivated by their potential and importance for applications in certain problems in probability, combinatorics, number theory, differential equations, numerical analysis, and other fields of mathematics and physics, several kinds of some special polynomials and numbers were recently studied by many authors (see [1][2][3][4][5][6][7]). Many mathematicians have studied the area of the degenerate Bernoulli polynomials, degenerate Euler polynomials, and degenerate tangent polynomials (see [8][9][10][11]). One of the important aspect of the study of any degenerate polynomials is to find their definition. Recently, Haroon and Khan [12] proposed the degenerate Hermite-Bernoulli polynomials, which are formulated in terms of p-adic invariant integral on Z p : Mathematicians have studied the differential equations arising from the generating functions of special polynomials (see [13][14][15][16][17][18]). Based on the results so far, in the present work, a new class of degenerate Hermite polynomials is constructed. We can derive the differential equations generated from the generating function of two variable degenerate Hermite polynomials. By using the coefficients of this differential equation, we have explicit identities for the two variable degenerate Hermite polynomials.
We remember that the classical Stirling numbers of the first kind S 1 (n, k) and the second kind S 2 (n, k) are defined by the relations (see [8][9][10][11][12]): respectively. Here, (x) n = x(x − 1) · · · (x − n + 1) denotes the falling factorial polynomial of order n. We also have: The rest of the paper is organized as follows. In Section 2, we introduce the two variable degenerate Hermite polynomials and obtain the basic properties of these polynomials. In Section 3, we give some symmetric identities for two variable degenerate Hermite polynomials. In Section 4, we derive the differential equations generated from the generating function of two variable degenerate Hermite polynomials. Using the coefficients of this differential equation, we have explicit identities for the two variable degenerate Hermite polynomials. In Section 5, we investigate the zeros of the two variable degenerate Hermite equations by using a computer. Further, we observe the pattern of scattering phenomenon for the zeros of two variable degenerate Hermite equations. Our paper ends with Section 6, where the conclusions and future directions of this work are presented.

Basic Properties for the Two Variable Degenerate Hermite Polynomials
In this section, a new class of the two variable degenerate Hermite polynomials is considered. Further, some properties of these polynomials are also obtained.
We define the two variable degenerate Hermite polynomials H n (x, y, λ) by means of the generating function: Since (1 + λ) t λ → e t as λ → 0, it is evident that Equation (4) reduces to Equation (1). Observe that Khan's degenerate Hermite-Bernoulli polynomials H B n (x, y|λ) vary from the two variable degenerate Hermite polynomials H n (x, y, λ). Another application of the differential equation for H n (x, y, λ) is as follows: Note that: Substitute the series in Equation (4) for F(t, x, y, λ) to get: Thus, the two variable degenerate Hermite polynomials H n (x, y, λ) in the generating function (4) are the solution of equation: Since λ log(1 + λ) → 1 as λ approaches zero, it is apparent that Equation (5) descends to Equation (2). By Equation (3), we also need the binomial theorem: for a variable x, The generating function (4) is useful for deriving several properties of the two variable degenerate Hermite polynomials H n (x, y, λ). For example, we have the following expression for these polynomials: Theorem 1. For any positive integer n, we have: where [ ] denotes taking the integer part.
Proof. By Equations (4) and (6), we have: On comparing the coefficients of t n n! , the expected result of Theorem 1 is achieved.
When λ → 0, then we have: The following basic properties of the two variable degenerate Hermite polynomials H n (x, y, λ) are derived from (4). We, therefore, chose to omit the details involved.

Theorem 2.
For any positive integer n, we have:

Symmetric Identities for the Two Variable Degenerate Hermite Polynomials
In this section, we give some new symmetric identities for the two variable degenerate Hermite polynomials. We also get some explicit formulas and properties for the two variable degenerate Hermite polynomials.
Let a, b > 0 and a = b. We start with: Then, the expression for F (t, λ) is symmetric in a and b: On similar lines, we can obtain that: Comparing the coefficients of t m m! in the last two equations, we have the following theorem.
Theorem 3. Let a, b > 0 and a = b. The following identity holds true: Again, we now use: For λ ∈ C, we define the degenerate Bernoulli polynomials given by the generating function: When x = 0 and B n (λ) = B n (0, λ) are called the degenerate Bernoulli numbers. The first few follow immediately from this generating function, For each integer k ≥ 0, T k (n) = 0 k + 1 k + 2 k + · · · + (n − 1) k is called the sum of integers. A generalized falling factorial sum σ k (n, λ) can be defined by the generation function: Note that lim λ→0 σ k (n, λ) = T k (n). From G(t, λ), we get the following result: In a similar fashion, we have: By comparing the coefficients of t m m! on the right-hand sides of the last two equations, we have the theorem below.
Theorem 4. Let a, b > 0 and a = b. The the following identity holds true: Making the N-times derivative for Equation (4) with respect to t, we have: By multiplying the exponential series e xt = ∑ ∞ m=0 x m t m m! in both sides of Equation (7) and the Cauchy product, we get: By Equation (8) and the Leibniz rule, we have: Hence, by Equations (8) and (9), comparing the coefficients of t m m! gives the below theorem.

Differential Equations Associated with Two Variable Degenerate Hermite Polynomials
In this section, we study the differential equations with coefficients a i (N, x, y, λ) arising from the generating functions of the two variable degenerate Hermite polynomials: By using the coefficients of this differential equation, we can derive explicit identities for the two variable degenerate Hermite polynomials H n (x, y, λ). Recall that: Then, by Equation (11), we have: 4xy tF(t, x, y, λ) Continuing this process as shown in Equation (13), we can guess that: Differentiating Equation (14) with respect to t, we have: Now, replacing N by N + 1 in Equation (14), we find: Comparing the coefficients on both sides of Equations (15) and (16), we obtain: and: In addition, by Equation (14), we have: By Equation (19), we get: a 0 (0, x, y, λ) = 1.

(20)
It is not difficult to show that: x, y, λ).
For i = 1 in Equation (18), we have: Continuing this process, we can deduce that, for 1 ≤ i ≤ N − 1,

Theorem 7.
For N = 0, 1, 2, . . . , the differential equation: where: Here is a plot of the surface for this solution.
Making the N-times derivative for Equation (4) with respect to t, we have: By Equation (28) and Theorem 7, we have: (a) (b) If we take m = 0 in Equation (29), then we have the corollary below.

Zeros of the Two Variable Degenerate Hermite Polynomials
This section shows the benefits of supporting theoretical prediction through numerical experiments and finding a new interesting pattern of the zeros of the two variable degenerate Hermite equations H n (x, y, λ) = 0. By using a computer, the two variable degenerate Hermite polynomials H n (x, y, λ) can be determined explicitly. We investigated the zeros of the two variable degenerate Hermite equations H n (x, y, λ) = 0. The zeros of the H n (x, y, λ) = 0 for n = 30, y = 2, −2, 2 + i, −2 − i, λ = 1/2, and x ∈ C are displayed in Figure 2.
In Figure 3a, we chose y = 2. In Figure 3b, we chose y = −2. In Figure 3c, we chose y = 2 + i. In Figure 3d, we chose y = −2 − i. Our numerical results for approximate solutions of real zeros of the two variable degenerate Hermite equations H n (x, y, λ) = 0 are displayed (Tables 1 and 2). We observed a remarkable regular structure of the complex roots of the two variable degenerate Hermite equations H n (x, y, λ) = 0 and also hoped to verify the same kind of regular structure of the complex roots of the two variable degenerate Hermite equations H n (x, y, λ) = 0 ( Table 1).
The plot of the real zeros of the two variable degenerate Hermite equations H n (x, y, λ) = 0 for the 1 ≤ n ≤ 30, λ = 1/2 structure is presented in Figure 4. In Figure 4a, we chose y = 2. In Figure 4b, we chose y = −2. In Figure 4c, we chose y = 2 + i. In Figure 4d, we chose y = −2 − i.

Conclusions and Future Directions
In this article, we defined the two variable degenerate Hermite polynomials and obtained some new symmetric identities for two variable degenerate Hermite polynomials. We constructed differential equations arising from the generating function of the two variable degenerate Hermite polynomials H n (x, y, λ). We also investigated the symmetry of the zeros of the two variable degenerate Hermite equations H n (x, y, λ) = 0 for various variables x and y. As a result, we found that the distribution of the zeros of two variable degenerate Hermite equations H n (x, y, λ) = 0 is a very regular pattern. Therefore, we made the following series of conjectures with numerical experiments: Let us use the following notations. R H n (x,y,λ) denotes the number of real zeros of H n (x, y, λ) = 0 lying on the real plane Im(x) = 0, and C H n (x,y,λ) denotes the number of complex zeros of H n (x, y, λ) = 0. Since n is the degree of the polynomial H n (x, y, λ), we have R H n (x,y,λ) = n − C H n (x,y,λ) .
We can see a good regular pattern of the complex roots of the two variable degenerate Hermite equations H n (x, y, λ) = 0 for y and λ. Therefore, the following conjecture is possible.
As a result of investigating more y and λ variables, it is still unknown whether the Conjectures 1 and 2 are true or false for all variables y and λ.
We observed that solutions of the two variable degenerate Hermite equations H n (x, y, λ) = 0 have no Re(x) = a reflection symmetry for a ∈ R. It is expected that solutions of the two variable degenerate Hermite equations H n (x, y, λ) = 0 have Re(x) = 0 reflection symmetry (see Figures 2-4).
Finally, we considered the more general problems. How many zeros does H n (x, y, λ) have? We were not able to decide if H n (x, y, λ) = 0 had n distinct solutions. We would like to know the number of complex zeros C H n (x,y,λ) of H n (x, y, λ) = 0.
As a result of investigating more n variables, it is still unknown whether the conjecture is true or false for all variables n (see Tables 1 and 2).
We expect that research in these directions will make a new approach using the numerical method related to the research of the two variable degenerate Hermite equations H n (x, y, λ) = 0, which appear in applied mathematics and mathematical physics.

Author Contributions:
We are all equally contributed to write this paper. All authors have read and agreed to the published version of the manuscript.
Funding: This work was supported by the Dong-A university research fund.

Conflicts of Interest:
The authors declare no conflict of interest.