Approximate Roots and Properties of Differential Equations for Degenerate q-Special Polynomials

Abstract

In this paper, we generate new degenerate quantum Euler polynomials (DQE polynomials), which are related to both degenerate Euler polynomials and q-Euler polynomials. We obtain several $(q,h)$-differential equations for DQE polynomials and find some relations of q-differential and h-differential equations. By varying the values of $q,\eta$, and h, we observe the values of DQE numbers and approximate roots of DQE polynomials to obtain some properties and conjectures.

1. Basic Concepts and Introduction

Before clarifying the objectives of this paper, we first introduce the necessary basic concepts. We identify several definitions and properties and present the goals of this paper based on them.

Let $\eta ,q\in \mathbb{R}$ with $q\ne 1$. The quantum number, q-number, discovered by Jackson is

$${\left[\eta \right]}_q=\frac{1-q^{\eta }}{1-q},$$

and we note that ${lim}_{q\to 1}{\left[\eta \right]}_q=\eta$. In particular, for $k\in \mathbb{Z}$, ${\left[k\right]}_q$ is called q-integer; see [1,2,3,4].

Many mathematicians in various fields have worked on the introduction of q-numbers, such as q-discrete distribution, q-differential equations, q-series, q-calculus, and so on; see [5,6,7,8].

The following equation,

$${\left[\begin{array}{c}m\\ r\end{array}\right]}_q=\frac{{\left[m\right]}_q!}{{[m-r]}_q!{\left[r\right]}_q!},$$

is the q-Gaussian binomial coefficient where m and r are non-negative integers; see [4,5]. For $r=0$, the value of q-Gaussian binomial coefficients is 1 since the numerator and the denominator are both empty products. One notes that ${\left[\eta \right]}_q!={\left[\eta \right]}_q{[\eta -1]}_q\cdots {\left[2\right]}_q{\left[1\right]}_q$ and ${\left[0\right]}_q!=1$.

In [9], a two-parameter time scale ${\mathbf{T}}_{q,h}$ was introduced as follows:

$${\mathbf{T}}_{q,h}:=\{q^{\eta }\psi +{\left[\eta \right]}_{q}h|\psi \in \mathbb{R},\eta \in \mathbb{Z},h,q\in {\mathbb{R}}^{+},q\ne 1\}\cup \left\{\frac{h}{1-q}\right\}.$$

Definition 1 

([9,10]). Let $f:{\mathbf{T}}_{q,h}\to \mathbb{R}$ be any function. Then, the delta $(q,h)$-derivative of f $D_{q,h}\left(f\right)$ is defined by

$$D_{q,h}f\left(\psi \right):=\frac{f(q\psi +h)-f\left(\psi \right)}{(q-1)\psi +h}.$$

From the above definition, we can see several properties as follows:

(i)

For $\psi \in {\mathbf{T}}_{q,h}$, $D_{q,h}f\left(\psi \right)=0$ if and only if $f\left(\psi \right)$ is a constant;

(ii)

$D_{q,h}f\left(\psi \right)=D_{q,h}g\left(\psi \right)$ for all $\psi \in {\mathbf{T}}_{q,h}$ if and only if $f\left(\psi \right)=g\left(\psi \right)+c$ with some constant c;

(iii)

For $\psi \in {\mathbf{T}}_{q,h}$, $D_{q,h}f\left(\psi \right)=c_1$ if and only if $f\left(\psi \right)=c_1\psi +c_2$, where $c_1$ and $c_2$ are constant.

In Definition 1, we can see that $D_{q,h}\left(f\right)$, the delta $(q,h)$-derivative of f reduces to $D_q\left(f\right)$, the q-derivative of f for $h\to 0$ and reduces to $D_h\left(f\right)$, the h-derivative of f for $q\to 1$.

In addition, we can find the product rule and quotient rule for the delta $(q,h)$-derivative.

Theorem 1 

([9,10]). Let $f,g$ be arbitrary functions.

$$\begin{gathered} \mathit{\left(}i\mathit{\right)}\mathit{Product} \mathit{rule} \\ {D}_{q,h}\left(f\left(\psi \right)g\left(\psi \right)\right)=g(q\psi +h)D_{q,h}f\left(\psi \right)+f\left(\psi \right)D_{q,h}g\left(\psi \right) \\ =f(q\psi +h)D_{q,h}g\left(\psi \right)+g\left(\psi \right)D_{q,h}f\left(\psi \right). \\ \mathit{\left(}\mathit{ii}\mathit{\right)}\mathit{Quotient} \mathit{rule} \\ {D}_{q,h}\left(\frac{f\left(\psi \right)}{g\left(\psi \right)}\right)=\frac{g\left(\psi \right)D_{q,h}f\left(\psi \right)-f\left(\psi \right)D_{q,h}g\left(\psi \right)}{g\left(\psi \right)g(q\psi +h)} \\ =\frac{g(q\psi +h)D_{q,h}f\left(\psi \right)-f(q\psi +h)D_{q,h}g\left(\psi \right)}{g\left(\psi \right)g(q\psi +h)}. \end{gathered}$$

Definition 2 

([10,11]). The generalized quantum binomial ${(\psi -{\psi }_0)}_{q,h}^{\eta }$ is defined by

$${(\psi -{\psi }_0)}_{q,h}^{\eta }:=\left\{\begin{array}{cc}1,& \mathrm{if}\eta =0,\hfill \\ \prod _{i=1}^{\eta }\left(\psi -(q^{i-1}{\psi }_0+{[i-1]}_{q}h)\right),& \mathrm{if}\eta >0,\hfill \end{array}\right.$$

where ${\psi }_0\in \mathbb{R}$.

The generalized quantum binomial reduces to q-binomial ${(\psi -{\psi }_0)}_q^{\eta }$ as $h\to 0$ and to h-binomial ${(\psi -{\psi }_0)}_h^{\eta }$ when $q\to 1$. Also, we note ${lim}_{(q,h)\to (1,0)}{(\psi -{\psi }_0)}_{q,h}^{\eta }={(\psi -{\psi }_0)}^{\eta }$.

Definition 3 

([10]). The generalized quantum exponential function ${\mathit{exp}}_{q,h}\left(\alpha \psi \right)$ is defined as

$${\mathit{exp}}_{q,h}\left(\alpha \psi \right):=\sum _{i=0}^{\infty }\frac{{\alpha }^i{(\psi -0)}_{q,h}^i}{{\left[i\right]}_q!},$$

where α is an arbitrary nonzero constant.

Clearly, we note ${\exp }_{q,h}\left(0\right)=1$. As $h\to 0$ and $\alpha =1$, the generalized quantum exponential function ${\exp }_{q,h}\left(\alpha \psi \right)$ becomes the so-called q-exponential function $e_q\left(\psi \right)$; see [4,5]. Also, as $q\to 1$ and $\alpha =1$, the generalized quantum exponential function ${\exp }_{q,h}\left(\alpha \psi \right)$ reduces to the so-called h-exponential function $e_h\left(\psi \right)={(1+h)}^{\frac{\psi }{h}}$; see [4].

Based on the above concepts, many mathematicians have studied q-special functions, q-differential equations, q-calculus, and so on; see [12,13,14,15,16]. For example, Duran, Acikgoz and Araci [16] considered different trigonometric functions and hyperbolic functions related to quantum numbers and looked for properties related to them. Mathematicians also discovered various theorems about basic concepts related to h-numbers. Benaoum [11] obtained Newton’s binomial formula in terms of $(q,h)$, and Cermak and Nechvatal [9] created a $(q,h)$ version of fractional calculus. In 2011, Rahmat [17] studied the $(q,h)$-Laplace transform. Silindir and Yantir [10] studied the quantum generalization of Taylor’s formula and the q-binomial coefficient in 2019; their work motivated the research reported in this paper. Mathematicians who study polynomials have already defined and characterized degenerate Euler polynomials. They also studied the definition and properties of Euler polynomials when combined with quantum numbers.

The main purpose of this paper is to construct degenerate quantum Euler polynomials (DQE polynomials) using properties of q-numbers and the $(q,h)$-derivative. The topic of this paper is a field of mathematics that can be expanded into subareas such as series methods or generalizations of existing series. Its results can also be applied in interdisciplinary areas such as nonlinear physics, as well as for solutions to nonlinear differential equations providing instrumental defects such as kinks, vortices, etc. The diagram below shows the relationship of Euler, q-Euler, and degenerate Euler polynomials to the degenerate quantum Euler polynomials (DQE polynomials) that we define here.

Definition 4 

([14,18]). q-Euler numbers and polynomials are defined as:

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\frac{2}{e_q\left(\tau \right)+1},\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}\left(\psi \right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\frac{2}{e_q\left(\tau \right)+1}{e}_q\left(\tau \psi \right).$$

Definition 5 

([19,20]). Degenerate Euler numbers and polynomials are defined as:

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta }\left(h\right)\frac{{\tau }^{\eta }}{\eta !}={\displaystyle \frac{2}{e_h(1:\tau )+1}},\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta }(\psi :h)\frac{{\tau }^{\eta }}{\eta !}={\displaystyle \frac{2}{e_h(1:\tau )+1}}{e}_h(\psi :\tau ).$$

The structure of this paper is as follows. In Section 2, we define DQE numbers and polynomials. We find several properties of these polynomials by using q-numbers and $(q,h)$-derivatives. In addition, we construct several higher-order differential equations whose solutions are DQE polynomials. Section 3 shows the structure of approximate roots of DQE polynomials, which are solutions of the higher-order differential equations obtained in Section 2. By observing various structures of these approximate roots, we can make several conjectures.

2. $(\mathit{q},\mathit{h})$-Differential Equations That are Related to DQE Polynomials

In this section, we define a degenerate q-exponential function and DQE numbers and polynomials. We find several $(q,h)$-differential equations that are related to DQE polynomials by the $(q,h)$-derivative. We also discuss how these DQE polynomials relate to both q-Euler polynomials and degenerate Euler polynomials. We introduce the following degenerate quantum exponential function:

$$e_{q,h}(\psi :\tau ):=\sum _{\eta =0}^{\infty }{\left(\psi \right)}_{q,h}^{\eta }\frac{{\tau }^{\eta }}{\eta !}=\sum _{\eta =0}^{\infty }\prod _{k=0}^{\eta }\left(\psi -{[k-1]}_{q}h\right)\frac{{\tau }^{\eta }}{\eta !}.$$

For example, substituting $\psi =1$ in the above equation, we have

$$e_{q,h}(1:\tau )=\sum _{\eta =0}^{\infty }{\left(1\right)}_{q,h}^{\eta }\frac{{\tau }^{\eta }}{\eta !},$$

where ${\left(1\right)}_{q,h}^n=1(1-h)\cdots (1-{[n-1]}_{q}h)$.

Definition 6. 

Let $\left|q\right|<1$ and let h be a non-negative integer. Then, we define the DQE polynomials ${\mathcal{E}}_{\eta ,q}(\psi :h)$ as:

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}(\psi :h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\frac{2}{e_{q,h}(1:\tau )+1}{e}_{q,h}(\psi :\tau ).$$

When $\psi =0$, we can note that

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}\left(h\right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\frac{2}{e_{q,h}(1:\tau )+1}.$$

We denote ${\mathcal{E}}_{\eta ,q}\left(h\right)$ as the DQE numbers. From Definition 6, we can see several relationships between Euler polynomials. Setting $h\to 0$ in Definition 6, we can find the q-Euler numbers ${\mathcal{E}}_{\eta ,q}$ and polynomials ${\mathcal{E}}_{\eta ,q}\left(\psi \right)$ as

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\frac{2}{e_q\left(\tau \right)+1},\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}\left(\psi \right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\frac{2}{e_q\left(\tau \right)+1}{e}_q\left(\tau \psi \right).$$

Let $h\to 0$ and $q\to 1$ in Definition 6. Then, we have the Euler numbers ${\mathcal{E}}_{\eta }$ and polynomials ${\mathcal{E}}_{\eta }\left(\psi \right)$ as

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta }\frac{{\tau }^{\eta }}{\eta !}=\frac{2}{e^{\tau }+1},\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta }\left(\psi \right)\frac{{\tau }^{\eta }}{\eta !}=\frac{2}{e^{\tau }+1}{e}^{\tau \psi }.$$

In addition, for $q\to 1$ in Definition 6, we can see the degenerate Euler numbers ${\mathcal{E}}_{\eta }\left(h\right)$ and polynomials ${\mathcal{E}}_{\eta }(\psi :h)$ as follows:

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta }\left(h\right)\frac{{\tau }^{\eta }}{\eta !}=\frac{2}{e_h(1:\tau )+1},\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta }(\psi :h)\frac{{\tau }^{\eta }}{\eta !}=\frac{2}{e_h(1:\tau )+1}{e}_h(\psi :\tau ),$$

where ${\mathcal{E}}_{\eta }\left(h\right)={\mathcal{E}}_{\eta }(0:h)$.

Theorem 2. 

Let $\left|q\right|<1$ and $h\in \mathbb{N}$. Then, we have

$$D_{q,h}{\mathcal{E}}_{\eta ,q}(\psi :h)={\left[\eta \right]}_q{\mathcal{E}}_{\eta -1,q}(\psi :h).$$

Proof. 

From the generating function of the DQE polynomials, we can find

$$\begin{gathered} \sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}(\psi :h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}\left(h\right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}\sum _{\eta =0}^{\infty }{\left(\psi \right)}_{q,h}^{\eta }\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!} \\ =\sum _{\eta =0}^{\infty }\left(\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q{\left(\psi \right)}_{q,h}^{\eta -k}{\mathcal{E}}_{k,q}\left(h\right)\right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}. \end{gathered}$$

Applying the coefficient comparison method to the above equation, we find a relation of DQE numbers and polynomials such as

$${\mathcal{E}}_{\eta ,q}(\psi :h)=\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q{\left(\psi \right)}_{q,h}^{\eta -k}{\mathcal{E}}_{k,q}\left(h\right).$$

(1)

Using the $(q,h)$-derivative in Equation (1), we can obtain the following equation:

$$D_{q,h}{\mathcal{E}}_{\eta ,q}(\psi :h)=\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q{[\eta -k]}_q{\left(\psi \right)}_{q,h}^{\eta -k-1}{\mathcal{E}}_{k,q}\left(h\right).$$

Considering Equation (1) and the above equation together, we have the required result. □

Corollary 1. 

Let k be a non-negative integer. From Theorem 2, one holds

$${\mathcal{E}}_{\eta -k,q}(\psi :h)=\frac{{[\eta -k]}_q!}{{\left[\eta \right]}_q!}{D}_{q,h}^{\left(k\right)}{\mathcal{E}}_{\eta ,q}(\psi :h).$$

Corollary 2. 

From Theorem 2, the following holds:

(i) Setting $q\to 1$ in Theorem 2, we have

$$D_h{\mathcal{E}}_{\eta }(\psi :h)=\eta {\mathcal{E}}_{\eta -1}(\psi :h),{\mathcal{E}}_{\eta -k}(\psi :h)=\frac{(\eta -k)!}{\eta !}{D}_h^{\left(k\right)}{\mathcal{E}}_{\eta }(\psi :h),$$

where $D_h$ is the h-derivative and ${\mathcal{E}}_{\eta }(\psi :h)$ is the degenerate Euler polynomial.

(ii) Putting $h\to 0$ in Theorem 2, we have

$$D_q{\mathcal{E}}_{\eta ,q}\left(\psi \right)={\left[\eta \right]}_q{\mathcal{E}}_{\eta -1,q}\left(\psi \right),{\mathcal{E}}_{\eta -k,q}\left(\psi \right)=\frac{{[\eta -k]}_q!}{{\left[\eta \right]}_q!}{D}_q^{\left(k\right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right),$$

where $D_q$ is the q-derivative and ${\mathcal{E}}_{\eta ,q}\left(\psi \right)$ is the q-Euler polynomial.

Theorem 3. 

The DQE polynomials represent a solution of the $(q,h)$-differential equation of higher order shown below:

$$\begin{gathered} \frac{{\left(1\right)}_{q,h}^{\eta }}{{\left[\eta \right]}_q!}{D}_{q,h}^{\left(\eta \right)}{\mathcal{E}}_{\eta ,q}(\psi :h)+\frac{{\left(1\right)}_{q,h}^{\eta -1}}{{[\eta -1]}_q!}{D}_{q,h}^{(\eta -1)}{\mathcal{E}}_{\eta ,q}(\psi :h)+\frac{{\left(1\right)}_{q,h}^{\eta -2}}{{[\eta -2]}_q!}{D}_{q,h}^{(\eta -2)}{\mathcal{E}}_{\eta ,q}(\psi :h) \\ +\cdots +\frac{{\left(1\right)}_{q,h}^2}{{\left[2\right]}_q!}{D}_{q,h}^{\left(2\right)}{\mathcal{E}}_{\eta ,q}(\psi :h)+{\left(1\right)}_{q,h}^1{D}_{q,h}^{\left(1\right)}{\mathcal{E}}_{\eta ,q}(\psi :h)+2{\mathcal{E}}_{\eta ,q}(\psi :h)-2{\left(\psi \right)}_{q,h}^{\eta }=0. \end{gathered}$$

Proof. 

Consider $e_{q,h}(1:\tau )\ne -1$ in the generating function of the DQE polynomials. Then, we have

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}(\psi :h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_{q,h}!}\left(e_{q,h}(1:\tau )+1\right)=2e_{q,h}(\psi :\tau ).$$

(2)

The left-hand side of Equation (2) is transformed as

$$\begin{gathered} \sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}(\psi :h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}\left(e_{q,h}(1:\tau )+1\right) \\ =\sum _{\eta =0}^{\infty }\left(\sum _{k=0}^{\eta }{\left[\begin{array}{c}n\\ k\end{array}\right]}_q{\left(1\right)}_{q,h}^k{\mathcal{E}}_{\eta -k,q,h}\left(\psi \right)+{\mathcal{E}}_{\eta ,q}(\psi :h)\right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}. \end{gathered}$$

The right-hand side of Equation (2) is changed as

$$2e_{q,h}(\psi :\tau )=2\sum _{\eta =0}^{\infty }{\left(\psi \right)}_{q,h}^{\eta }\frac{t^{\eta }}{{\left[\eta \right]}_q!}.$$

From the above equations, we can obtain

$$\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q{\left(1\right)}_{q,h}^k{\mathcal{E}}_{\eta -k,q}(\psi :h)+{\mathcal{E}}_{\eta ,q}(\psi :h)=2{\left(\psi \right)}_{q,h}^{\eta }.$$

(3)

Considering both Corollary 2 and Equation (3), we have

$$\sum _{k=0}^{\eta }\frac{{\left(1\right)}_{q,h}^k}{{\left[k\right]}_q!}{D}_{q,h}^{\left(k\right)}{\mathcal{E}}_{\eta ,q}(\psi :h)+{\mathcal{E}}_{\eta ,q}(\psi :h)-2{\left(\psi \right)}_{q,h}^{\eta }=0,$$

which is the desired result. □

Corollary 3. 

Setting $q\to 1$ in Theorem 3, it holds that:

$$\begin{gathered} \frac{{\left(1\right)}_h^{\eta }}{\eta !}{D}_h^{\left(\eta \right)}{\mathcal{E}}_{\eta }(\psi :h)+\frac{{\left(1\right)}_h^{\eta -1}}{(\eta -1)!}{D}_h^{(\eta -1)}{\mathcal{E}}_{\eta }(\psi :h)+\frac{{\left(1\right)}_h^{\eta -2}}{(\eta -1)!}{D}_h^{(\eta -2)}{\mathcal{E}}_{\eta }(\psi :h) \\ +\cdots +\frac{{\left(1\right)}_h^2}{2!}{D}_h^{\left(2\right)}{\mathcal{E}}_{\eta }(\psi :h)+{\left(1\right)}_h{D}_h^{\left(1\right)}{\mathcal{E}}_{\eta }(\psi :h)+2{\mathcal{E}}_{\eta }(\psi :h)-2{\left(\psi \right)}_h^{\eta }=0, \end{gathered}$$

where $D_h$ is the h-derivative and ${\mathcal{E}}_{\eta }(\psi :h)$ is the degenerate Euler polynomial.

Corollary 4. 

Let $h\to 0$ in Theorem 3. Then, it holds that:

$$\begin{array}{l}\frac{1}{{\left[\eta \right]}_q!}{D}_q^{\left(\eta \right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+\frac{1}{{[\eta -1]}_q!}{D}_q^{(\eta -1)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+\frac{1}{{[\eta -2]}_q!}{D}_q^{(\eta -2)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+\cdots \\ +\frac{1}{{\left[2\right]}_q!}{D}_q^{\left(2\right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+D_q^{\left(1\right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+2{\mathcal{E}}_{\eta ,q}\left(\psi \right)-2{\psi }^{\eta }=0,\end{array}$$

where $D_q$ is the q-derivative and ${\mathcal{E}}_{\eta ,q}\left(\psi \right)$ is the q-Euler polynomial.

Theorem 4. 

DQE polynomials are the solutions to the $(q,h)$-differential equation of higher order as

$$\begin{gathered} \frac{{\mathcal{E}}_{\eta ,q}(1:h)+{\mathcal{E}}_{\eta ,q}\left(h\right)}{{\left[\eta \right]}_q!}{D}_{q,h}^{\left(\eta \right)}{\mathcal{E}}_{\eta ,q}(\psi :h)+\frac{{\mathcal{E}}_{\eta -1,q}(1:h)+{\mathcal{E}}_{\eta -1,q}\left(h\right)}{{[\eta -1]}_q!}{D}_{q,h}^{(\eta -1)}{\mathcal{E}}_{\eta ,q}(\psi :h)+\cdots \\ +\frac{{\mathcal{E}}_{2,q}(1:h)+{\mathcal{E}}_{2,q}\left(h\right)}{{\left[2\right]}_q!}{D}_{q,h}^{\left(2\right)}{\mathcal{E}}_{\eta ,q}(\psi :h)+({\mathcal{E}}_{1,q}(1:h)+{\mathcal{E}}_{1,q}\left(h\right))D_{q,h}^{\left(1\right)}{\mathcal{E}}_{\eta ,q}(\psi :h) \\ +\left({\mathcal{E}}_{0,q}(1:h)+{\mathcal{E}}_{0,q}\left(h\right)-2\right){\mathcal{E}}_{\eta ,q}(\psi :h)=0. \end{gathered}$$

Proof. 

From Definition 6, we have

$$\begin{gathered} \sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}(\psi :h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!} \\ =\frac{2}{e_{q,h}(1:\tau )+1}{e}_{q,h}(\psi :\tau ) \\ =\frac{1}{2}\left(\frac{2}{e_{q,h}(1:\tau )+1}{e}_{q,h}(1:\tau )+\frac{2}{e_{q,h}(1:\tau )+1}\right)\frac{2}{e_{q,h}(1:\tau )+1}{e}_{q,h}(\psi :\tau ). \end{gathered}$$

Using the generating function of DQE polynomials, we find the following relation:

$$2\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}(\psi :h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\sum _{\eta =0}^{\infty }\left(\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q\left({\mathcal{E}}_{k,q,h}\left(1\right)+{\mathcal{E}}_{k,q,h}\right){\mathcal{E}}_{\eta -k,q,h}\left(\psi \right)\right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}.$$

Comparing the coefficients of both sides, we obtain

$$\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q\left({\mathcal{E}}_{k,q}(1:h)+{\mathcal{E}}_{k,q}\left(h\right)\right){\mathcal{E}}_{\eta -k,q}(\psi :h)-2{\mathcal{E}}_{\eta ,q}(\psi :h)=0.$$

(4)

Replacing $D_{q,h}^{\left(k\right)}{\mathcal{E}}_{\eta ,q}(\psi :h)$ for ${\mathcal{E}}_{\eta -k,q}(\psi :h)$ in Equation (4), we derive

$$\sum _{k=0}^{\eta }\frac{\left({\mathcal{E}}_{k,q}(1:h)+{\mathcal{E}}_{k,q}\left(h\right)\right)}{{\left[k\right]}_q!}{D}_{q,h}^{\left(k\right)}{\mathcal{E}}_{\eta ,q}(\psi :h)-2{\mathcal{E}}_{\eta ,q}(\psi :h)=0.$$

The equation above completes the proof of Theorem 4. □

Corollary 5. 

Setting $h\to 0$ in Theorem 4, it holds that:

$$\begin{array}{l}\frac{{\mathcal{E}}_{\eta ,q}\left(1\right)+{\mathcal{E}}_{\eta ,q}}{{\left[\eta \right]}_q!}{D}_q^{\left(\eta \right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+\frac{{\mathcal{E}}_{\eta -1,q}\left(1\right)+{\mathcal{E}}_{\eta -1,q}}{{[\eta -1]}_q!}{D}_{q,\psi }^{(\eta -1)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+\cdots \\ +\frac{{\mathcal{E}}_{2,q}\left(1\right)+{\mathcal{E}}_{2,q}}{{\left[2\right]}_q!}{D}_q^{\left(2\right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+({\mathcal{E}}_{1,q}\left(1\right)+{\mathcal{E}}_{1,q})D_q^{\left(1\right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)\\ +\left({\mathcal{E}}_{0,q,}\left(1\right)+{\mathcal{E}}_{0,q}-2\right){\mathcal{E}}_{\eta ,q}\left(\psi \right)=0,\end{array}$$

where $D_q$ is the q-derivative and ${\mathcal{E}}_{\eta ,q}\left(\psi \right)$ is the q-Euler polynomial.

Corollary 6. 

Putting $q\to 1$ in Theorem 4, the following holds:

$$\begin{gathered} \frac{{\mathcal{E}}_{\eta }(1:h)+{\mathcal{E}}_{\eta }\left(h\right)}{\eta !}{D}_h^{\left(\eta \right)}{\mathcal{E}}_{\eta }(\psi :h)+\frac{{\mathcal{E}}_{\eta -1}(1:h)+{\mathcal{E}}_{\eta -1}\left(h\right)}{(\eta -1)!}{D}_h^{(\eta -1)}{\mathcal{E}}_{\eta }(\psi :h)+\cdots \\ +\frac{{\mathcal{E}}_2(1:h)+{\mathcal{E}}_2\left(h\right)}{2!}{D}_h^{\left(2\right)}{\mathcal{E}}_{\eta }(\psi :h)+({\mathcal{E}}_1(1:h)+{\mathcal{E}}_1\left(h\right))D_h^{\left(1\right)}{\mathcal{E}}_{\eta }(\psi :h) \\ +\left({\mathcal{E}}_0(1:h)+{\mathcal{E}}_0\left(h\right)-2\right){\mathcal{E}}_{\eta }(\psi :h)=0, \end{gathered}$$

where $D_h$ is the h-derivative and ${\mathcal{E}}_{\eta }(\psi :h)$ is the degenerate Euler polynomial.

From a property of $e_{q,h}(\psi :\tau )$, we note a relation

$$\begin{gathered} e_{q,h}(q\psi :\tau ) \\ =\sum _{\eta =0}^{\infty }q\psi (q\psi -h)(q\psi -{\left[2\right]}_{q}h)(q\psi -{\left[3\right]}_{q}h)\cdots (q\psi -{[\eta -1]}_{q}h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!} \\ =e_{q,q^{-1}h}(\psi :q\tau ). \end{gathered}$$

(5)

Theorem 5. 

The following higher-order differential equation has DQE polynomials as the solution:

$$\begin{gathered} \frac{q^{\eta }\left({\mathcal{E}}_{\eta ,q}(1:q^{-1}h)+{\mathcal{E}}_{\eta ,q}\left(q^{-1}h\right)\right)}{{\left[\eta \right]}_q!}{D}_{q,h}^{\left(\eta \right)}{\mathcal{E}}_{\eta ,q}(\psi :h) \\ +\frac{q^{\eta -1}\left({\mathcal{E}}_{\eta -1,q}(1:q^{-1}h)+{\mathcal{E}}_{\eta -1,q}\left(q^{-1}h\right)\right)}{{[\eta -1]}_q!}{D}_{q,h}^{(\eta -1)}{\mathcal{E}}_{\eta ,q}(\psi :h) \\ +\cdots +\frac{q^2\left({\mathcal{E}}_{2,q}(1:q^{-1}h)+{\mathcal{E}}_{2,q}\left(q^{-1}h\right)\right)}{{\left[2\right]}_q!}{D}_{q,h}^{\left(2\right)}{\mathcal{E}}_{\eta ,q}(\psi :h) \\ +q({\mathcal{E}}_{1,q}(1:q^{-1}h)+{\mathcal{E}}_{1,q}\left(q^{-1}h\right))D_{q,h}^{\left(1\right)}{\mathcal{E}}_{\eta ,q}(\psi :h) \\ +\left({\mathcal{E}}_{0,q}(1:q^{-1}h)+{\mathcal{E}}_{0,q}\left(q^{-1}h\right)-2\right){\mathcal{E}}_{\eta ,q}(\psi :h)=0. \end{gathered}$$

Proof. 

Consider Equation (5) in the generating function of the DQE polynomials. Then, we find

$$\begin{gathered} \sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}(q\psi :h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=\frac{2}{e_{q,h}(1:\tau )+1}{e}_{q,h}(q\psi :\tau ) \\ =\frac{1}{2}\left(\frac{2}{e_{q,q^{-1}h}(1:q\tau )+1}{e}_{q,q^{-1}h}(1:q\tau )+\frac{2}{e_{q,q^{-1}h}(1:q\tau )+1}\right) \\ \times \frac{2}{e_{q,h}(1:\tau )+1}{e}_{q,h}(q\psi :\tau ). \end{gathered}$$

From ${\mathcal{E}}_{n,q}(\psi :h)$, we have the relation

$$\begin{gathered} 2\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}(q\psi :h)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!} \\ =\sum _{\eta =0}^{\infty }\left(\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q{q}^k\left({\mathcal{E}}_{k,q}(1:q^{-1}h)+{\mathcal{E}}_{k,q}\left(q^{-1}h\right)\right){\mathcal{E}}_{\eta -k,q}(q\psi :h)\right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}. \end{gathered}$$

From the above equation, we find

$$\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q{q}^k\left({\mathcal{E}}_{k,q}(1:q^{-1}h)+{\mathcal{E}}_{k,q}\left(q^{-1}h\right)\right){\mathcal{E}}_{\eta -k,q}(q\psi :h)-2{\mathcal{E}}_{\eta ,q}(q\psi :h)=0.$$

(6)

Substituting $q\psi$ for $\psi$ in Corollary 2, we note that

$${\mathcal{E}}_{\eta -k,q}(q\psi :h)=\frac{{[\eta -k]}_q!}{{\left[\eta \right]}_q!}{D}_{q,h}^{\left(k\right)}{\mathcal{E}}_{\eta ,q}(q\psi :h).$$

Applying the above equation in Equation (6), we have

$$\sum _{k=0}^{\eta }\frac{q^k\left({\mathcal{E}}_{k,q}(1:q^{-1}h)+{\mathcal{E}}_{k,q}\left(q^{-1}h\right)\right)}{{\left[k\right]}_q!}{D}_{q,h}^{\left(k\right)}{\mathcal{E}}_{\eta ,q}(q\psi :h)-2{\mathcal{E}}_{\eta ,q}(q\psi :h)=0,$$

which gives the required result. □

Corollary 7. 

Setting $h\to 0$ in Theorem 5, it holds that:

$$\begin{gathered} \frac{q^{\eta }\left({\mathcal{E}}_{\eta ,q}\left(1\right)+{\mathcal{E}}_{\eta ,q}\right)}{{\left[\eta \right]}_q!}{D}_q^{\left(\eta \right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+\frac{q^{\eta -1}\left({\mathcal{E}}_{\eta -1,q}\left(1\right)+{\mathcal{E}}_{\eta -1,q}\right)}{{[\eta -1]}_q!}{D}_q^{(\eta -1)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+\cdots \\ +\frac{q^2\left({\mathcal{E}}_{2,q}\left(1\right)+{\mathcal{E}}_{2,q}\right)}{{\left[2\right]}_q!}{D}_q^{\left(2\right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right)+q({\mathcal{E}}_{1,q}\left(1\right)+{\mathcal{E}}_{1,q})D_q^{\left(1\right)}{\mathcal{E}}_{\eta ,q}\left(\psi \right) \\ +\left({\mathcal{E}}_{0,q}\left(1\right)+{\mathcal{E}}_{0,q}-2\right){\mathcal{E}}_{\eta ,q}\left(\psi \right)=0, \end{gathered}$$

where $D_q$ is the q-derivative and ${\mathcal{E}}_{\eta ,q}\left(\psi \right)$ is the q-Euler polynomial.

3. Properties for Approximate Roots of DQE Polynomials

In this section, we use Mathematica to show the structures and shapes of approximate roots of DQE polynomials. These DQE polynomials share several properties with both degenerate Euler polynomials and q-Euler polynomials. Here, the purpose of changing the q-number is to explore the properties of q-Euler polynomials, while the reason for changing the value of h is to explore the properties of degenerate Euler polynomials.

Let $e_{q,h}(1:\tau )\ne -1$. From the generating function of DQE numbers, we have

$$\sum _{\eta =0}^{\infty }{\mathcal{E}}_{\eta ,q}\left(h\right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}\left(\sum _{\eta =0}^{\infty }{\left(1\right)}_{q,h}^{\eta }\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}+1\right)=2.$$

Using the Cauchy product in the above equation, we obtain

$$\sum _{\eta =0}^{\infty }\left(\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q{\left(1\right)}_{q,h}^k{\mathcal{E}}_{\eta -k,q}\left(h\right)+{\mathcal{E}}_{\eta ,q}\left(h\right)\right)\frac{{\tau }^{\eta }}{{\left[\eta \right]}_q!}=2.$$

Therefore, we derive

$$\sum _{k=0}^{\eta }{\left[\begin{array}{c}\eta \\ k\end{array}\right]}_q{\left(1\right)}_{q,h}^k{\mathcal{E}}_{\eta -k,q}\left(h\right)+{\mathcal{E}}_{\eta ,q}\left(h\right)=\left\{\begin{array}{cc}2,& \mathrm{if}\eta =0,\hfill \\ 0,& \mathrm{otherwise}.\hfill \end{array}\right.$$

From the above equation, we show several DQE numbers ${\mathcal{E}}_{\eta ,q,h}$ as follows:

$$\begin{array}{lll}{\mathcal{E}}_{0,q}\left(h\right)& =& 1,\\ {\mathcal{E}}_{1,q}\left(h\right)& =& \frac{1}{2},\\ {\mathcal{E}}_{2,q}\left(h\right)& =& \frac{1}{4}(1+2h-q),\\ {\mathcal{E}}_{3,q}\left(h\right)& =& \frac{1}{8}(-4h^2(1+q)+4h(2+q)+(1+q)(1+(-3+q)q)),\\ {\mathcal{E}}_{4,q}\left(h\right)& =& \frac{1}{16}(8h^3(1+q)(1+q+q^2)-(-1+q)(1+q)(1+(-4+q)q)(1+q+q^2)\\ & & +4h^2(1+q)(-5+q(-2+(-2+q)q))+2h(11+q(6+q-q^2(3+q(2+q))))),\\ & \cdots & .\end{array}$$

Figure 1 shows the values of ${\mathcal{E}}_{\eta ,q,h}$ obtained by varying the values of q, $\eta$, and h. Non-negative integers on the x-axis in Figure 1 represent the value of $\eta$, with the 0 mark on the x-axis corresponding to ${\mathcal{E}}_{0,q,h}$, the 1 mark indicating ${\mathcal{E}}_{1,q,h}$, and so on. The lines represent variations of the approximate values for DQE numbers. The approximate values of $(0.1,h)$-Euler numbers in Figure 1a are the blue dots, yellow squares, green rhombuses, and red triangles, respectively, for $\eta =0,2,4,6$. Here, we can think of the blue dots as approximate values of q-Euler numbers. In Figure 1b, the blue dots, yellow squares, green rhombuses, and red triangles are the approximate values of $(q,2)$-Euler numbers in the respective cases $q=0.3,0.5,0.7,0.9$. The red triangles can be thought of as approximate values of h-Euler numbers.

We next consider the approximate roots of DQE polynomials. In order to identify some of these approximate roots and their properties, we need the exact shapes of the polynomials. Several DQE polynomials ${\mathcal{E}}_{\eta ,q}(\psi :h)$ are shown in the following:

$$\begin{gathered} {\mathcal{E}}_{0,q}(\psi :h)=1, \\ {\mathcal{E}}_{1,q}(\psi :h)=-\frac{1}{2}+\psi , \\ {\mathcal{E}}_{2,q}(\psi :h)=-\frac{1}{4}(-1+2h+q-2(1+2h+q)\psi +4{\psi }^2), \\ {\mathcal{E}}_{3,q}(\psi :h)=\frac{1}{8}(-(1+q)(1+4h^2+4h(-1+q)+(-3+q)q)) \\ +\frac{1}{4}(-1+q^3+4h^2(1+q)+4h(1+q+q^2))\psi \\ -\frac{1}{2}(1+q+q^2+2h(2+q)){\psi }^2+{\psi }^3, \\ \cdots . \end{gathered}$$

Based on the DQE polynomials obtained above, we hope to find out about the various behaviors of their approximate roots according to the changes in the values of q and h. We will restrict the value of q to less than 1, since DQE polynomials become degenerate Euler polynomials when q approaches 1 Also, DQE polynomials become q-Euler polynomials when h goes to 0, so the value of h has to exclude 0. According to the above condition, we can check the structures of approximate roots of DQE polynomials. Consider the case when the value of h is changed and but the value of q is fixed.

Let us fix $0\le \eta \le 50$ and $q=0.9$. Then, Figure 2 illustrates the structures of the approximate roots of DQE polynomials obtained by varying h. The condition of the top-left panel Figure 2a shows the case where $h=10$; the condition of the top-right Figure 2b is shown when $h=5$; the condition of the bottom-left Figure 2c is shown when $h=1$; and the condition of the bottom-right Figure 2d is shown when $h=0$. From Figure 2d, it can be seen that, when the value of h becomes 0, the shape of the approximate roots reduces to the approximate roots of q-Euler polynomials. The blue dots are the positions of the approximate roots that appear when the value of n is small, and the red dots are the positions of the approximate roots that appear when $n=50$ in Figure 2, Figure 3, Figure 4 and Figure 5.

Looking at Figure 2 from above gives the features shown in Figure 3. There, we can see an interesting phenomenon: Figure 3 shows the change in h with $q=0.9$ fixed for $0\le n\le 50$. For $h=10$, it seems that all values of real numbers are approximate roots in Figure 3a. In Figure 3b with $h=1$, it can be seen that the approximate roots all have values of real numbers. In Figure 3c, where $h=0$, we see that the shapes of the approximate roots have symmetric properties.

Based on Figure 3, we can see real values shown in

Table 1. Numbers of approximate real zeros of ${\mathcal{E}}_{\eta ,0.9}(\psi :h)$.

	h	10	1	0
$\mathit{\eta }$
1		1	1	1
2		2	2	2
3		3	3	3
4		4	4	4
5		5	5	1
⋯		⋯	⋯	⋯
10		10	10	2
⋯		⋯	⋯	⋯
20		16	20	4
⋯		⋯	⋯	⋯
30		26	30	4
⋯		⋯	⋯	⋯
40		34	40	6
⋯		⋯	⋯	⋯
50		44	50	4
⋯		⋯	⋯	⋯

, which shows approximate values of real numbers that appear when h is changed.

Table 1 shows that the numbers of approximate real roots match the value of $\eta$ when h is 1. When h is 1, the numbers of approximated real roots are also equal to the value of n. However, when h is 0, we can see that the approximated roots come out with both real and imaginary values. Considering Figure 2 and Figure 3, and Table 1, we can make the following conjecture:

Conjecture 1. 

Let us fix $q=0.9$. If $0\le \eta \le 50$, $h=1$, then all values of approximate roots for DQE polynomials will be found on the real axis.

Figure 4 shows the shapes that appear when q is fixed at $q=0.1$. The conditions of panels (a), (b), and (c) in Figure 4 are as follows:

(a)

${\mathcal{E}}_{\eta ,0.1}(\psi :10)$ for $0\le \eta \le 50$,

(b)

${\mathcal{E}}_{\eta ,0.1}(\psi :5)$ for $0\le \eta \le 50$,

(c)

${\mathcal{E}}_{\eta ,0.1}(\psi :1)$ for $0\le \eta \le 50$.

In Figure 4, the blue dots correspond to $\eta =0$, and the red dots indicate $\eta =50$. In (a), (b), and (c) of Figure 4, it can be seen that the approximate roots continue to accumulate even if $\eta$ increases and the value of h changes at one position.

Table 2. Approximate real roots of ${\mathcal{E}}_{\eta ,0.1}(\psi :h)$.

	h	10	5	1
$\mathit{\eta }$
⋯		⋯	⋯	⋯
10		$0.303694$	$0.167572$	$0.000545164$
⋯		⋯	⋯	⋯
20		$0.165068$	$0.0512303$	$5.29824\times {10}^{-7}$
⋯		⋯	⋯	⋯
30		$0.0912093$	$0.0174946$	$5.17401\times {10}^{-10}$
⋯		⋯	⋯	⋯
40		$0.0520318$	$0.00635171$	$5.05275\times {10}^{-13}$
⋯		⋯	⋯	⋯
50		$0.030506$	$0.00237742$	$4.93432\times {10}^{-16}$

shows the approximate roots of one pillar location in Figure 4; we can see that the positions of approximate roots stacked are close to 0. Based on Figure 4 and Table 2, we can make the following conjecture:

Conjecture 2. 

One of the approximate roots of the DQE polynomial has a value close to zero under the conditions $q=0.1$, $0\le \eta \le 50$, and $1\le h\le 10$.

Next, we explore further properties of DQE polynomials by varying the q-number, given the conditions $h=1$ and $0\le \eta \le 50$. Figure 5 shows a view from above of the 3D shape that appears when the q-number is changed under these constraints. The blue dots represent the low values of n; the red dots indicate where $n=50$. In Figure 5, the q value of the shape on the left is $0.0001$, the q value in the middle is $0.001$, and the q value on the right is $0.5$. Figure 5 shows that the position of the red dots changes according to the value of q and, based on this, we can make the following conjecture:

Conjecture 3. 

Let us fix $h=1$. Then, all values of approximate roots for higher-order DQE polynomials will appear on the real axis as the q-number approaches 0.

4. Conclusions

In this paper, we have introduced DQE polynomials and found higher-order differential equations related to these polynomials. By separately varying the q-number and the h-number of these DQE polynomials, we have shown some properties of their approximate roots and the structure of those roots. We have proposed several further conjectures on questions of interest which we hope will lead to new theorems and properties. This result can be applied to nonlinear physics or problems of finding solutions to nonlinear differential equations. Furthermore, we think the study of quantum polynomials with two variables is an interesting topic.