Existence and stability of solution for multi-order nonlinear fractional differential equations

Leping Xie; Jueliang Zhou; Haiyun Deng; Yubo He; Leping Xie; Jueliang Zhou; Haiyun Deng; Yubo He

doi:10.3934/math.2022899

AIMS Mathematics

2022, Volume 7, Issue 9: 16440-16448. doi: 10.3934/math.2022899

Previous Article Next Article

Research article Special Issues

Existence and stability of solution for multi-order nonlinear fractional differential equations

1.
School of Mathematics and Computer Science, Huaihua University, Hunan 418000, China
2.
Key Laboratory of Intelligent Control Technology for Wuling-Mountain Ecological Agriculture in Hunan Province, Hunan 418000, China

Received: 13 April 2022 Revised: 20 June 2022 Accepted: 27 June 2022 Published: 07 July 2022
MSC : 26A33, 34A08, 34K20, 45D05

In this paper, relying on the Banach contraction mapping principle to discuss the existence of solution for a multi-order nonlinear fractional differential equations on the infinite interval $[0, +\infty)$ . Moreover the stability of Ulam-Hyers-Rassias and Ulam-Hyers to the initial value problem are obtained. An example that can illustrate the conclusions of this paper have been given at the end.

Keywords:

multi-order nonlinear fractional differential equation,
Banach contraction mapping principle,
existence,
stability

Citation: Leping Xie, Jueliang Zhou, Haiyun Deng, Yubo He. Existence and stability of solution for multi-order nonlinear fractional differential equations[J]. AIMS Mathematics, 2022, 7(9): 16440-16448. doi: 10.3934/math.2022899

Related Papers:

[1]	Irina Di Ruocco . A political concept for the Gragnano Valley of Mills (Valle dei Mulini). Urban redevelopment of cultural-industrial heritage. Urban Resilience and Sustainability, 2023, 1(4): 278-308. doi: 10.3934/urs.2023018
[2]	Aibin Yan, Dinghan Zheng . Restoration and integration of the Huang Family Garden within the contemporary urban fabric of Shanghai. Urban Resilience and Sustainability, 2024, 2(1): 27-44. doi: 10.3934/urs.2024003
[3]	Elham Zabetian Targhi, Niusha Fardnava, Saba Saghafi . Urban development criteria with a focus on resilience to pandemics: A case study of coronavirus (Covid-19). Urban Resilience and Sustainability, 2023, 1(1): 66-85. doi: 10.3934/urs.2023005
[4]	Jichao Wang, Xiaoning Sui, Jie Zhang, Wenjie Shi, Wayne L. Thompson . Relocating disaster-prone villages and improving villager well-being: Evidence from Beijing, China. Urban Resilience and Sustainability, 2024, 2(3): 236-255. doi: 10.3934/urs.2024012
[5]	João C. G. Lanzinha . Rehabilitation of existing building parks and its relationship with urban agglomerates–An urgent and demanding task for our common future. Urban Resilience and Sustainability, 2023, 1(2): 86-90. doi: 10.3934/urs.2023006
[6]	Adoyo Laji, Jeremiah N. Ayonga . Mainstreaming resilience to flood risk among households in informal settlements in Kisumu City, Kenya. Urban Resilience and Sustainability, 2024, 2(4): 326-347. doi: 10.3934/urs.2024017
[7]	Maria Helena Luengo-Duque . Erasing roots: The impact of urban development on historical memory and identity in San Juan. Urban Resilience and Sustainability, 2025, 3(1): 26-56. doi: 10.3934/urs.2025002
[8]	Erik Velasco . Circular economy in Singapore: waste management, food and agriculture, energy, and transportation. Urban Resilience and Sustainability, 2024, 2(2): 110-150. doi: 10.3934/urs.2024007
[9]	Cheikh Cisse . At the crossroads of datas: Using artificial intelligence for efficient and resilient planning in African cities. Urban Resilience and Sustainability, 2023, 1(4): 309-313. doi: 10.3934/urs.2023019
[10]	Belhaj Naoufel . Determinants of cities' financial governance for sustainable urban development in Morocco—Case of the regions of Casablanca-Settat and Rabat-Sale-Kenitra. Urban Resilience and Sustainability, 2023, 1(1): 1-19. doi: 10.3934/urs.2023001

Abstract

1. Definition and notations

The special polynomials and numbers play an important role in many branches of mathematics such as combinatorics, computer science, statistics, etc. The generating functions of the special polynomials are used in the investigation of numerous properties satisfied by polynomials. A large number of indispensable special polynomials such as the Bernoulli, Euler and monic Hermite polynomials belong to the Appell class of sequences. These polynomials play crucial roles in many aspects, most importantly, they provide solutions to certain specified problems of physics, engineering and mathematics (see for more details ^{[5,6,7,8,9,21,22,34,35]}).

Throughout this paper, we use the following standard notations: $\mathbb{N}: = {\{1, 2, 3, \cdots \}}, \mathbb{N}_{0}: = \mathbb{N}\cup\{0\}: = {\{0, 1, 2, 3, \cdots \}}, \mathbb{Z}^{-}: = {\{-1, -2, -3, \cdots \}}, \mathbb{Z}$ denotes the set of integers and $\mathbb{R}$ denotes the set of real numbers.

The generating function of the Euler polynomials $E_q(v)$ ^[23], which arise in various classical and numerical analysis (and also of analytic number theory) is defined as:

$\begin{equation} \frac{2}{e^w +1} e^{vw} = \sum\limits_{q = 0}^\infty E_q(v) \frac{w^q}{q!}, \; \; \quad (|w| < \pi). \end{equation}$

(1.1)

The corresponding Euler numbers (appear in the Taylor series expansions of the secant and the hyperbolic secant functions) $E_q : = E_q(0)$ are given as:

$\begin{equation} \frac{2}{e^w +1} = \sum\limits_{q = 0}^\infty E_q \frac{w^q}{q!}, \qquad (|w| < \pi). \end{equation}$

(1.2)

Recently, Kim et al., ^[31] introduced the Changhee polynomials and numbers and established certain explicit formulae and identities for these polynomials. The Changhee polynomials $\mathcal{C}h_q(v)$ are specified by the following generating function ^[31]:

$\begin{equation} \frac{2}{2+w} (1+w)^v = \sum\limits_{q = 0}^{\infty} \mathcal{C}h_q(v) \frac{w^q}{q!}, \end{equation}$

(1.3)

which for $v = 0$ gives the corresponding Changhee numbers $\mathcal{C}h_q : = \mathcal{C}h_q(0)$ . Thereafter, Lee et al., ^[32] introduced a new form of Changhee polynomials, called the Appell-type Changhee polynomials $\mathcal{C}h_q^* (v)$ , and defined by the following generating equation ^[32]:

$\begin{equation} \frac{2}{2+w} e^{vw} = \sum\limits_{q = 0}^\infty \mathcal{C}h_q^* (v) \frac{w^q}{q!}, \end{equation}$

(1.4)

and have the following explicit relation [32,P3,Theorem 1]:

$\begin{equation} \mathcal{C}h_q^* (v) = \sum\limits_{p = 0}^{q} \binom{q}{p} \mathcal{C}h_q^* \; v^{q-p}, \end{equation}$

(1.5)

where $\mathcal{C}h_q^* : = \mathcal{C}h_q^* (0)$ are the corresponding Appell-type Changhee numbers:

$\begin{equation} \frac{2}{2+w} = \sum\limits_{q = 0}^\infty \mathcal{C}h_q^* \frac{w^q}{q!}. \end{equation}$

(1.6)

The above Appell-type Changhee numbers $\mathcal{C}h_q^*$ satisfy the following relation [32,P4,Theorem 4]:

$\begin{equation} 2\, \mathcal{C}h_q^* + q\, \mathcal{C}h_{q-1}^* = 0, \qquad \, \text{in the case when} \, \, q \geq 1, \end{equation}$

(1.7)

with initial condition $\mathcal{C}h_0^* = 1$ .

Motivated by the work of Kim et al., ^{[24,26,27,28,29,30,31]} and Lee et al., ^[31,32], in this paper we introduce the generating function of the Appell-type Changhee-Euler polynomials (ATCEP) and derive its properties.

2. Appell-type Changhee-Euler polynomials

In this section, the Appell-type Changhee-Euler polynomials (ATCEP), denoted by ${\mathcal{C}h}E_q^*(v)$ are introduced and certain properties of these polynomials are established.

Definition 1. The Appell-type Changhee-Euler polynomials ${\mathcal{C}h}E_q^*(v)$ are defined by the following generating relation:

$\begin{equation} \frac{4}{(e^w +1)(2+w)} e^{vw} = \sum\limits_{q = 0}^\infty {\mathcal{C}h}E_q^*(v) \frac{w^q}{q!}. \end{equation}$

(2.1)

Utilizing Eqs (1.1) and (1.6) in the left hand side of generating function (2.1) and then using Cauchy-product rule in the left hand side of the obtained equation, we get

$\begin{equation*} \sum\limits_{q = 0}^\infty \sum\limits_{l = 0}^q \binom{q}{l} E_l(v)\; \mathcal{C}h_{q-l}^* \frac{w^q}{q!} = \sum\limits_{q = 0}^\infty {\mathcal{C}h}E_q^*(v) \frac{w^q}{q!}. \end{equation*}$

Comparing the coefficients of identical powers of $w$ in both sides of the foregoing equation, we get the following series expansion for Appell-type Changhee-Euler polynomials ${\mathcal{C}h}E_q^*(v)$ :

$\begin{equation} \sum\limits_{l = 0}^q \binom{q}{l} E_l(v)\; \mathcal{C}h^*_{q-l} = {\mathcal{C}h}E_q^*(v). \end{equation}$

(2.2)

Motivated by the importance of the Euler numbers ^[36,38] in various branches of mathematics specially in number theory and combinatorics, we introduce the numbers related to the polynomial families introduced here. Taking $v = 0$ in series definition (2.2) of the ATCEP ${\mathcal{C}h}E_q^*(v)$ and using notation:

$\begin{equation} {\mathcal{C}h}E_q^* : = {\mathcal{C}h}E_q^*(0), \end{equation}$

(2.3)

in the left hand side, we get the following series form for Appell-type Changhee-Euler number (ATCEN) ${\mathcal{C}h}E_q^*$ :

$\begin{equation} \sum\limits_{l = 0}^q \binom{q}{l} E_l\; \mathcal{C}h_{q-l}^* = {\mathcal{C}h}E_q^*. \end{equation}$

(2.4)

Also, we find the following generating functions for ATCEN ${\mathcal{C}h}E_q^*$ :

$\begin{equation} \frac{4}{(e^w +1) (2+w)} = \sum\limits_{q = 0}^\infty {\mathcal{C}h}E_q^* \frac{w^q}{q!}. \end{equation}$

(2.5)

Theorem 2.1. Let $u, v \in \mathbb{C}$ and $q \in \mathbb{N}_0$ , then the following addition formulas hold:

$\begin{equation} {\mathcal{C}h}E_q^*(u+v) = \sum\limits_{p = 0}^{q}\binom{q}{p} E_p(u)\; \mathcal{C}h_{q-p}^*(v). \end{equation}$

(2.6)

$\begin{equation} {\mathcal{C}h}E_q^*(u+v) = \sum\limits_{p = 0}^{q}\binom{q}{p} u^{q-p}\; {\mathcal{C}h}E_p^*(v). \end{equation}$

(2.7)

Proof. Replacing $v$ by $u +v$ in Eq (2.1), we get

$\begin{equation} \Bigg(\frac{2}{e^w +1}\Bigg) \Bigg(\frac{2}{2+w} \Bigg) e^{(v+u)w} = \sum\limits_{q = 0}^\infty {\mathcal{C}h}E_q^*(v+u) \frac{w^q}{q!}, \end{equation}$

(2.8)

which on using Eqs (1.4) and (1.1) in the l.h.s. and then on comparing the coefficients of the identical powers of $w$ in both sides of the obtained relation, we get assertion (2.6).

Equation (2.8) can be written as

$\begin{equation} \sum\limits_{q = 0}^\infty {\mathcal{C}h}E_q^*(v+u) \frac{w^q}{q!} = \frac{4}{(e^w +1) (2+w)} e^{vw}\; e^{uw}, \end{equation}$

(2.9)

which on using Eq (2.1) and expanding the second exponential of the right side gives

$\begin{equation} \sum\limits_{q = 0}^\infty {\mathcal{C}h}E_q^*(v+u) \frac{w^q}{q!} = \sum\limits_{q = 0}^\infty \frac{(uw)^q}{q!}\; \sum\limits_{p = 0}^\infty {\mathcal{C}h}E_p^*(v) \frac{w^p}{p!}. \end{equation}$

(2.10)

Using the Cauchy product rule on the right side of Eq (2.10) and on comparing the coefficients of $w$ , we are lead to assertion (2.7).

Theorem 2.2. Let $q \in \mathbb{N}_0$ and $v \in \mathbb{C}$ , then the following identities hold true:

$\begin{equation} 2\; E_q(v) = 2\; {\mathcal{C}h^*}E_q(v)+ q\; {\mathcal{C}h}E_{q-1}^*(v). \end{equation}$

(2.11)

$\begin{equation} 2\; \mathcal{C}h_q^*(v) = {\mathcal{C}h^*}E_q(v) + \sum\limits_{p = o}^{q} \binom{q}{p} {\mathcal{C}h}E_{q-p}^*(v). \end{equation}$

(2.12)

$\begin{equation} {\mathcal{C}h}E_q^*(v) = \sum\limits_{p = 0}^{q} \binom{q}{p} \frac{(-1)^p p!\; E_{q-p}(v)}{2^p}. \end{equation}$

(2.13)

$\begin{equation} {\mathcal{C}h}E_q^*(v) = \sum\limits_{p = 0}^{q} \binom{q}{p} {\mathcal{C}h}E_{q-p}^* (v)^p. \end{equation}$

(2.14)

$\begin{equation} {\mathcal{C}h}E_q^*(v) = \sum\limits_{p = 0}^{q} \binom{q}{p} {\mathcal{C}h^*}_{q-p}\left(\frac{v}{2}\right) E_p\left(\frac{v}{2}\right). \end{equation}$

(2.15)

$\begin{equation} \frac{d}{ dv} {\mathcal{C}h}E_q^*(v) = q\; {\mathcal{C}h}E_{q-1}^*(v). \end{equation}$

(2.16)

$\begin{equation} \frac{d^p}{dv^p} {\mathcal{C}h}E_q^*(v) = \frac{q!}{(q-p)!} {\mathcal{C}h}E_{q-p}^*(v). \end{equation}$

(2.17)

$\begin{equation} \int_{a}^{v} {\mathcal{C}h}E_{q-1}^*(u) du = \frac{1}{q}\; \left({\mathcal{C}h}E_{q}^*(v) - {\mathcal{C}h}E_{q}^*(a)\right). \end{equation}$

(2.18)

$\begin{equation} {\mathcal{C}h}E_q^* = \sum\limits_{p = 0}^{q} \binom{q}{p} (-v)^p\; {\mathcal{C}h}E_{q-p}^*(v). \end{equation}$

(2.19)

$\begin{equation} 2\; \mathcal{C}h_q^* = \sum\limits_{p = 0}^{q} \binom{q}{p} {\mathcal{C}h}E_{q-p}^*(v) (-v)^p \left( 1 + \left(1 -\frac{1}{v} \right)^p \right). \end{equation}$

(2.20)

Proof. These identities are easily derivable with the help of generating Eqs (1.1), (1.4) and (2.1) by manipulating double series. Now, rewriting the generating Eq (2.1) as

$\begin{equation} \frac{4}{(e^w +1)(2+w)} = \sum\limits_{q = 0}^\infty {\mathcal{C}h}E_q^*(v) \frac{w^q}{q!} \left(e^{-vw} \right), \end{equation}$

(2.21)

which on using generating relation (2.5) in the l.h.s. and expanding the exponential on the r.h.s. and then on equating the coefficients of identical powers of $w$ in both sides, we are led to assertion (2.19).

From Eq (2.21), we have

$\begin{equation} \frac{4}{2+w} = \sum\limits_{q = 0}^\infty {\mathcal{C}h}E_q^*(v) \frac{w^q}{q!} \left(e^{-vw}(e^w +1) \right), \end{equation}$

(2.22)

which on using generating relation (1.6) in l.h.s. and expanding the exponentials on the r.h.s. and then on equating the coefficients of identical powers of $w$ in both sides of the obtained result, yields assertion (2.20).

Many researchers have established the determinant form of certain special and $q$ -special polynomials ^{[4,12,16,18,19,20,33]}. Costabile ^[10,11,13] has given several approaches to Bernoulli polynomials. An important approach based on a determinant definition was given in ^[10]. This approach is further extended to provide determinant definitions of the Appell and Sheffer polynomial sequences by Costabile and Longo in ^[12] and ^[14] respectively. Motivated by the recent work on the determinant approaches, here we give the determinant representation of ATCEP ${\mathcal{C}h}E_q^*(v)$ . We recall the determinant form of the Euler polynomials $E_q(v)$ :

The Euler polynomials $E_q(v)$ of degree $q = 0, 1, 2, \cdots$ are defined by

$\begin{equation} \begin{array}{l} E_0(v) = 1, \\ E_q(v) = {(-1)}^q \left|\begin{array}{ccccccc} 1 & v & v^2 & \cdots & v^{q-1} & v^q\\ \\ 1 & \frac{1}{2} & \frac{1}{2} & \cdots & \frac{1}{2} & \frac{1}{2}\\ \\ 0 & 1 & \frac{1}{2}\binom{2}{1} & \cdots & \frac{1}{2} \binom{q-1}{1}& \frac{1}{2} \binom{q}{1}\\ \\ 0 & 0 & 1 & \cdots & \frac{1}{2} \binom{q-1}{2} & \frac{1}{2} \binom{q}{2}\\ . & . & . & \cdots & . & . \\ . & . & .& \cdots & . & . \\ 0 & 0 & 0 & \cdots & 1 & \frac{1}{2} \binom{q}{q-1} \end{array} \right|, \; \; q = 1, 2, .... \end{array} \end{equation}$

(2.23)

Definition 2. The Appell-type Changhee-Euler polynomials ${\mathcal{C}h}E_q^*(v)$ of degree $q$ are defined by

$\begin{equation} \begin{array}{l} {\mathcal{C}h}E_0^*(v) = 1, \\ {\mathcal{C}h}E_q^*(v) = {(-1)}^q \left|\begin{array}{ccccccc} 1 & \mathcal{C}h^*_1(v) & \mathcal{C}h^*_2(v) & \cdots & \mathcal{C}h^*_{q-1}(v) & \mathcal{C}h^*_q(v)\\ \\ 1 & \frac{1}{2} & \frac{1}{2} & \cdots & \frac{1}{2} & \frac{1}{2}\\ \\ 0 & 1 & \frac{1}{2}\binom{2}{1} & \cdots & \frac{1}{2} \binom{q-1}{1}& \frac{1}{2} \binom{q}{1}\\ \\ 0 & 0 & 1 & \cdots & \frac{1}{2} \binom{q-1}{2} & \frac{1}{2} \binom{q}{2}\\ . & . & . & \cdots & . & . \\ . & . & .& \cdots & . & . \\ 0 & 0 & 0 & \cdots & 1 & \frac{1}{2} \binom{q}{q-1} \end{array} \right|, \; \; q = 1, 2, .... \end{array} \end{equation}$

(2.24)

where the Appell-type Changhee polynomials $\mathcal{C}h^*_q(v)\; (q = 0, 1, 2, ...)$ are given by Eq (1.4).

In the forthcoming section, a family of linear differential equations arising from the generating function of ATCEN ${\mathcal{C}h}E_q^*$ is derived by following the approach presented in ^[26].

3. Differential equations

Differential equations, besides playing important role in pure mathematics, constitute fundamental part of mathematical description of physical processes. Thus, obtaining the solutions for differential equations is of paramount importance. Few types of differential equations allow explicit and straightforward analytical solutions. Many researchers obtained the differential equations for special polynomials and numbers utilizing their generating functions ^[17,25,30].

Theorem 3.1. For $n \in \mathbb{N}$ the following family of differential equations for Appell-type Changhee-Euler polynomials ${\mathcal{C}h}E_q^*(v)$ holds true:

$\begin{equation} \sum\limits_{j = 1}^n \beta_j(n) e^w G^{(j-1)}(w) + (e^w +1) G^{(n)}(w) = (-1)^{n}\; \frac{\beta_0(n)}{(2+w)^{(n+1)}}, \end{equation}$

(3.1)

where $G^{(n)}(w) : = \big(\frac{d}{dw}\big)^n G(w),$

$\begin{equation} \; \; \; \; \; \; G = G(w) = \frac{4}{(e^w +1) (2+w)}, \end{equation}$

(3.2)

$\begin{equation} \beta_0(n) = 4 \, n!, \; \; \; \; \; \; \beta_1(n) = 1, \; \; \; \; \; \; \beta_n(n) = n, \end{equation}$

(3.3)

and

$\begin{equation} \; \; \; \; \; \; \; \beta_i(n+1) = 1+ \sum\limits_{j_1 = 0}^{n-i+1} \beta_{i-1}(n-j_1), \; \; \; \; \; \quad (2 \leq i\leq n). \end{equation}$

(3.4)

Proof. Let

$\begin{equation*} G = G(w) = \frac{4}{(e^w +1) (2+w)}, \end{equation*}$

which can be rewritten as:

$\begin{equation} (e^w +1) G(w) = \frac{4}{2+w}. \end{equation}$

(3.5)

Taking the derivative of Eq (3.5) with respect to $w,$ it follows that:

$\begin{equation} (e^w +1) G^{(1)}(w) + e^w G(w) = (-1)\; \frac{4}{(2+w)^{2}}, \end{equation}$

(3.6)

Again, by taking the derivative with respect to $w$ of Eq (3.6), we find

$\begin{equation} (e^w +1) G^{(2)}(w) + 2 e^w G^{(1)}(w) + e^w G(w) = (-1)^2\; \frac{8}{(2+w)^{3}}, \end{equation}$

(3.7)

and consequently

$\begin{equation} (e^w+1) G^{(3)}(w) + 3 e^w G^{(2)}(w) + 3 e^w G^{(1)}(w) + e^w G(w) = (-1)^3\; \frac{24}{(2 +w)^{4}}. \end{equation}$

(3.8)

Continuing this process, we find

$\begin{eqnarray} \sum\limits_{j = 1}^n \beta_j(n) e^w G^{(j-1)}(w) + (e^w +1) G^{(n)}(w) = (-1)^{n}\; \frac{\beta_0(n)}{(2+w)^{(n+1)}}. \end{eqnarray}$

(3.9)

Further, taking the derivative with respect to $w$ of Eq (3.9), we get

$\begin{eqnarray} \begin{split} &\beta_1(n) e^w G(w) + \sum\limits_{j = 2}^n e^w (\beta_j(n) + \beta_{j-1}(n)) G^{(j-1)}(w) + (\beta_n(n) +1) e^w G^{(n)}(w)\\ &\quad + (e^w +1) G^{(n+1)}(w) = (-1)^{n+1} (n+1) \beta_0(n) (2+w)^{-(n+2)}. \end{split} \end{eqnarray}$

(3.10)

Replacing $n$ by $n+1$ , Eq (3.9) may be rewritten as:

$\begin{eqnarray} \begin{split} \sum\limits_{j = 1}^{n+1} \beta_j(n+1) e^w G^{(j-1)}(w) &+ (e^w +1) G^{(n+1)}(w)\\ &\quad = (-1)^{n+1}\; \frac{\beta_0(n+1)}{(2+w)^{(n+2)}}. \end{split} \end{eqnarray}$

(3.11)

Comparison of the coefficients of $G^{(i)}$ 's on both sides of Eqs (3.10) and (3.11), yields the following recursive formulas:

$\begin{eqnarray} \beta_0(n+1) = (n+1) \beta_0(n), \end{eqnarray}$

(3.12)

$\begin{equation} \beta_1(n+1) = \beta_1(n), \; \; \; \; \; \; \; \; \; \; \end{equation}$

(3.13)

$\begin{equation} \; \; \; \; \; \beta_j(n+1) = \beta_j(n) + \beta_{j-1}(n), \; \; \; \; \; \; \quad (2 \leq j \leq n) \end{equation}$

(3.14)

and

$\begin{equation} \beta_{n+1} (n+1) = \beta_n(n) + 1. \end{equation}$

(3.15)

From Eqs (3.6) and (3.9), we have

$\begin{equation} \beta_0(1) = 4, \; \; \; \; \; \quad \beta_1(1) = 1. \end{equation}$

(3.16)

Again, from Eq (3.12), we note that

$\begin{equation*} \beta_0(n+1) = (n+1) \beta_0(n) = (n+1) n \beta_0(n) = .\; .\; . = (n+1) n (n-1) .\; .\; . 2.1 \, \beta_0(1), \end{equation*}$

that is

$\begin{equation} \; \; \beta_0(n +1) = 4\, (n+1)!\; . \; \; \; \; \; \; \; \; \; \; \; \end{equation}$

(3.17)

In view of Eq (3.13), we have

$\begin{equation} \beta_1(n+1) = \beta_1(n) = \beta_1(n-1) = .\; .\; . = \beta_1(1) = 1. \end{equation}$

(3.18)

Further, from Eq (3.15), it follows that:

$\begin{equation*} \beta_{n+1}(n+1) = \beta_n(n) +1 = \beta_{n-1}(n-1) + 2 = .\; .\; . = \beta_1(1) + n, \end{equation*}$

that is

$\begin{equation} \beta_{n+1} (n+1) = 1+ n.\; \; \; \; \; \; \; \; \; \; \; \; \end{equation}$

(3.19)

Taking $j = 2$ in Eq (3.14), we find

$\begin{equation*} \; \; \; \; \; \; \; \; \; \; \; \; \; \; \beta_2(n+1) = \beta_2(n) + \beta_1(n) = \beta_2(n-1) + \beta_1(n-1) + \beta_1(n) \end{equation*}$

$\begin{equation*} \; \; \; \; \; \; \; = .\; .\; . \; \; \; \; = \sum\limits_{j = 0}^{n-2} \beta_1(n-j) + \beta_2(2), \end{equation*}$

that is

$\begin{equation} \; \; \; \beta_2(n +1) = 1 + \sum\limits_{j = 0}^{n-1} \beta_1(n-j).\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{equation}$

(3.20)

Similarly, for $j = 3$ Eq (3.14) gives

$\begin{equation} \beta_3(n+1) = 1 + \sum\limits_{j = 0}^{n-2} \beta_2(n-j), \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{equation}$

(3.21)

consequently for $j = 4$ , we have

$\begin{equation} \beta_4(n+1) = 1+ \sum\limits_{j = 0}^{n-3} \beta_3(n-j).\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{equation}$

(3.22)

Continuing this process, we deduce that

$\begin{equation} \; \; \; \; \; \; \; \; \; \; \; \beta_i(n+1) = 1+ \sum\limits_{j_1 = 0}^{n-i+1} \beta_{i-1}(n-j_1), \; \; \; \; \; \; \; \quad (2 \leq j \leq n). \end{equation}$

(3.23)

Equation (3.9) together with Eqs (3.17), (3.18) and (3.23) completes the proof of Theorem 3.1.

Theorem 3.2. For $i, n \in \mathbb{N}, p \in \mathbb{N}_{0}$ the following identity for the Appell-type Changhee-Euler numbers ${\mathcal{C}h}E_q^*$ holds true:

$\begin{equation} \begin{split} \sum\limits_{j = 1}^n \sum\limits_{p = 0}^i \binom{i}{p} \beta_j(n)\; {\mathcal{C}h}E_{p+j-1}^* &+ \sum\limits_{p = 0}^{i} \binom{i}{p}\; {\mathcal{C}h}E_{p+n}^* + {\mathcal{C}h}E_{i+n}^*\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \\ &\quad = (-1)^{n+i} (n+1)_i \left( \frac{\beta_0(n)}{2^{n+i+1}}\right) , \end{split} \end{equation}$

(3.24)

where $(n+1)_i$ is the pochhammer symbol [,P21] and $\beta_j(n), \; \; (0\leq j \leq n)$ defined in Eqs (3.3) and (3.4).

Proof. Making use of Eq (2.5) in the first term of the left hand side of Eq (3.1), we find

$\begin{equation*} \sum\limits_{j = 1}^n \beta_j(n) e^w G^{(j-1)}(w) = \sum\limits_{j = 1}^n \beta_j(n) \sum\limits_{l = 0}^\infty \frac{w^l}{l!} \sum\limits_{p = 0}^\infty {\mathcal{C}h}E_{p+j-1}^* \frac{w^p}{p!}, \end{equation*}$

or equivalently

$\begin{equation} \; \; \; \; \; \; \sum\limits_{j = 1}^n \beta_j(n) e^w G^{(j-1)}(w) = \sum\limits_{i = 0}^\infty \Bigg( \sum\limits_{j = 1}^n \sum\limits_{p = 0}^i \binom{i}{p} \beta_j(n)\; {\mathcal{C}h}E_{p+j-1}^* \Bigg) \frac{w^i}{i!}. \end{equation}$

(3.25)

Again, using Eq (2.5) in the second term of the left hand side of Eq (3.1), we find

$\begin{equation*} (e^w +1) G^{(n)}(w) = \Bigg( \sum\limits_{p = 0}^\infty \frac{w^p}{p!} +1 \Bigg) \Bigg( \Big( \frac{d}{dw}\Big)^n \sum\limits_{l = 0}^\infty {\mathcal{C}h}E_l^* \frac{w^l}{l!}\Bigg), \end{equation*}$

that is

$\begin{equation} (e^w +1) G^{(n)}(w) = \sum\limits_{i = 0}^\infty \Bigg( \sum\limits_{p = 0}^{i} \binom{i}{p} {\mathcal{C}h}E_{p+n}^* + {\mathcal{C}h}E_{i+n}^* \Bigg) \frac{w^i}{i!}.\; \; \; \; \; \end{equation}$

(3.26)

Finally, from the right hand side of Eq (3.1), we have

$\begin{equation*} \beta_0(n)(-1)^n(2+w)^{-(n+1)} = \beta_0(n)(-1)^n 2^{-(n+1)}(1+ w/2)^{-(n+1)}, \end{equation*}$

which on using the following identity [37,P34]:

$\begin{equation*} (1-z)^{-a} = \sum\limits_{n = 0}^{\infty} (a)_n\; \frac{z^n}{n!}, \end{equation*}$

gives

$\begin{equation} \beta_0(n) (-1)^n (2+w)^{-(n+1)} = \sum\limits_{i = 0}^\infty \Bigg(\frac{(-1)^{n+i} \beta_0(n) (n+1)_i }{2^{n+i+1}} \Bigg) \frac{w^i}{i!}. \end{equation}$

(3.27)

Use of Eqs (3.25)–(3.27) in Eq (3.1), yields assertion (3.24).

In the next section, the graphical representations of the ATCEP ${\mathcal{C}h}E_q^*(v)$ and distribution of their zeros are discussed.

4. Graphical approach

Over the years, there has been increasing interest in solving mathematical problems with the aid of computers. The software Mathematica is used to show the behaviour of these newly introduced polynomials using the graphs plotted for the special values of indices. The manual computation of the zeros is too complicated, therefore, we use Mathematica to investigate the zeros of the ATCEP ${\mathcal{C}h}E_q^*(v)$ . By using numerical investigations and computer experiments, we find the real and complex zeros for certain values of index $q$ . The investigation in this direction will lead to a new approach employing numerical methods in the field of the hybrid special polynomials to appear in mathematics and physics, see for example ^[18,19].

With the help of Mathematica and using Eqs (1.5) and (1.7) in relation (2.2), and are drawn for odd values of index $q$ :

Figure 1. Curve of

${\mathcal{C}h}E_{15}^*(v)$ .

DownLoad: Full-Size Img PowerPoint

Figure 2. Curve of

${\mathcal{C}h}E_{25}^*(v)$ .

DownLoad: Full-Size Img PowerPoint

Similarly, for even values of index $q$ , Figures 3 and 4 are drawn:

Figure 3. Curve of

${\mathcal{C}h}E_{20}^*(v)$ .

DownLoad: Full-Size Img PowerPoint

Figure 4. Curve of

${\mathcal{C}h}E_{30}^*(v)$ .

DownLoad: Full-Size Img PowerPoint

Next, a remarkably regular structure of the zeros of ATCEP ${\mathcal{C}h}E_q^*(v)$ are explored with the help of Mathematica in Figures 5–8.

Figure 5. Zeros of

${\mathcal{C}h}E_{15}^*(v)$ .

DownLoad: Full-Size Img PowerPoint

Figure 6. Zeros of

${\mathcal{C}h}E_{20}^*(v)$ .

DownLoad: Full-Size Img PowerPoint

Figure 7. Zeros of

${\mathcal{C}h}E_{25}^*(v)$ .

DownLoad: Full-Size Img PowerPoint

Figure 8. Zeros of

${\mathcal{C}h}E_{30}^*(v)$ .

DownLoad: Full-Size Img PowerPoint

Further properties, applications and theoretical investigation of scattering of zeros of the ATCEP ${\mathcal{C}h}E_{q}^*(v)$ introduced here are left to the authors and the interested researchers, for future study.

5. Concluding remarks

Exclusive role have been played by special polynomials in applied mathematics. It is not astonishing when new classes of special polynomials are established as the issues associated with special polynomials are too immense. In this paper, by incorporating the Appell-type Changhee polynomials $\mathcal{C}h_q^* (v)$ and the Euler polynomials $E_q(v)$ in a natural way, the so-called Appell-type Changhee-Euler polynomials ${\mathcal{C}h}E_{q}^*(v)$ are introduced. Then, we investigate certain properties and identities for these new polynomials such as series representation, determinant form and some novel identities. We also present non-homogeneous linear differential equations for the Appell-type Changhee-Euler numbers ${\mathcal{C}h}E_{q}^*$ . Further we discuss the shape and zero distributions of these polynomials by observing their graphs drawn by Mathematica. The ATCEP ${\mathcal{C}h}E_{q}^*(v)$ were found to be the member of well-known Appell family and therefore many important properties and identities for this new polynomials were easily established by employing those in the well-developed theory of Appell polynomials. A unifying tool for studying polynomial sequences, namely the representation of Appell polynomials in matrix form has been studied in ^[1] and extended in ^[3] to the hypercomplex case. After that, Aceto and Caçāo ^[2] extends this approach to find the matrix representation of the Sheffer polynomials. The matrix which represents the Sheffer polynomial coefficients can be factorized into two matrices, one associated to Appell polynomials and the other linked to the binomial type polynomial sequences. This approach can be further extended to find the matrix representations of mixed special polynomials related to the Appell polynomials.

Posing a problem

If we consider the Genocchi polynomials ^[15] in place of Euler polynomials in Eq (2.1), we get the following generating function for the Appell-type Changhee-Genocchi polynomials (ATCGP) ${\mathcal{C}h}G_{q}^*(v)$ :

$\begin{equation} \frac{4w}{(e^w +1)(2+w)} e^{vw} = \sum\limits_{q = 0}^\infty {\mathcal{C}h}G_q^*(v) \frac{w^q}{q!}. \end{equation}$

(5.1)

We are believed to be able to establish the corresponding results for the Appell-type Changhee-Genocchi polynomials ${\mathcal{C}h}G_{q}^*(v)$ as those established in this paper. This posed problem is left to the interested researchers and authors for further investigation.

Acknowledgments

The authors are thankful to the reviewers for several useful comments and suggestions towards the improvement of this paper.

Conflict of interest

No conflict of interest was declared by the authors.

References

[1]	W. Tan, W. Pan, M. Xu, A note on unsteady flows of a viscoelastic fluid with fractional Maxwell model between two parallel plates, Int. J. Nonlin. Mech., 38 (2003), 645–650. http://dx.doi.org/10.1016/S0020-7462(01)00121-4 doi: 10.1016/S0020-7462(01)00121-4
[2]	A. Kilbas, H. Srivastava, J. Trujillo, Theory and applications of fractional differential equations, Boston: Elsevier, 2006.
[3]	S. Ray, Fractional calculus with applications for nuclear reactor dynamics, Boca Raton: CRC Press, 2015. http://dx.doi.org/10.1201/b18684
[4]	L. Liu, Q. Dong, G. Li, Exact solutions and Hyers-Ulam stability for fractional oscillation equations with pure delay, Appl. Math. Lett., 112 (2021), 106666. http://dx.doi.org/10.1016/j.aml.2020.106666 doi: 10.1016/j.aml.2020.106666
[5]	C. Chen, M. Bohner, B. Jia, Ulam-Hyers stability of Caputo fractional difference equations, Math. Method. Appl. Sci., 42 (2019), 7461–7470. http://dx.doi.org/10.1002/mma.5869 doi: 10.1002/mma.5869
[6]	S. Peng, J. Wang, Existence and Ulam-Hyers stability of ODEs involving two Caputo fractional derivatives, Electron. J. Qual. Theory Differ. Equ., 52 (2015), 1–16. http://dx.doi.org/10.14232/ejqtde.2015.1.52 doi: 10.14232/ejqtde.2015.1.52
[7]	Z. Akbar, W. Ali, S. Farina, Hyers-Ulam stability of nonlinear differential equations with fractional integrable impulses, Math. Method. Appl. Sci., 40 (2017), 5502–5514. http://dx.doi.org/10.1002/mma.4405 doi: 10.1002/mma.4405
[8]	Z. Yang, T. Xu, M. Qi, Ulam-Hyers stability for fractional differential equations in quaternionic analysis, Adv. Appl. Clifford Algebras, 26 (2016), 469–478. http://dx.doi.org/10.1007/s00006-015-0576-3 doi: 10.1007/s00006-015-0576-3
[9]	M. Akkouchi, Hyers-Ulam-Rassias stability of nonlinear volterra integral equations via a fixed point approach, Acta Universitatis Apulensis, 26 (2011), 257–266.
[10]	J. Wang, Y. Zhang, Ulam-Hyers-Mittag-Leffler stability of fractional-order delay differential equations, Optimization, 63 (2014), 1181–1190. http://dx.doi.org/10.1080/02331934.2014.906597 doi: 10.1080/02331934.2014.906597
[11]	M. Abdoab, S. Panchala, H. Wahash, Ulam-Hyers-Mittag-Leffler stability for a psi-Hilfer problem with fractional order and infinite delay, Results in Applied Mathematics, 7 (2020), 100115. http://dx.doi.org/10.1016/j.rinam.2020.100115 doi: 10.1016/j.rinam.2020.100115
[12]	M. Abdo, T. Abdeljawad, K. Shah, S. Ali, On nonlinear coupled evolution system with nonlocal subsidiary conditions under fractal-fractional order derivative, Math. Method. Appl. Sci., 44 (2021), 6581–6600. http://dx.doi.org/10.1002/mma.7210 doi: 10.1002/mma.7210
[13]	S. Ali, M. Abdo, Qualitative analysis for multiterm Langevin systems with generalized Caputo fractional operators of different orders, Math. Probl. Eng., 2022 (2022), 1879152. http://dx.doi.org/10.1155/2022/1879152 doi: 10.1155/2022/1879152
[14]	M. Abdo, S. Panchal, H. Shafei Hussien, Fractional integro-differential equations with nonlocal conditions and $\psi$ -Hilfer fractional derivative, Math. Model. Anal., 24 (2019), 564–584. http://dx.doi.org/10.3846/mma.2019.034 doi: 10.3846/mma.2019.034
[15]	H. Wahash, M. Abdo, S. Panchal, Fractional integrodifferential equations with nonlocal conditions and generalized Hilfer fractional derivative, Ufa Math. J., 11 (2019), 151–170. http://dx.doi.org/10.13108/2019-11-4-151 doi: 10.13108/2019-11-4-151
[16]	J. Sun, Y. Zhao, Multiplicity of positive solutions of a class of nonlinear fractional differential equations, Comput. Math. Appl., 49 (2005), 73–80. http://dx.doi.org/10.1016/j.camwa.2005.01.006 doi: 10.1016/j.camwa.2005.01.006
[17]	J. Zhou, S. Zhang, Y. He, Existence and stability of solution for nonlinear differential equations with $\psi$ -Hilfer fractional derivative, Appl. Math. Lett., 121 (2021), 107457. http://dx.doi.org/10.1016/j.aml.2021.107457 doi: 10.1016/j.aml.2021.107457
[18]	J. Zhou, S. Zhang, Y. He, Existence and stability of solution for a nonlinear fractional differential equation, J. Math. Anal. Appl., 498 (2021), 124921. http://dx.doi.org/10.1016/j.jmaa.2020.124921 doi: 10.1016/j.jmaa.2020.124921
[19]	X. Su, S. Zhang, Unbounded solutions to a boundary value problem of fractional order on the half-line, Comput. Math. Appl., 61 (2011), 1079–1087. http://dx.doi.org/10.1016/j.camwa.2010.12.058 doi: 10.1016/j.camwa.2010.12.058
[20]	X. Su, Solutions to boundary value problem of fractional order on unbounded domains in a Banach space, Nonlinear Anal.-Theor., 74 (2011), 2844–2852. http://dx.doi.org/10.1016/j.na.2011.01.006 doi: 10.1016/j.na.2011.01.006
[21]	C. Kou, H. Zhou, Y. Yan, Existence of solutions of initial value problems for nonlinear fractional differential equations on the half-axis, Nonlinear Anal.-Theor., 74 (2011), 5975–5986. http://dx.doi.org/10.1016/j.na.2011.05.074 doi: 10.1016/j.na.2011.05.074
[22]	J. Diaz, B. Margolis, A fixed point theorem of the alternative for contractions on a generalized complete metric space, Bull. Amer. Math. Soc., 74 (1968), 305–309. http://dx.doi.org/10.1090/S0002-9904-1968-11933-0 doi: 10.1090/S0002-9904-1968-11933-0
[23]	L. C $\breve{a}$ dariu, L. G $\breve{a}$ vruţa, P. G $\breve{a}$ vruţa, Weighted space method for the stability of some nonlinear equations, Appl. Anal. Discr. Math., 6 (2012), 126–139. http://dx.doi.org/10.2298/AADM120309007C doi: 10.2298/AADM120309007C

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Mathematics

1.8 3.4

Metrics

Article views(2086) PDF downloads(121) Cited by(9)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Mathematics

Existence and stability of solution for multi-order nonlinear fractional differential equations

Related Papers:

Abstract

1. Definition and notations

2. Appell-type Changhee-Euler polynomials

3. Differential equations

4. Graphical approach

5. Concluding remarks

Acknowledgments

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Catalog

AIMS Mathematics

Existence and stability of solution for multi-order nonlinear fractional differential equations

Related Papers:

Abstract

1. Definition and notations

2. Appell-type Changhee-Euler polynomials

3. Differential equations

4. Graphical approach

5. Concluding remarks

Acknowledgments

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog