In this paper, a fractional order HIV/HTLV co-infection model with HIV-specific antibody immune response is established. Two cases are considered: constant control and optimal control. For the constant control system, the existence and uniqueness of the positive solutions are proved, and then the sufficient conditions for the existence and stability of five equilibriums are obtained. For the second case, the Pontryagin's Maximum Principle is used to analyze the optimal control, and the formula of the optimal solution are derived. After that, some numerical simulations are performed to validate the theoretical prediction. Numerical simulations indicate that in the case of HIV/HTLV co-infection, the concentration of CD4+T cells is no longer suitable as an effective reference data for understanding the development process of the disease. On the contrary, the number of HIV virus particles should be used as an important indicator for reference.
Citation: Ruiqing Shi, Yihong Zhang. Dynamic analysis and optimal control of a fractional order HIV/HTLV co-infection model with HIV-specific antibody immune response[J]. AIMS Mathematics, 2024, 9(4): 9455-9493. doi: 10.3934/math.2024462
[1] | Weaam Alhejaili, Mohammed. K. Elboree, Abdelraheem M. Aly . A symbolic computation approach and its application to the Kadomtsev-Petviashvili equation in two (3+1)-dimensional extensions. AIMS Mathematics, 2022, 7(11): 20085-20104. doi: 10.3934/math.20221099 |
[2] | Wafaa B. Rabie, Hamdy M. Ahmed, Taher A. Nofal, Soliman Alkhatib . Wave solutions for the (3+1)-dimensional fractional Boussinesq-KP-type equation using the modified extended direct algebraic method. AIMS Mathematics, 2024, 9(11): 31882-31897. doi: 10.3934/math.20241532 |
[3] | Junjie Li, Gurpreet Singh, Onur Alp İlhan, Jalil Manafian, Yusif S. Gasimov . Modulational instability, multiple Exp-function method, SIVP, solitary and cross-kink solutions for the generalized KP equation. AIMS Mathematics, 2021, 6(7): 7555-7584. doi: 10.3934/math.2021441 |
[4] | Abeer S. Khalifa, Hamdy M. Ahmed, Niveen M. Badra, Jalil Manafian, Khaled H. Mahmoud, Kottakkaran Sooppy Nisar, Wafaa B. Rabie . Derivation of some solitary wave solutions for the (3+1)- dimensional pKP-BKP equation via the IME tanh function method. AIMS Mathematics, 2024, 9(10): 27704-27720. doi: 10.3934/math.20241345 |
[5] | Boyu Wang . A splitting lattice Boltzmann scheme for (2+1)-dimensional soliton solutions of the Kadomtsev-Petviashvili equation. AIMS Mathematics, 2023, 8(11): 28071-28089. doi: 10.3934/math.20231436 |
[6] | Jamilu Sabi'u, Sekson Sirisubtawee, Surattana Sungnul, Mustafa Inc . Wave dynamics for the new generalized (3+1)-D Painlevé-type nonlinear evolution equation using efficient techniques. AIMS Mathematics, 2024, 9(11): 32366-32398. doi: 10.3934/math.20241552 |
[7] | Zheng Dou, Kexin Luo . Global weak solutions of nonlinear rotation-Camassa-Holm model. AIMS Mathematics, 2023, 8(7): 15285-15298. doi: 10.3934/math.2023781 |
[8] | Yaya Wang, Md Nurul Raihen, Esin Ilhan, Haci Mehmet Baskonus . On the new sine-Gordon solitons of the generalized Korteweg-de Vries and modified Korteweg-de Vries models via beta operator. AIMS Mathematics, 2025, 10(3): 5456-5479. doi: 10.3934/math.2025252 |
[9] | Ninghe Yang . Exact wave patterns and chaotic dynamical behaviors of the extended (3+1)-dimensional NLSE. AIMS Mathematics, 2024, 9(11): 31274-31294. doi: 10.3934/math.20241508 |
[10] | Yunxi Guo, Ying Wang . The Cauchy problem to a gkCH equation with peakon solutions. AIMS Mathematics, 2022, 7(7): 12781-12801. doi: 10.3934/math.2022707 |
In this paper, a fractional order HIV/HTLV co-infection model with HIV-specific antibody immune response is established. Two cases are considered: constant control and optimal control. For the constant control system, the existence and uniqueness of the positive solutions are proved, and then the sufficient conditions for the existence and stability of five equilibriums are obtained. For the second case, the Pontryagin's Maximum Principle is used to analyze the optimal control, and the formula of the optimal solution are derived. After that, some numerical simulations are performed to validate the theoretical prediction. Numerical simulations indicate that in the case of HIV/HTLV co-infection, the concentration of CD4+T cells is no longer suitable as an effective reference data for understanding the development process of the disease. On the contrary, the number of HIV virus particles should be used as an important indicator for reference.
The prominent Banach contraction in metric spaces has laid a solid foundation for fixed point theory in metric space. The basic idea of the contraction mapping principle has been fine-tuned in several domains (see, e.g. [12,13,15]). The applications of fixed point range across inequalities, approximation theory, optimization and so on. Researchers in this area have introduced several new concepts in metric space and obtained a great deal of fixed point results for linear and nonlinear contractions. Recently, Karapınar et al. [7] introduced a new notion of hybrid contraction which is a unification of some existing linear and nonlinear contractions in metric space.
On the other hand, Mustafa [8] pioneered an extension of metric space by the name, generalized metric space (or more precisely, G-metric space) and proved some fixed point results for Banach-type contraction mappings. This new generalization was brought to spotlight by Mustafa and Sims [9]. Subsequently, Mustafa et al. [11] obtained some engrossing fixed point results for Lipschitzian-type mappings on G-metric space. However, Jleli and Samet [5], as well as Samet et al. [16] noted that most of the fixed point results in G-metric space are direct consequences of existence results in corresponding metric space. Jleli and Samet [5] further observed that if a G-metric is consolidated into a quasi-metric, then the resultant fixed point results become the known fixed point results in the setting of quasi-metric space. Motivated by the latter observation, many investigators (see for instance, [3,6]) have established techniques of obtaining fixed point results in symmetric G-metric space that are not deducible from their ditto ones in metric space or quasi-metric space.
Following the existing literature, we realize that hybrid fixed point results in G-metric space are not adequately investigated. Hence, motivated by the ideas in [3,6,7], we introduce a new concept of hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction in G-metric space and prove some related fixed point theorems. An example is constructed to demonstrate that our result is valid, an improvement of existing result and the main ideas obtained herein do not reduce to any existence result in metric space. A corollary is presented to show that the concept proposed herein is a generalization and improvement of well-known fixed point result in metric space. Finally, one of our obtained corollaries is applied to establish novel existence conditions for solution of a class of integral equations.
In this section, we will present some fundamental notations and results that will be deployed subsequently.
All through, every set Φ is considered non-empty, N is the set of natural numbers, R represents the set of real numbers and R+, the set of non-negative real numbers.
Definition 2.1. [9] Let Φ be a non-empty set and let G:Φ×Φ×Φ⟶R+ be a function satisfying:
(G1) G(r,s,t)=0 if r=s=t;
(G2) 0<G(x,r,s) for all r,s∈Φ with r≠s;
(G3) G(r,r,s)≤G(r,s,t), for all r,s,t∈Φ with t≠s;
(G4) G(r,s,t)=G(r,t,s)=G(s,r,t)=... (symmetry in all variables);
(G5) G(r,s,t)≤G(r,u,u)+G(u,s,t), for all r,s,t,u∈Φ (rectangular inequality).
Then the function G is called a generalised metric, or more precisely, a G-metric on Φ, and the pair (Φ,G) is called a G-metric space.
Example 2.2. [11] Let (Φ,d) be a usual metric space, then (Φ,Gp) and (Φ,Gm) are G-metric space, where
Gp(r,s,t)=d(r,s)+d(s,t)+d(r,t)∀r,s,t∈Φ. | (2.1) |
Gm(r,s,t)=max{d(r,s),d(s,t),d(r,t)}∀r,s,t∈Φ. | (2.2) |
Definition 2.3. [11] Let (Φ,G) be a G-metric space and let {rn} be a sequence of points of Φ. Then {rn} is said to be G-convergent to r if limn,m→∞G(r,rn,rm)=0; that is, for any ϵ>0, there exists n0∈N such that G(r,rn,rm)<ϵ, ∀n,m≥n0. We refer to r as the limit of the sequence {rn}.
Proposition 1. [11] Let (Φ,G) be a G-metric space. Then the following are equivalent:
(i) {rn} is G-convergent to r.
(ii) G(r,rn,rm)⟶0, as n,m→∞.
(iii) G(rn,r,r)⟶0, as n→∞.
(iv) G(rn,rn,r)⟶0, as n→∞.
Definition 2.4. [11] Let (Φ,G) be a G-metric space. A sequence {rn} is called G-Cauchy if for any ϵ>0, we can find n0∈N such that G(rn,rm,rl)<ϵ, ∀n,m,l≥n0, that is, G(rn,rm,rl)⟶0, as n,m,l→∞.
Proposition 2. [11] If (Φ,G) is a G-metric space, the following statements are equivalent:
(i) The sequence {rn} is G-Cauchy.
(ii) For every ϵ>0, there exists n0∈N such that G(rn,rm,rm)<ϵ, ∀n,m≥n0.
Definition 2.5. [11] Let (Φ,G) and (Φ′,G′) be two G-metric spaces and f:(Φ,G)⟶(Φ′,G′) be a function. Then f is G-continuous at a point u∈Φ if and only if for any ϵ>0, there exists δ>0 such that r,s∈Φ; and G(u,r,s)<δ implies G′(f(u),f(r),f(s))<ϵ. A function f is G-continuous on Φ if and only if it is G-continuous at all u∈Φ.
Proposition 3. [11] Let (Φ,G) and (Φ′,G′) be two G-metric spaces. Then a function f:(Φ,G)⟶(Φ′,G′) is said to be G-continuous at a point r∈Φ if and only if it is G-sequentially continuous at r, that is, whenever {rn} is G-convergent to r, {frn} is G-convergent to fr.
Definition 2.6. [11] A G-metric space (Φ,G) is called symmetric G-metric space if
G(r,r,s)=G(s,r,r)∀r,s∈Φ. |
Proposition 4. [11] Let (Φ,G) be a G-metric space. Then the function G(r,s,t) is jointly continuous in all variables.
Proposition 5. [11] Every G-metric space (Φ,G) defines a metric space (Φ,dG) by
dG(r,s)=G(r,s,s)+G(s,r,r)∀r,s∈Φ. | (2.3) |
Note that for a symmetric G-metric space (Φ,G),
(Φ,dG)=2G(r,s,s)∀r,s∈Φ. | (2.4) |
On the other hand, if (Φ,G) is not symmetric, then by the G-metric properties,
32G(r,s,s)≤dG(r,s)≤3G(r,s,s)∀r,s∈Φ, | (2.5) |
and that in general, these inequalities are sharp.
Definition 2.7. [11] A G-metric space (Φ,G) is referred to as G-complete (or complete G-metric) if every G-Cauchy sequence in (Φ,G) is G-convergent in (Φ,G).
Proposition 6. [11] A G-metric space (Φ,G) is G-complete if and only if (Φ,dG) is a complete metric space.
Popescu [14] gave the following definition in the setting of metric space.
Definition 2.8. [14] Let α:Φ×Φ⟶R+ be a function. A self-mapping Γ:Φ⟶Φ is referred to as α-orbital admissible if for all r∈Φ,
α(r,Γr)≥1 ⇒ α(Γr,Γ2r)≥1.
We modify the above definitions in the framework of G-metric space as follows:
Definition 2.9. Let α:Φ×Φ×Φ⟶R+ be a function. A self-mapping Γ:Φ⟶Φ is called (G-α)-orbital admissible if for all r∈Φ,
α(r,Γr,Γ2r)≥1 implies α(Γr,Γ2r,Γ3r)≥1.
Definition 2.10. [2] Let α:Φ×Φ×Φ⟶R+ be a mapping. The set Φ is called regular with respect to α if and only if for a sequence {rn} in Φ such that α(rn,rn+1,rn+2)≥1, for all n and rn→r∈Φ as n→∞, we have α(rn,r,r)≥1 for all n.
The following Propositions 7 and 8 were studied in [7], where the constant C∈R+. However, we noticed that the arguments in their proofs are only valid if we restrict R+ to [0,1). Hence, we reexamine them in the latter interval.
Proposition 7. Given C∈[0,1), let {ρn}⊂R+ be a sequence such that
ρn+2≤Cmax{ρn,ρn+1}∀n∈N. | (2.6) |
Let K=max{ρ0,ρ1}. Then
ρ2n≤CnK,ρ2n+1≤CnK∀n≥1. | (2.7) |
Proof. The proof is by induction.
For n=0, we have from (2.6) that
ρ2≤Cmax{ρ0,ρ1}=CK. |
For n=1, (2.6) becomes
ρ3≤Cmax{ρ1,ρ2}≤Cmax{ρ1,Cmax{ρ0,ρ1}}≤Cmax{ρ1,CK}≤CK. |
Suppose that (2.7) holds for some n∈N. Then
ρ2n+2≤Cmax{ρ2n,ρ2n+1}≤Cmax{CnK,CnK}=Cn+1K;ρ2n+3≤Cmax{ρ2n+1,ρ2n+2}≤Cmax{CnK,Cn+1K}=Cn+1K. |
This completes the induction. Hence, the proof.
Lemma 2.11. Let {rn} be a sequence in a G-metric space (Φ,G). Suppose that there exists C∈[0,1) such that
G(rn+2,rn+3,rn+4)≤Cmax{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}∀n∈N. |
Then {rn} is a G-Cauchy sequence in (Φ,G).
Proof. Consider the sequence {ρn} in Φ defined by
ρn=G(rn,rn+1,rn+2)∀n∈N. |
Then by the hypothesis, {ρn} satisfies (2.6). Hence, by Proposition 7, we obtain
G(r2n,r2n+1,r2n+2)=ρ2n≤CnK;G(r2n+1,r2n+2,r2n+3)=ρ2n+1≤CnK, |
for all n∈N, where K=max{ρ0,ρ1}. In particular,
G(r2n,r2n+1,r2n+2)+G(r2n+1,r2n+2,r2n+3)≤2CnK∀n∈N. | (2.8) |
If C=0 or K=0, then {rn}n≥2 is constant, hence G-Cauchy. Assume that C>0 and K>0. To see that {rn}n∈N is G-Cauchy in (Φ,G), take arbitrary ϵ>0. Since ϵ2K>0 and 0<C<1, then we can find n0∈N such that
∞∑i=n0Ci<ϵ2K. |
In particular,
2Kp∑i=n0Ci<2K∞∑i=n0Ci<ϵ∀p∈N;p≥n0. |
Let l,m,n∈N, where 2n0≤n<m<l. Let p∈N be such that p≥n0+2 and 2p≥l. Then by rectangular inequality and (2.8), we have
G(rn,rm,rl)≤l−2∑j=nG(rj,rj+1,rj+2)≤2p−2∑j=2n0G(rj,rj+1,rj+2)=p−2∑q=n0[G(r2q,r2q+1,r2q+2)+G(r2q+1,r2q+2,r2q+3)]≤p−2∑q=n02CqK=2Kp∑q=n0Cq<2K∞∑q=n0Cq<ϵ. |
This completes the proof.
Proposition 8. Given C∈[0,1) and σ∈(0,1), let {ρn}⊂(0,1) be a sequence such that
ρn+2≤Cmax{ρn,ρn+1}σ∀n∈N. | (2.9) |
Let K=max{ρ0,ρ1}. Then
ρ2n≤C1+σ+σ2+...+σn−1Kσn,ρ2n+1≤C1+σ+σ2+...+σn−1Kσn∀n≥1. | (2.10) |
Therefore,
lim supn→∞ρn≤C11−σ. |
Proof. The result holds if C=0. Assume C>0. Since σ∈(0,1), then obviously,
σm+1<σm<σm−1<...<σ2<σ∀m∈N;m≥1. |
Therefore, since K∈(0,1), we have
Kσm+1<Kσm<Kσm−1<...<Kσ2<Kσ∀m∈N;m≥1. |
Consider (2.9) and let n=0. Then
ρ2≤Cmax{ρ0,ρ1}σ≤CKσ. |
For n=1, we obtain
ρ3≤Cmax{ρ1,ρ2}σ≤Cmax{ρ1,CKσ}σ≤Cmax{K,CKσ}σ=Cmax{Kσ,CσKσ2}≤Cmax{Kσ,CσKσ}=CKσmax{1,Cσ}=CKσ. |
Now, assume that (2.10) holds for some n∈N. We then prove that it holds for n+1. Therefore,
ρ2n+2≤Cmax{ρ2n,ρ2n+1}σ≤Cmax{C1+σ+σ2+...+σn−1Kσn,C1+σ+σ2+...+σn−1Kσn}σ=C(C1+σ+σ2+...+σn−1Kσn)σ=C(Cσ+σ2+...+σnKσn+1)=C1+σ+σ2+...+σnKσn+1. |
Similarly,
ρ2n+3≤Cmax{ρ2n+1,ρ2n+2}σ≤Cmax{C1+σ+σ2+...+σn−1Kσn,C1+σ+σ2+...+σnKσn+1}σ=Cmax{Cσ+σ2+...+σnKσn+1,Cσ+σ2+...+σn+1Kσn+2}≤Cmax{Cσ+σ2+...+σnKσn+1,Cσ+σ2+...+σn+1Kσn+1}=CKσn+1max{Cσ+σ2+...+σn,Cσ+σ2+...+σn+1}=CKσn+1(Cσ+σ2+...+σn)=C1+σ+σ2+...+σnKσn+1. |
This completes the induction, therefore, verifying (2.10). Noting that {σn}n∈N→0 as n→∞, it is obvious that {Kσn}n∈N→K0=1. Also,
limn→∞(1+σ+σ2+...+σn)≤∞∑i=0σi=11−σ. |
Therefore,
limn→∞C1+σ+σ2+...+σn=C(11−σ). |
Corollary 1. Let {rn} be a sequence in a G-metric space (Φ,G). Assume that there exist C∈[0,1) and θ,η∈[0,1] with θ+η=1 such that
G(rn+2,rn+3,rn+4)≤C[G(rn,rn+1,rn+2)θ⋅G(rn+1,rn+2,rn+3)η]∀n∈N. |
Then {rn} is a G-Cauchy sequence in (Φ,G).
Proof. Notice that
G(rn+2,rn+3,rn+4)≤C[G(rn,rn+1,rn+2)θ⋅G(rn+1,rn+2,rn+3)η]≤C[max{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}θ⋅max{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}η]=Cmax{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}θ+η=Cmax{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}. |
Then the result follows from Lemma 2.11.
Karapınar et al. [7] gave the following definition of hybrid-interpolative Reich-Istr˘aţescu-type contraction in metric space.
Definition 2.12. [7] Let (Φ,d) be a metric space and let α:Φ×Φ⟶R+ be a function. A self-mapping Γ:Φ⟶Φ is called hybrid-interpolative Reich-Istr˘aţescu-type contraction if for some q∈R+, there exist constants μ∈(0,1), δ≥0 and λi≥0 with i=1,2,...,5 such that for all distinct r,s∈Φ∖Fix(Γ),
α(r,s)d(Γ2r,Γ2s)≤μM(r,s), | (2.11) |
where
M(r,s)={[λ1d(r,s)q+λ2d(r,Γr)q+λ3d(s,Γs)q+λ4d(Γr,Γs)q+λ5d(Γr,Γ2r)q+δd(Γs,Γ2s)q]1q,forq>0,with∑5i=1λi+δ≤1;[d(r,s)]λ1⋅[d(r,Γr)]λ2⋅[d(s,Γs)]λ3⋅[d(Γr,Γs)]λ4⋅[d(Γr,Γ2r)]λ5⋅[d(Γs,Γ2s)]δ,forq=0,with∑5i=1λi+δ=1, | (2.12) |
and Fix(Γ)={r∈Φ:Γr=r}.
In this section, we introduce a new concept of hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction in G-metric space.
Definition 3.1. Let (Φ,G) be a G-metric space and let α:Φ×Φ⟶R+ be a function. A self-mapping Γ:Φ⟶Φ is called hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction if for some q∈R+, there exist constants μ∈(0,1), δ≥0 and λi≥0 with i=1,2,...,5 such that for all r,s∈Φ∖Fix(Γ),
α(r,s,Γs)G(Γ2r,Γ2s,Γ3s)≤μM(r,s,Γs), | (3.1) |
where
M(r,s,Γs)={[λ1G(r,s,Γs)q+λ2G(r,Γr,Γ2r)q+λ3G(s,Γs,Γ2s)q+λ4G(Γr,Γs,Γ2s)q+λ5G(Γr,Γ2r,Γ3r)q+δG(Γs,Γ2s,Γ3s)q]1q,forq>0,with∑5i=1λi+δ≤1;[G(r,s,Γs)]λ1⋅[G(r,Γr,Γ2r)]λ2⋅[G(s,Γs,Γ2s)]λ3⋅[G(Γr,Γs,Γ2s)]λ4⋅[G(Γr,Γ2r,Γ3r)]λ5⋅[G(Γs,Γ2s,Γ3s)]δ,forq=0,with∑5i=1λi+δ=1, | (3.2) |
and Fix(Γ)={r∈Φ:Γr=r}.
Our main result is the following.
Theorem 3.2. Let (Φ,G) be a complete G-metric space and let Γ:Φ⟶Φ be a hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction satisfying the following conditions:
(i) Γ is (G-α)-orbital admissible;
(ii) Γ is continuous;
(iii) there exists r0∈Φ such that α(r0,Γr0,Γ2r0)≥1.
Then Γ has at least a fixed point in Φ.
Proof. Let r0∈Φ be such that α(r0,Γr0,Γ2r0)≥1. Since Γ is (G-α)-orbital admissible, then α(Γr0,Γ2r0,Γ3r0)≥1 and by induction, we have α(Γnr0,Γn+1r0,Γn+2r0)≥1 for any n∈N. Let {rn} be a sequence in Φ defined by rn=Γnr0 for all n∈N. If there exists some m∈N such that Γrm=rm+1=rm, then clearly, rm is a fixed point of Γ. Assume now that rn≠rn+1 for any n∈N. Since Γ is hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction, then we have from (3.1) that
G(rn+2,rn+3,rn+4)≤α(rn,rn+1,Γrn+1)G(Γ2rn,Γ2rn+1,Γ3rn+1)≤μM(rn,rn+1,Γrn+1). | (3.3) |
We now consider the following cases:
Case 1: For q>0, we have
M(rn,rn+1,Γrn+1)=[λ1G(rn,rn+1,Γrn+1)q+λ2G(rn,Γrn,Γ2rn)q+λ3G(rn+1,Γrn+1,Γ2rn+1)q+λ4G(Γrn,Γrn+1,Γ2rn+1)q+λ5G(Γrn,Γ2rn,Γ3rn)q+δG(Γrn+1,Γ2rn+1,Γ3rn+1)q]1q=[λ1G(rn,rn+1,rn+2)q+λ2G(rn,rn+1,rn+2)q+λ3G(rn+1,rn+2,rn+3)q+λ4G(rn+1,rn+2,rn+3)q+λ5G(rn+1,rn+2,rn+3)q+δG(rn+2,rn+3,rn+4)q]1q=[(λ1+λ2)G(rn,rn+1,rn+2)q+(λ3+λ4+λ5)G(rn+1,rn+2,rn+3)q+δG(rn+2,rn+3,rn+4)q]1q≤[(λ1+λ2+λ3+λ4+λ5)max{G(rn,rn+1,rn+2)q,G(rn+1,rn+2,rn+3)q}+δG(rn+2,rn+3,rn+4)q]1q. | (3.4) |
Therefore, (3.3) becomes
G(rn+2,rn+3,rn+4)q≤μq[(λ1+λ2+λ3+λ4+λ5)max{G(rn,rn+1,rn+2)q,G(rn+1,rn+2,rn+3)q}+δG(rn+2,rn+3,rn+4)q] |
so that
(1−μqδ)G(rn+2,rn+3,rn+4)q≤μq(λ1+λ2+λ3+λ4+λ5)max{G(rn,rn+1,rn+2)q,G(rn+1,rn+2,rn+3)q}≤μq(1−δ)max{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}q, |
implying that for all n∈N,
G(rn+2,rn+3,rn+4)q≤(μq(1−δ)1−μqδ)max{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}q, |
which is equivalent to
G(rn+2,rn+3,rn+4)≤Cmax{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}, |
where
C=(μq(1−δ)1−μqδ)∈(0,1). |
Hence, by Lemma 2.11, {rn}n∈N is a G-Cauchy sequence in (Φ,G).
Case 2: For q=0, we have
M(rn,rn+1,Γrn+1)=G(rn,rn+1,Γrn+1)λ1⋅G(rn,Γrn,Γ2rn)λ2⋅G(rn+1,Γrn+1,Γ2rn+1)λ3⋅G(Γrn,Γrn+1,Γ2rn+1)λ4⋅G(Γrn,Γ2rn,Γ3rn)λ5⋅G(Γrn+1,Γ2rn+1,Γ3rn+1)δ=G(rn,rn+1,rn+2)λ1⋅G(rn,rn+1,rn+2)λ2⋅G(rn+1,rn+2,rn+3)λ3⋅G(rn+1,rn+2,rn+3)λ4⋅G(rn+1,rn+2,rn+3)λ5⋅G(rn+2,rn+3,rn+4)δ=G(rn,rn+1,rn+2)(λ1+λ2)⋅G(rn+1,rn+2,rn+3)(λ3+λ4+λ5)⋅G(rn+2,rn+3,rn+4)δ. |
Therefore, (3.3) becomes
G(rn+2,rn+3,rn+4)≤μG(rn,rn+1,rn+2)(λ1+λ2)⋅G(rn+1,rn+2,rn+3)(λ3+λ4+λ5)⋅G(rn+2,rn+3,rn+4)δ. | (3.5) |
If δ=1, then λ1=λ2=λ3=λ4=λ5=0, and so we obtain
0<G(rn+2,rn+3,rn+4)≤μG(rn+2,rn+3,rn+4), |
which is a contradiction. Therefore, δ<1, so that ∑5i=1λi=1−δ>0, implying that
θ=λ1+λ21−δ,η=λ3+λ4+λ51−δ |
satisfying θ+η=1. Hence, (3.5) becomes
G(rn+2,rn+3,rn+4)(1−δ)≤μG(rn,rn+1,rn+2)(λ1+λ2)⋅G(rn+1,rn+2,rn+3)(λ3+λ4+λ5) |
so that
G(rn+2,rn+3,rn+4)≤μ(11−δ)G(rn,rn+1,rn+2)(λ1+λ21−δ)⋅G(rn+1,rn+2,rn+3)(λ3+λ4+λ51−δ)=μ(11−δ)G(rn,rn+1,rn+2)θ⋅G(rn+1,rn+2,rn+3)η |
for all n in N. Noting that μ∈(0,1), we have
0<1−δ≤1⇒1≤11−δ⇒μ(11−δ)≤μ<1. |
Hence, since 0<μ(11−δ)<1 and θ+η=1, then by Corollary 1, we can conclude that {rn}n∈N is G-Cauchy in (Φ,G).
Therefore, for all q≥0, we have established that {rn}n∈N is a G-Cauchy sequence in (Φ,G) and so by the completeness of (Φ,G), there exists a point c in Φ such that limn→∞rn=c. Moreover, since Γ is continuous, then we can conclude that Γc=c, that is, c is a fixed point of Γ.
Theorem 3.3. Let (Φ,G) be a complete G-metric space and let Γ:Φ⟶Φ be a hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction satisfying the following conditions:
(i) Γ is (G-α)-orbital admissible;
(ii) there exists r0∈Φ such that α(r0,Γr0,Γ2r0)≥1;
(iii) Γ3 is continuous and α(r,Γr,Γ2r)≥1 for any r∈Fix(Γ3).
Then Γ has at least a fixed point in Φ.
Proof. Let r0∈Φ be arbitrary and define a sequence {rn}n∈N in Φ by rn=Γnr0. We have shown in Theorem 3.2 that there exists c∈Φ such that rn→c. Since Γ3 is continuous, then
Γ3c=limn→∞Γ3rn=c, | (3.6) |
that is, c is a fixed point of Γ3. This implies that Γ3 has at least one fixed point in Φ, that is, Fix(Γ3) is nonempty. Moreover,
Γ4c=Γc. | (3.7) |
To see that c is a fixed point of Γ, assume contrary that Γc≠c. In this case, Γc is not a fixed point of Γ3 either, since Γc=Γ3c=c, which is a contradiction. Also, by (3.1), we have
G(c,Γc,Γ2c)≤α(c,Γc,Γ2c)G(c,Γc,Γ2c)=α(c,Γc,Γ2c)G(Γ3c,Γ4c,Γ2c)=α(c,Γc,Γ2c)G(Γ3c,Γc,Γ2c)=α(c,Γc,Γ2c)G(Γc,Γ2c,Γ3c)≤μM(c,Γc,Γ2c). | (3.8) |
Considering Case 1, we obtain
M(c,Γc,Γ2c)=[λ1G(c,Γc,Γ2c)q+λ2G(c,Γc,Γ2c)q+λ3G(Γc,Γ2c,Γ3c)q+λ4G(Γc,Γ2c,Γ3c)q+λ5G(Γc,Γ2c,Γ3c)q+δG(Γ2c,Γ3c,Γ4c)q]1q=[λ1G(c,Γc,Γ2c)q+λ2G(c,Γc,Γ2c)q+λ3G(c,Γc,Γ2c)q+λ4G(c,Γc,Γ2c)q+λ5G(c,Γc,Γ2c)q+δG(c,Γc,Γ2c)q]1q=[(λ1+λ2+λ3+λ4+λ5+δ)G(c,Γc,Γ2c)q]1q=(λ1+λ2+λ3+λ4+λ5+δ)1qG(c,Γc,Γ2c)≤G(c,Γc,Γ2c) |
implying that
0<G(c,Γc,Γ2c)≤μG(c,Γc,Γ2c), |
which is a contradiction. Hence, Γc=c.
Similarly, for Case 2, we obtain
M(c,Γc,Γ2c)=G(c,Γc,Γ2c)λ1⋅G(c,Γc,Γ2c)λ2⋅G(Γc,Γ2c,Γ3c)λ3⋅G(Γc,Γ2c,Γ3c)λ4⋅G(Γc,Γ2c,Γ3c)λ5⋅G(Γ2c,Γ3c,Γ4c)δ=G(c,Γc,Γ2c)λ1⋅G(c,Γc,Γ2c)λ2⋅G(c,Γc,Γ2c)λ3⋅G(c,Γc,Γ2c)λ4⋅G(c,Γc,Γ2c)λ5⋅G(c,Γc,Γ2c)δ=G(c,Γc,Γ2c)(λ1+λ2+λ3+λ4+λ5+δ)=G(c,Γc,Γ2c), |
so that (3.8) becomes
0<G(c,Γc,Γ2c)≤μG(c,Γc,Γ2c), |
which is also a contradiction. Hence, Γc=c.
Therefore, c is a fixed point of Γ in Φ.
Theorem 3.4. If in addition to the hypotheses of Theorem 3.3, we assume supplementary that α(r,s,Γs)≥1 for any r,s∈Fix(Γ), then the fixed point of Γ is unique.
Proof. Let v,c∈Φ be any two fixed point of Γ with v≠c. By replacing this in (3.1) and noting the additional hypothesis, we have:
G(c,v,Γv)≤α(c,v,Γv)G(Γ2c,Γ2v,Γ3v)≤μM(c,v,Γv). | (3.9) |
By Case 1, we obtain
M(c,v,Γv)=[λ1G(c,v,Γv)q+λ2G(c,Γc,Γ2c)q+λ3G(v,Γv,Γ2v)q+λ4G(Γc,Γv,Γ2v)q+λ5G(Γc,Γ2c,Γ3c)q+δG(Γv,Γ2v,Γ3v)q]1q=[λ1G(c,v,Γv)q+λ2G(c,c,c)q+λ3G(v,v,v)q+λ4G(c,v,Γv)q+λ5G(c,c,c)q+δG(v,v,v)q]1q=[(λ1+λ4)G(c,v,Γv)q]1q=(λ1+λ4)1qG(c,v,Γv)≤G(c,v,Γv). |
Hence, (3.9) becomes
0<G(c,v,Γv)≤μG(c,v,Γv), |
which is a contradiction. Therefore, v=c, and so the fixed point of Γ is unique. In the following result, we examine the existence of fixed point of Γ when the G-metric space (Φ,G) is regular.
Theorem 3.5. Let (Φ,G) be a complete G-metric space and let Γ:Φ⟶Φ be a hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction for q=0 such that λ2>0 and λ5>0. Suppose further that:
(i) Γ is (G-α)-orbital admissible;
(ii) there exists r0∈Φ such that α(r0,Γr0,Γ2r0)≥1;
(iii) (Φ,G) is regular with respect to α.
Then Γ has a fixed point.
Proof. In Theorem 3.2, we have established that for any q≥0, the sequence {rn}n∈N is G-Cauchy and by the completeness of the G-metric space (Φ,G), there exists a point c in Φ such that rn→c. To prove that c is a fixed point of Γ, suppose contrary that Γc≠c.
Assume {rn}n∈N is such that rn≠rm whenever n≠m, for all n,m∈N. Then there exists n0∈N such that rn and c are distinct and not in Fix(Γ) for all n≥n0. We will verify that c is a fixed point of Γ3. Indeed, for all n≥n0,
G(rn+2,Γ2c,Γ3c)≤α(rn,c,Γc)G(Γ2rn,Γ2c,Γ3c)≤μM(rn,c,Γc)=μG(rn,c,Γc)λ1⋅G(rn,Γrn,Γ2rn)λ2⋅G(c,Γc,Γ2c)λ3⋅G(Γrn,Γc,Γ2c)λ4⋅G(Γrn,Γ2rn,Γ3rn)λ5⋅G(Γc,Γ2c,Γ3c)δ=μG(rn,c,Γc)λ1⋅G(rn,rn+1,rn+2)λ2⋅G(c,Γc,Γ2c)λ3⋅G(rn+1,Γc,Γ2c)λ4⋅G(rn+1,rn+2,rn+3)λ5⋅G(Γc,Γ2c,Γ3c)δ. |
Since λ2>0 and λ5>0, then letting n→∞ and noting Proposition 1, it is verified that c is a fixed point of Γ3. Hence, by (3.6) and (3.7), we obtain a contradiction. Therefore Γc=c, implying that c is a fixed point of Γ.
Example 3.6. Let Φ=[−1,1] and let Γ:Φ⟶Φ be a self-mapping on Φ defined by Γr=r2 for all r∈Φ. Define G:Φ×Φ×Φ⟶R+ by
G(r,s,t)=|r−s|+|r−t|+|s−t| ∀ r,s,t∈Φ.
Then (Φ,G) is a complete G-metric space. Define α:Φ×Φ×Φ⟶R+ by
α(r,s,t)={1,ifr,s,t∈{−1}∪[0,1];0,otherwise. | (3.10) |
Then obviously, Γ is a (G-α)-orbital admissible and Γ is continuous for all r∈Φ. Also, there exists r0=13∈Φ such that α(13,Γ(13),Γ2(13))=α(13,16,112)≥1. Hence, conditions (i)-(iii) of Theorems 3.2 and 3.3 are satisfied.
To see that Γ is a hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction, let μ=12. Notice that α(r,s,Γs)=0 for all r,s∈(−1,0). Hence, inequality (3.1) holds for all r,s∈(−1,0).
Now for r,s∈{−1,1}, if r=s, then letting λ1=15, λ2=λ3=λ4=0, λ5=δ=25 and q=2, we obtain
α(r,s,Γs)G(Γ2r,Γ2s,Γ3s)=α(1,1,12)G(14,14,18)=α(−1,−1,−12)G(−14,−14,−18)=14<25=12(45)=12(M(1,1,12))=12(M(−1,−1,−12))=μM(r,s,Γs). |
Also, if q=0, we have
α(r,s,Γs)G(Γ2r,Γ2s,Γ3s)=14<12(45)=μM(r,s,Γs). |
If r≠s, then letting λ1=35, λ2=λ3=λ5=0, λ4=δ=15 and q=2, we obtain
α(r,s,Γs)G(Γ2r,Γ2s,Γ3s)=α(−1,1,12)G(−14,14,18)=α(1,−1,−12)G(14,−14,−18)=1<85=12(165)=12(M(−1,1,12))=12(M(1,−1,−12))=μM(r,s,Γs). |
Also, for q=0, we take λ1=35, λ2=λ5=δ=0, λ3=λ4=15. Then
α(r,s,Γs)G(Γ2r,Γ2s,Γ3s)=α(−1,1,12)G(−14,14,18)=α(1,−1,−12)G(14,−14,−18)=1<75=12(145)=12(M(−1,1,12))=12(M(1,−1,−12))=μM(r,s,Γs). |
Finally, for all r,s∈(0,1), we take λ1=1, λ2=λ3=λ4=λ5=δ=0. Then
α(r,s,Γs)G(Γ2r,Γ2s,Γ3s)=G(Γ2r,Γ2s,Γ3s)=|r4−s4|+|r4−s8|+|s4−s8|=14(|r−s|+|r−s2|+|s−s2|)=14G(r,s,Γs)<12G(r,s,Γs)=μM(r,s,Γs) |
for all q≥0.
Hence, inequality (3.1) is satisfied for all r,s∈Φ∖Fix(Γ). Therefore, Γ is a hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction which satisfies all the assumptions of Theorems 3.2 and 3.3. The point r=0 is the fixed point of Γ in Φ.
We now demonstrate that our result is independent and an improvement of the results of Karapınar et al. [7]. Let α:Φ×Φ⟶R+ be as given by Definition (2.12), r0∈Φ be such that α(r0,Γr0)≥1 and d:Φ×Φ⟶R+ be defined by
d(r,s)=|r−s| ∀ r,s∈Φ.
Consider r,s∈{−1,1} and take for Case 1, r≠s, λ1=15, λ2=λ3=λ5=0, λ4=12, δ=310 and q=1. Then inequality (3.1) becomes
α(r,s,Γs)G(Γ2r,Γ2s,Γ3s)=α(−1,1,12)G(−14,14,18)=α(1,−1,−12)G(14,−14,−18)=1≤101100=12(202100)=12(M(−1,1,12))=12(M(1,−1,−12))=μM(r,s,Γs), |
while inequality (2.11) due to Karapınar et al. [7] yields
α(r,s)d(Γ2r,Γ2s)=α(−1,1)d(−14,14)=α(1,−1)d(14,−14)=12>49100=12(98100)=12(M(−1,1))=12(M(1,−1))=μM(r,s). |
For Case 2, Karapınar et al. [7] have noted that their result is indeterminate for r=s, if either λ1=0 or λ4=0, since
[d(r,s)]λ1=[d(Γr,Γs)]λ4=00. |
Hence, they declared that r and s are distinct and that 00=1, which in contrast, is unconventional. But our result is valid for all r,s∈Φ∖Fix(Γ). Therefore, hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction is not hybrid-interpolative Reich-Istr˘aţescu-type contraction defined by Karapınar et al. [7], and so Theorems 12 and 17 due to Karapınar et al. [7] are not applicable to this example.
The following is an Istr˘aţescu-type (see [7]) consequence of our result.
Corollary 2. Let (Φ,G) be a complete G-metric space and Γ:Φ⟶Φ be a continuous self-mapping such that there exist η,λ∈(0,1) with η+λ<1, satisfying
G(Γ2r,Γ2s,Γ3s)≤ηG(r,s,Γs)+λG(Γr,Γs,Γ2s) |
for all r,s∈Φ. Then Γ has a fixed point in Φ.
Proof. Let μ=η+λ∈(0,1). Consider Definition (3.1) and let α(r,s,Γs)=1, q=1, λ1=ημ, λ4=λμ and λ2=λ3=λ5=δ=0. Then for all r,s∈Φ, we have
G(Γ2r,Γ2s,Γ3s)≤μM(r,s,Γs)=μ[ημG(r,s,Γs)+λμG(Γr,Γs,Γ2s)]=ηG(r,s,Γs)+λG(Γr,Γs,Γ2s), | (3.11) |
implying that inequality (3.1) holds for all r,s∈Φ∖Fix(Γ).
Let r0∈Φ be arbitrary and define a sequence {rn}n∈N in Φ by rn=Γnr0. Then by (3.11), we have
G(rn+2,rn+3,rn+4)≤ηG(rn,rn+1,rn+2)+λG(rn+1,rn+2,rn+3)≤ηmax{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}+λmax{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}=(η+λ)max{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}=μmax{G(rn,rn+1,rn+2),G(rn+1,rn+2,rn+3)}. |
Hence, by Lemma 2.11, {rn}n∈N is a G-Cauchy sequence in (Φ,G) and since (Φ,G) is complete, then there exists c∈Φ such that rn→c. Since Γ is continuous, then we can conclude that c is a fixed point of Γ, that is, Γc=c.
In Definition (3.1), if we specialize the parameters λi (i=1,2,...,5), δ and q, as well as let α(r,s,Γs)=1 for all r,s∈Φ and Γs=t, we obtain the following corollary, which is a consequence of Theorem 3.2.
Corollary 3. Let (Φ,G) be a complete G-metric space and Γ:Φ⟶Φ be a continuous self-mapping such that there exists μ∈(0,1), satisfying
G(Γ2r,Γ2s,Γ2t)≤μG(r,s,t) | (3.12) |
for all r,s,t∈Φ. Then Γ has a fixed point in Φ.
In this section, an existence theorem for a solution of a class of integral equations is provided using Corollary 3. For similar results, we refer to [1,4,10,17].
Consider the integral equation
r(y)=∫baL(y,x)f(x,r(x))dx,y∈[a,b]. | (4.1) |
Let Φ=C([a,b],R) be the set of all continuous real-valued functions. Define G:Φ×Φ×Φ⟶R+ by
G(r,s,t)=maxy∈[a,b]|r(y)−s(y)|+maxy∈[a,b]|r(y)−t(y)|+maxy∈[a,b]|s(y)−t(y)|∀r,s,t∈Φ,y∈[a,b]. | (4.2) |
Then, (Φ,G) is a complete G-metric space.
Define a function Γ:Φ⟶Φ as follows:
Γr(y)=∫baL(y,x)f(x,r(x))dx,y∈[a,b]. | (4.3) |
Then a point u∗ is said to be a fixed point of Γ if and only if u∗ is a solution to (4.1).
Now, we study existence conditions of the integral equation (4.1) under the following hypotheses.
Theorem 4.1. Assume that the following conditions are satisfied:
(C1) L:[a,b]×[a,b]⟶R+ and f:[a,b]×R⟶R are continuous;
(C2) for all r,s∈Φ, x∈[a,b], we have |f(x,r(x))−f(x,s(x))|≤|r(x)−s(x)|;
(C3) maxy∈[a,b]∫baL(y,x)dx≤λ for some λ<1.
Then, the integral equation (4.1) has a solution u∗ in Φ.
Proof. Observe that for any r,s∈Φ, using (4.3) and the above hypotheses, we obtain
|Γr(y)−Γs(y)|=|∫ba[L(y,x)f(x,r(x))−L(y,x)f(x,s(x))]dx|≤∫baL(y,x)|f(x,r(x))−f(x,s(x))|dx≤∫baL(y,x)|r(x)−s(x)|dx≤∫baL(y,x)maxx∈[a,b]|r(x)−y(x)|dx≤λmaxy∈[a,b]|r(y)−s(y)|, |
so that
|Γ2r(y)−Γ2s(y)|≤λmaxy∈[a,b]|Γr(y)−Γs(y)|≤λmaxy∈[a,b][λmaxy∈[a,b]|r(y)−s(y)|]≤λ2maxy∈[a,b]|r(y)−s(y)|. |
Using this in (4.2), we have
G(Γ2r,Γ2s,Γ2t)=maxy∈[a,b]|Γ2r−Γ2s|+maxy∈[a,b]|Γ2r−Γ2t|+maxy∈[a,b]|Γ2s−Γ2t|≤λ2maxy∈[a,b]|r−s|+λ2maxy∈[a,b]|r−t|+λ2maxy∈[a,b]|s−t|=λ2(maxy∈[a,b]|r−s|+maxy∈[a,b]|r−t|+maxy∈[a,b]|s−t|)=μG(r,s,t), |
where μ=λ2<1.
Hence, all the hypotheses of Corollary 3 are verified, implying that there exists a solution u∗ in Φ of the integral equation (4.1).
Conversely, if u∗ is a solution of (4.1), then u∗ is also a solution of (4.3), so that Γu∗=u∗, that is, u∗ is a fixed point of Γ.
Remark 1.
(i) We can deduce many other corollaries by particularizing some of the parameters in Definition (3.1).
(ii) None of the results presented in this work is expressible in the form G(s,r,r) or G(s,s,r). Hence, they cannot be obtained from their corresponding versions in metric spaces.
A generalization of metric space was introduced by Mustafa and Sims [9], namely G-metric space and several fixed point results were studied in that space. However, Jleli and Samet [5] as well as Samet et al. [16] established that most fixed point theorems obtained in G-metric space are direct consequences of their analogues in metric space. Contrary to the above observation, a new family of contraction, called hybrid-interpolative Reich-Istr˘aţescu-type (G-α-μ)-contraction is introduced in this manuscript and some fixed point theorems that cannot be deduced from their corresponding ones in metric space are proved. The main distinction of this class of contraction is that its contractive inequality is expressible in a number of ways with respect to multiple parameters. Consequently, a few corollaries, including some recently announced results in the literature are highlighted and analyzed. Nontrivial comparative examples are constructed to validate the assumptions of our obtained theorems. Furthermore, one of our obtained corollaries is applied to set up novel existence conditions for solution of a class of integral equations.
The authors declare that they have no competing interests.
[1] |
S. M. Salman, Memory and media coverage effect on an HIV/AIDS epidemic model with treatment, J. Comput. Appl. Math., 385 (2021), 113203. http://dx.doi.org/10.1016/j.cam.2020.113203 doi: 10.1016/j.cam.2020.113203
![]() |
[2] | HIV, World Health Organization, 2022. |
[3] | A. S. Perelson, Modeling the interaction of the immune system with HIV, In: Mathematical and statistical approaches to AIDS epidemiology, Berlin, Heidelberg: Springer, 1989. http://dx.doi.org/10.1007/978-3-642-93454-4_17 |
[4] | H. Ye, Modeling and analyzing of the dynamics of HIV infections based on fractional differential equations, Doctoral thesis, Donghua University, 2009. |
[5] |
R. Xu, C. Song, Dynamics of an HIV infection model with virus diffusion and latently infected cell activation, Nonlinear Anal. Real World Appl., 67 (2022), 103618. https://doi.org/10.1016/j.nonrwa.2022.103618 doi: 10.1016/j.nonrwa.2022.103618
![]() |
[6] |
P. Wu, S. Zheng, Z. He, Evolution dynamics of a time-delayed reaction-diffusion HIV latent infection model with two strains and periodic therapies, Nonlinear Anal. Real World Appl., 67 (2022), 103559. https://doi.org/10.1016/j.nonrwa.2022.103559 doi: 10.1016/j.nonrwa.2022.103559
![]() |
[7] |
P. Wu, H. Zhao, Mathematical analysis of an age-structured HIV/AIDS epidemic model with HAART and spatial diffusion, Nonlinear Anal. Real World Appl., 60 (2021), 103289. https://doi.org/10.1016/j.nonrwa.2021.103289 doi: 10.1016/j.nonrwa.2021.103289
![]() |
[8] |
B. J. Nath, K. Dehingia, K. Sadri, H. K. Sarmah, K. Hosseini, C. Park, Optimal control of combined antiretroviral therapies in an HIV infection model with cure rate and fusion effect, Int. J. Biomath, 16 (2023), 2250062. https://doi.org/10.1142/S1793524522500620 doi: 10.1142/S1793524522500620
![]() |
[9] |
Q. Liu, D. Jiang, T. Hayat, A. Alsaedi, Stationary distribution and extinction of a stochastic HIV-1 infection model with distributed delay and logistic growth, J. Nonlinear Sci., 30 (2020), 369–395. https://doi.org/10.1007/s00332-019-09576-x doi: 10.1007/s00332-019-09576-x
![]() |
[10] |
K. Qi, D. Jiang, The impact of virus carrier screening and actively seeking treatment on dynamical behavior of a stochastic HIV/AIDS infection model, Appl. Math. Model., 85 (2020), 378–404. https://doi.org/10.1016/j.apm.2020.03.027 doi: 10.1016/j.apm.2020.03.027
![]() |
[11] |
Q. Liu, Dynamics of a stochastic SICA epidemic model for HIV transmission with higher-order perturbation, Stoch. Anal. Appl., 40 (2022), 209–235. https://doi.org/10.1080/07362994.2021.1898979 doi: 10.1080/07362994.2021.1898979
![]() |
[12] |
J. Ren, Q. Zhang, X. Li, F. Cao, M. Ye, A stochastic age-structured HIV/AIDS model based on parameters estimation and its numerical calculation, Math. Comput. Simulat., 190 (2021), 159–180. https://doi.org/10.1016/j.matcom.2021.04.024 doi: 10.1016/j.matcom.2021.04.024
![]() |
[13] |
Y. Tan, Y. Cai, X. Sun, K. Wang, R. Yao, W. Wang, et al., A stochastic SICA model for HIV/AIDS transmission, Chaos Soliton Fract., 165 (2022), 112768. https://doi.org/10.1016/j.chaos.2022.112768 doi: 10.1016/j.chaos.2022.112768
![]() |
[14] |
F. Rao, J. Luo, Stochastic effects on an HIV/AIDS infection model with incomplete diagnosis, Chaos Soliton Fract., 152 (2021), 111344. https://doi.org/10.1016/j.chaos.2021.111344 doi: 10.1016/j.chaos.2021.111344
![]() |
[15] |
R. Shi, T. Lu, C. Wang, Dynamic analysis of a fractional-order model for HIV with drug-resistance and CTL immune response, Math. Comput. Simulat., 188 (2021), 509–536. https://doi.org/10.1016/j.matcom.2021.04.022 doi: 10.1016/j.matcom.2021.04.022
![]() |
[16] |
H. Singh, Analysis of drug treatment of the fractional HIV infection model of CD4+T-cells, Chaos Soliton Fract., 146 (2021), 11068. https://doi.org/10.1016/j.chaos.2021.110868 doi: 10.1016/j.chaos.2021.110868
![]() |
[17] |
Y. Zhao, E. E. Elattar, M. A. Khan, Fatmawati, M. Asiri, P. Sunthrayuth, The dynamics of the HIV/AIDS infection in the framework of piecewise fractional differential equation, Results Phys., 40 (2022), 105842. https://doi.org/10.1016/j.rinp.2022.105842 doi: 10.1016/j.rinp.2022.105842
![]() |
[18] |
M. Jafari, H. Kheiri, Free terminal time optimal control of a fractional-order model for the HIV/AIDS epidemic, Int. J. Biomath., 15 (2022), 2250022. https://doi.org/10.1142/S179352452250022X doi: 10.1142/S179352452250022X
![]() |
[19] |
B. Asquith, C. Bangham, Quantifying HTLV-I dynamics, Immunol. Cell Biol., 85 (2007), 280–286. https://doi.org/10.1038/sj.icb.7100050 doi: 10.1038/sj.icb.7100050
![]() |
[20] |
L. M. Mansky, In vivo analysis of human T-cell leukemia virus type Ⅰ reverse transcription accuracy, J. Virol., 74 (2000), 9525–9531. https://doi.org/10.1128/JVI.74.20.9525-9531.2000 doi: 10.1128/JVI.74.20.9525-9531.2000
![]() |
[21] |
Q. Kai, D. Jiang, Threshold behavior in a stochastic HTLV-I infection model with CTL immune response and regime switching, Math. Method. Appl. Sci., 41 (2018), 6866–6882. https://doi.org/10.1002/mma.5198 doi: 10.1002/mma.5198
![]() |
[22] |
Z. Shi, D. Jiang, Dynamical behaviors of a stochastic HTLV-I infection model with general infection form and Ornstein-Uhlenbeck process, Chaos Soliton Fract., 165 (2022), 112789. https://doi.org/10.1016/j.chaos.2022.112789 doi: 10.1016/j.chaos.2022.112789
![]() |
[23] |
S. Bera, S. Khajanchi, T. K. Roy, Dynamics of an HTLV-I infection model with delayed CTLs immune response, Appl. Math. Comput., 430 (2022), 127206. https://doi.org/10.1016/j.amc.2022.127206 doi: 10.1016/j.amc.2022.127206
![]() |
[24] |
A. M. Elaiw, A. S. Shflot, A. D. Hobiny, Global stability of a general HTLV-I infection model with Cytotoxic T-Lymphocyte immune response and mitotic transmission, Alexandria Eng., 67 (2023), 77–91. https://doi.org/10.1016/j.aej.2022.08.021 doi: 10.1016/j.aej.2022.08.021
![]() |
[25] |
S. Khajanchi, S. Bera, T. K. Roy, Mathematical analysis of the global dynamics of a HTLV-I infection model, considering the role of cytotoxic T-lymphocytes, Math. Comput. Simulat., 180 (2021), 354–378. https://doi.org/10.1016/j.matcom.2020.09.009 doi: 10.1016/j.matcom.2020.09.009
![]() |
[26] |
N. Kobayashi, Y. Hamamoto, N. Yamamoto, Production of tumor necrosis factors by human T cell lines infected with HTLV-1 may cause their high susceptibility to human immunodeficiency virus infection, Med. Microbiol. Immunol., 179 (1990), 115–122. https://doi.org/10.1007/BF00198532 doi: 10.1007/BF00198532
![]() |
[27] |
C. D. Mendoza, E. Caballero, A. Aguilera, R. Benito, D. Maciá, J. García-Costa, et al., HIV co-infection in HTLV-1 carriers in Spain, Virus Res., 266 (2019), 48–51. https://doi.org/10.1016/j.virusres.2019.04.004 doi: 10.1016/j.virusres.2019.04.004
![]() |
[28] |
M. A. Alshaikh, N. H. Alshamrani, A. M. Elaiw, Stability of HIV/HTLV co-infection model with effective HIV-specific antibody immune response, Results Phys., 27 (2021), 104448. https://doi.org/10.1016/j.rinp.2021.104448 doi: 10.1016/j.rinp.2021.104448
![]() |
[29] |
A. M. Elaiw, N. H. Alshamrani, Analysis of a within-host HIV/HTLV-I co-infection model with immunity, Virus Res., 295 (2021), 198204. https://doi.org/10.1016/j.virusres.2020.198204 doi: 10.1016/j.virusres.2020.198204
![]() |
[30] |
A. M. Elaiw, N. H. Alshamrani, E. Dahy, A. A. Abdellatif, Stability of within host HTLV-I/HIV-1 co-infection in the presence of macrophages, Int. J. Biomath, 16 (2023), 2250066. https://doi.org/10.1142/S1793524522500668 doi: 10.1142/S1793524522500668
![]() |
[31] |
Z. Guo, H. Huo, H. Xiang, Optimal control of TB transmission based on an age structured HIV-TB co-infection model, J. Franklin Inst., 359 (2022), 4116–4137. https://doi.org/10.1016/j.jfranklin.2022.04.005 doi: 10.1016/j.jfranklin.2022.04.005
![]() |
[32] |
A. Mallela, S. Lenhart, N. K. Vaidya, HIV-TB co-infection treatment: Modeling and optimal control theory perspectives, J. Comput. Appl. Math., 307 (2016), 143–161. https://doi.org/10.1016/j.cam.2016.02.051 doi: 10.1016/j.cam.2016.02.051
![]() |
[33] |
Tanvi, R. Aggarwal, Stability analysis of a delayed HIV-TB co-infection model in resource limitation settings, Chaos Soliton Fract., 140 (2020), 110138. https://doi.org/10.1016/j.chaos.2020.110138 doi: 10.1016/j.chaos.2020.110138
![]() |
[34] |
L. Zhang, M. U. Rahman, M. Arfan, A. Ali, Investigation of mathematical model of transmission co-infection TB in HIV community with a non-singular kernel, Results Phys., 28 (2021), 104559. https://doi.org/10.1016/j.rinp.2021.104559 doi: 10.1016/j.rinp.2021.104559
![]() |
[35] |
N. H. Shah, N. Sheoran, Y. Shah, Dynamics of HIV-TB co-infection model, Axioms, 9 (2020), 29. https://doi.org/10.3390/axioms9010029 doi: 10.3390/axioms9010029
![]() |
[36] |
I. Ahmed, E. F. D. Goufo, A. Yusuf, P. Kumam, K. Nonlaopon, An epidemic prediction from analysis of a combined HIV-COVID-19 co-infection model via ABC-fractional operator, Alexandria Eng. J., 60 (2021), 2979–2995. https://doi.org/10.1016/j.aej.2021.01.041 doi: 10.1016/j.aej.2021.01.041
![]() |
[37] |
N. Ringa, M. L. Diagne, H. Rwezaura, A. Omame, S. Y. Tchoumi, J. M. Tchuenche, HIV and COVID-19 co-infection: A mathematical model and optimal control, Inform. Med. Unlocked, 31 (2022), 100978. https://doi.org/10.1016/j.imu.2022.100978 doi: 10.1016/j.imu.2022.100978
![]() |
[38] |
A. Omame, M. E. Isah, M. Abbas, A. H. Abdel-Aty, C. P. Onyenegecha, A fractional order model for Dual Variants of COVID-19 and HIV co-infection via Atangana-Baleanu derivative, Alexandria Eng. J., 61 (2022), 9715–9731. https://doi.org/10.1016/j.aej.2022.03.013 doi: 10.1016/j.aej.2022.03.013
![]() |
[39] |
V. Lakshmikantham, A. S. Vatsala, Basic theory of fractional differential equations, Nonlinear Anal. Theor., 69 (2008), 2677–2682. https://doi.org/10.1016/j.na.2007.08.042 doi: 10.1016/j.na.2007.08.042
![]() |
[40] | I. Podlubny, Fractional differential equations: An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications, San Diego: Academic Press, 1999. |
[41] |
K. Hattaf, A new generalized definition of fractional derivative with non-singular kernel, Computation, 8 (2020), 49. https://doi.org/10.3390/computation8020049 doi: 10.3390/computation8020049
![]() |
[42] |
M. Bachraoui, K. Hattaf, N. Yousfi, Analysis of a fractional reaction-diffusion HBV model with cure of infected cells, Discrete Dyn. Nat. Soc., 2020 (2020), 3140275. https://doi.org/10.1155/2020/3140275 doi: 10.1155/2020/3140275
![]() |
[43] |
C. Huang, L. Cai, J. Cao, Linear control for synchronization of a fractional-order time-delayed chaotic financial system, Chaos Soliton Fract., 113 (2018), 326–332. https://doi.org/10.1016/j.chaos.2018.05.022 doi: 10.1016/j.chaos.2018.05.022
![]() |
[44] |
R. Rakkiyappan, G. Velmurugan, J. Cao, Stability analysis of fractional-order complex-valued neural networks with time delays, Chaos Soliton Fract., 78 (2015), 297–316. https://doi.org/10.1016/j.chaos.2015.08.003 doi: 10.1016/j.chaos.2015.08.003
![]() |
[45] |
H. Li, Y. Shen, Y. Han, J. Dong, J. Li, Determining Lyapunov exponents of fractional-order systems: A general method based on memory principle, Chaos Soliton Fract., 168 (2023), 113167. https://doi.org/10.1016/j.chaos.2023.113167 doi: 10.1016/j.chaos.2023.113167
![]() |
[46] |
Q. Gao, J. Cai, Y. Liu, Y. Chen, L. Shi, W. Xu, Power mapping-based stability analysis and order adjustment control for fractional-order multiple delayed systems, ISA Trans., 138 (2023), 10–19. https://doi.org/10.1016/j.isatra.2023.02.019 doi: 10.1016/j.isatra.2023.02.019
![]() |
[47] |
C. Pinto, A. Carvalho, The role of synaptic transmission in a HIV model with memory, Appl. Math. Comput., 292 (2017), 76–95. https://doi.org/10.1016/j.amc.2016.07.031 doi: 10.1016/j.amc.2016.07.031
![]() |
[48] |
T. Sardar, S. Rana, S. Bhattacharya, K. Al-Khaled, J. Chattopadhyay, A generic model for a single strain mosquito-transmitted disease with memory on the host and the vector, Math. Biosci., 263 (2015), 18–36. https://doi.org/10.1016/j.mbs.2015.01.009 doi: 10.1016/j.mbs.2015.01.009
![]() |
[49] |
K. Diethelm, Monotonicity of functions and sign changes of their Caputo derivatives, Fract. Calc. Appl. Anal., 19 (2016), 561–566. https://doi.org/10.1515/FCA-2016-0029 doi: 10.1515/FCA-2016-0029
![]() |
[50] | C. Kou, Y. Yan, J. Liu, Stability analysis for fractional differential equations and their applications in the models of HIV-1 infection, Comput. Model. Eng. Sci., 39 (2009), 301–317. |
[51] |
W. Lin, Global existence theory and chaos control of fractional differential equations, J. Math. Anal. Appl., 332 (2007), 709–726. https://doi.org/10.1016/j.jmaa.2006.10.040 doi: 10.1016/j.jmaa.2006.10.040
![]() |
[52] |
P. Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental systems of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
![]() |
[53] |
E. Ahmed, A. M. A. El-Sayed, H. A. A. El-Saka, On some Routh-Hurwitz conditions for fractional order differential equations and their applications in Lorenz, Rssler, Chua and Chen systems, Phys. Lett. A, 358 (2006), 1–4. https://doi.org/10.1016/j.physleta.2006.04.087 doi: 10.1016/j.physleta.2006.04.087
![]() |
[54] |
J. P. LaSalle, Stability theory for ordinary differential equations, J. Differ. Equ., 4 (1968), 57–65. https://doi.org/10.1016/0022-0396(68)90048-X doi: 10.1016/0022-0396(68)90048-X
![]() |
[55] |
R. Shi, T. Lu, Dynamic analysis and optimal control of a fractional order model for hand-foot-mouth Disease, J. Appl. Math. Comput., 64 (2020), 565–590. https://doi.org/10.1007/s12190-020-01369-w doi: 10.1007/s12190-020-01369-w
![]() |
[56] |
E. Roxin, Differential equations: Classical to controlled, Am. Math. Mon., 92 (1985), 223–225. https://doi.org/10.1080/00029890.1985.11971586 doi: 10.1080/00029890.1985.11971586
![]() |
[57] | L. S. Pontryagin, V. G. Boltyanskii, R. V. Gramkrelidze, E. F. Mischenko, The mathematical theory of optimal processes, New York: Interscience Publishers, 1962. |
[58] |
N. H. Sweilam, S. M. Al-Mekhlafi, On the optimal control for fractional multi-strain TB model, Optim. Contr. Appl. Met., 37 (2016), 1355–1374. https://doi.org/10.1002/oca.2247 doi: 10.1002/oca.2247
![]() |
[59] | L. Zhang, HIV viral load and CD4+T lymphocyte count in HIV-1/HTLV-1 co-infected patients, Foreign Med. Sci. Sect. Virol., 5 (1998), 27–29. |
1. | Chander Bhan, Ravi Karwasra, Sandeep Malik, Sachin Kumar, Ahmed H. Arnous, Nehad Ali Shah, Jae Dong Chung, Bifurcation, chaotic behavior and soliton solutions to the KP-BBM equation through new Kudryashov and generalized Arnous methods, 2024, 9, 2473-6988, 8749, 10.3934/math.2024424 | |
2. | Chong-Dong Cheng, Bo Tian, Tian-Yu Zhou, Yuan Shen, Nonlinear localized waves and their interactions for a (2+1)-dimensional extended Bogoyavlenskii-Kadomtsev-Petviashvili equation in a fluid, 2024, 125, 01652125, 103246, 10.1016/j.wavemoti.2023.103246 |