
We focused on the analytical solution of strong nonlinearity and complicated time-fractional evolution equations, including the Sawada-Kotera equation, Ito equation, and Kaup-Kupershmidt equation, using an effective and accurate method known as the Aboodh residual power series method (ARPSM) in the framework of the Caputo operator. Therefore, the Caputo operator and the ARPSM are practical for figuring out a linear or nonlinear system with a fractional derivative. This technique was effectively proposed to obtain a set of analytical solutions for various types of fractional differential equations. The derived solutions enabled us to understand the mechanisms behind the propagation and generation of numerous nonlinear phenomena observed in diverse scientific domains, including plasma physics, fluid physics, and optical fibers. The fractional property also revealed some ambiguity that may be observed in many natural phenomena, and this is one of the most important distinguishing factors between fractional differential equations and non-fractional ones. We also helped clarify fractional calculus in nonlinear dynamics, motivating researchers to work in mathematical physics.
Citation: Humaira Yasmin, Aljawhara H. Almuqrin. Efficient solutions for time fractional Sawada-Kotera, Ito, and Kaup-Kupershmidt equations using an analytical technique[J]. AIMS Mathematics, 2024, 9(8): 20441-20466. doi: 10.3934/math.2024994
[1] | D. Ajay, P. Chellamani, G. Rajchakit, N. Boonsatit, P. Hammachukiattikul . Regularity of Pythagorean neutrosophic graphs with an illustration in MCDM. AIMS Mathematics, 2022, 7(5): 9424-9442. doi: 10.3934/math.2022523 |
[2] | Supriya Bhunia, Ganesh Ghorai, Qin Xin . On the fuzzification of Lagrange's theorem in (α,β)-Pythagorean fuzzy environment. AIMS Mathematics, 2021, 6(9): 9290-9308. doi: 10.3934/math.2021540 |
[3] | Akarachai Satirad, Ronnason Chinram, Aiyared Iampan . Pythagorean fuzzy sets in UP-algebras and approximations. AIMS Mathematics, 2021, 6(6): 6002-6032. doi: 10.3934/math.2021354 |
[4] | Supriya Bhunia, Ganesh Ghorai, Qin Xin . On the characterization of Pythagorean fuzzy subgroups. AIMS Mathematics, 2021, 6(1): 962-978. doi: 10.3934/math.2021058 |
[5] | Mehwish Shehzadi, Hanan Alolaiyan, Umer Shuaib, Abdul Razaq, Qin Xin . Crafting optimal cardiovascular treatment strategy in Pythagorean fuzzy dynamic settings. AIMS Mathematics, 2024, 9(11): 31495-31531. doi: 10.3934/math.20241516 |
[6] | Muhammad Zeeshan Hanif, Naveed Yaqoob, Muhammad Riaz, Muhammad Aslam . Linear Diophantine fuzzy graphs with new decision-making approach. AIMS Mathematics, 2022, 7(8): 14532-14556. doi: 10.3934/math.2022801 |
[7] | Abdul Razaq, Ghaliah Alhamzi . On Pythagorean fuzzy ideals of a classical ring. AIMS Mathematics, 2023, 8(2): 4280-4303. doi: 10.3934/math.2023213 |
[8] | Damchaa Adiyanyam, Enkhbayar Azjargal, Lkhagva Buyantogtokh . Bond incident degree indices of stepwise irregular graphs. AIMS Mathematics, 2022, 7(5): 8685-8700. doi: 10.3934/math.2022485 |
[9] | Hifza, Muhammad Gulistan, Zahid Khan, Mohammed M. Al-Shamiri, Muhammad Azhar, Asad Ali, Joseph David Madasi . A new fuzzy decision support system approach; analysis and applications. AIMS Mathematics, 2022, 7(8): 14785-14825. doi: 10.3934/math.2022812 |
[10] | Saleem Abdullah, Muhammad Qiyas, Muhammad Naeem, Mamona, Yi Liu . Pythagorean Cubic fuzzy Hamacher aggregation operators and their application in green supply selection problem. AIMS Mathematics, 2022, 7(3): 4735-4766. doi: 10.3934/math.2022263 |
We focused on the analytical solution of strong nonlinearity and complicated time-fractional evolution equations, including the Sawada-Kotera equation, Ito equation, and Kaup-Kupershmidt equation, using an effective and accurate method known as the Aboodh residual power series method (ARPSM) in the framework of the Caputo operator. Therefore, the Caputo operator and the ARPSM are practical for figuring out a linear or nonlinear system with a fractional derivative. This technique was effectively proposed to obtain a set of analytical solutions for various types of fractional differential equations. The derived solutions enabled us to understand the mechanisms behind the propagation and generation of numerous nonlinear phenomena observed in diverse scientific domains, including plasma physics, fluid physics, and optical fibers. The fractional property also revealed some ambiguity that may be observed in many natural phenomena, and this is one of the most important distinguishing factors between fractional differential equations and non-fractional ones. We also helped clarify fractional calculus in nonlinear dynamics, motivating researchers to work in mathematical physics.
Fractional calculus has garnered significant attention across various fields of study because of its ability to capture non-local and long-range dependencies in data analysis and modeling. Fractional calculus provides a more flexible framework for capturing intricate dynamics in real-world systems. By incorporating non-integer orders, it becomes possible to model processes that exhibit fractal or self-similar behavior, which are prevalent in many natural and man-made systems [2,18]. This flexibility allows for a more accurate representation of complex phenomena, leading to improved understanding and enhanced predictive capabilities. Moreover, the application of fractional calculus can unveil hidden patterns and long-range dependencies that may go unnoticed when using integer-order operators [3,4,12,13]. In physics, fractional derivatives have been used to model anomalous diffusion [41], viscoelastic materials [17], and complex transport phenomena [20,43].
A wide range of real-world processes exhibit memory effects or non-local interactions, where the behavior of the system is significantly influenced by distant events or past observations [18]. For a detailed discussion, we refer to the monographs by Kilbas et al. [27], Lakshmikanthem et al. [28], Miller and Ross [29], Podlubny [31], Almeida et al. [11], and the references therein. Kilbas et al. [27] provided generalized Riemann-Liouville fractional operators to define the derivative with respect to a function. Several fractional derivatives concerning a function have been introduced, and their properties have been studied [5,8,9,14,24,34,37,38,39]. Katugampola expands upon the established framework of Erd'elyi-Kober operators, which are commonly used in the field of fractional derivatives. However, Katugampola introduced a specific modification and variation to these operators, which is referred to as the Hilfer-Katugampola fractional derivative [26]. These properties provide insights into the behavior of such fractional derivatives and their interactions with other mathematical operations.
Terminal value problems have applications in various fields where the ability to specify a desired final state is crucial to the problem at hand. For instance, in optimal control, terminal value problems are employed to determine the optimal control that minimizes a cost function over a given time frame. The terminal value serves as a constraint, allowing the specification of a particular final position or velocity for the system. Researchers have extensively investigated the existence and uniqueness of solutions to the terminal value problem, exploring the conditions under which a unique value can be determined for the system. For example, Benchohra et al. [16] investigated the existence and uniqueness of solutions for the terminal value problem in the Hilfer-Katugampola operator by using different types of classical fixed point theory. Almalahi et al. [6,7], developed sufficient conditions for the existence and uniqueness of solutions as well as stability results for the new terminal problems with respect to another function. Shah et al. [35] introduced some important notes on terminal value problems for fractional differential equations (FDEs) on infinite intervals. Abdo et al. [1] investigated the novel properties for the solution of a coupled system for terminal value problems involving a generalized Hilfer fractional operator. On the other hand, the symmetric solutions offer valuable insights into the behavior, stability, and simplification of mathematical and physical systems. They play a crucial role in simplifying analysis, facilitating physical interpretation, aiding in stability analysis, enhancing numerical computations, and contributing to the development of mathematical theories [42,44].
All of the aforementioned studies solely examined the properties of solutions within a single Banach space. However, our study seeks to expand upon this by investigating how the behavior of solutions is influenced by different weighted Banach spaces. Also, we aim to explore the impact of altering the function w on the solutions, while also examining multiple symmetric cases to observe the resulting variations in their behavior.
In this study, we investigate the existence and uniqueness of solutions for a system of terminal FDEs with w-Hilfer fractional derivatives and study their Ulam-Hyers (UH) stability [19,33,40]. Specifically, we study the fractional system
{Dpi,αa+,w(ι)ςi(ι)=Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)), ι∈(a,b],a>0,ςi(b)=υi∈R, i=1,2,....,n, | (1.1) |
where
● w∈C1([a,b],R) is an increasing function such that w′(ι)≠0, for all ι∈(a,b).
● Dpi,αa+,w(ι) denotes the w -Hilfer fractional derivatives of order pi,(0<pi<1) and type α,(0≤α≤1), for each i=1,2,....,n.
● The functions Fi:J×Rn→R,i=1,2,....,n, are continuous and satisfy some conditions that will be described later.
One of the primary motivations for studying terminal value problems with fractional derivatives is to analyze the stability properties of the corresponding dynamical systems. The w-Hilfer fractional derivative can capture complex behaviors such as subdiffusion or superdiffusion, which are not possible with classical integer-order derivatives. By formulating the problem in a weighted space, stability analysis can be performed in a manner that accounts for the importance or significance of different components or variables within the system. The innovative aspects of this work compared to the existing literature can be summarized as follows:
● We extend the analysis from single Banach spaces to different weighted Banach spaces, exploring how the change of weights affects solution properties. This allows us to capture additional structural information and understand the impact of weight functions on solutions.
● We investigate the sensitivity of solutions to changes in the weight function, providing insights into the relationship between the weighted functions and solution behavior by considering multiple symmetric cases to observe variations in solution behavior under different symmetry conditions.
● Furthermore, our study examines the system's properties in various Banach weighted spaces. We establish the existence, uniqueness, and stability of solutions for the system of w-Hilfer FDEs with minimal hypotheses, contributing to the theoretical understanding and providing a foundation for further analysis and applications.
● To the best of our knowledge, this is the first study to establish the existence, uniqueness, and stability of solutions for a system of w-Hilfer FDEs in different Banach-weighted spaces. Our approach expands the understanding of the system's behavior and broadens the range of potential applications.
● Additionally, the study addresses various parallel problems by considering specific function instances w. System (1.1) encompasses a fractional derivative that integrates several traditional fractional derivatives. This broadens the scope of the analysis and allows for a comprehensive investigation of the system's properties.
Notations: To enhance the readability, we fix the above notations and the following ones, we use these without any further mention.
(1) For 0<a<b, let J:=(a,b]⊂R.
(2)Γ(p) is the standard Gamma function given in [15].
(3)C(ˉJ,R) is the Banach space for all real-valued continuous functions defined on ˉJ, equipped with the norm ‖ς‖=maxι∈J|ς(ι)|, where ˉJ=[a,b]. The space ∏ni=1C(ˉJ,R) is the product space of n copies of C(ˉJ,R).
(4) For each i=1,2,...,n, let γi=pi+α(1−pi), such that pi∈(0,1) and α∈[0,1]. Then, we define the following spaces:
C1−γi,w(J,R)={ς:J→R;(w(ι)−w(a))1−γiςi(ι)∈C(ˉJ,R)}, |
and
Cγi1−γi,w(J,R)={ς∈C1−γi,w(J,R):Dγia+,wς∈C1−γi,w(J,R)}. |
The spaces C1−γi,w(J,R) are called weighted spaces. Clearly, these are Banach spaces equipped with the norm
‖ςi‖1−γi,w=maxι∈J|(w(ι)−w(a))1−γiςi(ι)|. |
(5)B=∏ni=1C1−γi,w(J,R) is the product space of C1−γi,w(J,R),i=1,2,....,n. The product spaces B are Banach spaces equipped with the norm
‖(ς1,ς2,.......,ςn)‖B=n∑i=1‖ςi‖1−γi. |
Hypotheses: To obtain our results, the following hypotheses must be satisfied:
(Hy1) For each (ι,ς1,ς2,....,ςn),(ι,ς∗1,ς∗2,....,ς∗n)∈J×Rn, there exist constants κjFi>0,i,j=1,2,....,n such that
|Fi(ι,ς1(ι),ς2(ι),....,ςn(ι))−Fi(ι,ς∗1(ι),ς∗2(ι),....,ς∗n(ι))|≤n∑j=1κjFi|ςj(ι)−ς∗j(ι)|. |
(Hy2) For each (ι,ς1,ς2,....ςn)∈J×Rn, there exist φjFi>0 with
|Fi(ι,ς1(ι),ς2(ι),....ςn(ι))|≤n∑j=1φjFi|ςj(ι)|. |
This paper is organized as follows. In Section 2, we state the needed fundamental notions of fractional derivatives, lemmas, and theorems that will be used in the conversion of the system into a fixed-point problem. In Section 3, we study the existence and uniqueness results of system (1.1), and discuss some properties of the symmetric cases of system (1.1). In Section 4, we investigate the UH stability of the solutions of system (1.1) and its symmetric cases. To illustrate our findings, Section 5 exhibits numerical examples with four cases and provides graphs of the solutions for all cases. We end the paper with concluding remarks.
In this section, the background definitions, lemmas, and theorems are given [8,27,30,37].
Definition 2.1. Let 0<p<1, ς be a continuous function on [a,b], and w be an increasing function on (a,b) such that w′(ι)≠0 for all ι∈(a,b). Then
(1) The w-Riemann-Liouville fractional integral for ς of order p is given as
Ipa+,w(ι)ς(ι)=1Γ(p)∫ιaw′(s)(w(ι)−w(s))p−1ς(s)ds. |
(2) The Riemann-Liouville fractional derivative for ς of order p with respect to w is given as
Dpa+,w(ι)ς(ι)=(1w′(ι)ddι)I1−pa+,w(ι)ς(ι). |
Definition 2.2. Let 0<p<1, 0≤α≤1, ι>a, and w be an increasing function on (a,b) such that w′(ι)≠0 for all ι∈(a,b). Then, the w−Hilfer fractional derivative of ς of order p and type α is defined by
Dp,αa+,w(ι)ς(ι)=Iα(1−p)a+,w(ι)(1w′(ι)ddι) I(1−α)(1−p)a+,w(ι) ς(ι). | (2.1) |
Clearly,
Dp,αa+,w(ι)ς(ι)=Iγ−pa+,w(ι) Dγa+,w(ι)ς(ι), |
where γ=p+α(1−p).
Lemma 2.3. For p,η,δ>0, we have
(1) Ipa+,w(ι) Iηa+,w(ι)ς(ι)=Ip+ηa+,w(ι)ς(ι).
(2) Ipa+,w(ι)(w(ι)−w(a))δ−1=Γ(δ)Γ(p+δ)(w(ι)−w(a))p+δ−1.
Note that Dγa+,w(ι)(w(ι)−w(a))γ−1=0, where γ=p+α(1−p).
Lemma 2.4. Let ς∈C(J,R), p∈(0,1) and α∈[0,1]. For γ=p+α(1−p), we have
Ipa+,w(ι) Dpa+,w(ι)ς(ι)=ς(ι)−(w(ι)−w(a))γ−1Γ(γ) I1−γa+,w(ι)ς(ι). |
Lemma 2.5. For each 0<p<1,0≤α≤1,γ=p+α(1−p), and ς∈Cγ1−γ;w(ι)(J,R), we have
Dp,αa+,w(ι)Ipa+,w(ι)ς(ι)=ς(ι), Iγa+,w(ι)Dγa+,w(ι)ς(ι)=Ipa+,w(ι) Dp,αa+,w(ι)ς(ι), |
and
Dγa+,w(ι)Ipa+,w(ι)ς(ι)=Dα(1−p)a+,w(ι)ς(ι). |
Definition 2.6. (Contraction mapping) Let (X,d) be a metric space, and Q:X→X be a mapping. Then, Q is a contraction mapping if there exists a constant 0≤L<1 such that for all ς,y∈X, the following inequality holds:
d(Q(ς),Q(y))≤L d(ς,y). |
The constant L is referred to as the contraction constant. If L=0, the mapping is called a strict contraction.
Theorem 2.7. (Banach fixed point theorem) Let (X,d) be a complete metric space, and Q:X→X be a contraction mapping. Then, the mapping Q has a unique fixed point in X.
Definition 2.8. Let a<b be real numbers, and X be the space of all continuous function h on [a,b]. The set Ω⊂X is said to be
(1) a uniformly bounded set if there exists M>0 with ‖h‖=supι∈[a,b]|h(ι)|<M, for every h∈Ω.
(2) an equicontinuous set if for every ε>0 there exists δ>0 such that
|ι−ι∗|<δ⟹|h(ι)−h(ι∗)|<ε, |
for all h∈Ω and for any ι,ι∗∈[a,b].
Theorem 2.9. (Arzelá-Ascoli's theorem) Let a<b be real numbers, and X be the space of all continuous function on [a,b]. The set Ω⊂X is relatively compact if it is a uniformly bounded and equicontinuous set.
Theorem 2.10. (Schauder's fixed-point theorem) Let X be a Banach space and Ω⊂X be a non-empty, compact, convex subset. If Q:Ω→Ω is a continuous mapping, and Q(Ω) is a relatively compact subset of X, then Q has at least one fixed point in Ω.
In this section, we establish and develop sufficient conditions for the existence and uniqueness as well as UH stability results for the system (1.1). We begin by stating an equivalent integrals for the considered system.
In the following theorem, we convert system (1.1) into equivalent integral equations.
Theorem 3.1. Let 0<pi<1, 0≤α≤1, and γi=pi+α(1−pi),i=1,2,....,n. If Fi(⋅,ς1(⋅),ς2(⋅),......,ςn(⋅))∈C1−γ,w(J,R), and ςi∈Cγi1−γi,w(J,R), then ςi satisfies system (1.1) if and only if ςi is given by:
ςi(ι)=(w(ι)−w(a)w(b)−w(a))γi−1[υi−Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)]+Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)), | (3.1) |
where
Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)=1Γ(pi)∫baw′(s)(w(b)−w(s))pi−1Fi(s,ς1(s),..,ςn(s))ds, |
and
Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))=1Γ(pi)∫ιaw′(s)(w(ι)−w(s))pi−1Fi(s,ς1(s),....,ςn(s))ds. |
Proof. Let ςi∈Cγi1−γi,w(J,R) be a solution of system (1.1). We prove that ςi is also a solution of Eq (3.1). From the definition of Cγi1−γi,w(J,R), and using the definition 2.2, we have
I1−γia+,w(ι)ςi(ι)∈C1−γi,w(J,R) and Dγia+,w(ι)ςi(ι)=D1a+,w(ι)I1−γia+,w(ι)ςi(ι)∈C1−γi,w(J,R). |
By the definition of the space Cn1−γi,w(J,R), it follows that
I1−γia+,w(ι)ςi(ι)∈C11−γi,w(J,R). | (3.2) |
Using Lemma 2.4, with p=γi, we obtain
(Iγia+,w(ι) Dγia+,w(ι)ςi)(ι)=ςi(ι)−(w(ι)−w(a))γi−1Γ(γi) limι→a(I1−γia+,w(ι)ςi)(ι). | (3.3) |
By hypothesis (ςi∈Cγi1−γi,w(J,R)) and Lemma 2.5, we have
(Iγia+,w(ι) Dγia+,w(ι)ςi)(ι)=Ipia+,w(ι) Dpi,αa+,w(ι)ςi(ι). |
Since Dpi,αa+,w(ι)ςi(ι)=Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)), we have
(Iγia+,w(ι) Dγia+,w(ι)ςi)(ι)=Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). | (3.4) |
Comparing Eqs (3.3) and (3.4), we see that
ςi(ι)=(w(ι)−w(a))γi−1Γ(γi)limι→a(I1−γia+,w(ι)ςi)(ι)+Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). | (3.5) |
Replacing ι with b in Eq (3.5), we get
ςi(b)=(w(b)−w(a))γi−1Γ(γi)limι→a(I1−γia+,w(ι)ςi)(ι)+Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). |
Using the terminal condition ςi(b)=υi, we get
(w(b)−w(a))γi−1Γ(γi)limι→a(I1−γia+,w(ι)ςi)(ι)+Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))=υi. |
Thus,
limι→a(I1−γia+,w(ι)ςi)(ι)=Γ(γi)(w(b)−w(a))γi−1[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))]. | (3.6) |
Putting Eq (3.6) into Eq (3.5), we get
ςi(ι)=((w(ι)−w(a))(w(b)−w(a)))γi−1[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))]+Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). |
Hence, ςi(ι) satisfies Eq (3.1).
Conversely, let ςi∈Cγi1−γi,w(J,R) be functions satisfying (3.1). We prove that ςi is also a solution of the system (1.1). Applying the operator Dγia+,w(ι) on both sides of Eq (3.1), we get
Dγia+,w(ι)ςi(ι)=Dγia+,w(ι)((w(ι)−w(a))(w(b)−w(a)))γi−1[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))]+Dγia+,w(ι)Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). | (3.7) |
By Lemma 2.3, we have Dγia+,w(ι)(w(ι)−w(a))γi−1=0, and hence Eq (3.7) becomes
Dγia+,w(ι)ςi(ι)=Dγia+,w(ι)Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). | (3.8) |
Then, from Lemma 2.5, the right side of Eq (3.8) becomes
Dγia+,w(ι)Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))=Dα(1−pi)a+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). |
Thus, Eq (3.8) becomes
Dγia+,w(ι)ςi(ι)=Dα(1−pi)a+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). |
From Eq (3.2), we have Dγia+,w(ι)ςi(ι)∈C1−γi,w(J,R), and hence, Eq (3.8) implies
Dγia+,w(ι)ςi(ι)=Dα(1−pi)a+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))∈C1−γi,w(J,R). | (3.9) |
As Fi(⋅,ς1(⋅),ς2(⋅),......,ςn(⋅))∈C1−γi,w(J,R), it follows that
I1−α(1−pi)a+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))∈C1−γi,w(J,R). | (3.10) |
From Eqs (3.9) and (3.10) and the definition of the space Cn1−γi,w(J,R), we get
I1−α(1−pi)a+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))∈C11−γi,w(J,R). |
Now, by applying the operator Iα(1−pi)a+,w(ι) on both sides of Eq (3.9) and using Lemma 2.4, we have
Iα(1−pi)a+,w(ι)Dγia+,w(ι)ςi(ι)=ςi(ι)−I1−α(1−pi)a+,w(ι)Fi(a,ς1(a),ς2(a),......,ςn(a))Γ(α(1−pi))(w(ι)−w(a))α(1−pi)−1=Dpi,αa+,w(ι)ςi(ι)=Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). | (3.11) |
From Eq (2.1), Eq (3.11) reduces to
Dpi,αa+,w(ι)ςi(ι)=Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). |
Thus, Eq (1.1) holds.
Let Q:B→B be an operator defined by Q=(Q1,Q2,.....,Qn), where
Qi(ς1,ς2,....,ςn)(ι)=(w(ι)−w(a)w(b)−w(a))γi−1[υi−Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)]+Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)). |
Theorem 3.2. Assume both (Hy1) and (Hy2) hold. For r=1,2,.....,n, let
Pr=n∑i=12(w(b)−w(a))1−γiφrFiΓ(γr)Γ(pi+γr)(w(b)−w(a))pi+γr−1, |
and
ζr=n∑i=12(w(b)−w(a))1−γiκrFiΓ(γr)Γ(pi+γr)(w(b)−w(a))pi+γr−1. |
If ζ=maxr{ζr}<1 and P=maxr{Pr}<1, then system (1.1) has at least one solution.
Proof. Let β≥11−P∑ni=1|υi|(w(b)−w(a))γi−1, and define a closed ball Sβ as
Sβ={(ς1,ς2,....ςn)∈B:‖(ς1,ς2,....ςn)‖B≤β}. |
To apply Schauder's fixed point theorem, we divide the proof into the following steps.
Step 1: Q(Sβ)⊂Sβ. For any (ς1,ς2,....ςn)∈Sβ, we have
‖Q(ς1,ς2,....ςn)‖B≤n∑i=1‖Qi(ς1,ς2,....ςn)‖C1−γi,w. | (3.12) |
From (Hy2), we get
|Qi(ς1,ς2,....ςn)(ι)|≤(w(ι)−w(a)w(b)−w(a))γi−1[|υi|+Ipia+,w(b)|Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|]+Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))≤(w(ι)−w(a)w(b)−w(a))γi−1[|υi|+ Ipia+,w(b)(n∑j=1φjFi|ςj(b)|)]+Ipia+,w(ι)(n∑j=1φjFi|ςj(ι)|). | (3.13) |
Taking into consideration that (ς1,ς2,....ςn)∈B, we have
Ipia+,w(ι)(n∑j=1φFji|ςj(ι)|)=n∑j=1Ipia+,w(ι)φjFi(w(ι)−w(a))γj−1|(w(ι)−w(a))1−γj(ςj(ι))|≤n∑j=1φjF1Γ(γj)Γ(pi+γj)(w(b)−w(a))pi+γj−1‖ςj‖C1−γj,w. | (3.14) |
From (3.13) and (3.14), it follows that
|(w(ι)−w(a))1−γiQi(ς1,ς2,....ςn)(ι)|≤(w(b)−w(a))1−γi[|υi|+n∑j=1φjFiΓ(γj)Γ(pi+γj)(w(b)−w(a))pi+γj−1‖ςj‖C1−γj,w]+(w(ι)−w(a))1−γin∑j=1φjFiΓ(γj)Γ(pi+γj)(w(b)−w(a))pi+γj−1‖ςj‖C1−γj,w. |
Since γi<1, then 1−γi>0, and hence (w(ι)−w(a))1−γi<(w(b)−w(a))1−γi. Then, we get
‖Qi(ς1,ς2,....ςn)‖C1−γi,w≤|υi|(w(b)−w(a))1−γi+2(w(b)−w(a))1−γi×n∑j=1φjFiΓ(γj)Γ(pi+γj)(w(b)−w(a))pi+γj−1‖ςj‖C1−γj,w. |
Thus,
‖Q(ς1,ς2,....ςn)‖B≤n∑i=1‖Qi((ς1,ς2,....ςn)‖C1−γi,w≤n∑i=1|υi|(w(b)−w(a))1−γi+n∑i=12(w(b)−w(a))1−γi×n∑j=1φjFiΓ(γj)Γ(pi+γj)(w(b)−w(a))pi+γj−1‖ςj‖C1−γj,w≤n∑i=1|υi|(w(b)−w(a))1−γi+P‖(ς1,ς2,....ςn)‖B≤β. |
Hence, Q(Sβ)⊂Sβ.
Step 2: Q is continuous. Let (ς1k,ς2k,....,ςnk) be a sequence in Sβ such that (ς1k,ς2k,....,ςnk)→(ς1,ς2,....,ςn) in Sβ as k→∞. Then, we have
‖Q(ς1k,ς2k,....,ςnk)(ι)−Q(ς1,ς2,....,ςn)(ι)‖B≤n∑i=1‖(Qi(ς1k,ς2k,....,ςnk)−Qi(ς1,ς2,....,ςn))(ι)‖C1−γi,w≤n∑i=1(w(b)−w(a))1−γiIpia+,w(ι)|Fi(ι,ς1k(ι),ς2k(ι),....,ςnk(ι))(b)−Fi(ι,ς1(ι),ς2(ι),....,ςn(ι))(b)|+n∑i=1(w(ι)−w(a))1−γiIpia+,w(ι)|Fi(ι,ς1k(ι),ς2k(ι),....,ςnk(ι))−Fi(ι,ς1(ι),ς2(ι),....,ςn(ι))|≤n∑i=12(w(b)−w(a))1−γin∑j=1κiFiΓ(γj)Γ(pi+γj)(w(b)−w(a))pi+γj−1‖ςjk−ςj‖C1−γj,w≤ζ‖(ς1k,ς2k,....,ςnk)−(ς1,ς2,....ςn)‖B. |
This implies that ‖Q(ς1k,ς2k,....,ςnk)−Q(ς1,ς2,....,ςn)‖B→0 as k→∞. So, Q is continuous.
Step 3: Q(Sβ)is equicontinuous.
For ι1,ι2∈J with ι1<ι2, and for any (ς1,ς2,....,ςn)∈Sβ, we have
|Q(ς1,ς2,....,ςn)(ι2)−Q(ς1,ς2,....,ςn)(ι1)|≤n∑i=1|Qi(ς1,ς2,....,ςn)(ι2)−Qi(ς1,ς2,....,ςn)(ι1)|≤n∑i=1(w(ι2)−w(a))γi−1−(w(ι1)−w(a))γi−1(w(b)−w(a))γi−1×[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)]+n∑i=1|Ipia+,w(ι2)Fi(ι2,ς1,ς2,....,ςn)−Ip1a+,w(ι1)Fi(ι1,ς1,ς2,....,ςn)|. |
Since Fi(⋅,ς1(⋅),....,ςn(⋅)) are continuous on J, there exist ξFi∈R such that
|Fi(⋅,ς1(⋅),....,ς1(⋅))|≤ξFi. |
Hence,
|Q(ς1,ς2,....,ςn)(ι2)−Q(ς1,ς2,....,ςn)(ι1)|≤n∑i=1(w(ι2)−w(a))γi−1−(w(ι1)−w(a))γi−1(w(b)−w(a))γi−1[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)]+n∑i=1|Ipia+,w(ι2)Fi(ι2,ς1(ι2),ς2(ι2),....,ςn(ι2))−Ipia+,w(ι1)Fi(ι1,ς1(ι1),ς2(ι1),....,ςn(ι1))|≤n∑i=1(w(ι2)−w(a))γi−1−(w(ι1)−w(a))γi−1(w(b)−w(a))γi−1[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)]+n∑i=11Γ(pi)∫ι1aw′(ι)[(w(ι2)−w(ι))pi−1−(w(ι1)−w(ι))pi−1]|Fi(ι,ς1(ι),ς2(ι),....,ςn(ι))|dι+n∑i=11Γ(pi)∫ι2ι1w′(ι)(w(ι2)−w(ι))pi−1|Fi(ι,ς1(ι),ς2(ι),....,ςn(ι))|dι≤n∑i=1(w(ι2)−w(a))γi−1−(w(ι1)−w(a))γi−1(w(b)−w(a))γi−1[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)]+n∑i=1ξFiΓ(pi)∫ι1aw′(ι)[(w(ι2)−w(ι))pi−1−(w(ι1)−w(ι))pi−1]dι+n∑i=1ξFiΓ(pi)∫ι2ι1w′(ι)(w(ι2)−w(ι))pi−1dι≤n∑i=1(w(ι2)−w(a))γi−1−(w(ι1)−w(a))γi−1(w(b)−w(a))γi−1[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)]+n∑i=1ξFiΓ(pi+1)[(w(ι2)−w(a))pi−(w(ι1)−w(a))pi]. | (3.15) |
Thus, we get
|Q(ς1,ς2,....,ςn)(ι2)−Q(ς1,ς2,....,ςn)(ι1)|≤n∑i=1(w(ι2)−w(a))γi−1−(w(ι1)−w(a))γi−1(w(b)−w(a))γi−1[υi−Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))(b)]+n∑i=12ξFiΓ(pi+1)[(w(ι2)−w(a))pi−(w(ι1)−w(a))pi]. |
Since w(ι) is an increasing function, then
|Q(ς1,ς2,....,ςn)(ι2)−Q(ς1,ς2,....,ςn)(ι1)|→0, as ι2→ι1. |
This implies that Q(Sβ) is equicontinuous. Since Q(Sβ) is also uniformly bounded, then by Arzelá-Ascoli's theorem, it is relatively compact. By Theorem 3.1, Q has a fixed point. Consequently, system (1.1) possesses at least one solution.
Theorem 3.3. Suppose that hypothesis (Hy1) holds. For r=1,2,...,n, let
ζr=n∑i=12(w(b)−w(a))1−γiκrFiΓ(γr)Γ(pi+γr)(w(b)−w(a))pi+γr−1. |
If ζ=maxr{ζr}<1, then system (1.1) has a unique solution.
Proof. To show that system (1.1) has a unique solution, we show that the operator Q is a contraction, where Φ as defined in Theorem 3.1.
For all ι∈J and (ς1,ς2,....,ςn),(ς∗1,ς∗2,....,ς∗n)∈B, we have
‖Q(ς1,ς2,....,ςn)−Q(ς∗1,ς∗2,....,ς∗n)‖B≤n∑i=1‖Qi(ς1,ς2,....,ςn)−Qi(ς∗1,ς∗2,....,ς∗n)‖C1−γi,w≤n∑i=1(w(b)−w(a))1−γiIpia+,w(b)|Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))−Fi(ι,ς∗1(ι),ς∗2(ι),....,ς∗n(ι))|+n∑i=1(w(ι)−w(a))1−γiIpia+,w(ι)|Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))−Fi(ι,ς∗1(ι),ς∗2(ι),....,ς∗n(ι))|≤n∑i=12(w(b)−w(a))1−γi(n∑j=1κjFiΓ(γj)Γ(pi+γj)(w(b)−w(a))pi+γj−1‖ςj−ς∗j‖C1−γj,w)≤ζ‖(ς1,ς2,....,ςn)−(ς∗1,ς∗2,....,ς∗n)‖B. |
Since ζ<1, the operator Q is a contraction map. Thus, by Theorem 2.7, system (1.1) has a unique solution.
In this subsection, we consider some symmetric cases of system (1.1):
Case1: w(ι)=ι. In this case, system (1.1) reduces to a system of Hilfer fractional equations in the following form
{Dpi,αa+,ιςi(ι)=Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)), ι∈J,ςi(b)=υi∈R,i=1,2,....,n, | (3.16) |
where Dpi,αa+,ι represents the Hilfer fractional derivative of order pi [23]. Let B1=∏ni=1C1−γi,ι(J,R), where C1−γi,ι is the weighted space defined as in the notation part and equipped with the norm
‖ςi‖1−γi=maxι∈J|(ι−a)1−γiςi(ι)|. | (3.17) |
Theorems 3.2 and 3.3 imply the following corollary.
Corollary 3.4. Assume that both (Hy1) and (Hy2) hold. For 1≤r≤n, let
P⋆r=n∑i=12(b−a)γi−1φrFiΓ(γr)Γ(pi+γr)(b−a)pi+γr−1, | (3.18) |
and
ζ⋆r=n∑i=12(b−a)γi−1κrFiΓ(γr)Γ(pi+γr)(b−a)pi+γr−1. | (3.19) |
If P⋆=maxr{P⋆r}<1, and ζ⋆=maxr{ζ⋆r}<1, then system (3.16) has at least one solution (ς1,ς2,....,ςn)∈B1, where ςi is given by
ςi(ι)=(ι−a)γi−1(b−a)γi−1[υi−Ipia+,bFi(ι,ς1(ι),ς2(ι),......,ςn(ι))]+Ipia+,ιFi(ι,ς1(ι),ς2(ι),......,ςn(ι)), |
where
Ipia+,bFi(ι,ς1(ι),ς2(ι),....,ςn(ι))=1Γ(pi)∫ba(b−s)pi−1Fi(s,ς1(s),ς2(s),....,ςn(s))ds, |
and
Ipia+,ιFi(ι,ς1(ι),ς2(ι),....,ςn(ι))=1Γ(pi)∫ιa(ι−s)pi−1Fi(s,ς1(s),ς2(s),....,ςn(s))ds. |
Case 2: w(ι)=logι. In this case, system (1.1) has the form
{Dpi,α1+,logιςi(ι)=Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)), ι∈J,a=1,ςi(b)=υi∈R,i=1,2,....,n, | (3.20) |
where, Dpi,α1+,logι is the Hilfer-Hadamard fractional derivative of order pi [22,32]. Let B2=∏ni=1C1−γi,logι(J,R), where C1−γi,logι is equipped with the norm
‖ςi‖1−γi,logι=maxι∈J|(logι)1−γiςi(ι)|. | (3.21) |
Theorems 3.2 and 3.3 then imply the following corollary.
Corollary 3.5. Assume that both (Hy1) and (Hy2) hold. For 1≤r≤n, let
P∗∗r=n∑i=12(logι)γi−1φrFiΓ(γr)Γ(pi+γr)(logι)pi+γr−1, | (3.22) |
and
ζ∗∗j=2(logι)γi−1n∑j=1κjFiΓ(γj)Γ(pi+γj)(logι)pi+γj−1. | (3.23) |
If P∗∗=maxr{P∗∗r}<1, and ζ∗∗=maxr{ζ∗∗r}<1, then, system (3.20) has at least one solution (ς1,ς2,....,ςn)∈B2, where ςi is given by
ςi(ι)=(logι)γi−1(logb)γi−1[υi−Ipi1+,logbFi(ι,ς1(ι),ς2(ι),......,ςn(ι))]+Ipi1+,logιFi(ι,ς1(ι),ς2(ι),......,ςn(ι)), |
where
Ipi1+,logbFi(ι,ς1(ι),ς2(ι),....,ςn(ι))=1Γ(pi)∫b11s(logbs)pi−1Fi(s,ς1(s),ς2(s),....,ςn(s))ds, |
and
Ipi1+,logιFi(ι,ς1(ι),ς2(ι),....,ςn(ι))=1Γ(pi)∫ι11s(logιs)pi−1Fi(s,ς1(s),ς2(s),....,ςn(s))ds. |
Case 3: w(ι)=ιq, where q>0. In this case, system (1.1) becomes a system of Hilfer-Katugampola in the following form.
{Dpi,αa+,ιqςi(ι)=Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)), ι∈J,a>0,ςi(b)=υi∈R,i=1,2,....,n, | (3.24) |
where Dpi,αa+,ιp is the Hilfer-Katugampola fractional derivative of order pi [26,30]. Let B3=∏ni=1C1−γi,ιq(J,R), where C1−γi,ιq is the weighted space equipped with the norm
‖ςi‖1−γi,ιq=maxι∈J|(ιq−aq)1−γiςi(ι)|. | (3.25) |
Theorems 3.2 and 3.3 imply the following corollary.
Corollary 3.6. Assume that (Hy1) and (Hy2) hold. For 1≤r≤n, let
P∗∗∗r=n∑i=12(bq−aq)γi−1φrFiΓ(γr)Γ(pi+γr)(bq−aq)pi+γr−1, | (3.26) |
and
ζ∗∗∗r=n∑i=12(bq−aq)γi−1κrFiΓ(γr)Γ(pi+γr)(bq−aq)pi+γr−1. | (3.27) |
If P∗∗∗=maxr{P∗∗∗r}<1, and ζ∗∗∗=maxr{ζ∗∗∗r}<1, then system (3.20) has at least one solution (ς1,ς2,....,ςn)∈B3, where ςi is given by
ςi(ι)=(ιq−aqq)γi−1(bq−aqq)γi−1[υi−Ipia+,bqFi(ι,ς1(ι),ς2(ι),......,ςn(ι))]+Ipia+,ιqFi(ι,ς1(ι),ς2(ι),......,ςn(ι)), |
where
Ipia+,bqFi(ι,ς1(ι),ς2(ι),....,ςn(ι))=1Γ(pi)∫baqsq−1(bq−sq)pi−1Fi(s,ς1(s),ς2(s),....,ςn(s))ds, |
and
Ipia+,ιqFi(ι,ς1(ι),ς2(ι),....,ςn(ι))=1Γ(pi)∫ιaqsq−1(ιq−sq)pi−1Fi(s,ς1(s),ς2(s),....,ςn(s))ds. |
Remark 3.7. There are many symmetric cases of the function w and parameter α that result in symmetric systems, some of which have already been discussed in the literature. For example.
1) For w(ι)=ι and α=0, system (1.1) is equivalent to the Riemann-Liouville system [27].
2) For w(ι)=ι and α=1, system (1.1) is equivalent to the Caputo system [27].
3) For w(ι)=ιq and α=1, system (1.1) is equivalent to the Caputo-Katugampola system [10].
4) For w(ι)=ιq and α=0, system (1.1) is equivalent to the Katugampola system [26].
5) For w(ι)=ι and α∈(0,1), system (1.1) is equivalent to the Hilfer system [23].
6) For w(ι)=ιq,q>0 and α∈(0,1) system (1.1) is equivalent to the Hilfer-Katugampola system [30].
7) For w(ι)=logι and α∈(0,1), system (1.1) is equivalent to the Hilfer-Hadamard system [32].
8) For w(ι)=logι and α=0), system (1.1) is equivalent to the Hadamard system [22].
9) For w(ι)=logι and α=1, system (1.1) is equivalent to the Caputo-Hadamard system [25].
The UH stability of system (1.1) is discussed in this section. Following the notation of this paper, we state the following definition for UH for system (1.1) [19,33,40].
Definition 4.1. System (1.1) is UH stable if there exists a positive constant Υ=max{Υ1,Υ2,...,Υn} such that for some ε=max{ε1,ε2,...,εn}>0 and for each solution (z1,z2,....,zn) of system (1.1) satisfying the following inequality
|Dpi,αa+,w(ι)zi(ι)−Fi(ι,z1(ι),z2(ι),......,zn(ι))|≤εi, | (4.1) |
there exists a unique solution (ς1,ς2,....,ςn) for system (1.1) such that
‖(z1,z2,....,zn)−(ς1,ς2,....,ςn)‖B≤Υε. | (4.2) |
Remark 4.2. The functions (z1,z2,....,zn) are a solution of inequality (4.1) if and only if there exist functions σi(ι)∈C1−γi,w,i=1,2,....,n (which depends on solution zi) such that
● |σi(ι)|≤εi.
● Dpi,αa+,w(ι)zi(ι)=Fi(ι,z1(ι),z2(ι),......,zn(ι))+σi(ι).
Theorem 4.3. Under hypothesis (Hy1), system (1.1) is UH stable, provided that
2κiFiΓ(γi)Γ(pi+γi)(w(b)−w(a))pi<1, |
for all 1≤i≤n.
Proof. Let ε=max{ε1,ε2,....,εn}>0 and let (z1,z2,....,zn)∈B be a solution of system (1.1) satisfying inequality (4.1), and let (ς1,ς2,....,ςn)∈B be the unique solution of system (1.1). By virtue of Theorem 3.1, we obtain
ςi(ι)=Aςi−((w(ι)−w(a))(w(b)−w(a)))γi−1Ipia+,w(b)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))+Ipia+,w(ι)Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)), | (4.3) |
where
Aςi=υi((w(ι)−w(a))(w(b)−w(a)))γi−1, | (4.4) |
since, we have assumed that (z1,z2,....,zn)∈B is a solution of system (1.1) satisfying the inequality (4.1). Then, by Remark 4.2, we conclude that (z1,z2,....,zn) satisfies the following system
{Dpi,αa+,w(ι)zi(ι)=Fi(ι,z1(ι),z2(ι),......,zn(ι))+σi(ι), ι∈J,zi(b)=υi∈R,i=1,2,....,n. | (4.5) |
By Theorem 3.1, the solution zi(ι) of system (4.5) is given as
zi(ι)=Azi−((w(ι)−w(a))(w(b)−w(a)))γi−1Ipia+,w(b)[Fi(ι,z1(ι),z2(ι),...,zn(ι))+σi(ι)]+Ipia+,w(ι)[Fi(ι,z1(ι),z2(ι),...,zn(ι))+σi(ι)], |
where
Azi=υi((w(ι)−w(a))(w(b)−w(a)))γi−1. | (4.6) |
Hence, we have
|zi(ι)−ςi(ι)|≤|Azi−Aςi|+((w(ι)−w(a))(w(b)−w(a)))γi−1Ipia+,w(b)|σi(ι)|+Ipia+,w(ι)|σi(ι)|+((w(ι)−w(a))(w(b)−w(a)))γi−1Ipia+,w(b)|Fi(ι,z1(ι),z2(ι),......,zn(ι))−Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|+Ipia+,w(ι)|Fi(ι,z1(ι),z2(ι),......,zn(ι))−Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|. | (4.7) |
By (4.4) and (4.6), we have Azi=Aςi. Since the function w is increasing for each ι∈J and γi<1 for i=1,2,...n, we get
((w(ι)−w(a))(w(b)−w(a)))γi−1<1. |
Hence, Eq (4.7) becomes
|zi(ι)−ςi(ι)|≤Ipia+,w(b)|σi(ι)|+Ipia+,w(ι)|σi(ι)|+Ipia+,w(b)|Fi(ι,z1(ι),z2(ι),......,zn(ι))−Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|+Ipia+,w(ι)|Fi(ι,z1(ι),z2(ι),......,zn(ι))−Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|. |
By the first part of Remark 4.2 and (Hy1), we get
|zi(ι)−ςi(ι)|≤εi2(w(b)−w(a))piΓ(pi+1)+2Ipia+,w(b)n∑j=1κjFi|zj(ι)−ςj(ι)|. |
Thus, for i=1, we have
|z1(ι)−ς1(ι)|≤ε12(w(b)−w(a))p1Γ(p1+1)+2Ip1a+,w(b)n∑j=1κjF1|zj(ι)−ςj(ι)|. |
Since ‖ςi‖1−γi,w=maxι∈J|(w(ι)−w(a))1−γiςi(ι)|, and w is increasing function, then
‖z1−ς1‖1−γ1,w≤ε12(w(b)−w(a))p1Γ(p1+1)+2n∑j=1(w(b)−w(a))1−γ1Ip1a+,w(ι)κjF1|zj(ι)−ςj(ι)|≤ε12(w(b)−w(a))p1Γ(p1+1)+2κ1F1Γ(γ1)Γ(p1+γ1)(w(b)−w(a))p1‖z1−ς1‖C1−γ1,w+2n∑j=2κjF1Γ(γj)Γ(p1+γj)(w(b)−w(a))p1+γj−γ1‖zj−ςj‖C1−γj,w. |
Hence, we have
(1−2κ1F1Γ(γ1)Γ(p1+γ1)(w(b)−w(a))p1)‖z1−ς1‖1−γ1,w≤ε12(w(b)−w(a))p1+1−γ1Γ(p1+1)+2n∑j=2κjF1Γ(γj)Γ(p1+γj)(w(b)−w(a))p1+γj−γ1‖zj−ςj‖C1−γj,w≤ε12(w(b)−w(a))p1+1−γ1Γ(p1+1)+2ε1n∑j=2κjF1Γ(γj)Γ(p1+γj)(w(b)−w(a))p1+γj−γ1‖zj−ςj‖C1−γj,w. |
Note that the positive values of the Gamma function on (0,1] are necessary to obtain the last inequality.
Let c1=1−2κ1F1Γ(γ1)Γ(p1+γ1)(w(b)−w(a))p1, and
Υ1=2(w(b)−w(a))p1Γ(p1+1)+2n∑j=2κjF1Γ(γj)Γ(p1+γj)(w(b)−w(a))p1+γj−γ1‖zj−ςj‖C1−γj,w. |
If 2κ1F1Γ(γ1)Γ(p1+γ1)(w(b)−w(a))p1<1, then c1>0. Then, we obtain
‖z1−ς1‖1−γ1,w≤ε1c1Υ1. |
Similarly, for each i=2,3,..,n, provided that 2κiFiΓ(γi)Γ(pi+γi)(w(b)−w(a))pi<1, we have
‖zi−ςi‖1−γi,w≤εiciΥi, |
where ci=1−2κiFiΓ(γi)Γ(pi+γi)(w(b)−w(a))pi, and
Υi=2(w(b)−w(a))piΓ(pi+1)+n∑j≠iκjFiΓ(γj)Γ(pi+γj)(w(b)−w(a))pi+γj−γi‖zj−ςj‖C1−γj,w. |
Therefore, if 2κiFiΓ(γi)Γ(pi+γi)(w(b)−w(a))pi<1,i=1,2,..n, then
‖(z1,z2,....,zn)−(ς1,ς2,....,ςn)‖B≤n∑i=1‖zi−ςi‖1−γi,w≤n∑i=1εiciΥi≤εΥ, |
where ε=maxi{εi} and Υ=∑ni=1Υici. By Definition 4.1, the solution of system (1.1) is UH stable.
According to Theorem 4.3, we can easily prove the UH stability of symmetric systems as follows.
(1) Under hypothesis (Hy1), system (3.16) is UH stable, provided that
2κiFiΓ(γi)Γ(pi+γi)(b−a)pi<1, |
for all 1≤i≤n.
(2) Under hypothesis (Hy1), system (3.20) is UH stable, provided that
2κiFiΓ(γi)Γ(pi+γi)(logba)pi<1, |
for all 1≤i≤n.
(3) Under the hypothesis (Hy1), the system (3.24) is UH stable, provided that
2κiFiΓ(γi)Γ(pi+γi)(bq−aq)pi<1, |
for all 1≤i≤n.
In this section, we consider an example of a system of w-Hilfer FDEs with different functions w.
Example 5.1. We study the system
{D12,120,w(ι)ςi(ι)=118(52ι)[1+∑nj=1|ςj(ι)|], ι∈(0,1],ςi(1)=3i2, i=1,2,...,n, | (5.1) |
for the cases w(ι)=eι,w(ι)=ι,w(ι)=logι, and w(ι)=ιq.
Here, pi=12,αi=12,γi=34,a=0,b=1,υi=3i2, and
Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))=118(52ι)[1+∑nj=1|ςj(ι)|]. |
Clearly, Fi(ι,ς1(ι),ς2(ι),......,ςn(ι)) are continuous functions on [0,1]×Rn and satisfy
|Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|≤118, |
and
|Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))−Fi(ι,z1(ι),z1(ι),......,zn(ι))|≤118n∑j=1|ςj−zj|, |
for all (ς1,ς2,....ςn),(z1,z2,....zn), i.e., φjFi=κjFi=118.
We fix the number of equations to be n=5, and investigate the following cases:
Case 1: w(ι)=eι. In this case, we have P=max{P1,...,P5}≃0.75<1, and ζ=max{ζ1,...,ζ5}=0.75<1, where Pr and ζr are given in Theorems 3.2 and 3.3, respectively.
Hence, all conditions in Theorem 3.2 and Theorem 3.3 are satisfied. Thus, the system has a unique solution (ς1,ς2,...,ς5). This solution has the form
ςi(ι)=(eι−1)−14(e−1)−14[3i2−1Γ(12)∫10es(e−es)−1218(52s)[1+∑5j=1|ςj(s)|]ds]+1Γ(12)∫ι0es(eι−es)−1218(52s)[1+∑5j=1|ςj(s)|]ds. |
On the other hand, let ι∈(0,1] and εi>0,i=1,2,...,n satisfying the following inequality
|D12,120+,eιςi(ι)−Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|≤εi. |
Then, system (1.1) is UH stable with
‖(z1,z2,....,zn)−(ς1,ς2,....,ςn)‖B≤n∑i=1εiciΥi, |
where
Υi=2(eι−1)12Γ(12+1)+n∑j≠iκjFiΓ(γj)Γ(12+γj)(e−1)12+γj−γi>0, |
and
ci=1−2118Γ(34)Γ(12+34)(e−1)12=0.8>0. |
Symmetric cases:
Case 2: w(ι)=ι. In this case, we have P∗=max{P∗1,...,P∗5}≃0.37<1, and ζ∗=maxι∈J{ζ∗1,...,ζ∗5}=0.37<1, where P∗r and ζ∗r are given in (3.18) and (3.19), respectively.
Hence, all conditions in Theorem 3.2 and Theorem 3.3 are satisfied. Thus, the system has at least one solution (ς1,ς2,...,ς5) given by
ςi(ι)=ι−14[3i2−1Γ(12)∫10(1−s)−1218(52s)[1+∑5j=1|ςj(s)|]ds]+1Γ(12)∫ι0(ι−s)−1218(52s)[1+∑5j=1|ςj(s)|]ds. | (5.2) |
On the other hand, let ι∈(0,1] and εi>0,i=1,2,...,n, satisfying the following inequality,
|D12,120+,ιςi(ι)−Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|≤εi. |
Then, system (3.16) is UH stable with
‖(z1,z2,....,zn)−(ς1,ς2,....,ςn)‖B≤n∑i=1εiciΥi, |
where
Υi=2Γ(12+1)+n∑j≠iκjFiΓ(γj)Γ(12+γj)>0, |
and
ci=1−2118Γ(34)Γ(12+34)=0.9>0. |
Case 3: w(ι)=logι. For this case, we need to change the interval J, as the function w is not defined on (0,1]. So, in this case, only a and b will be changed to be a=1,b=2, P∗∗=max{P∗∗1,...,P∗∗5}≃0.75<1, and ζ∗∗=max{ζ∗∗1,...,ζ∗∗5}=0.76<1, where P∗∗r and ζ∗∗r are defined in (3.22) and (3.23), respectively. Hence, all conditions in Theorem 3.2 and Theorem 3.3 are satisfied. Thus, the system has at least one solution (ς1,ς2,...,ς5) given by
ςi(ι)=(logι)−14(log2)−14[3i2−1Γ(12)∫21s−1(log2−logs)−1218(52s)[1+∑5j=1|ςj(s)|]ds]+1Γ(12)∫ι1s−1(logι−logs)−1218(52s)[1+∑5j=1|ςj(s)|]ds. | (5.3) |
On the other hand, let ι∈(1,2] and εi>0,i=1,2,...,n, satisfying the following inequality.
|D12,120+,logιςi(ι)−Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|≤εi. |
Then, system (3.20) is UH stable with
‖(z1,z2,....,zn)−(ς1,ς2,....,ςn)‖B≤n∑i=1εiciΥi, |
where
Υi=2(log2−log1)12Γ(12+1)+n∑j≠iκjFiΓ(γj)Γ(12+γj)(log2−log1)12+γj−γi>0, |
and
ci=1−2118Γ(34)Γ(12+34)(log2−log1)12=0.9>0. |
Case 4: w(ι)=ιq,q>0. We also consider this case on (1,2], even though we can consider the system on [0,1], where it will be generalizing case 2. In this case, we have P=max{P∗∗∗1,..,P∗∗∗5}≃0.75<1, and ζ=max{ζ∗∗∗1,..,ζ∗∗∗5}=0.75<1, where P∗∗∗r and ζ∗∗∗r are defined in (3.26) and (3.27), respectively. Hence, all conditions in Theorem 3.2 are satisfied. Thus, the system has at least one solution (ς1,ς2,...,ςn) given by
ςi(ι)=(ιq−1)−14(2q−1)−14[3i2−1Γ(12)∫21qsq−1(2q−sq)−1218(52s)[1+∑5j=1|ςj(s)|]ds]+1Γ(12)∫ι1qsq−1(ιq−sq)−1218(52s)[1+∑5j=1|ςj(s)|]ds. | (5.4) |
Let q=2 and εi>0,i=1,2,...,n, satisfying the following inequality.
|D12,120+,ι2ςi(ι)−Fi(ι,ς1(ι),ς2(ι),......,ςn(ι))|≤εi. |
Then, system (3.24) is UH stable with
‖(z1,z2,....,zn)−(ς1,ς2,....,ςn)‖B≤n∑i=1εiciΥi, |
where
Υi=2(22−1)12Γ(12+1)+n∑j≠iκjFiΓ(γj)Γ(12+γj)(22−1)12+γj−γi>0, |
and
ci=1−2118Γ(34)Γ(12+34)(22−1)12=0.7>0. |
We fix i=2 and present the graphical presentations of solutions for the above four cases in Figures 1, 2, 3, and 4 respectively as follows:
In Figure 1, we have presented the graphical presentation of Case I using i=2 to investigate the properties of the system under the Banach weighted spaces and provide insights into the system's response by using the weighted measurable function w(ι)=eι over the interval [0,1]. From Figure 1, the UH stability is obvious and the symmetry in the approximation solution also arises. Moreover, the uniqueness of the solutions here clearly corresponds to the coupled system of two equations. Remember, we use for graphical presentation (1.1) using i=2.
In Figure 2, we have presented the approximate solution for the system using i=2 to investigate the properties of the system under the Banach weighted spaces and provide insights into the system's response by using the weighted measurable function w(ι)=ι over the interval [0,1]. The graphical illustration for the approximate solution demonstrates the symmetry in both solutions for the given weighted measurable function. Here, we have also used (5.2) to illustrate graphs for the approximate solution using i=2.
In Figure 3, the approximate solutions using the weighted measurable function w(ι)=logι in Banach weighted spaces are displayed by selecting the interval [1,2]. We considered again (5.3) using i=2 and demonstrated the approximate solutions to elaborate the behaviors of both graphs. From the graphs, we see that, as logι is considered, we have displayed the approximate solutions over the interval [1,2]. The UH stability can be observed in the graphs and the symmetry also is understandable. There is a little bit of difference in both graphs' behaviors is due to the reason that we do not know the exact solution as the problem is very highly nonlinear.
In Figure 4, we have displayed the approximate solution of (5.4) using i=2 by utilizing the weighted measurable function w(ι)=ιq in Banach weighted spaces. The solutions for the coupled system have been plotted over [1,2]. From the graphs, we see the symmetric behaviors and UH stability. From these discussions, we conclude that symmetry concepts are significant when it comes to natural laws. Further, it is highly interesting that symmetry is the key to nature, but symmetry-breaking processes are mostly responsible for the texture of the universe. There are many different ways that nature's symmetry can be obscured or disturbed. All these discussions about the symmetry in solutions of type as mentioned above integrals will be helpful for those researchers working in physical sciences including quantum physics, mechanics, etc.
In our study, we focused on exploring the existence and uniqueness of symmetric solutions for a system of terminal FDEs in the context of weighted spaces. The specific operator we utilized was the w-Hilfer FD, which is known for its importance in fractional calculus. To examine the existence and uniqueness of symmetric solutions, we employed two fixed-point theorems, the Schauder and Banach fixed-point theorems. These theorems provide the mathematical tools for establishing the existence and uniqueness of solutions in various settings. By applying Schauder and Banach fixed-point theorems to our system, we were able to obtain two key theoretical findings and four symmetry cases. These findings clarify the behavior and properties of symmetric solutions in the considered Banach weighted spaces. By mathematical analysis, we discussed two types of stability in the Ulam sense. The investigation of terminal FDEs using the w-Hilfer FD operator in Banach weighted spaces represents a recent and cutting-edge research direction. We believe this study contributes to the advancement of knowledge in this particular area. Overall, our study not only addresses the challenge of examining the existence and uniqueness of symmetric solutions, but also provides valuable insights into the dynamics and applications of such solutions in fractal mediums. By offering a comprehensive explanation of these processes, our research contributes to the existing literature and expands the understanding of the real-world application of fractional calculus. In future work, we will extend these results with a piecewise fractional operator with the concept of Atangana-Baleanu.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (grant no. 6057). This study is supported via funding from Prince Sattam bin Abdulaziz University, project number (PSAU/2024/R/1445). The researchers wish to extend their sincere gratitude to the Deanship of Scientific Research at the Islamic University of Madinah for the support provided to the Post-Publishing Program.
The author declares that the research was conducted in the absence of any conflict of interest.
[1] |
R. Khalil, M. Al Horani, A. Yousef, M. Sababheh, A new definition of fractional derivative, J. Comput. Appl. Math., 264 (2014), 65–70. http://dx.doi.org/10.1016/j.cam.2014.01.002 doi: 10.1016/j.cam.2014.01.002
![]() |
[2] |
R. Mabel Lizzy, K. Balachandran, J. Trujillo, Controllability of nonlinear stochastic fractional neutral systems with multiple time varying delays in control, Chaos Soliton. Fract., 102 (2017), 162–167. http://dx.doi.org/10.1016/j.chaos.2017.04.024 doi: 10.1016/j.chaos.2017.04.024
![]() |
[3] | K. Diethelm, The analysis of fractional differential equations, Berlin: Springer, 2010. http://dx.doi.org/10.1007/978-3-642-14574-2 |
[4] |
H. Yasmin A. Alderremy, R. Shah, A. Hamid Ganie, S. Aly, Iterative solution of the fractional Wu-Zhang equation under Caputo derivative operator, Front. Phys., 12 (2024), 1333990. http://dx.doi.org/10.3389/fphy.2024.1333990 doi: 10.3389/fphy.2024.1333990
![]() |
[5] |
M. Kbiri Alaoui, K. Nonlaopon, A. Zidan, A. Khan, R. Shah, Analytical investigation of fractional-order Cahn-Hilliard and gardner equations using two novel techniques, Mathematics, 10 (2022), 1643. http://dx.doi.org/10.3390/math10101643 doi: 10.3390/math10101643
![]() |
[6] |
T. Botmart, R. Agarwal, M. Naeem, A. Khan, R. Shah, On the solution of fractional modified Boussinesq and approximate long wave equations with non-singular kernel operators, AIMS Mathematics, 7 (2022), 12483–12513. http://dx.doi.org/10.3934/math.2022693 doi: 10.3934/math.2022693
![]() |
[7] |
H. Yasmin, N. Aljahdaly, A. Saeed, R. Shah, Probing families of optical soliton solutions in fractional perturbed radhakrishnan-kundu-akshmanan model with improved versions of extended direct algebraic method, Fractal Fract., 7 (2023), 512. http://dx.doi.org/10.3390/fractalfract7070512 doi: 10.3390/fractalfract7070512
![]() |
[8] |
H. Yasmin, N. Aljahdaly, A. Saeed, R. Shah, Investigating families of soliton solutions for the complex structured coupled fractional biswas-arshed model in birefringent fibers using a novel analytical technique, Fractal Fract., 7 (2023), 491. http://dx.doi.org/10.3390/fractalfract7070491 doi: 10.3390/fractalfract7070491
![]() |
[9] |
L. Barros, M. Lopes, F. Pedro, E. Esmi, J. Santos, D. Sanchez, The memory effect on fractional calculus: an application in the spread of COVID-19, Comp. Appl. Math., 40 (2021), 72. http://dx.doi.org/10.1007/s40314-021-01456-z doi: 10.1007/s40314-021-01456-z
![]() |
[10] |
R. Pakhira, U. Ghosh, S. Sarkar, Study of memory effects in an inventory model using fractional calculus, Applied Mathematical Sciences, 12 (2018), 797–824. http://dx.doi.org/10.12988/ams.2018.8578 doi: 10.12988/ams.2018.8578
![]() |
[11] |
S. Alkhateeb, S. Hussain, W. Albalawi, S. El-Tantawy, E. El-Awady, Dissipative Kawahara ion-acoustic solitary and cnoidal waves in a degenerate magnetorotating plasma, J. Taibah Univ. Sci., 17 (2023), 2187606. http://dx.doi.org/10.1080/16583655.2023.2187606 doi: 10.1080/16583655.2023.2187606
![]() |
[12] |
R. Alharbey, W. Alrefae, H. Malaikah, E. Tag-Eldin, S. El-Tantawy, Novel approximate analytical solutions to the nonplanar modified Kawahara equation and modeling nonlinear structures in electronegative plasmas, Symmetry, 15 (2023), 97. http://dx.doi.org/10.3390/sym15010097 doi: 10.3390/sym15010097
![]() |
[13] |
S. El-Tantawy, A. Salas, H. Alyouse, M. Alharthi, Novel exact and approximate solutions to the family of the forced damped Kawahara equation and modeling strong nonlinear waves in a plasma, Chinese J. Phys., 77 (2022), 2454–2471. http://dx.doi.org/10.1016/j.cjph.2022.04.009 doi: 10.1016/j.cjph.2022.04.009
![]() |
[14] |
H. Alyousef, A. Salas, M. Alharthi, S. El-tantawy, New periodic and localized traveling wave solutions to a Kawahara-type equation: applications to plasma physics, Complexity, 2022 (2022), 9942267. http://dx.doi.org/10.1155/2022/9942267 doi: 10.1155/2022/9942267
![]() |
[15] |
M. Alharthi, R. Alharbey, S. El-Tantawy, Novel analytical approximations to the nonplanar Kawahara equation and its plasma applications, Eur. Phys. J. Plus, 137 (2022), 1172. http://dx.doi.org/10.1140/epjp/s13360-022-03355-6 doi: 10.1140/epjp/s13360-022-03355-6
![]() |
[16] |
S. El-Tantawy, L. El-Sherif, A. Bakry, W. Alhejaili, A. Wazwaz, On the analytical approximations to the nonplanar damped Kawahara equation: cnoidal and solitary waves and their energy, Phys. Fluids, 34 (2022), 113103. http://dx.doi.org/10.1063/5.0119630 doi: 10.1063/5.0119630
![]() |
[17] |
R. Shah, H. Khan, P. Kumam, M. Arif, An analytical technique to solve the system of nonlinear fractional partial differential equations, Mathematics, 7 (2019), 505. http://dx.doi.org/10.3390/math7060505 doi: 10.3390/math7060505
![]() |
[18] |
H. Khan, R. Shah, D. Baleanu, P. Kumam, M. Arif, Analytical solution of fractional-order hyperbolic telegraph equation, using natural transform decomposition method, Electronics, 8 (2019), 1015. http://dx.doi.org/10.3390/electronics8091015 doi: 10.3390/electronics8091015
![]() |
[19] |
X. Li, Y. Sun, Application of RBF neural network optimal segmentation algorithm in credit rating, Neural Comput. Appl., 33 (2021), 8227–8235. http://dx.doi.org/10.1007/s00521-020-04958-9 doi: 10.1007/s00521-020-04958-9
![]() |
[20] |
T. Ali, Z. Xiao, H. Jiang, B. Li, A class of digital integrators based on trigonometric quadrature rules, IEEE T. Ind. Electron., 71 (2024), 6128–6138. http://dx.doi.org/10.1109/TIE.2023.3290247 doi: 10.1109/TIE.2023.3290247
![]() |
[21] |
B. Chen, J. Hu, B. Ghoso, Finite-time observer based tracking control of heterogeneous multi-AUV systems with partial measurements and intermittent communication, Sci. China Inform. Sci., 67 (2024), 152202. http://dx.doi.org/10.1007/s11432-023-3903-6 doi: 10.1007/s11432-023-3903-6
![]() |
[22] |
B. Chen, J. Hu, Y. Zhao, B. Ghosh, Finite-time observer based tracking control of uncertain heterogeneous underwater vehicles using adaptive sliding mode approach, Neurocomputing, 481 (2022), 322–332. http://dx.doi.org/10.1016/j.neucom.2022.01.038 doi: 10.1016/j.neucom.2022.01.038
![]() |
[23] |
C. Guo, J. Hu, Time base generator based practical predefined-time stabilization of high-order systems with unknown disturbance, IEEE T. Circuits-II, 70 (2023), 2670–2674. http://dx.doi.org/10.1109/TCSII.2023.3242856 doi: 10.1109/TCSII.2023.3242856
![]() |
[24] |
S. Lin, J. Zhang, C. Qiu, Asymptotic analysis for one-stage stochastic linear complementarity problems and applications, Mathematics, 11 (2023), 482. http://dx.doi.org/10.3390/math11020482 doi: 10.3390/math11020482
![]() |
[25] |
L. Liu, S. Zhang, L. Zhang, G. Pan, J. Yu, Multi-UUV maneuvering counter-game for dynamic target scenario based on fractional-order recurrent neural network, IEEE T. Cybernetics, 53 (2023), 4015–4028. http://dx.doi.org/10.1109/TCYB.2022.3225106 doi: 10.1109/TCYB.2022.3225106
![]() |
[26] | Y. Kai, S. Chen, K. Zhang, Z. Yin, Exact solutions and dynamic properties of a nonlinear fourth-order time-fractional partial differential equation, Wave. Random Complex, in press. http://dx.doi.org/10.1080/17455030.2022.2044541 |
[27] |
D. Kaup, On the inverse scattering problem for cubic eigenvalue problems of the class ψxxx+6Qψx+6Rψ=λψ, Stud. Appl. Math., 62 (1980), 189–216. http://dx.doi.org/10.1002/sapm1980623189 doi: 10.1002/sapm1980623189
![]() |
[28] |
B. Kupershmidt, A super Korteweg-de Vries equation: an integrable system, Phys. Lett. A, 102 (1984), 213–215. http://dx.doi.org/10.1016/0375-9601(84)90693-5 doi: 10.1016/0375-9601(84)90693-5
![]() |
[29] | O. Abdulaziz, I. Hashim, M. Chowdhury, A. Zulkifle, Assessment of decomposition method for linear and nonlinear fractional differential equations, Far East Journal of Applied Mathematics, 28 (2007), 95–112. |
[30] |
J. He, X. Wu, Construction of solitary solution and compacton-like solution by variational iteration method, Chaos Soliton. Fract., 29 (2006), 108–113. http://dx.doi.org/10.1016/j.chaos.2005.10.100 doi: 10.1016/j.chaos.2005.10.100
![]() |
[31] |
Z. Odibat, S. Momani, Application of variational iteration method to nonlinear differential equations of fractional order, Int. J. Nonlin. Sci. Num., 7 (2006), 27–34. http://dx.doi.org/10.1515/IJNSNS.2006.7.1.27 doi: 10.1515/IJNSNS.2006.7.1.27
![]() |
[32] |
J. He, Application of homotopy perturbation method to nonlinear wave equations, Chaos Soliton. Fract., 26 (2005), 695–700. http://dx.doi.org/10.1016/j.chaos.2005.03.006 doi: 10.1016/j.chaos.2005.03.006
![]() |
[33] |
Z. Odibat, Exact solitary solutions for variants of the KdV equations with fractional time derivatives, Chaos Soliton. Fract., 40 (2009), 1264–1270. http://dx.doi.org/10.1016/j.chaos.2007.08.080 doi: 10.1016/j.chaos.2007.08.080
![]() |
[34] |
Q. Wang, Homotopy perturbation method for fractional KdV-Burgers equation, Chaos Soliton. Fract., 35 (2008), 843–850. http://dx.doi.org/10.1016/j.chaos.2006.05.074 doi: 10.1016/j.chaos.2006.05.074
![]() |
[35] |
S. Liao, An approximate solution technique not depending on small parameters: a special example, Int. J. NonLin. Mech., 30 (1995), 371–380. http://dx.doi.org/10.1016/0020-7462(94)00054-E doi: 10.1016/0020-7462(94)00054-E
![]() |
[36] | M. El-Tawil, S. Huseen, The q-homotopy analysis method (q-HAM), Int. J. Appl. Math. Mech., 8 (2012), 51–75. |
[37] | J. Biazar, K. Hosseini, P. Gholamin, Homotopy perturbation method for solving KdV and Sawada-Kotera equations, Journal of Operational Research in its Applications, 6 (2009), 11–16. |
[38] | S. Dinarvand, S. Khosravi, A. Doosthoseini, M. Rashidi, The homotopy analysis method for solving the Sawada-Kotera and Lax's fifth-order KdV equations, Adv. Theor. Appl. Mech., 1 (2008), 327–335. |
[39] |
O. Abu Arqub, Series solution of fuzzy differential equations under strongly generalized differentiability, J. Adv. Res. Appl. Math., 5 (2013), 31–52. http://dx.doi.org/10.5373/jaram.1447.051912 doi: 10.5373/jaram.1447.051912
![]() |
[40] |
O. Abu Arqub, Z. Abo-Hammour, R. Al-Badarneh, S. Momani, A reliable analytical method for solving higher-order initial value problems, Discrete Dyn. Nat. Soc., 2013 (2013), 673829. http://dx.doi.org/10.1155/2013/673829 doi: 10.1155/2013/673829
![]() |
[41] |
O. Abu Arqub, A. El-Ajou, Z. Zhour, S. Momani, Multiple solutions of nonlinear boundary value problems of fractional order: a new analytic iterative technique, Entropy, 16 (2014), 471–493. http://dx.doi.org/10.3390/e16010471 doi: 10.3390/e16010471
![]() |
[42] |
A. El-Ajou, O. Abu Arqub, S. Momani, Approximate analytical solution of the nonlinear fractional KdV-Burgers equation: a new iterative algorithm, J. Comput. Phys., 293 (2015), 81–95. http://dx.doi.org/10.1016/j.jcp.2014.08.004 doi: 10.1016/j.jcp.2014.08.004
![]() |
[43] |
F. Xu, Y. Gao, X. Yang, H. Zhang, Construction of fractional power series solutions to fractional Boussinesq equations using residual power series method, Math. Probl. Eng., 2016 (2016), 5492535. http://dx.doi.org/10.1155/2016/5492535 doi: 10.1155/2016/5492535
![]() |
[44] |
J. Zhang, Z. Wei, L. Li, C. Zhou, Least-squares residual power series method for the time-fractional differential equations, Complexity, 2019 (2019), 6159024. http://dx.doi.org/10.1155/2019/6159024 doi: 10.1155/2019/6159024
![]() |
[45] |
I. Jaradat, M. Alquran, R. Abdel-Muhsen, An analytical framework of 2D diffusion, wave-like, telegraph, and Burgers' models with twofold Caputo derivatives ordering, Nonlinear Dyn., 93 (2018), 1911–1922. http://dx.doi.org/10.1007/s11071-018-4297-8 doi: 10.1007/s11071-018-4297-8
![]() |
[46] |
I. Jaradat, M. Alquran, K. Al-Khaled, An analytical study of physical models with inherited temporal and spatial memory, Eur. Phys. J. Plus, 133 (2018), 162. http://dx.doi.org/10.1140/epjp/i2018-12007-1 doi: 10.1140/epjp/i2018-12007-1
![]() |
[47] | M. Alquran, K. Al-Khaled, S. Sivasundaram, H. Jaradat, Mathematical and numerical study of existence of bifurcations of the generalized fractional Burgers-Huxley equation, Nonlinear Stud., 24 (2017), 235–244. |
[48] |
M. Alquran, M. Alsukhour, M. Ali, I. Jaradat, Combination of Laplace transform and residual power series techniques to solve autonomous n-dimensional fractional nonlinear systems, Nonlinear Engineering, 10 (2021), 282–292. http://dx.doi.org/10.1515/nleng-2021-0022 doi: 10.1515/nleng-2021-0022
![]() |
[49] | A. Khan, M. Junaid, I. Khan, F. Ali, K. Shah, D. Khan, Application of homotopy analysis natural transform method to the solution of nonlinear partial differential equations, Sci. Int. (Lahore), 29 (2017), 297–303. |
[50] | M. Zhang, Y. Liu, X. Zhou, Efficient homotopy perturbation method for fractional non-linear equations using Sumudu transform, Therm. Sci., 19 (2015), 1167–1171. |
[51] | R. Al-Deiakeh, M. Ali, M. Alquran, T. Sulaiman, S. Momani, M. Al-Smadi, On finding closed-form solutions to some nonlinear fractional systems via the combination of multi-Laplace transform and the Adomian decomposition method, Rom. Rep. Phys., 74 (2022), 111. |
[52] |
H. Eltayeb, A. Kilicman, A note on double Laplace transform and telegraphic equations, Abstr. Appl. Anal., 2013 (2013), 932578. http://dx.doi.org/10.1155/2013/932578 doi: 10.1155/2013/932578
![]() |
[53] | M. Alquran, K. Al-Khaled, M. Ali, A. Ta'any, The combined Laplace transform-differential transform method for solving linear non-homogeneous PDEs, J. Math. Comput. Sci., 2 (2012), 690–701. |
[54] | K. Aboodh, The new integral transform "Aboodh transform'', Global Journal of Pure and Applied Mathematics, 9 (2013), 35–43. |
[55] | S. Aggarwal, R. Chauhan, A comparative study of Mohand and Aboodh transforms, International Journal of Research in Advent Technology, 7 (2019), 520–529. |
[56] |
M. Benattia, K. Belghaba, Application of the Aboodh transform for solving fractional delay differential equations, Universal Journal of Mathematics and Applications, 3 (2020), 93–101. http://dx.doi.org/10.32323/ujma.702033 doi: 10.32323/ujma.702033
![]() |
[57] |
B. Delgado, J. Macias-Diaz, On the general solutions of some non-homogeneous Div-curl systems with Riemann-Liouville and Caputo fractional derivatives, Fractal Fract., 5 (2021), 117. http://dx.doi.org/10.3390/fractalfract5030117 doi: 10.3390/fractalfract5030117
![]() |
[58] |
S. Alshammari, M. Al-Smadi, I. Hashim, M. Alias, Residual power series technique for simulating fractional Bagley-Torvik problems emerging in applied physics, Appl. Sci., 9 (2019), 5029. http://dx.doi.org/10.3390/app9235029 doi: 10.3390/app9235029
![]() |
[59] |
S. Almutlak, S. Parveen, S. Mahmood, A. Qamar, B. Alotaibi, S. El-Tantawy, On the propagation of cnoidal wave and overtaking collision of slow shear Alfvén solitons in low β−magnetized plasmas, Phys. Fluids, 35 (2023), 075130. http://dx.doi.org/10.1063/5.0158292 doi: 10.1063/5.0158292
![]() |
[60] |
W. Albalawi, S. El-Tantawy, A. Salas, On the rogue wave solution in the framework of a Korteweg-de Vries equation, Results Phys., 30 (2021), 104847. http://dx.doi.org/10.1016/j.rinp.2021.104847 doi: 10.1016/j.rinp.2021.104847
![]() |
[61] |
T. Hashmi, R. Jahangir, W. Masood, B. Alotaibi, S. Ismaeel, S. El-Tantawy, Head-on collision of ion-acoustic (modified) Korteweg-de Vries solitons in Saturn's magnetosphere plasmas with two temperature superthermal electrons, Phys. Fluids, 35 (2023), 103104. http://dx.doi.org/10.1063/5.0171220 doi: 10.1063/5.0171220
![]() |
[62] |
A. Wazwaz, W. Alhejaili, S. El-Tantawy, Physical multiple shock solutions to the integrability of linear structures of Burgers hierarchy, Phys. Fluids, 35 (2023), 123101. http://dx.doi.org/10.1063/5.0177366 doi: 10.1063/5.0177366
![]() |
[63] |
S. El-Tantawy, R. Matoog, R. Shah, A. Alrowaily, S. Ismaeel, On the shock wave approximation to fractional generalized Burger-Fisher equations using the residual power series transform method, Phys. Fluids, 36 (2024), 023105. http://dx.doi.org/10.1063/5.0187127 doi: 10.1063/5.0187127
![]() |
[64] |
S. El-Tantawy, A. Salas, H. Alyousef, M. Alharthi, Novel approximations to a nonplanar nonlinear Schrodinger equation and modeling nonplanar rogue waves/breathers in a complex plasma, Chaos Soliton. Fract., 163 (2022), 112612. http://dx.doi.org/10.1016/j.chaos.2022.112612 doi: 10.1016/j.chaos.2022.112612
![]() |
[65] |
S. El-Tantawy, A. Wazwaz, R. Schlickeiser, Solitons collision and freak waves in a plasma with Cairns-Tsallis particle distributions, Plasma Phys. Control. Fusion, 57 (2015), 125012. http://dx.doi.org/10.1088/0741-3335/57/12/125012 doi: 10.1088/0741-3335/57/12/125012
![]() |
1. | Muhammad Akram, Ayesha Shareef, Ahmad N. Al-Kenani, Pythagorean fuzzy incidence graphs with application in one-way toll road network, 2024, 9, 2364-4966, 10.1007/s41066-024-00455-9 | |
2. | Muhammad Akram, Uzma Ahmad, Mohammed M. Ali Al-Shamiri, Ayesha Shareef, Algorithms for computing Pythagorean fuzzy average edge connectivity of Pythagorean fuzzy graphs, 2024, 1598-5865, 10.1007/s12190-023-01970-9 |