Next Article in Journal
On Z-Intuitionistic Fuzzy Fractional Valuations for Medical Diagnosis: An Intuitionistic Fuzzy Knowledge-Based Expert System
Previous Article in Journal
Fractional Modeling Applied to the Dynamics of the Action Potential in Cardiac Tissue
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study of Generalized Hybrid Discrete Pantograph Equation via Hilfer Fractional Operator

by
Wafa Shammakh
1,
A. George Maria Selvam
2,
Vignesh Dhakshinamoorthy
2 and
Jehad Alzabut
3,4,*
1
Department of Mathematics, Faculty of Science, University of Jeddah, Jeddah 21493, Saudi Arabia
2
Department of Mathematics, Sacred Heart College (Autonomous), Tirupattur 635601, Tamil Nadu, India
3
Department of Mathematics and Sciences, Prince Sultan University, Riyadh 11586, Saudi Arabia
4
Department of Industrial Engineering, OSTİM Technical University, Ankara 06374, Turkey
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(3), 152; https://doi.org/10.3390/fractalfract6030152
Submission received: 12 February 2022 / Revised: 4 March 2022 / Accepted: 7 March 2022 / Published: 10 March 2022

Abstract

:
Pantograph, a device in which an electric current is collected from overhead contact wires, is introduced to increase the speed of trains or trams. The work aims to study the stability properties of the nonlinear fractional order generalized pantograph equation with discrete time, using the Hilfer operator. Hybrid fixed point theorem is considered to study the existence of solutions, and the uniqueness of the solution is proved using Banach contraction theorem. Stability results in the sense of Ulam and Hyers, and its generalized form of stability for the considered initial value problem are established and we depict numerical simulations to demonstrate the impact of the fractional order on stability.

1. Introduction

Differential equations with delay have a wide range of applications in science and engineering. Pantograph equations, in general referring to differential equations with proportional delays, have their origin in 1971 from the work of Ockendon and Tayler [1]. Since then, continuous and discrete versions of the equation have attracted a large number of physicists and mathematicians to study and analyze the importance of the equations in science, engineering and technology [2,3,4,5,6,7,8,9]. The generalized pantograph equation
υ ( ζ ) = A υ ( ζ ) + B υ ( δ ζ ) + C υ ( δ ζ ) ,
where 0 < ϑ < 1 , was proposed by Iserles [10]. Some notable fields of application of pantograph equations include number theory, cell growth, economy, quantum mechanics, biology, electrodynamics, chemical kinetics, physics, medicine and so on.
Fractional calculus was known to mathematicians and scientists for more than three centuries, but the application of fractional calculus in various fields of engineering, science and economics was understood only during last few decades. The reality of nature is better translated by fractional order calculus due its non-local distribution and memory effects. Though different fractional order derivatives and integrals like the Hadamard fractional derivative, Riesz derivative, and so on, are defined, only few derivatives like the Caputo fractional derivative and Riemann Liouville derivative are very widely used to model real life phenomena [11]. Generalization of the Riemann-Liouville and Caputo type derivatives into a fractional order derivative of type 0 μ 1 was proposed by R. Hilfer in [12]. The fractional derivative proposed by Hilfer yields Caputo ( μ = 1 ) and Riemann Liouville derivative ( μ = 0 ) as particular cases. Some recent works on fractional order Hilfer derivatives can be found in [13,14,15,16,17,18]. Perturbation techniques are powerful tools in nonlinear analysis for studying diverse aspects of the solution of nonlinear dynamical systems. They are useful in describing, predicting and demonstrating the nonlinear effects caused in vibrating systems. Hybrid fixed point theory is a common approach to tackle perturbed nonlinear equations. Hybrid differential equations are nonlinear differential equations with perturbation of the equation involving multiplication or division by a term (quadratic perturbation) [19]. Hybrid equations of fractional order have attracted researchers in recent times to the extent that they embrace various dynamic systems as particular cases [20,21,22,23,24,25]. The existence of results for the solution of the hybrid pantograph equation with fractional order, given by
D 0 + η υ ( ζ ) f ( ζ , υ ( ζ ) , υ ( δ ζ ) ) = Θ ( ζ , υ ( ζ ) , υ ( σ ζ ) ) , 0 < ζ < 1 , υ ( 0 ) = 0 ,
was studied in [2]. The generalized hybrid fractional pantograph equation
D 0 + η υ ( ζ ) f ( ζ , υ ( ζ ) , υ ( ϕ ( ζ ) ) ) = Θ ( ζ , υ ( ζ ) , υ ( ρ ( ζ ) ) ) , 0 < ζ < 1 , υ ( 0 ) = 0 ,
was considered to study the existence of solution in [26].
Discrete time fractional order calculus was enriched by the contributions of Atici et al. [27,28,29], Anastassiou [30], Goodrich [31,32], Holm [33], and so on. The definition of Hilfer fractional sum and differences are proposed in [34]. As for hybrid equations with discrete fractional operators, the authors in [35] considered the hybrid fractional sum-difference equations and investigated the existence of solutions. The qualitative properties of the discrete fractional hybrid equations are yet to be explored. The stability analysis of Hilfer type hybrid fractional equations has not been studied, to fill this gap, we consider an application of Hilfer fractional sum and difference to generalized hybrid pantograph equation and perform stability analysis.
The paper is formatted as follows: Essential definitions and lemmas are provided in Section 2 and the mathematical representation of the discrete time hybrid fractional pantograph equation is presented in Section 3. The existence of a unique solution for the Hilfer type discrete fractional generalized hybrid pantograph equation is illustrated in Section 4. Section 5 and Section 6 present the stability results and application of the main result is demonstrated with a numerical example, respectively.

2. Prerequisites

Definition 1
([31]). ( η t h Fractional Sum ) Let η > 0 and p : N λ R . Then the delta fractional sum of p is
Δ λ η p ( ζ ) = = λ ζ η h η 1 ( ζ , σ ( ) ) p ( ) ,
where h η ( ζ , ) = ( ζ ) ( η ) Γ ( η + 1 ) is the η t h fractional Taylor monomial, σ ( ) = + 1 and ζ ( η ) = Γ ( ζ + 1 ) Γ ( ζ η + 1 ) is the falling factorial function.
Definition 2
([28]). The η t h Riemann–Liouville type fractional difference of function p is defined by
Δ λ η p ( ζ ) = Δ m Δ λ ( m η ) p ( ζ ) , = = λ ζ + η h m η 1 ( ζ , σ ( ) ) p ( ) , ζ N λ + m η ,
where η > 0 ,   p : N λ R and m 1 < η < m for m N 1 .
Theorem 1
([36]). Let ρ , η 1 > 0 and p : N R . Then the inequalities
(1).
Δ η 1 [ Δ ρ p ( ζ ) ] = Δ ( ρ + η 1 ) p ( ζ ) = Δ ρ [ Δ η 1 p ( ζ ) ] ,
(2).
Δ η 1 Δ p ( ζ ) = Δ Δ η 1 p ( ζ ) ( ζ a ) ( η 1 1 ) Γ ( η 1 ) p ( a ) ,
hold.
Definition 3
([37]). The η t h Caputo type fractional difference of function p is defined by
x C Δ λ η p ( ζ ) = Δ λ ( m η ) Δ m p ( ζ ) , = = λ ζ ( m η ) h m η 1 ( ζ , σ ( ) ) Δ m p ( ) ζ N λ + m η ,
where η > 0 ,   p : N λ R and m 1 < η < m for m N 1 .
Definition 4
([34]). The Hilfer fractional difference of order m 1 < η < m and type 0 μ 1 of function p is defined by
Δ λ η , μ p ( ζ ) = Δ λ + ( 1 μ ) ( m η ) μ ( m η ) Δ Δ λ ( 1 μ ) ( m η ) p ( ζ ) , f o r ζ N λ + 1 η .
where p : N λ R .
Lemma 1
([34]). Assume 0 < η < 1 ,   0 μ 1 , and function p : N λ R , then for ζ N λ + 1 the composition properties are
(I).
Δ λ + 1 η η Δ λ η , μ p ( ζ ) = Δ λ + ( 1 μ ) ( 1 η ) ( η + μ η μ ) Δ Δ λ ( 1 μ ) ( 1 η ) p ( ζ ) .
(II).
Δ λ + 1 η η Δ λ η , μ p ( ζ ) = Δ λ + ( 1 μ ) ( 1 η ) ( η + μ η μ ) Δ λ ( η + μ η μ ) p ( ζ ) .
(III).
Δ λ + η η , μ Δ λ η p ( ζ ) = Δ λ + ( 1 μ + η μ ) μ ( 1 η ) Δ λ μ ( 1 η ) p ( ζ ) .
(IV).
Δ λ + η η , μ Δ λ η p ( ζ ) = p ( ζ ) Δ λ ( 1 μ ( 1 η ) ) p ( λ + 1 μ ( 1 η ) ) h μ ( 1 η ) 1 ( ζ , λ + 1 μ ( 1 η ) ) .
Theorem 2
([38]). (Banach Contraction Mapping Principle)
A contraction mapping on a complete metric space has exactly one fixed point.
Theorem 3
([20]). (Hybrid Fixed Point Theorem)
Let Y be the nonempty, closed, bounded and convex subset of Banach algebra B . Let the operators be P 1 : B B , P 2 : Y B such that
(i).
P 1 is Lipschitz continuous with constants θ.
(ii).
P 2 is completely continuous.
(iii).
υ = P 1 υ P 2 y υ Y y Y .
(iv).
θ A < 1 , where A = P 2 ( Y ) .
then, P 1 υ P 2 υ = υ has a solution.
Theorem 4
([39]). (Arzela-Ascoli Theorem) A set of functions in C ( [ a , b ] ) with supremum norm is relatively compact if, and only if, it is uniformly bounded and equicontinuous on [ a , b ] .

3. Discrete Fractional Hybrid Pantograph Equation

This section is devoted to the description of the hybrid pantograph model and approximate solution. Stability analysis of fractional hybrid equations by [40] and works on the hybrid pantograph equation by [2,26] have inspired and motivated us to investigate the stability of initial value Hilfer type discrete fractional generalized hybrid pantograph equation (HDFGHPE). Let us denote D = [ λ , T ] N λ . Let ϑ = η + μ η μ with 0 < η 1 ,   0 μ 1 . We have
Δ λ η , μ υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) = w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) , Δ λ ( 1 ϑ ) υ ( λ + 1 ϑ ) Θ ( λ + 1 ϑ , υ ( λ + 1 ϑ ) , υ ( ψ ( λ + 1 ϑ ) ) ) = Λ ,
where ζ [ λ + 1 η , T ] N λ + 1 η with λ R and T N where N j = j , j + 1 , , j R . Here Δ η , μ is Hilfer type fractional difference operator of order η and type μ , Λ R , ψ , φ : D [ 0 , 1 ] ,   Θ : D × R × R R 0 and w : D × R × R R are continuous functions.
Lemma 2.
Let 0 < η < 1 ,   0 μ 1 ,   Θ : D × R × R R 0 and w : D × R × R R . Then HDFGHPE (6) with initial condition has an unique solution
υ ( ζ ) = Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) [ Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) ] ,
for ζ D .
Proof. 
Applying Δ λ + 1 η η for (6), we have
Δ λ + 1 η η Δ λ η , μ υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) = Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) .
Apply Lemma (1) to left-hand side of (8) to obtain
Δ λ + 1 η η Δ λ η , μ υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) = Δ λ + ( 1 η ) ( 1 μ ) ( η + μ η μ ) Δ λ ( η + μ η μ ) υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) .
Now, consider
Δ λ + ( 1 η ) ( 1 μ ) ( η + μ η μ ) Δ λ ( η + μ η μ ) υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) = Δ λ + ( 1 ϑ ) ϑ Δ Δ λ ( 1 ϑ ) υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) .
Using the Theorem (1) yields
Δ λ + ( 1 η ) ( 1 μ ) ( η + μ η μ ) Δ λ ( η + μ η μ ) υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) = υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Λ h ϑ 1 ( ζ , λ + 1 ϑ ) .
Thus, Equation (7) holds. This completes the proof. □

4. Fixed Point Operators of HDFGHPE (6)

This section defines a fixed point operator and establishes the existence of an unique solution of (6). Let B with norm υ = sup ζ D υ ( ζ ) be the space of all functions υ ( ζ ) . Clearly, B is a Banach space.
Let us define the operator Υ : B B by
Υ υ ( ζ ) = Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) [ Λ h ϑ 1 ( ζ , λ + 1 ϑ ) +   Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) ] ,
for ζ D .
For all υ , υ * R and ζ D , we make the following assumptions
(J1):
There exists χ w ( ζ ) C ( D , R + ) such that
| w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) | χ w ( ζ ) .
(J2):
There exist ξ 1 ( ζ ) C ( D , R + ) with bound ξ 1 such that
| Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Θ ( ζ , υ * ( ζ ) , υ * ( ψ ζ ) ) | ξ 1 ( ζ ) max | υ ( ζ ) υ * ( ζ ) | , | υ ( ψ ( ζ ) ) υ * ( ψ ( ζ ) ) | .
(J3):
There exist ξ 2 ( ζ ) C ( D , R + ) with bound ξ 2 such that
| w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) w ( ζ + η 1 , υ * ( ζ + η 1 ) , υ * ( φ ( ζ + η 1 ) ) ) | ξ 2 ( ζ ) max | υ ( ζ ) υ * ( ζ ) | , | υ ( φ ( ζ ) ) υ * ( φ ( ζ ) ) | .

4.1. Uniqueness of Solution of HDFGHPE (6)

Theorem 5.
Assume that ( J 1 ) , ( J 2 ) , ( J 3 ) hold and there exists χ Θ ( ζ ) C ( D , R + ) such that
| Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) | χ Θ ( ζ ) ,
υ R and ζ D . Then υ ( ζ ) is a unique solution of hybrid fractional difference Equation (6) if
Ω = ξ 1 Λ h ϑ 1 ( T , λ + 1 ϑ ) + χ w ( T + η λ 1 ) ( η ) Γ ( η + 1 ) + χ Θ ξ 2 ( T + η λ 1 ) ( η ) Γ ( η + 1 ) < 1 .
Proof. 
Let χ w = max ζ D χ w ( ζ ) and χ Θ = max ζ D χ Θ ( ζ ) .
We aim to prove that the mapping Υ defined in (9) is a contraction. For υ , υ * R and ζ D , we have
| Υ υ ( ζ ) Υ υ * ( ζ ) | = | Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) [ Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) ] Θ ( ζ , υ * ( ζ ) , υ * ( ψ ( ζ ) ) ) [ Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ * ( ζ + η 1 ) , υ * ( φ ( ζ + η 1 ) ) ) ] | , | ξ 1 ( ζ ) | | h ϑ 1 ( ζ , λ + 1 ϑ ) | Λ + | Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Θ ( ζ , υ * ( ζ ) , υ * ( ψ ( ζ ) ) ) | Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) + | Θ ( ζ , υ * ( ζ ) , υ * ( ψ ( ζ ) ) ) | Δ λ + 1 η η | w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) w ( ζ + η 1 , υ * ( ζ + η 1 ) , υ * ( φ ( ζ + η 1 ) ) ) | , ξ 1 Λ h ϑ 1 ( T , λ + 1 ϑ ) + ( χ w ξ 1 + χ Θ ξ 2 ) ( T + η λ 1 ) ( η ) Γ ( η + 1 ) υ υ * , Υ υ Υ υ *   Ω υ υ * .
Thus, the mapping Υ is a contraction mapping with Ω < 1 . Therefore, the initial value discrete time Hilfer type fractional generalized pantograph Equation (6) has a unique solution. □

4.2. Existence Results for HDFGHPE (6)

The existence of the solution of HDFGHPE (6) is established using hybrid fixed point theorem (3) due to Dhage [20].
Theorem 6.
Assume that ( J 1 ) , ( J 2 ) hold. Then HDFGHPE (6) has a solution υ ( ζ ) for ζ D if
ξ 1 Λ ( T λ 1 + ϑ ) ( ϑ 1 ) Γ ( η ) + χ w Γ ( ϑ ) ( T λ 1 + η ) ( η ) < Γ ( ϑ ) Γ ( η ) .
Proof. 
Let Y =   υ B : υ S , where S is a real number such that
S Θ 0 | Ξ | 1 ξ 1 | Ξ | ,
where Θ 0 = max ζ D Θ ( ζ , 0 , 0 ) and Ξ = Λ | h ϑ 1 ( T , λ + 1 ϑ ) | + χ w ( T λ 1 + η ) ( η ) Γ ( η ) .
Clearly, Y is a closed, bounded and convex subset of B . By Lemma (2), let the operators P 1 : B B and P 2 : Y B be
P 1 υ ( ζ ) = Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) , P 2 υ ( ζ ) = Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) ,
where ζ D . Thus, solution of (6) is equivalent to
P 1 υ ( ζ ) P 2 υ ( ζ ) = υ ( ζ ) , ζ D .
We shall now prove that the operators P 1 and P 2 satisfy the conditions of the Theorem (3) in following steps.
Step 1:
We shall show that the operator P 1 is Lipschitz continuous on B .
From the condition J 1 , we have,
P 1 υ ( ζ ) P 1 υ * ( ζ ) = | Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Θ ( ζ υ * ( ζ ) , υ * ( ψ ( ζ ) ) ) | ,   | ξ 1 ( ζ ) | max | υ ( ζ ) υ * ( ζ ) | , | υ ( ψ ( ζ ) ) υ * ( ψ ( ζ ) ) | , P 1 υ P 1 υ *   ξ 1 υ υ * .
Thus, P 1 is lipschitz continuous with constants ξ 1 .
Step 2:
We proceed to prove that the operator P 2 is completely continuous on Y .
The continuity of w implies the continuity of operator P 2 on Y .
First, we shall prove the uniform boundedness of the operator P 2 in Y .
P 2 υ ( ζ ) =   Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) , Λ | h ϑ 1 ( ζ , λ + 1 ϑ ) | + Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) ,   | h ϑ 1 ( T , λ + 1 ϑ ) | Λ + χ w Γ ( T + η λ 1 ) η Γ ( η ) = ϖ , P 2 υ ϖ .
Therefore, P 2 is uniformly bounded on Y .
We prove the equicontinuity of the the operator P 2 . For any ε > 0 , let there exist ζ a , ζ b D ( ζ a < ζ b ) such that
Γ ( ζ b + η ) Γ ( ζ a ) Γ ( ζ b ) Γ ( ζ a + η ) Γ ( ζ b ) Γ ( ζ a + η )   <   Γ ( η + 1 ) Γ ( T ) χ w Γ ( T + η ) ε 2 ,
and
Γ ( ζ b λ + ϑ ) Γ ( ζ a λ + 1 ) Γ ( ζ b λ + 1 ) Γ ( ζ a λ ϑ ) 1   <   Γ ( T λ + 1 ) Λ Γ ( T λ ϑ ) ε 2 ,
then,
P 2 υ ( ζ b ) P 2 υ ( ζ a )   = | Λ h ϑ 1 ( ζ b , λ + 1 ϑ ) Λ h ϑ 1 ( ζ a , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ b + η 1 , υ ( ζ b + η 1 ) , υ ( φ ( ζ b + η 1 ) ) ) Δ λ + 1 η η w ( ζ a + η 1 , υ ( ζ a + η 1 ) , υ ( φ ( ζ a + η 1 ) ) ) | , Λ ( ζ b λ 1 + ϑ ) ( ϑ 1 ) Γ ( ϑ ) ( ζ b λ 1 + ϑ ) ( ϑ 1 ) Γ ( ϑ ) + 1 Γ ( η ) | = λ + 1 η ζ b η ( ζ b 1 ) ( η 1 ) 1 Γ ( η ) = λ + 1 η ζ a η ( ζ a 1 ) ( η 1 ) w ( + η , υ ( + η ) , υ ( σ ( + η ) ) ) | , Λ ( ζ b λ 1 + ϑ ) ( ϑ 1 ) Γ ( ϑ ) ( ζ b λ 1 + ϑ ) ( ϑ 1 ) Γ ( ϑ ) + χ w Γ ( η ) Γ ( ζ b + η ) Γ ( ζ b ) Γ ( ζ a + η ) Γ ( ζ a ) , Λ Γ ( T λ ϑ ) Γ ( T λ + 1 ) Γ ( ζ b λ + ϑ ) Γ ( ζ a λ + 1 ) Γ ( ζ b λ + 1 ) Γ ( ζ a λ ϑ ) 1 + χ w Γ ( T + η ) Γ ( T ) Γ ( η + 1 ) Γ ( ζ b + η ) Γ ( ζ a ) Γ ( ζ a + η ) Γ ( ζ b ) Γ ( ζ a + η ) Γ ( ζ b ) , P 2 υ ( ζ b ) P 2 υ ( ζ a ) < ε .
which implies the equicontinuity of P 2 in B . By Arzela Ascoli’s theorem, the operator P 2 is completely continuous.
Step 3:
We prove υ = P 1 υ P 2 y υ Y y Y .
Let υ B , y Y be arbitrary such that υ = P 1 υ P 2 y .
υ ( ζ )   | P 1 υ ( ζ ) | | P 2 y ( ζ ) | , Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) | Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , y ( ζ + η 1 ) , y ( φ ( ζ + η 1 ) ) ) | , Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Θ ( ζ , 0 , 0 ) + Θ ( ζ , 0 , 0 ) | Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , y ( ζ + η 1 ) , y ( φ ( ζ + η 1 ) ) ) | , | ξ 1 ( ζ ) | Λ | h ϑ 1 ( T , λ + 1 ϑ ) | + χ w ( T λ 1 + η ) ( η ) Γ ( η ) max | υ ( ζ ) | , | υ ( ψ ( ζ ) ) | + Θ 0 Λ | h ϑ 1 ( T , λ + 1 ϑ ) | + χ w ( T λ 1 + η ) ( η ) Γ ( η ) ,
Θ 0 | Ξ | 1 ξ 1 | Ξ | , υ   S .
Therefore υ ( ζ ) Y .
Step 4:
We show that θ A < 1 . Here,
A = sup ζ D P 2 ( Y ) , Λ | h ϑ 1 ( T , λ + 1 ϑ ) | + χ w ( T λ 1 + η ) ( η ) Γ ( η ) .
With θ = ξ 1 and (12), the condition θ A < 1 is satisfied.
Evidently, P 1 υ ( ζ ) P 2 y ( ζ ) = υ ( ζ ) , ζ D has a solution in Y which implies that HDFGHPE (6) has a solution υ ( ζ ) for ζ D . This completes the proof. □

5. Stability of HDFGHPE (6)

Stability is a condition in which trajectories of the system would not exhibit any significant changes under small disturbance. The asymptotic stability analysis of nonlinear discrete fractional equations were studied by Fulai Chen in [36,41]. Several authors have contributed on the stability analysis of various applications of fractional order discrete time equations as in [6,42,43,44,45,46,47]. We devote this section to study the stability of the HDFGHPE (6). Consider the discrete time Hilfer fractional initial value problem (6) and the following inequality
Δ λ η , μ υ 1 ( ζ ) Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) ) ε , ζ D ,
where ε > 0 and υ 1 C ( D , R ) .
Definition 5
([48] Ulam–Hyers Stability). The discrete time Hilfer fractional initial value problem (6) is Ulam–Hyers stable if there exists a real constant V > 0 such that for each ε > 0 and for every solution υ 1 C ( D , R ) of inequality (16), there exists a solution υ C ( D , R ) of (6) with
υ 1 ( ζ ) υ ( ζ ) < V ε , ζ D .
Definition 6
([48] Generalized Ulam–Hyers Stability). The HDFGHPE (6) is generalized Ulam–Hyers stable if σ C ( R + , R + ) with σ ( 0 ) = 0 exists such that for every solution υ 1 C ( D , R ) of inequality (16), there exists υ C ( D , R ) of (6) with
υ 1 ( ζ ) υ ( ζ )   < σ ( ε ) , ζ D .
The following remark is essential for proving the stability results.
Remark 1.
A function υ 1 ( ζ ) B solves (16) iff a function q : D × R R exists satisfying
(H1):
q ( ζ + η 1 , υ 1 ( ζ + η 1 ) )   ε , ζ D .
(H2):
Δ λ η , μ υ 1 ( ζ ) Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) )   = q ( ζ + η 1 , υ 1 ( ζ + η 1 ) ) .
Lemma 3.
If the function υ 1 ( ζ ) is the solution of the inequality (16) then
| υ 1 ( ζ ) Λ h ϑ 1 ( ζ , λ + 1 ϑ ) Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) Δ λ + 1 η η w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) ) | ε χ Θ ( T + η λ 1 ) ( η ) Γ ( η + 1 ) ,
for ζ D .
Proof. 
If υ 1 ( ζ ) satisfy (16), using Remark (1) and Lemma (1) the solution of ( H 2 ) is
υ 1 ( ζ ) = Θ ( ζ , υ ( ζ ) , υ 1 ( ψ ( ζ ) ) ) ( Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) ) + q ( ζ + η 1 , υ 1 ( ζ + η 1 ) ) )
where ζ D .
Hence,
| υ 1 ( ζ ) Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) [ Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) ) ] | = ( Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) ) Δ λ + 1 η η q ( ζ + η 1 , υ 1 ( ζ + η 1 ) ) , ε | χ Θ ( ζ ) | Γ ( η ) = λ + 1 η ζ η ( ζ 1 ) ( η 1 ) , ε χ Θ ( T + η λ 1 ) ( η ) Γ ( η + 1 ) .
This completes the proof. □
Theorem 7.
Assume that conditions ( J 1 ) , ( J 2 ) , ( J 3 ) and (10) hold. Let υ 1 B solve (16) and let υ B be the solution of
Δ λ η , μ υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) = 0 , 0 < η 1 , Δ λ ( 1 ϑ ) υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) ζ = λ + 1 ϑ = Δ λ ( 1 ϑ ) υ 1 ( ζ ) Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) ζ = λ + 1 ϑ .
Then (6) is Ulam-Hyers stable provided
Γ ( η + 1 ) ξ 1 Λ ( T λ 1 + ϑ ) ( ϑ 1 ) + Γ ( ϑ ) ( T + η λ 1 ) ( η ) [ ξ 1 χ w + χ Θ ξ 2 ] < Γ ( η + 1 ) Γ ( ϑ ) .
Proof. 
Using Lemma (2), we have
υ ( ζ ) = Δ λ ( 1 ϑ ) υ 1 ( λ + 1 ϑ ) Θ ( λ + 1 ϑ , υ 1 ( λ + 1 ϑ ) , υ 1 ( ψ ( λ + 1 ϑ ) ) ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) + Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) , ζ D .
Thus,
υ 1 ( ζ ) υ ( ζ )   = | υ 1 ( ζ ) Λ h ϑ 1 ( ζ , λ + 1 ϑ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) | , = | υ 1 ( ζ ) Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) [ Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) ) ] υ 1 ( ζ ) υ ( ζ ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) [ Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) ] + Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) ) | , | υ 1 ( ζ ) [ Λ h ϑ 1 ( ζ , λ + 1 ϑ ) + Δ λ + 1 η η ( w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) ) ) ] Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) | + Λ | h ϑ 1 ( ζ , λ + 1 ϑ ) | | Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) | + | Θ ( ζ , υ 1 ( ζ ) , υ 1 ( ψ ( ζ ) ) ) Δ η w ( ζ + η 1 , υ 1 ( ζ + η 1 ) , υ 1 ( φ ( ζ + η 1 ) ) ) Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Δ η w ( ζ + η 1 , υ ( ζ + η 1 ) , υ ( φ ( ζ + η 1 ) ) ) | ,
υ 1 υ   ε χ Θ ( T + η λ 1 ) ( η ) Γ ( η + 1 ) + [ ξ 1 Λ h ϑ 1 ( T , λ + 1 ϑ ) + χ w ( T + η λ 1 ) ( η ) Γ ( η + 1 ) + χ Θ ξ 2 ( T + η λ 1 ) ( η ) Γ ( η + 1 ) ] υ 1 υ , V ε ,
where V = χ Θ ( T + η λ 1 ) ( η ) Γ ( η + 1 ) ( 1 Ω ) with Ω defined in (11) is the stability constant. This completes the proof. The generalized Ulam–Hyers stability of the HDFGHPE (6) is established by replacing σ ( ε ) = V ε with σ ( 0 ) = 0 .

6. Numerical Examples

This section presents a numerical example to illustrate the results obtained in previous sections.
Example 1.
Let us consider the Hilfer type discrete fractional generalized hybrid pantograph equation of the form
Δ 0.3 0.6 , 0.5 10 υ ( ζ ) 3 + cos ( υ ( ζ ) ) = 1 30 1 2 + ( υ ( ζ 0.4 ) ) 3 Δ 0.2 10 υ ( 0.5 ) 3 + cos ( υ ( 0.5 ) ) = 0.9
where ζ [ 0 , 20 ] N 0.7 . We shall now establish the Ulam- Hyers stability of (21). Comparing (21) and (6), we have η = 0.6 , μ = 0.5 , ϑ = 0.8 , λ = 0.3 , Λ = 0.9 , T = 20 , D = [ 0.3 , 20 ] N 0.3 and
Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) = 0.3 + 1 10 cos ( υ ( ζ ) ) , w ( ζ 0.4 , υ ( ζ 0.4 ) , υ ( φ ( ζ 0.4 ) ) ) = 1 30 1 2 + ( υ ( ζ 0.4 ) ) 3 .
Let max ζ D υ ( ζ ) M 2 , where M R + . Using the conditions ( J 1 ) , ( J 2 ) , ( J 3 ) and (10), we yield
| w ( ζ 0.4 , υ ( ζ 0.4 ) , υ ( φ ( ζ 0.4 ) ) ) | 1 30 1 2 + ( M ) 3 8 , | Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) Θ ( ζ , υ * ( ζ ) , υ * ( ψ ( ζ ) ) ) | 0.1 | υ ( ζ ) υ * ( ζ ) | , | w ( ζ 0.4 , υ ( ζ 0.4 ) , υ ( φ ( ζ 0.4 ) ) ) w ( ζ 0.4 υ * ( ζ 0.4 ) , υ * ( φ ( ζ 0.4 ) ) ) | M 2 40 | υ ( ζ ) υ * ( ζ ) | , | Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) | 0.4
That is χ w = 1 30 1 2 + ( M ) 3 8 , χ Θ = 0.4 , ξ 1 = 0.1 , ξ 2 = M 2 40 , where M = 1 . From Theorem (5), we get Ω = 0.122794 < 1 .
Using (13), we have S 0.08044 with Θ 0 = 0.4 . Thus, (21) has a unique solution.
Ulam–Hyers stability of HDFGHPE (21) is evident from Theorem (7).
The stability condition (11) thus obtained for HDFGHPE (21) is tabulated for different fractional order η ( 0 , 1 ) in Table 1 and represented in Figure 1. For an increase in time, the stability condition ( Ω ) increases gradually for all values of η in ( 0 , 1 ) , and the Ω is clearly less than 1 satisfying the condition obtained in Theorem (5). An important observation to be made is when order ( η ) is small the value of Ω decreases with increasing time. As the order ( η ) increases, this trend changes with value of Ω increasing for increase in time.
The values in the Table 1 and Table 2 are obtained by substitution of the numerically calculated values that satisfy the conditions ( J 1 ) , ( J 2 ) , ( J 3 ) and (10) and employing the definition of falling factorial function in the inequality (11). As we already know, the Hilfer operator generalizes both the Riemann– Liouville and Caputo type operator for particular cases of μ . The analysis of stability condition on varying μ [ 0 , 1 ] is carried out and tabulated in Table 2 with presentation in Figure 2. The change in type μ of the operator between [ 0 , 1 ] results in increase of stability condition ( Ω ) as in Figure 2.
Example 2.
This example establishes the Hyers–Ulam stability of Hilfer type discrete fractional equation of the form
Δ 0.4 0.5 , 0.4 υ ( ζ ) = 1 10 υ ( ζ 0.5 ) + 0.01 sin 2 ( υ ( ζ 0.5 ) ) Δ 0.3 υ ( 0.7 ) = 0.6
where ζ [ 0 , 15 ] N 0.9 . Making an analogy between (22) and (6), we have η = 0.5 , μ = 0.4 , ϑ = 0.7 ,   λ = 0.4 , Λ = 0.6 , T = 15 , D = [ 0.4 , 15 ] N 0.4 and
Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) = 1 , w ( ζ 0.5 , υ ( ζ 0.5 ) , υ ( φ ( ζ 0.5 ) ) ) = 1 10 υ ( ζ 0.5 ) + 0.01 sin 2 ( υ ( ζ 0.5 ) ) .
Let max ζ D υ ( ζ ) M , where M R + . Using the conditions ( J 1 ) , ( J 2 ) , ( J 3 ) and (10), we yield
| w ( ζ 0.4 , υ ( ζ 0.4 ) , υ ( φ ( ζ 0.4 ) ) ) | M 10 + 0.01 , | w ( ζ 0.4 , υ ( ζ 0.4 ) , υ ( φ ( ζ 0.4 ) ) ) w ( ζ 0.4 υ * ( ζ 0.4 ) , υ * ( φ ( ζ 0.4 ) ) ) | 0.11 | υ ( ζ ) υ * ( ζ ) | , | Θ ( ζ , υ ( ζ ) , υ ( ψ ( ζ ) ) ) | 1
That is χ w = 0.06 , χ Θ = 1 , ξ 1 = 0 , ξ 2 = M 10 + 0.01 , where M = 0.5 . From Theorem (5), we get Ω = 0.470225 < 1 .
Using (13), we have S 0.256486 with Θ 0 = 1 . Thus, existence of unique solution for (22) is evident from the numerical calculations. Theorem (7) ensures the stability of Hilfer type discrete fractional problem (22) in the sense of Hyers and Ulam.

7. Conclusions

The article concentrated on obtaining the stability results generalizing Riemann–Liouville and Caputo type fractional derivatives in discrete time with the help of the Hilfer type discrete fractional operator. The application of the generalized discrete fractional operator (Hilfer) to a hybrid pantograph equation was considered in this work. Hybrid fixed point theory was used to develop the existence of a solution, and Banach contraction theorem was used to prove the uniqueness of the solution of the generalized fractional hybrid pantograph equation with discrete time. Stability analysis in the sense of Ulam–Hyers is performed, and the impact of the fractional order ( η ) and type ( μ ) are carried out as an example, with simulations present supporting the results obtained.

Author Contributions

Formal analysis and writing—review & editing, W.S.; Supervision and validation, A.G.M.S.; Writing—original draft and resources, V.D.; Project administration and methodology, J.A. All authors have read and agreed to the published version of the manuscript.

Funding

J. Alzabut is thankful to Prince Sultan University and OSTIM Technical University for their endless support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was funded by the University of Jeddah, Jeddah, Saudi Arabia, under grant No. (UJ-21-DR-37). The authors, therefore, acknowledge with the University of Jeddah for the technical and financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ockendon, J.; Tayler, A.B. The dynamics of a current collection system for an electric locomotive. Proc. R. Soc. Math. Phys. Eng. Sci. 1971, 322, 447–468. [Google Scholar]
  2. Darwish, M.A.; Sadarangani, K. Existence of solutions for hybrid fractional pantograph equations. Appl. Anal. Discret. Math. 2015, 9, 150–167. [Google Scholar] [CrossRef] [Green Version]
  3. Wongcharoen, A.; Ntouyas, S.K.; Tariboon, J. Nonlocal boundary value problems for Hilfer type pantograph fractional differential equations and inclusions. Adv. Differ. Equ. 2020, 2020, 279. [Google Scholar] [CrossRef]
  4. Derfel, G.; Iserles, A. The pantograph equation in the complex plane. J. Math. Anal. Appl. 1997, 213, 117–132. [Google Scholar] [CrossRef] [Green Version]
  5. Li, D.; Liu, M.Z. RungeKutta methods for the multi-pantograph delay equation. Appl. Math. Comput. 2005, 163, 383–395. [Google Scholar]
  6. Alzabut, J.; Selvam, A.G.M.; El-Nabulsi, R.A.; Vignesh, D.; Samei, M.E. Asymptotic Stability of Nonlinear Discrete Fractional Pantograph Equations with Non-Local Initial Conditions. Symmetry 2021, 13, 473. [Google Scholar] [CrossRef]
  7. Harikrishnan, S.; Elsayed, E.M.; Kanagarajan, K. Existence theory and Stability analysis of nonlinear neutral pantograph equations via Hilfer-Katugampola fractional derivative. J. Adv. Appl. Comput. Math. 2020, 7, 1–7. [Google Scholar] [CrossRef]
  8. Othman, M.I.A.; Said, S.; Marin, M. A novel model of plane waves of two-temperature fiber-reinforced thermoelastic medium under the effect of gravity with three-phase lag model. Int. J. Numer. Methods Heat Fluid Flow 2019, 29, 4788–4806. [Google Scholar] [CrossRef]
  9. Doha, E.H.; Bhrawy, A.H.; Baleanu, D.; Hafez, R.M. A new Jacobi rational-Gauss collocation method for numerical solution of generalized pantograph equations. Appl. Numer. Math. 2014, 77, 43–54. [Google Scholar] [CrossRef]
  10. Iserles, A. On the generalized pantograph functional differential equation. J. Appl. Math. 1993, 4, 1–38. [Google Scholar] [CrossRef]
  11. Podlubny, I. Geometric and physical interpretation of fractional integration and differentiation. Fract. Calc. Appl. Anal. 2002, 5, 367–386. [Google Scholar]
  12. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  13. Gu, H.; Trujillo, J.J. Existence of mild solution for evolution equation with Hilfer fractional derivative. Appl. Math. Comput. 2015, 257, 344–354. [Google Scholar] [CrossRef]
  14. Wang, J.R.; Zhang, Y. Nonlocal initial value problems for differential equations with Hilfer fractional derivative. Appl. Math. Comput. 2015, 266, 850–859. [Google Scholar] [CrossRef]
  15. Furati, K.M.; Kassim, N.D.; Tatar, N.E. Existence and uniqueness for a problem involving Hilfer fractional derivative. Comput. Math. Appl. 2012, 64, 1616–1626. [Google Scholar] [CrossRef] [Green Version]
  16. Ahmad, M.; Zada, A.; Alzabut, J. Hyers-Ulam stability of a coupled system of fractional differential equations of Hilfer—Hadamard type. Demonstr. Math. 2019, 52, 283–295. [Google Scholar] [CrossRef]
  17. Thaiprayoon, C.; Sudsutad, W.; Alzabut, J.; Etamed, S.; Rezapour, S. On the qualitative analysis of the fractional boundary value problem describing thermostat control model via ψ-Hilfer fractional operator. Adv. Differ. Equ. 2021, 2021, 201. [Google Scholar] [CrossRef]
  18. Rezazadeh, H.; Aminikhah, H.; Refahi Sheikhani, A. Stability Analysis of Hilfer fractional systems. Math. Commun. 2015, 21, 45–64. [Google Scholar]
  19. Dhage, B.C. Quadratic perturbation of periodic boundary value problems of second order ordinary differential equations. Differ. Equ. Appl. 2010, 2, 465–486. [Google Scholar] [CrossRef] [Green Version]
  20. Dhage, B.C. A fixed point theorem in Banach algebras with applications to fractional integral equations. Kyungpook Math. J. 2004, 44, 145–155. [Google Scholar]
  21. Ahmad, B.; Ntouyas, S.K. An existence theorem for fractional hybrid differential inclusions of Hadamard type with Dirichlet boundary conditions. Abstr. Appl. Anal. 2014, 2014, 705809. [Google Scholar] [CrossRef]
  22. Dhage, B.C.; Ntouyas, S.K. Existence results for boundary value problems for fractional hybrid differential inclusions. Topol. Methods Nonlinear Anal. 2014, 44, 229–238. [Google Scholar] [CrossRef]
  23. Sitho, S.; Ntouyas, S.K.; Tariboon, J. Existence results for hybrid fractional integro-differential equations. Bound. Value Probl. 2015, 2015, 1–13. [Google Scholar] [CrossRef] [Green Version]
  24. Bashiri, T.; Vaezpour, S.M.; Park, C. Existence results for fractional hybrid differential systems in Banach algebras. Adv. Differ. Equ. 2016, 2016, 1–13. [Google Scholar] [CrossRef] [Green Version]
  25. Herzallah, M.A.; Baleanu, D. On Fractional Order Hybrid Differential Equations. Abstr. Appl. Anal. 2014, 2014, 389386. [Google Scholar] [CrossRef]
  26. Karimov, E.T.; Lopez, B.; Sadarangani, K. About the existence of solutions for a hybrid nonlinear generalized fractional pantograph equation. Fract. Differ. Calc. 2016, 6, 95–110. [Google Scholar] [CrossRef] [Green Version]
  27. Atici, F.M.; Eloe, P.W. A transform method in discrete fractional calculus. Int. J. Differ. Equ. 2007, 2, 165–176. [Google Scholar]
  28. Atici, F.M.; Eloe, P.W. Initial value problems in discrete fractional calculus. Proc. Am. Math. Soc. 2009, 137, 981–989. [Google Scholar] [CrossRef]
  29. Atici, F.M.; Eloe, P.W. Two-point boundary value problems for finite fractional difference equations. J. Differ. Equ. Appl. 2011, 17, 445–456. [Google Scholar] [CrossRef]
  30. Anastassiou, G.A. Nabla discrete fractional calculus and nabla inequalities. Math. Comput. Model. 2010, 51, 562–571. [Google Scholar] [CrossRef] [Green Version]
  31. Goodrich, C.S.; Peterson, A.C. Discrete Fractional Calculus, 1st ed.; Springer: New York, NY, USA, 2015. [Google Scholar]
  32. Goodrich, C.S. Existence and uniqueness of solutions to a fractional difference equation with nonlocal conditions. Comput. Math. Appl. 2011, 61, 191–202. [Google Scholar] [CrossRef]
  33. Holm, M. Sum and difference compositions in discrete fractional calculus. Cubo 2011, 13, 153–184. [Google Scholar] [CrossRef] [Green Version]
  34. Haider, S.S.; Rehman, M.U.; Abdeljawad, T. On Hilfer fractional difference operator. Adv. Differ. Equ. 2020, 2020, 122. [Google Scholar] [CrossRef]
  35. Chasreechai, S.; Sitthiwirattham, T. Existence Results of Initial Value Problems for Hybrid Fractional Sum-Difference Equations. Discret. Dyn. Nat. Soc. 2018, 2018, 5268528. [Google Scholar] [CrossRef] [Green Version]
  36. Chen, F.; Liu, Z. Asymptotic stability results for nonlinear fractional difference equations. J. Appl. Math. 2012, 2012, 879657. [Google Scholar] [CrossRef]
  37. Abdeljawad, T. On Riemann and Caputo fractional differences. Comput. Math. Appl. 2011, 62, 1602–1611. [Google Scholar] [CrossRef] [Green Version]
  38. Kolmogorov, A.N.; Fomin, S.V. Elements of Function Theory and Functional Analysis; Nauka: Moscow, Russia, 1981. [Google Scholar]
  39. Griffel, D.H. Applied Functional Analysis; Ellis Horwood Publishers: Chichester, UK, 1981. [Google Scholar]
  40. Lachouri, A.; Ardjouni, A.; Djoudi, A. Existence and Ulam stability results for nonlinear hybrid implicit Caputo fractional differential equations. Math. Moravica 2020, 24, 109–122. [Google Scholar] [CrossRef]
  41. Chen, F. Fixed points and asymptotic stability of nonlinear fractional difference equations. Electron. J. Qual. Theory Differ. Equ. 2011, 39, 1–18. [Google Scholar] [CrossRef]
  42. Chen, C.; Bohner, M.; Jia, B. Ulam-Hyers stability of Caputo fractional difference equations. Math. Methods Appl. Sci. 2019, 42, 7461–7470. [Google Scholar] [CrossRef]
  43. Chen, F.; Zhou, Y. Existence of Ulam Stability of Solutions for Discrete Fractional Boundary Value problem. Discret. Dyn. Nat. Soc. 2013, 2013, 459161. [Google Scholar] [CrossRef]
  44. Selvam, G.M.; Alzabut, J.; Dhakshinamoorthy, V.; Jonnalagadda, J.M.; Abodayeh, K. Existence and stability of nonlinear discrete fractional initial value problems with application to vibrating eardrum. Math. Biosci. Eng. 2021, 18, 3907–3921. [Google Scholar] [CrossRef] [PubMed]
  45. Alzabut, J.; Selvam, A.G.M.; Vignesh, D.; Gholami, Y. Solvability and stability of nonlinear hybrid Δ-difference equations of fractional order. Int. J. Nonlinear Sci. Numer. Simul. 2021. [Google Scholar] [CrossRef]
  46. Tabouche, N.; Berhail, A.; Matar, M.M.; Alzabut, J.; Selvam, A.G.M.; Vignesh, D. Existence and stability analysis of solution for Mathieu fractional differential equations with applications on some physical phenomena. Iran. J. Sci. Technol. Transcation Sci. 2021, 45, 973–982. [Google Scholar] [CrossRef]
  47. Selvam, A.G.M.; Baleanu, D.; Alzabut, J.; Vignesh, D.; Abbas, S. On Hyers-Ulam Mittag-Leffler stability of discrete fractional duffing equation with applications on inverted pendulum. Adv. Differ. Equ. 2020, 2020, 1–15. [Google Scholar] [CrossRef]
  48. Rus, I.A. Ulam stabilities of ordinary differential equations in a Banach space. Carpathian J. Math. 2010, 26, 103–107. [Google Scholar]
Figure 1. Representation of impact of fractional order ( η ) on the stability condition Ω .
Figure 1. Representation of impact of fractional order ( η ) on the stability condition Ω .
Fractalfract 06 00152 g001
Figure 2. Representation of effect of change in μ [ 0 , 1 ] for different fractional order η ( 0 , 1 ) on stability condition Ω .
Figure 2. Representation of effect of change in μ [ 0 , 1 ] for different fractional order η ( 0 , 1 ) on stability condition Ω .
Fractalfract 06 00152 g002
Table 1. Representation of impact of fractional order η on stability condition Ω .
Table 1. Representation of impact of fractional order η on stability condition Ω .
Time η = 0.15 η = 0.30 η = 0.45 η = 0.60 η = 0.75 η = 0.90
Ω
0.30.090000.090000.089990.090000.090000.09000
1.30.063830.070580.077330.084080.090830.09758
2.30.054640.063970.073790.084130.094970.10632
3.30.049910.060690.072580.085610.099820.11526
4.30.046940.058770.072280.087610.104880.12422
5.30.044870.057530.072440.089840.109980.13315
6.30.043340.056710.072850.092160.115060.14202
7.30.042140.056130.073410.094520.120100.15084
8.30.041180.055740.074060.096890.125080.15959
9.30.040400.055460.074770.099250.129990.16829
10.30.039740.055280.075520.101580.134850.17693
11.30.039180.055160.076290.103900.139640.18552
12.30.038690.055100.077070.106180.144380.19406
13.30.038270.055080.077860.108430.149050.20255
14.30.037900.055090.078660.110660.153680.21100
15.30.037580.055120.079460.112850.158250.21940
16.30.037290.055180.080250.115020.162770.22776
17.30.037030.055250.081040.117160.167240.23609
18.30.036790.055340.081830.119270.171670.24437
19.30.036580.055430.082600.121350.176060.25262
Table 2. Effect of change in μ [ 0 , 1 ] for different fractional order η ( 0 , 1 ) on stability condition Ω .
Table 2. Effect of change in μ [ 0 , 1 ] for different fractional order η ( 0 , 1 ) on stability condition Ω .
μ η = 0.15 η = 0.30 η = 0.45 η = 0.60 η = 0.75 η = 0.90
Ω
00.021330.036460.060660.098600.157050.24566
0.20.024400.041110.067010.106280.164770.25050
0.40.031000.063970.076660.116510.173920.25567
0.60.043960.063410.091010.129980.184720.26117
0.80.067900.086330.111910.147580.197440.26704
1.00.110180.122740.141840.170370.212350.27329
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shammakh, W.; Selvam, A.G.M.; Dhakshinamoorthy, V.; Alzabut, J. A Study of Generalized Hybrid Discrete Pantograph Equation via Hilfer Fractional Operator. Fractal Fract. 2022, 6, 152. https://doi.org/10.3390/fractalfract6030152

AMA Style

Shammakh W, Selvam AGM, Dhakshinamoorthy V, Alzabut J. A Study of Generalized Hybrid Discrete Pantograph Equation via Hilfer Fractional Operator. Fractal and Fractional. 2022; 6(3):152. https://doi.org/10.3390/fractalfract6030152

Chicago/Turabian Style

Shammakh, Wafa, A. George Maria Selvam, Vignesh Dhakshinamoorthy, and Jehad Alzabut. 2022. "A Study of Generalized Hybrid Discrete Pantograph Equation via Hilfer Fractional Operator" Fractal and Fractional 6, no. 3: 152. https://doi.org/10.3390/fractalfract6030152

APA Style

Shammakh, W., Selvam, A. G. M., Dhakshinamoorthy, V., & Alzabut, J. (2022). A Study of Generalized Hybrid Discrete Pantograph Equation via Hilfer Fractional Operator. Fractal and Fractional, 6(3), 152. https://doi.org/10.3390/fractalfract6030152

Article Metrics

Back to TopTop