Next Article in Journal
Survival Risk Prediction of Esophageal Cancer Based on the Kohonen Network Clustering Algorithm and Kernel Extreme Learning Machine
Next Article in Special Issue
Numerical Studies of Channel Management Strategies for Nonstationary Immersion Environments: EURUSD Case Study
Previous Article in Journal
Robust State Estimation for Uncertain Discrete Linear Systems with Delayed Measurements
Previous Article in Special Issue
Finite-Time Stability Analysis of Linear Differential Systems with Pure Delay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on Dynamics of CD4+ T-Cells under the Effect of HIV-1 Infection Based on a Mathematical Fractal-Fractional Model via the Adams-Bashforth Scheme and Newton Polynomials

by
Hashem Najafi
1,†,
Sina Etemad
2,†,
Nichaphat Patanarapeelert
3,†,
Joshua Kiddy K. Asamoah
4,†,
Shahram Rezapour
2,5,*,† and
Thanin Sitthiwirattham
6,*,†
1
Department of Mathematics, College of Sciences, Shiraz University, Shiraz 7187919556, Iran
2
Department of Mathematics, Azarbaijan Shahid Madani University, Tabriz 3751-71379, Iran
3
Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
4
Department of Mathematics, Kwame Nkrumah University of Science and Technology, Kumasi, Ghana
5
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
6
Mathematics Department, Faculty of Science and Technology, Suan Dusit University, Bangkok 10300, Thailand
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2022, 10(9), 1366; https://doi.org/10.3390/math10091366
Submission received: 19 February 2022 / Revised: 3 April 2022 / Accepted: 5 April 2022 / Published: 19 April 2022

Abstract

:
In recent decades, AIDS has been one of the main challenges facing the medical community around the world. Due to the large human deaths of this disease, researchers have tried to study the dynamic behaviors of the infectious factor of this disease in the form of mathematical models in addition to clinical trials. In this paper, we study a new mathematical model in which the dynamics of CD 4 + T-cells under the effect of HIV-1 infection are investigated in the context of a generalized fractal-fractional structure for the first time. The kernel of these new fractal-fractional operators is of the generalized Mittag-Leffler type. From an analytical point of view, we first derive some results on the existence theory and then the uniqueness criterion. After that, the stability of the given fractal-fractional system is reviewed under four different cases. Next, from a numerical point of view, we obtain two numerical algorithms for approximating the solutions of the system via the Adams-Bashforth method and Newton polynomials method. We simulate our results via these two algorithms and compare both of them. The numerical results reveal some stability and a situation of lacking a visible order in the early days of the disease dynamics when one uses the Newton polynomial.

1. Introduction

According to medical definitions and clinical findings and virology, human immunodeficiency virus (HIV) is a type of retrovirus that leads to acquired immunodeficiency syndrome (AIDS) in humans [1]. In fact, CD 4 + T-lymphocytes are the largest number of white blood cells in the human immune system that are attacked by HIV viruses, which attack CD 4 + T-cells and infect them, reducing their number and efficiency. They disrupt cells. Therefore, this process reduces the resistance of the immune system in the human body and weakens it [1]. Although the HIV virus infects other cells, it causes the most damage to T cells by causing the degradation and destruction of CD 4 + T-cells. As a result, the affected person’s body gradually becomes sensitive to various types of infections and contaminants and completely loses its immunity. Usually, according to laboratory results, a healthy person has white blood cells with a normal CD 4 + T-lymphocyte count of 800 to 1200/mm3. The decrease in CD 4 + T-lymphocyte count may be due to thymus insufficiency or bone defects. One way to determine whether a person is infected with the virus is recognized through very small number of CD 4 + T-lymphocytes (less than 200/mm3).
Currently, AIDS is one of the most unknown and dangerous diseases of our time. According to UNAIDS 2017, which is published annually, about 36.7 million people worldwide are living with HIV, and nearly 1.8 million people have recently been infected with HIV, and almost one million people worldwide are affected by various complication, and died of AIDS in 2016. There are even several countries, especially in Africa and some less developed countries, where up to thirty-five percent of the population between the ages of fifteen and fifty are infected with the HIV virus [2]. Despite significant advances in the control of the disease, no vaccine for HIV has yet been found.
In recent decades, many attempts have been made to design, analyze, simulate and solve mathematical dynamic models, which include the basic rules in the analysis of control and prevention of the spread of various diseases and infections. In this regard, biology and engineering are among the most widely used fields for various types of mathematical modeling, which are particularly popular among researchers. For such a purpose, derivation operators and differential equations play an important role. To demonstrate such high efficiency, we can refer to a sample of articles in which various types of mathematical models are observed, including time delying model of COVID-19 [3], anthrax model for animals [4], Hepatitis C model [5], memristor chaotic model of circuit [6], Lorenz-Stenflo hyperchaotic model [7], dynamics of environmental persistence of infections [8], Langevin model [9], Mump virus [10], Zika virus [11], mosaic disease [12], Computer viruses [13], thermostat control model [14,15], pantograph model [16,17], canine distemper virus [18], Lassa fever [19], hybrid model of p-Laplacian operators [20], co-dynamics of COVID-19 and diabetes [21], checmical modeling of cyclohexane [22,23], Navier systems [24], etc.
Designing mathematical models to analyze the dynamics of HIV infection is also known as a valuable and efficient measure in this regard [25,26,27]. Usually in almost all of these dynamic mathematical models, we can see the association between HIV viruses and non-infected CD 4 + T-cells and the effects of a variety of drug therapies and clinical therapies on reducing and controlling infected cells. Perlson [28] designed and presented a simple model for primary HIV infection in 1989. This model, which is related to HIV infection, is important, and many other models have been developed later, which are inspired by this model. Four years later, Perlson et al. generalized the model and examined some of the behaviors and characteristics of the new model [29]. They considered the model by defining four different state functions: virus population, non-infected CD 4 + T-cells, productive infected CD 4 + T-cells, and latent infected CD 4 + T-cells.
In recent years, some other mathematicians developed the basic integer order models of HIV infection in the context of the fractional order systems. In [30], Ding et al. introduced a fractional version of the infection of HIV for CD 4 + T-cells and analyzed the non-negative solutions of this system. Arafa et al. [31] used the generalized Euler technique for his model to obtain solutions and investigated the impact of changes of viral particles on the blood. After that, Bulut et al. [32] combined the homotopy and sumudu techniques to staudy the dynamics of their fractional model of the infection of HIV on CD 4 + T-cells. Lichae et al. [33] extended their fractional HIV model based on the effect of drug therapy and solved three-compartmental model of CD4+ with the help of LADM (Laplace Adomian decomposition method). In 2020, Nazir et al. [34] investigated not only existenece theory, but also stability criteria for the Caputo-Fabrizo HIV model of CD 4 + T-cells. Wang et al. [35] studied the time periodic reaction-diffusion equations for modeling 2-LTR dynamics in HIV-infected patients.
In 2017, Atangana [36] presented a generalized kind of operators entitled fractal-fractional operators. This definition used the existing notions of fractional calculus and fractal calculus with together. These operators are the convolution of the power-law, exponential law and generalized Mittag–Leffler law with fractal derivatives. There exist two components for such fractal-fractional operators: the fractional order and fractal dimension (order). Fractal-fractional differential equations transfer the order and dimension of every dynamical system into a rational order one. In fact, we are able to generalize each standard differential equation to the generalized systems having arbitrary order and dimension of derivatives. The main goal of such a combination is to analyze a vast range of nonlocal BVPs or IVPs that possess fractal behaviors. In this direction, a limited researchers obtained some results in which we see that the generalized fractal-fractional operators give accurate and more exact simulations for describing mathematical models of real-world phenomena. Some of new works in this regard are [37,38,39,40,41].
Due to the novelty and efficacy of these new fractal-fractional operators, in this paper, we aim to design a mathematical model of CD 4 + T-cells under the effect of HIV-1 infection in which derivatives are fractal-fractional operators in the sense of Atangana-Baleanu. It is notable that in 2021, Ahmad et al. [42] studied the dynamics of HIV primary infection in the context of a model designed by the fractal-fractional operators. The main contribution and novelty of our work in comparison to their paper [42] is that we analyze all qualitative behaviors of such a fractal-fractional system. In other words, we first review the existence and uniqueness of solutions and further, we complete our study by giving new results about the stability (Ulam-Hyers-Rassias) of solutions. Also, for the first time, in this paper, we derive numerical schemes for the fractal-fractional CD 4 + -HIV-1 model with the help of the Newton polynomials and by applying some real data, we compare our results with the Adams-Bashforth simulations. In this direction, we can see some dynamical behaviors of the solutions in our simulations.
The arrangement of the paper is as follows: we introduce our fractal-fractional model in Section 2 and describe its parameters and coefficients. The existence results are given in Section 3 and further, in Section 4, we investigate the uniqueness. The stability criterion are implemented in Section 5. In the sequel, we derive two numerical schemes. In other words, in Section 6, the Adams-Bashforth method are done and then, we derive another algorithm in Section 7 by using the Newton polynomials. We present some simulations and discussions about both numerical methods in Section 8. We end the paper by giving conclusions in Section 9.

2. The Structure of the Model for CD 4 + T-Cells and HIV-1

In 2006, Wang and Li [1] formulated an integer-order classical mathematical structure of dynamics of CD 4 + T-cells under the HIV-1 infection in three-compartmental model as
T ( s ) = θ q V T ρ T + ζ U , U ( s ) = q V T ( ζ + κ ) U , V ( s ) = κ N U ϑ V ,
via the initial values T ( 0 ) = T 0 , U ( 0 ) = U 0 , and V ( 0 ) = V 0 , and also the state functions T ( s ) , U ( s ) , V ( s ) are an amount of susceptible CD 4 + T-cells, an amount of infectious CD 4 + T-cells, and the free particles of the infection of the HIV virus in the blood at the time s J : = [ 0 , S ] , ( S > 0 ) , respectively. Moreover, the parameter N stands for the average number of infected particles by an existing infected cell, ϑ is the natural rate of death for the virus, ζ is the return rate of infected cells to susceptible compartment, κ is the rate of death for infected T-cells, q is the rate of infection T-cells, ρ is the natural rate of death, and θ shows the supply rate for new T-cells. They considered all parameters as positive values and T 0 , U 0 , V 0 0 .
To upgrade and improve the exact results, inspired by the standard model (1), we present a mathematical fractal-fractional model on dynamics of CD 4 + T-cells under the effect of HIV-1 infection via the generalized Mittag-Leffler-type kernel (fractal-fractional CD 4 + -HIV-1-model) as
FFML D 0 , s ( δ , σ ) T ( s ) = θ q V ( s ) T ( s ) ρ T ( s ) + ζ U ( s ) , FFML D 0 , s ( δ , σ ) U ( s ) = q V ( s ) T ( s ) ( ζ + κ ) U ( s ) , FFML D 0 , s ( δ , σ ) V ( s ) = κ N U ( s ) ϑ V ( s ) ,
subject to
T ( 0 ) = T 0 0 , U ( 0 ) = U 0 0 , V ( 0 ) = V 0 0 ,
where all assumptions and parameters are similar to above classical model (1). Also, FFML D 0 , s ( δ , σ ) is the ( δ , σ ) -fractal-fractional derivative with the fractional order δ ( 0 , 1 ] and the fractal order σ ( 0 , 1 ] via the Mittag-Leffler-type kernel.
In other words, let a continuous map Ψ : ( a , b ) [ 0 , ) be fractal differentiable of dimension σ . In this case, the ( δ , σ ) -fractal-fractional derivative of Ψ of the generalized Mittag-Leffler-type kernel of order δ in the Riemann-Liouville sense is defined as
FFML D a , s ( δ , σ ) Ψ ( s ) = AB ( δ ) 1 δ d d s σ a s E δ δ 1 δ ( s w ) δ Ψ ( w ) d w , 0 < δ , σ 1 ,
where
d Ψ ( w ) d w σ = lim s w Ψ ( s ) Ψ ( w ) s σ w σ ,
is the fractal derivative and AB ( δ ) = 1 δ + δ Γ ( δ ) and AB ( 0 ) = AB ( 1 ) = 1 [36].
Moreover, accordingly, for such a function Ψ , the ( δ , σ ) -fractal-fractional integral via the Mittag-Leffler-type kernel is given by
FFML I a , s ( δ , σ ) Ψ ( s ) = δ σ AB ( δ ) Γ ( δ ) a s w σ 1 ( s w ) δ 1 Ψ ( w ) d w + ( 1 δ ) σ s σ 1 AB ( δ ) Ψ ( s ) ,
if it exists, where δ , σ > 0 [36].

3. Existence Property

In this section, the existence property is investigated based on fixed point theory. For the qualitative analysis, make the Banach space X = M 3 , where M = C ( J , R ) with
K X = T , U , V X = max | W ( s ) | : s J ,
for | W | : = | T | + | U | + | V | . We reformulate the R.H.S. of the fractal-fractional CD 4 + -HIV-1-model (2) as:
Q 1 s , T ( s ) , U ( s ) , V ( s ) = θ q V ( s ) T ( s ) ρ T ( s ) + ζ U ( s ) , Q 2 s , T ( s ) , U ( s ) , V ( s ) = q V ( s ) T ( s ) ( ζ + κ ) U ( s ) , Q 3 s , T ( s ) , U ( s ) , V ( s ) = κ N U ( s ) ϑ V ( s ) .
In this case, the fractal-fractional CD 4 + -HIV-1-model (2) is transformed into the following system
ABR D 0 , s δ T ( s ) = σ s σ 1 Q 1 s , T ( s ) , U ( s ) , V ( s ) , ABR D 0 , s δ U ( s ) = σ s σ 1 Q 2 s , T ( s ) , U ( s ) , V ( s ) , ABR D 0 , s δ V ( s ) = σ s σ 1 Q 3 s , T ( s ) , U ( s ) , V ( s ) .
In view of (6), we rewrite the developed tripled-system with the compact IVP which takes the form
ABR D 0 , s δ K ( s ) = σ s σ 1 Q s , K ( s ) , K ( 0 ) = K 0 ,
where
K ( s ) = T ( s ) , U ( s ) , V ( s ) T , K 0 = T 0 , U 0 , V 0 T , δ , σ ( 0 , 1 ] ,
and
Q s , K ( s ) = Q 1 s , T ( s ) , U ( s ) , V ( s ) , Q 2 s , T ( s ) , U ( s ) , V ( s ) , Q 3 s , T ( s ) , U ( s ) , V ( s ) , s J .
By the definition and by (7), we have
AB ( δ ) 1 δ d d s 0 s E δ δ 1 δ ( s w ) δ K ( w ) d w = σ s σ 1 Q s , K ( s ) .
In the sequel, operating the fractal-fractional Atangana-Baleanu integral on (10), we get
K ( s ) = K ( 0 ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q ( w , K ( w ) ) d w + ( 1 δ ) σ s σ 1 AB ( δ ) Q ( s , K ( s ) ) .
Due to the above compact form of the fractal-fractional integral equation, the extended representation of it is illustrated as
T ( s ) = T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) d w , U ( s ) = U 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 2 ( s , T ( s ) , U ( s ) , V ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 2 ( w , T ( w ) , U ( w ) , V ( w ) ) d w , V ( s ) = V 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 3 ( s , T ( s ) , U ( s ) , V ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 3 ( w , T ( w ) , U ( w ) , V ( w ) ) d w .
We consider a self-map to derive a fixed-point problem, by defining F : X X as
F ( K ( s ) ) = K ( 0 ) + ( 1 δ ) σ s σ 1 AB ( δ ) Q ( s , K ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q ( w , K ( w ) ) d w .
The Leray-Schauder theorem is utilized to prove the existence property in relation to the fractal-fractional CD 4 + -HIV-1-model (2).
Theorem 1.
[43] Let X be a Banach space, E X a convex closed bounded set, O E an open set, and 0 O . Then for the continuous and compact map F : O ¯ E , either:
(HY1)
y O ¯ s . t . y = F ( y ) , or
(HY2)
y O and 0 < μ < 1 s.t. y = μ F ( y ) .
Theorem 2.
Let Q C ( J × X , X ) . Assume:
(P1) 
φ L 1 ( J , R + ) and A C ( [ 0 , ) , ( 0 , ) ) (A is non-decreasing) s.t. s J and K X ,
| Q ( s , K ( s ) ) | φ ( s ) A ( | K ( s ) | ) ;
(P2) 
ω > 0 s.t.
ω K 0 + ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) φ 0 * A ( ω ) > 1 ,
with φ 0 * = sup s J | φ ( s ) | .
Then there is a solution for the fractal-fractional problem (7) and accordingly, for the fractal-fractional CD4 + -HIV-1-model (2) on J .
Proof. 
First, consider F : X X which is defined by (13) and assume
N r = K X : K X r ,
for some r > 0 . Clearly, as Q is continuous, thus F is so. From (P1), we get
| F ( K ( s ) ) | | K ( 0 ) | + ( 1 δ ) σ s σ 1 AB ( δ ) | Q ( s , K ( s ) ) | + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 | Q ( w , K ( w ) ) | d w K 0 + ( 1 δ ) σ s σ 1 AB ( δ ) φ ( s ) A ( | K ( s ) | ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 φ ( w ) A ( | K ( w ) | ) d w K 0 + ( 1 δ ) σ S σ 1 AB ( δ ) φ 0 * A ( r ) + δ σ S δ + σ 1 B ( δ , σ ) AB ( δ ) Γ ( δ ) φ 0 * A ( r ) = K 0 + ( 1 δ ) σ S σ 1 AB ( δ ) φ 0 * A ( r ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) φ 0 * A ( r ) ,
for K N r . Hence
F K X K 0 + ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) φ 0 * A ( r ) < .
Thus F is uniformly bounded on X . Now, take s , v [ 0 , S ] s.t. s < v and K N r . By denoting
sup ( s , K ) J × N r | Q ( s , K ( s ) ) | = Q * < ,
we estimate
| F ( K ( v ) ) F ( K ( s ) ) | | ( 1 δ ) σ v σ 1 AB ( δ ) Q ( v , K ( v ) ) ( 1 δ ) σ s σ 1 AB ( δ ) Q ( s , K ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 v w σ 1 ( v w ) δ 1 Q ( w , K ( w ) ) d w δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q ( w , K ( w ) ) d w | ( 1 δ ) σ Q * AB ( δ ) ( v σ 1 s σ 1 ) + δ σ Q * AB ( δ ) Γ ( δ ) | 0 v w σ 1 ( v w ) δ 1 d w 0 s w σ 1 ( s w ) δ 1 d w | ( 1 δ ) σ Q * AB ( δ ) ( v σ 1 s σ 1 ) + δ σ Q * B ( δ , σ ) AB ( δ ) Γ ( δ ) v δ + σ 1 s δ + σ 1 = ( 1 δ ) σ Q * AB ( δ ) ( v σ 1 s σ 1 ) + δ σ Q * Γ ( σ ) AB ( δ ) Γ ( δ + σ ) v δ + σ 1 s δ + σ 1 .
We see that (16) approaches to 0 independent of K , as v s . Consequently
F ( K ( v ) ) F ( K ( s ) ) X 0 ,
when v s . This gives the equicontinuity of F, and accordingly the compactness of F on N r by the Arzelá–Ascoli thoerem. Since all conditions of Theorem 1 are fulfilled on F, we have one of (HY1) or (HY2). From (P2), we set
Φ : = K X : K X < ω ,
for some ω > 0 s.t.
K 0 + ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) φ 0 * A ( ω ) < ω .
From (P1) and (15), we have
F K X K 0 + ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) φ 0 * A ( K X ) .
Suppose that there are K Φ and 0 < μ < 1 s.t. K = μ F ( K ) . Then by (17), we write
ω = K X = μ F K X < K 0 + ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) φ 0 * A ( K X ) < K 0 + ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) φ 0 * A ( ω ) < ω ,
which cannot be held. Therefore, (HY2) does not hold and F possesses a fixed–point in Φ ¯ by Theorem 1. This confirms the existence of a solution to the fractal-fractional CD 4 + -HIV-1-model (2). □

4. Uniqueness Property

Here, on the fractal-fractional CD 4 + -HIV-1-model (2), we investigate the Lipschitz property in the first step and further, the uniqueness property.
Lemma 1.
Consider T , U , V , T * , U * , V * M : = C ( J , R ) , and let
(C1) 
T β 1 , U β 2 , V β 3 for some constants β 1 , β 2 , β 3 > 0 .
Then the kernels Q 1 , Q 2 , Q 3 defined in (5) are fulfilled the Lipschitz property with constants α 1 , α 2 , α 3 > 0 w.r.t. the relevant components, where
α 1 = ρ + q β 3 , α 2 = ζ + κ , α 3 = ϑ .
Proof. 
For Q 1 , we take T , T * M : = C ( J , R ) arbitrarily, and we have
Q 1 s , T ( s ) , U ( s ) , V ( s ) Q 1 s , T * ( s ) , U ( s ) , V ( s ) = θ q V ( s ) T ( s ) ρ T ( s ) + ζ U ( s ) θ q V ( s ) T * ( s ) ρ T * ( s ) + ζ U ( s ) ρ + q V ( s ) T ( s ) T * ( s ) ρ + q β 3 T ( s ) T * ( s ) = α 1 T ( s ) T * ( s ) .
From the above, we find out that Q 1 is Lipschitz w.r.t. T under the constant α 1 > 0 . For Q 2 , we choose two arbitrary elements U , U * M : = C ( J , R ) , and estimate
Q 2 s , T ( s ) , U ( s ) , V ( s ) Q 2 s , T ( s ) , U * ( s ) , V ( s ) = q V ( s ) T ( s ) ( ζ + κ ) U ( s ) q V ( s ) T ( s ) ( ζ + κ ) U * ( s ) ζ + κ U ( s ) U * ( s ) = α 2 U ( s ) U * ( s ) .
The above inequality means that Q 2 is Lipschitz w.r.t. U under the constant α 2 > 0 . Finally, for both arbitrary elemets V , V * M : = C ( J , R ) , we have
Q 3 s , T ( s ) , U ( s ) , V ( s ) Q 3 s , T ( s ) , U ( s ) , V * ( s ) = κ N U ( s ) ϑ V ( s ) κ N U ( s ) ϑ V * ( s ) ϑ V ( s ) V * ( s ) = α 3 V ( s ) V * ( s ) .
This confirms that Q 3 is Lipschitz w.r.t. V under the constant α 3 > 0 . Therefore kernel functions Q 1 , Q 2 , Q 3 are Lipschitze w.r.t. under the constants α 1 , α 2 , α 3 > 0 , respectively. □
Now, by invoking the above lemma, we are able to prove the uniqueness property for solutions of the fractal-fractional system (2).
Theorem 3.
Consider (C1) to be held. Then the fractal-fractional CD 4 + -HIV-1-model (2) possesses exactly one solution if
( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α j < 1 ,
for j { 1 , 2 , 3 } , and the Lipschitz constants α j > 0 introduced by (18).
Proof. 
To prove the desired result, we consider the contrary of the conclusion of theorem. That is, consider the existence of another solution for the fractal-fractional CD 4 + -HIV-1-model (2), namely T * ( s ) , U * ( s ) , V * ( s ) under the initial conditions
T * ( 0 ) = T 0 , U * ( 0 ) = U 0 , V * ( 0 ) = V 0 .
From (12), we have
T * ( s ) = T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ,
U * ( s ) = U 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 2 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 2 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ,
and
V * ( s ) = V 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 3 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 3 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w .
In this case, we estimate
| T ( s ) T * ( s ) | ( 1 δ ) σ s σ 1 AB ( δ ) | Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) | + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 × | Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) | d w ( 1 δ ) σ s σ 1 AB ( δ ) α 1 T T * + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 α 1 T T * d w ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 T T * ,
and so
1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 T T * 0 .
From (19), we know that the above inequality holds if T T * = 0 , or T = T * . In the similar manner, from
U U * ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 2 U U * ,
we obtain
1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 2 U U * 0 ,
and this gives U U * = 0 or U = U * . Further,
V V * ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 3 V V * ,
yields
1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 3 V V * 0 .
Hence V = V * . In consequence,
T ( s ) , U ( s ) , V ( s ) = T * ( s ) , U * ( s ) , V * ( s ) .
The latter equality confirms that the fractal-fractional CD 4 + -HIV-1-model (2) possesses a unique solution. □

5. Ulam-Hyers-Rassias Stability

In this part, the stability result of solutions in the context of four types of the Ulam–Hyers, Ulam–Hyers–Rassias and their generalizations are proved on the tripled system of the fractal-fractional CD 4 + -HIV-1-model (2).
Definition 1.
The fractal-fractional CD 4 + -HIV-1-model (2) is Ulam–Hyers stable if a Q 1 , a Q 2 , a Q 3 R + such that r j > 0 , j = 1 , 2 , 3 , and T * , U * , V * X satisfying
| FFML D 0 , s ( δ , σ ) T * ( s ) Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) | < r 1 , | FFML D 0 , s ( δ , σ ) U * ( s ) Q 2 s , T * ( s ) , U * ( s ) , V * ( s ) | < r 2 , | FFML D 0 , s ( δ , σ ) V * ( s ) Q 3 s , T * ( s ) , U * ( s ) , V * ( s ) | < r 3 ,
T , U , V X satisfying the fractal-fractional CD 4 + -HIV-1-model (2) s.t.
| T * ( s ) T ( s ) | a Q 1 r 1 , | U * ( s ) U ( s ) | a Q 2 r 2 , | V * ( s ) V ( s ) | a Q 3 r 3 .
Definition 2.
The fractal-fractional CD 4 + -HIV-1-model (2) is generalized Ulam–Hyers stable if a Q j C ( R + , R + ) , ( j { 1 , 2 , 3 } ) with a Q j ( 0 ) = 0 s.t. r j > 0 and T * , U * , V * X fulfilling (20), there is T , U , V X as a solution of the given fractal-fractional CD 4 + -HIV-1-model (2) s.t.
| T * ( s ) T ( s ) | a Q 1 ( r 1 ) , | U * ( s ) U ( s ) | a Q 2 ( r 2 ) , | V * ( s ) V ( s ) | a Q 3 ( r 3 ) .
Remark 1.
T * , U * , V * X is a solution for (20) iff z 1 , z 2 , z 3 C ( [ 0 , S ] , R ) (each of them depend on T * , U * , V * , respectively) s.t. s J ,
(i) 
| z j ( s ) | < r j ;
(ii) 
We have
FFML D 0 , s ( δ , σ ) T * ( s ) = Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) + z 1 ( s ) , FFML D 0 , s ( δ , σ ) U * ( s ) = Q 2 s , T * ( s ) , U * ( s ) , V * ( s ) + z 2 ( s ) , FFML D 0 , s ( δ , σ ) V * ( s ) = Q 3 s , T * ( s ) , U * ( s ) , V * ( s ) + z 3 ( s ) .
Definition 3.
The fractal-fractional CD 4 + -HIV-1-model (2) is Ulam–Hyers–Rassias stable w.r.t. functions ħ j , ( j { 1 , 2 , 3 } ) if 0 < a ( Q j , ħ j ) R s.t. r j > 0 and T * , U * , V * X fulfilling
| FFML D 0 , s ( δ , σ ) T * ( s ) Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) | < r 1 ħ 1 ( s ) , | FFML D 0 , s ( δ , σ ) U * ( s ) Q 2 s , T * ( s ) , U * ( s ) , V * ( s ) | < r 2 ħ 2 ( s ) , | FFML D 0 , s ( δ , σ ) V * ( s ) Q 3 s , T * ( s ) , U * ( s ) , V * ( s ) | < r 3 ħ 3 ( s ) ,
T , U , V X as a solution of the fractal-fractional CD 4 + -HIV-1-model (2) s.t.
| T * ( s ) T ( s ) | r 1 a ( Q 1 , ħ 1 ) ħ 1 ( s ) , s J , | U * ( s ) U ( s ) | r 2 a ( Q 2 , ħ 2 ) ħ 2 ( s ) , s J , | V * ( s ) V ( s ) | r 3 a ( Q 3 , ħ 3 ) ħ 3 ( s ) , s J ,
in which ħ 1 , ħ 2 , ħ 3 C ( [ 0 , S ] , R + ) .
Definition 4.
The fractal-fractional CD 4 + -HIV-1-model (2) is generalized Ulam–Hyers–Rassias stable w.r.t. functions ħ j if 0 < a ( Q j , ħ j ) R s.t. T * , U * , V * X satisfying
| FFML D 0 , s ( δ , σ ) T * ( s ) Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) | < ħ 1 ( s ) , | FFML D 0 , s ( δ , σ ) U * ( s ) Q 2 s , T * ( s ) , U * ( s ) , V * ( s ) | < ħ 2 ( s ) , | FFML D 0 , s ( δ , σ ) V * ( s ) Q 3 s , T * ( s ) , U * ( s ) , V * ( s ) | < ħ 3 ( s ) ,
T , U , V X as a solution of the fractal-fractional CD 4 + -HIV-1-model (2) s.t.
| T * ( s ) T ( s ) | a ( Q 1 , ħ 1 ) ħ 1 ( s ) , | U * ( s ) U ( s ) | a ( Q 2 , ħ 2 ) ħ 2 ( s ) , | V * ( s ) V ( s ) | a ( Q 3 , ħ 3 ) ħ 3 ( s ) .
If we take ħ j ( s ) = 1 , in this case Definition 3 yields the Ulam-Hyers criterion for the stability of solutions.
Remark 2.
T * , U * , V * X is a solution for (21) iff z 1 , z 2 , z 3 C ( [ 0 , S ] , R ) (each of them depend on T * , U * , V * , respectively) s.t. s J ,
(i) 
| z j ( s ) | < r j ħ j ( s ) ;
(ii) 
We have
FFML D 0 , s ( δ , σ ) T * ( s ) = Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) + z 1 ( s ) , FFML D 0 , s ( δ , σ ) U * ( s ) = Q 2 s , T * ( s ) , U * ( s ) , V * ( s ) + z 2 ( s ) , FFML D 0 , s ( δ , σ ) V * ( s ) = Q 3 s , T * ( s ) , U * ( s ) , V * ( s ) + z 3 ( s ) .
The following lemmas are useful for our main theorems.
Lemma 2.
For each r 1 , r 2 , r 3 > 0 , suppose that ( T * , U * , V * ) X is considered as a solution of (20). Then the functions T * , U * , V * M fulfill the inequalities
| T * ( s ) ( T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 × Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ) | ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r 1 ,
and
| U * ( s ) ( U 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 2 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 × Q 2 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ) | ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r 2 ,
and
| V * ( s ) ( V 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 3 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 × Q 3 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ) | ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r 3 .
Proof. 
Let r 1 > 0 be arbitrary. Since T * M satisfies
| FFML D 0 , s ( δ , σ ) T * ( s ) Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) | < r 1 ,
so, by Remark 1, we are allowed to take a function z 1 ( s ) s.t.
FFML D 0 , s ( δ , σ ) T * ( s ) = Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) + z 1 ( s ) ,
and | z 1 ( s ) | r 1 . Clearly,
T * ( s ) = T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + z 1 ( s ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) + z 1 ( w ) d w .
In this case, we estimate
| T * ( s ) ( T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ) | ( 1 δ ) σ s σ 1 AB ( δ ) | z 1 ( s ) | + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 | z 1 ( w ) | d w ( 1 δ ) σ S σ 1 AB ( δ ) r 1 + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r 1 = ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r 1 .
This states that the inequality (22) is fulfilled. Similarly, we can obtain the inequalities (23) and (24). □
To prove the next result, we regard the following:
(C2)
increasing mappings ħ j C ( [ 0 , S ] , R + ) , ( j { 1 , 2 , 3 } ) and Δ ħ j > 0 provided that
FFML I 0 , s ( δ , σ ) ħ j ( s ) < Δ ħ j ħ j ( s ) , ( j { 1 , 2 , 3 } ) , s J .
Lemma 3.
Let ( C 2 ) to be held. For each r 1 , r 2 , r 3 > 0 , suppose that ( T * , U * , V * ) X is considered as a solution of (21). Then the functions T * , U * , V * M fulfill the inequalities
| T * ( s ) T 0 + FFML I 0 , s ( δ , σ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) | r 1 Δ ħ 1 ħ 1 ( s ) ,
| U * ( s ) U 0 + FFML I 0 , s ( δ , σ ) Q 2 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) | r 2 Δ ħ 2 ħ 2 ( s ) ,
| V * ( s ) V 0 + FFML I 0 , s ( δ , σ ) Q 3 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) | r 3 Δ ħ 3 ħ 3 ( s ) .
Proof. 
Let r 1 > 0 . Since T * M satisfies
| FFML D 0 , s ( δ , σ ) T * ( s ) Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) | < r 1 ħ 1 ( s ) ,
thus, from Remark 2, we are allowed to take a function z 1 ( s ) s.t.
FFML D 0 , s ( δ , σ ) T * ( s ) = Q 1 s , T * ( s ) , U * ( s ) , V * ( s ) + z 1 ( s ) ,
and | z 1 ( s ) | r 1 ħ 1 ( s ) . Evidently,
T * ( s ) = T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + z 1 ( s ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) + z 1 ( w ) d w .
Then, we estimate
| T * ( s ) T 0 + FFML I 0 , s ( δ , σ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) | = | T * ( s ) ( T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ) | ( 1 δ ) σ s σ 1 AB ( δ ) | z 1 ( s ) | + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 | z 1 ( w ) | d w = FFML I 0 , s ( δ , σ ) | z 1 ( s ) | FFML I 0 , s ( δ , σ ) r 1 ħ 1 ( s ) r 1 Δ ħ 1 ħ 1 ( s ) .
Similarly, we can obtain the remaining inequalities. □
The Ulam–Hyers stability is checked about the fractal-fractional CD 4 + -HIV-1-model (2).
Theorem 4.
Let (C1) be fulfilled. Then the fractal-fractional CD 4 + -HIV-1-model (2) is Ulam–Hyers stable on J : = [ 0 , S ] and also is generalized Ulam–Hyers stable s.t.
( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α j < 1 , j { 1 , 2 , 3 } ,
in which α 1 , α 2 , α 3 > 0 are given by (18).
Proof. 
Let r 1 > 0 and T * M be an arbitrary solution of (20). Also, from Theorem 3, we assume T M as a unique solution of the fractal-fractional CD 4 + -HIV-1-model (2). Then T ( s ) is defined as
T ( s ) = T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) d w .
Therefore, by Lemma 2 and with the help of the triangle inequality, we estimate
| T * ( s ) T ( s ) | | T * ( s ) T 0 ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) d w | | T * ( s ) ( T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ) | + ( 1 δ ) σ s σ 1 AB ( δ ) | Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) | + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 | Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) | d w ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r 1 + ( 1 δ ) σ S σ 1 AB ( δ ) α 1 T * T + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 T * T ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r 1 + ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 T T * .
Hence,
T * T ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r 1 1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 .
We set a Q 1 = ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) 1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 . In this case, T * T a Q 1 r 1 .
Similarly, we obtain
U * U a Q 2 r 2 ,
V * V a Q 3 r 3 ,
where
a Q j = ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) 1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α j , ( j { 2 , 3 } ) .
Hence, the Ulam–Hyers stability of the fractal-fractional CD 4 + -HIV-1-model (2) is fulfilled. On the other hand, if we take
a Q j ( r j ) = ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) r j 1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α j , ( j { 1 , 2 , 3 } ) ,
then a Q j ( 0 ) = 0 , and the generalized Ulam–Hyers stability is simply proved. □
The Ulam–Hyers-Rassias stability is checked about the fractal-fractional CD 4 + -HIV-1-model (2) in the next theorem.
Theorem 5.
The hypotheses (C1) and (C2) are considered to be held. Then the given fractal-fractional CD 4 + -HIV-1-model (2) is the Ulam–Hyers–Rassias and generalized Ulam–Hyers–Rassias stable.
Proof. 
Let r 1 > 0 , and T * M satisfying (21). By Theorem 3, let T M be the unique solution of the given fractal-fractional CD 4 + -HIV-1-model (2). Then T ( s ) becomes
T ( s ) = T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) d w .
Therefore, by Lemma 3 and with the help of the triangle inequality, we estimate
| T * ( s ) T ( s ) | | T * ( s ) T 0 ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) d w | | T * ( s ) ( T 0 + ( 1 δ ) σ s σ 1 AB ( δ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) d w ) | + ( 1 δ ) σ s σ 1 AB ( δ ) | Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) | + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 | Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) | d w | T * ( s ) T 0 + FFML I 0 , s ( δ , σ ) Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) | + ( 1 δ ) σ s σ 1 AB ( δ ) | Q 1 ( s , T * ( s ) , U * ( s ) , V * ( s ) ) Q 1 ( s , T ( s ) , U ( s ) , V ( s ) ) | + δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 | Q 1 ( w , T * ( w ) , U * ( w ) , V * ( w ) ) Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) | d w r 1 Δ ħ 1 ħ 1 ( s ) + ( 1 δ ) σ S σ 1 AB ( δ ) α 1 T * T + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 T * T r 1 Δ ħ 1 ħ 1 ( s ) + ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 T T * .
Accordingly, it gives
T * T r 1 Δ ħ 1 ħ 1 ( s ) 1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 .
Set
a ( Q 1 , ħ 1 ) = Δ ħ 1 1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α 1 .
Then T * T r 1 a ( Q 1 , ħ 1 ) ħ 1 ( s ) . Similarly,
U * U r 2 a ( Q 2 , ħ 2 ) ħ 2 ( s ) ,
V * V r 3 a ( Q 3 , ħ 3 ) ħ 3 ( s ) ,
where
a ( Q j , ħ j ) = Δ ħ j 1 ( 1 δ ) σ S σ 1 AB ( δ ) + δ σ S δ + σ 1 Γ ( σ ) AB ( δ ) Γ ( δ + σ ) α j , ( j { 2 , 3 } ) .
In consequence, the fractal-fractional CD 4 + -HIV-1-model (2) is stable in the context of the Ulam–Hyers–Rassias criterion. Along with this, by defining r j = 1 , ( j { 1 , 2 , 3 } ) , the mentioned fractal-fractional CD 4 + -HIV-1-model (2) is generalized Ulam–Hyers–Rassias stable. □

6. Numerical Scheme via Adams-Bashforth Method

In the present section of the paper, we aim to derive a numerical algorithm for the fractal-fractional CD 4 + -HIV-1-model (2). To do this, we apply a technique based on two-step Lagrange polynomials entitled fractional Adams-Bashforth method. We re-define fractal-fractional integral equations (12) at s k + 1 . In fact, we discretize these integral equations (12) for s = s k + 1 as follows
T ( s k + 1 ) = T 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 1 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s k + 1 w σ 1 ( s k + 1 w ) δ 1 Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) d w , U ( s k + 1 ) = U 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 2 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s k + 1 w σ 1 ( s k + 1 w ) δ 1 Q 2 ( w , T ( w ) , U ( w ) , V ( w ) ) d w , V ( s k + 1 ) = V 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 3 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ AB ( δ ) Γ ( δ ) 0 s k + 1 w σ 1 ( s k + 1 w ) δ 1 Q 3 ( w , T ( w ) , U ( w ) , V ( w ) ) d w .
The approximation of above integrals are formulated by
T ( s k + 1 ) = T 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 1 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ AB ( δ ) Γ ( δ ) = 1 k s s + 1 w σ 1 ( s k + 1 w ) δ 1 Q 1 ( w , T ( w ) , U ( w ) , V ( w ) ) d w , U ( s k + 1 ) = U 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 2 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ AB ( δ ) Γ ( δ ) = 1 k s s + 1 w σ 1 ( s k + 1 w ) δ 1 Q 2 ( w , T ( w ) , U ( w ) , V ( w ) ) d w , V ( s k + 1 ) = V 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 3 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ AB ( δ ) Γ ( δ ) = 1 k s s + 1 w σ 1 ( s k + 1 w ) δ 1 Q 3 ( w , T ( w ) , U ( w ) , V ( w ) ) d w .
Next, we approximate three functions w σ 1 Q j ( w , T ( w ) , U ( w ) , V ( w ) ) , j = 1 , 2 , 3 , on the interval [ s , s + 1 ] by applying two-step Lagrange interpolation polynomials under the step size h = s s 1 . By straighforward computations, we obtain algorithms which yield the numerical solutions to the fractal-fractional CD 4 + -HIV-1-model (2) as
T k + 1 = T 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 1 ( s k , T k , U k , V k ) + σ h δ AB ( δ ) Γ ( δ + 2 ) × = 1 k s σ 1 Q 1 s , T , U , V Y ^ 1 ( k , ) s 1 σ 1 Q 1 s 1 , T 1 , U 1 , V 1 Y ^ 2 ( k , ) ,
U k + 1 = U 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 2 ( s k , T k , U k , V k ) + σ h δ AB ( δ ) Γ ( δ + 2 ) × = 1 k s σ 1 Q 2 s , T , U , V Y ^ 1 ( k , ) s 1 σ 1 Q 2 s 1 , T 1 , U 1 , V 1 Y ^ 2 ( k , ) ,
V k + 1 = V 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 3 ( s k , T k , U k , V k ) + σ h δ AB ( δ ) Γ ( δ + 2 ) × = 1 k s σ 1 Q 3 s , T , U , V Y ^ 1 ( k , ) s 1 σ 1 Q 3 s 1 , T 1 , U 1 , V 1 Y ^ 2 ( k , ) ,
where
Y ^ 1 ( k , ) = ( k + 1 ) δ ( k + 2 + δ ) ( k ) δ ( k + 2 + 2 δ ) ,
and
Y ^ 2 ( k , ) = ( k + 1 ) δ + 1 ( k ) δ ( k + 1 + δ ) .

7. Numerical Scheme via Newton Polynomials Method

In this section, we produce a new numerical scheme for solutions of our fractal-fractional CD 4 + -HIV-1-model (2) which was introduced by Atangana and Araz in the book [44] in 2021. To do this, we again use the compact form of IVP (7) with the conditions (8) and (9). In this case, we have
K ( s ) K ( 0 ) = δ σ AB ( δ ) Γ ( δ ) 0 s w σ 1 ( s w ) δ 1 Q ( w , K ( w ) ) d w + ( 1 δ ) σ s σ 1 AB ( δ ) Q ( s , K ( s ) ) .
Set Q * ( s , K ( s ) ) = σ s σ 1 Q ( s , K ( s ) ) . Then
K ( s ) K ( 0 ) = δ AB ( δ ) Γ ( δ ) 0 s ( s w ) δ 1 Q * ( w , K ( w ) ) d w + ( 1 δ ) AB ( δ ) Q * ( s , K ( s ) ) .
By discretizing the above equation at s = s k + 1 = ( k + 1 ) h , we get
K ( s k + 1 ) K ( 0 ) = δ AB ( δ ) Γ ( δ ) 0 s k + 1 ( s k + 1 w ) δ 1 Q * ( w , K ( w ) ) d w + ( 1 δ ) AB ( δ ) Q * ( s k , K ( s k ) ) .
If we approximate the above integral, then it becomes
K ( s k + 1 ) = K 0 + ( 1 δ ) AB ( δ ) Q * ( s k , K ( s k ) ) + δ AB ( δ ) Γ ( δ ) = 2 k s s + 1 ( s k + 1 w ) δ 1 Q * ( w , K ( w ) ) d w .
In this step, the function Q * ( s , K ( s ) ) is approximated by the Newton polynomial as
P k * ( w ) = Q * ( s k 2 , K ( s k 2 ) ) + Q * ( s k 1 , K ( s k 1 ) ) Q * ( s k 2 , K ( s k 2 ) ) h ( w s k 2 ) + Q * ( s k , K ( s k ) ) 2 Q * ( s k 1 , K ( s k 1 ) ) + Q * ( s k 2 , K ( s k 2 ) ) 2 h 2 ( w s k 2 ) ( w s k 1 ) .
Substitute (33) into (32):
K k + 1 = K 0 + ( 1 δ ) AB ( δ ) Q * ( s k , K ( s k ) ) + δ AB ( δ ) Γ ( δ ) = 2 k s s + 1 [ Q * ( s 2 , K 2 ) + Q * ( s 1 , K 1 ) Q * ( s 2 , K 2 ) h ( w s 2 ) + Q * ( s , K ) 2 Q * ( s 1 , K 1 ) + Q * ( s 2 , K 2 ) 2 h 2 ( w s 2 ) ( w s 1 ) ] × ( s k + 1 w ) δ 1 d w .
We simplify the above relations, and we get
K k + 1 = K 0 + ( 1 δ ) AB ( δ ) Q * ( s k , K ( s k ) ) + δ AB ( δ ) Γ ( δ ) = 2 k [ s s + 1 Q * ( s 2 , K 2 ) ( s k + 1 w ) δ 1 d w + s s + 1 Q * ( s 1 , K 1 ) Q * ( s 2 , K 2 ) h ( w s 2 ) ( s k + 1 w ) δ 1 d w + s s + 1 Q * ( s , K ) 2 Q * ( s 1 , K 1 ) + Q * ( s 2 , K 2 ) 2 h 2 ( w s 2 ) ( w s 1 ) × ( s k + 1 w ) δ 1 d w ] .
In consequence,
K k + 1 = K 0 + ( 1 δ ) AB ( δ ) Q * ( s k , K ( s k ) ) + δ AB ( δ ) Γ ( δ ) = 2 k Q * ( s 2 , K 2 ) s s + 1 ( s k + 1 w ) δ 1 d w + δ AB ( δ ) Γ ( δ ) = 2 k Q * ( s 1 , K 1 ) Q * ( s 2 , K 2 ) h s s + 1 ( w s 2 ) ( s k + 1 w ) δ 1 d w + δ AB ( δ ) Γ ( δ ) = 2 k Q * ( s , K ) 2 Q * ( s 1 , K 1 ) + Q * ( s 2 , K 2 ) 2 h 2 × s s + 1 ( w s 2 ) ( w s 1 ) × ( s k + 1 w ) δ 1 d w ] .
On the other hand, we compute above three integrals separately, and we get
s s + 1 ( s k + 1 w ) δ 1 d w = h δ δ ( k + 1 ) δ ( k ) δ ,
and
s s + 1 ( w s 2 ) ( s k + 1 w ) δ 1 d w = h δ + 1 δ ( δ + 1 ) [ ( k + 1 ) δ ( k + 3 + 2 δ ) ( k + 1 ) δ ( k + 3 + 3 δ ) ] ,
and
s s + 1 ( w s 2 ) ( w s 1 ) ( s k + 1 w ) δ 1 d w = h δ + 2 δ ( δ + 1 ) ( δ + 2 ) ( ( k + 1 ) δ [ 2 ( k ) 2 + ( 3 δ + 10 ) ( k ) + 2 δ 2 + 9 δ + 12 ] ( k ) δ [ 2 ( k ) 2 + ( 5 δ + 10 ) ( k ) + 6 δ 2 + 18 δ + 12 ] ) .
By putting (35)–(37) in (34), we obtain
K k + 1 = K 0 + ( 1 δ ) AB ( δ ) Q * ( s k , K ( s k ) ) + δ h δ AB ( δ ) Γ ( δ + 1 ) = 2 k Q * ( s 2 , K 2 ) ( k + 1 ) δ ( k ) δ + δ h δ AB ( δ ) Γ ( δ + 2 ) = 2 k Q * ( s 1 , K 1 ) Q * ( s 2 , K 2 ) × ( k + 1 ) δ ( k + 3 + 2 δ ) ( k + 1 ) δ ( k + 3 + 3 δ ) + δ h δ 2 AB ( δ ) Γ ( δ + 3 ) = 2 k Q * ( s , K ) 2 Q * ( s 1 , K 1 ) + Q * ( s 2 , K 2 ) × [ ( k + 1 ) δ 2 ( k ) 2 + ( 3 δ + 10 ) ( k ) + 2 δ 2 + 9 δ + 12 ( k ) δ [ 2 ( k ) 2 + ( 5 δ + 10 ) ( k ) + 6 δ 2 + 18 δ + 12 ] ] .
Lastly, we replace Q * ( s , K ( s ) ) = σ s σ 1 Q ( s , K ( s ) ) in (38), and we get
K k + 1 = K 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q ( s k , K ( s k ) ) + δ σ h δ AB ( δ ) Γ ( δ + 1 ) = 2 k s 2 σ 1 Q ( s 2 , K 2 ) Ψ ^ 1 ( k , , δ ) + δ σ h δ AB ( δ ) Γ ( δ + 2 ) = 2 k s 1 σ 1 Q ( s 1 , K 1 ) s 2 σ 1 Q ( s 2 , K 2 ) Ψ ^ 2 ( k , , δ ) + δ σ h δ 2 AB ( δ ) Γ ( δ + 3 ) = 2 k s σ 1 Q ( s , K ) 2 s 1 σ 1 Q ( s 1 , K 1 ) + s 2 σ 1 Q ( s 2 , K 2 ) Ψ ^ 3 ( k , , δ ) ,
where
Ψ ^ 1 ( k , , δ ) = ( k + 1 ) δ ( k ) δ , Ψ ^ 2 ( k , , δ ) = ( k + 1 ) δ ( k + 3 + 2 δ ) ( k + 1 ) δ ( k + 3 + 3 δ ) , Ψ ^ 3 ( k , , δ ) = ( k + 1 ) δ 2 ( k ) 2 + ( 3 δ + 10 ) ( k ) + 2 δ 2 + 9 δ + 12 ( k ) δ 2 ( k ) 2 + ( 5 δ + 10 ) ( k ) + 6 δ 2 + 18 δ + 12 .
Based on the numerical scheme obtained in (39), numerical solutions of the fractal-fractional CD 4 + -HIV-1-model (2) are given by
T k + 1 = T 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 1 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ h δ AB ( δ ) Γ ( δ + 1 ) = 2 k s 2 σ 1 Q 1 ( s 2 , T 2 , U 2 , V 2 ) Ψ ^ 1 ( k , , δ ) + δ σ h δ AB ( δ ) Γ ( δ + 2 ) = 2 k s 1 σ 1 Q 1 ( s 1 , T 1 , U 1 , V 1 ) s 2 σ 1 Q 1 ( s 2 , T 2 , U 2 , V 2 ) Ψ ^ 2 ( k , , δ ) + δ σ h δ 2 AB ( δ ) Γ ( δ + 3 ) = 2 k [ s σ 1 Q 1 ( s , T , U , V ) 2 s 1 σ 1 Q 1 ( s 1 , T 1 , U 1 , V 1 ) + s 2 σ 1 Q 1 ( s 2 , T 2 , U 2 , V 2 ) ] Ψ ^ 3 ( k , , δ ) ,
and
U k + 1 = U 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 2 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ h δ AB ( δ ) Γ ( δ + 1 ) = 2 k s 2 σ 1 Q 2 ( s 2 , T 2 , U 2 , V 2 ) Ψ ^ 1 ( k , , δ ) + δ σ h δ AB ( δ ) Γ ( δ + 2 ) = 2 k s 1 σ 1 Q 2 ( s 1 , T 1 , U 1 , V 1 ) s 2 σ 1 Q 2 ( s 2 , T 2 , U 2 , V 2 ) Ψ ^ 2 ( k , , δ ) + δ σ h δ 2 AB ( δ ) Γ ( δ + 3 ) = 2 k [ s σ 1 Q 2 ( s , T , U , V ) 2 s 1 σ 1 Q 2 ( s 1 , T 1 , U 1 , V 1 ) + s 2 σ 1 Q 2 ( s 2 , T 2 , U 2 , V 2 ) ] Ψ ^ 3 ( k , , δ ) ,
and
V k + 1 = V 0 + ( 1 δ ) σ s k σ 1 AB ( δ ) Q 3 ( s k , T ( s k ) , U ( s k ) , V ( s k ) ) + δ σ h δ AB ( δ ) Γ ( δ + 1 ) = 2 k s 2 σ 1 Q 3 ( s 2 , T 2 , U 2 , V 2 ) Ψ ^ 1 ( k , , δ ) + δ σ h δ AB ( δ ) Γ ( δ + 2 ) = 2 k s 1 σ 1 Q 3 ( s 1 , T 1 , U 1 , V 1 ) s 2 σ 1 Q 3 ( s 2 , T 2 , U 2 , V 2 ) Ψ ^ 2 ( k , , δ ) + δ σ h δ 2 AB ( δ ) Γ ( δ + 3 ) = 2 k [ s σ 1 Q 3 ( s , T , U , V ) 2 s 1 σ 1 Q 3 ( s 1 , T 1 , U 1 , V 1 ) + s 2 σ 1 Q 3 ( s 2 , T 2 , U 2 , V 2 ) ] Ψ ^ 3 ( k , , δ ) ,
where the constants Ψ ^ j ( k , , δ ) are introduced in (40) for j = 1 , 2 , 3 .

8. Simulations and Discussions

In this section, we simulate and discuss dynamical behaviors of the fractal-fractional CD 4 + -HIV-1-model (2) via real data for parameters and initial values assumed by [45]. Based on the given data in this source, we let θ = 10 , q = 0.000024 , ρ = 0.01 , ζ = 0.2 , κ = 0.16 , ϑ = 3.4 , and N = 1000 . Moreover, the initial values for three state functions are T ( 0 ) = T 0 = 1000 , U ( 0 ) = U 0 = 0 , V ( 0 ) = V 0 = 0.001 . It is assumed that all parameters are in days.
In the provided figures, we show the behaviors of three state functions T , U , V under the effect of different values for the fractal-fractional orders δ = σ = 1.00 , 0.98 , 0.96 , 0.94 , 0.92 , 0.90 simultaneously. We also compare the results in the graphs with respect to two numerical algorithms given by (29)–(31) and (41)–(43) as shown in Table 1, Table 2 and Table 3.
In these graphs, the quantity of the time s interprets the number of days. Note that Figure 1 and Figure 2 demonstrate the dynamics of the CD 4 + T-cells via the Adams-Bashforth method (Section 6). Figure 1a indicates that when the fractal dimensions and fractional orders decrease, the number of susceptible CD 4 + T-cells steadily increases from the 42nd day to the 64th day, and steadily decreases from the 65th day to the 102nd day. Further it increases again from the 103rd day and finally converges to the integer-order at the 250th day to the end of the simulation period.
Figure 1b,c show that when the fractal dimensions and fractional orders move from the integer order, the peak of the amount of infectious CD 4 + T-cells and the free particles of the infection of the HIV-1 virus in the blood decrease respectively and also give slight differences in their asymptotic stabilities.
In Figure 2a,c, we show a 0.01 -variation in fractal-fractional order. In this case, the numerical trajectories show that a slight change in the fractal-fractional order produces a slight changes in the asymptotic behavior of the HIV-1 virus on CD 4 + T-cells.
Figure 3 and Figure 4 demonstrate the dynamics of the CD 4 + T-cells via the Newton polynomials method (Section 7). Figure 3a also indicates that when the fractal dimensions and fractional orders decrease, the number of susceptible CD 4 + T-cells steadily increases from the 42nd day to the 64th day, and steadily decreases from the 65th day to the 102nd day. Also it increases again from the 103rd day and finally converges to the integer-order at the 250th day to the end of the simulation period, but it shows a disorganized behavior from the 11th day to the 24th day, which indicates the early stages of HIV in the CD 4 + T-cells.
Figure 3b,c show that when the fractal dimensions and fractional orders move from the integer order, the peak of the amount of infectious CD 4 + T-cells and the free particles of the infection of the HIV-1 virus in the blood decrease, respectively and also give slight differences in their asymptotic stabilities, but again, they show a disorganized behavior from the 11th day to the 24th day.
In Figure 4a,c, we show a 0.01 -variation in fractal-fractional orders, and the numerical trajectories show that a slight change in the fractal-fractional orders produce slight changes in the asymptotic behavior of the HIV-1 virus on CD 4 + T-cells.
In Figure 5, based on the Adams-Bashforth method, we see the numerical trajectories of three state functions by varying the average number of infected particles N for the values N = 500 , 600 , 700 , 800 , 900 , 1000 under the fractal-fractional order δ = σ = 0.95 . Thus, in Figure 5a, it shows that the susceptible CD 4 + T-cells increase as the average number of infected particles by an existing infected cell reduces by 10 % , and that of the peak of the amount of infectious CD 4 + T-cells and the free particles of the infection of the HIV-1 virus in the blood decrease, respectively in Figure 5b,c.
In Figure 6, based on the Adams-Bashforth method, we see the numerical trajectories by varying the supply rate for new T-cells for θ = 10 , 15 , 20 , 25 , 30 , 35 under the fractal order and fractional order δ = σ = 0.95 . Thus, in Figure 6a, it shows that the susceptible CD 4 + T-cells increase for a period of time and gradually decrease as the number of the supply rate θ increases before converging at the 150th day.
In Figure 6b,c, we noticed that the amount of infectious CD 4 + T-cells and the free particles of the infection of the HIV-1 virus in the blood increase, respectively. Figure 7a–c display the numerical trajectories when we compare the numerical solutions of the Adams-Bashforth method with the Newton polynomial method under the fractal dimension and fractional order δ = σ = 0.90 .
In Table 1, Table 2 and Table 3, we see some results of the two numerical schemes (the Adams-Bashforth method and Newton polynomials method) for all of three state functions under the fractal dimension and fractional order δ = σ = 0.95 with step size h = 0.1 .

9. Conclusions

In this paper, we designed a fractal-fractional CD 4 + -HIV-1-model and analyzed the dynamics of CD 4 + T-cells under the infection of HIV-1 virus. We considered three compartments for this model by defining three state functions T , U , and V for the amount of susceptible CD 4 + T-cells, amount of infectious CD 4 + T-cells, and the free particles of the infection of the HIV virus in the blood. We derived three fractal-fractional integral equations and proved that their kernels are Lipschitz. In this direction, we could prove the existence and uniqueness criteria for solutions of the fractal-fractional CD 4 + -HIV-1-model. In the sequel, we investigated four stability results with the help of two auxiliary inequality. We extracted two algorithms via the Adams-Bashforth method and also via the Newton polynomials and simulated our real data in relation to the given fractal-fractional CD 4 + -HIV-1-model. The numerical and graphical results showed that these two numerical algorithms give the same outcomes and differences are small. Also, we investigated the effect of frcatal dimensions and fractional orders on these simulations. Also, the effect of different values for the average number of infected particles and the supply rate of new T-cells were simulated in some graphs under the Adams-Bashforth method. This stude showed that we can predict the next behavior of the fractal-fractional CD 4 + -HIV-1-model via the two mentioned numerical methods and their results are more accurate and identical. This shows the power of simulation of the frcatal-fractional models in comparison to the fractional models. In the next researches, we can develope our numerical methods on different fractal-fractional models of diseases.

Author Contributions

Conceptualization, H.N. and S.E.; Formal analysis, H.N., S.E., N.P., J.K.K.A., S.R. and T.S.; Funding acquisition, N.P. and T.S.; Methodology, H.N., S.E., N.P., J.K.K.A., S.R. and T.S.; Software, S.E. and J.K.K.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Mongkut’s University of Technology North Bangkok. Contract no. KMUTNB-62-KNOW-27.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.

Acknowledgments

The first author would like to thank Shiraz University. Also, the second and fifth authors would like to thank Azarbaijan Shahid Madani University. Also, we would like to thank dear reviewers for their constructive and helpful comments and remarks to improve the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, L.; Li, M.Y. Mathematical analysis of the global dynamics of a model for HIV infection of CD4+ T cells. Math. Biosci. 2006, 200, 44–57. [Google Scholar] [CrossRef] [PubMed]
  2. UNAIDS. Report on the Global AIDS Epidemic. 2016. Available online: http://www.unaids.org (accessed on 18 February 2022).
  3. Kumar, P.; Erturk, V.S. The analysis of a time delay fractional COVID-19 model via Caputo type fractional derivative. Math. Methods Appl. Sci. 2022. [Google Scholar] [CrossRef] [PubMed]
  4. Rezapour, S.; Etemad, S.; Mohammadi, H. A mathematical analysis of a system of Caputo-Fabrizio fractional differential equationsfor the anthrax disease model in animals. Adv. Differ. Equ. 2020, 2020, 481. [Google Scholar] [CrossRef]
  5. Alshehri, H.M.; Khan, A. A fractional order Hepatitis C mathematical model with Mittag-Leffler kernel. J. Funct. Spaces 2021, 2021, 2524027. [Google Scholar] [CrossRef]
  6. Deressa, C.T.; Etemad, S.; Rezapour, S. On a new four-dimensional model of memristor-based chaotic circuit in the context of nonsingular Atangana-Baleanu-Caputo operators. Adv. Differ. Equ. 2021, 2021, 444. [Google Scholar] [CrossRef]
  7. Deressa, C.T.; Etemad, S.; Kaabar, M.K.A.; Rezapour, S. Qualitative analysis of a hyperchaotic Lorenz-Stenflo mathematical model via the Caputo fractional operator. J. Funct. Spaces 2022, 2022, 4975104. [Google Scholar] [CrossRef]
  8. Kumar, P.; Erturk, V.S. Environmental persistence influences infection dynamics for a butterfly pathogen via new generalised Caputo type fractional derivative. Chaos Solitons Fractals 2021, 144, 110672. [Google Scholar] [CrossRef]
  9. Devi, A.; Kumar, A.; Abdeljawad, T.; Khan, A. Stability analysis of solutions and existence theory of fractional Lagevin equation. Alex. Eng. J. 2021, 60, 3641–3647. [Google Scholar] [CrossRef]
  10. Mohammadi, H.; Kumar, S.; Rezapour, S.; Etemad, S. A theoretical study of the Caputo-Fabrizio fractional modeling for hearing loss due to Mumps virus with optimal control. Chaos Solitons Fractals 2021, 144, 110668. [Google Scholar] [CrossRef]
  11. Begum, R.; Tunc, O.; Khan, H.; Gulzar, H.; Khan, A. A fractional order Zika virus model with Mittag-Leffler kernel. Chaos Solitons Fractals 2021, 146, 110898. [Google Scholar] [CrossRef]
  12. Kumar, P.; Erturk, V.S.; Almusawa, H. Mathematical structure of mosaic disease using microbial biostimulants via Caputo and Atangana-Baleanu derivatives. Results Phys. 2021, 24, 104186. [Google Scholar] [CrossRef]
  13. Zarin, R.; Khaliq, H.; Khan, A.; Khan, D.; Akgul, A.; Humphries, U.W. Deterministic and fractional modeling of a computer virus propagation. Results Phys. 2022, 33, 105130. [Google Scholar] [CrossRef]
  14. Baleanu, D.; Etemad, S.; Rezapour, S. A hybrid Caputo fractional modeling for thermostat with hybrid boundary value conditions. Bound. Value Probl. 2020, 2020, 64. [Google Scholar] [CrossRef] [Green Version]
  15. Thaiprayoon, C.; Sudsutad, W.; Alzabut, J.; Etemad, S.; Rezapour, S. On the qualitative analysis of the fractional boundary valueproblem describing thermostat control model via ψ-Hilfer fractional operator. Adv. Differ. Equ. 2021, 2021, 201. [Google Scholar] [CrossRef]
  16. Alzabut, J.; Selvam, G.M.; El-Nabulsi, R.A.; Vignesh, D.; Samei, M.E. Asymptotic stability of nonlinear discrete fractional pantograph equations with non-local initial conditions. Symmetry 2021, 13, 473. [Google Scholar] [CrossRef]
  17. Wongcharoen, A.; Ntouyas, S.K.; Tariboon, J. Nonlocal boundary value problems for Hilfer type pantograph fractional differentialequations and inclusions. Adv. Differ. Equ. 2020, 2020, 279. [Google Scholar] [CrossRef]
  18. Kumar, P.; Erturk, V.S.; Yusuf, A.; Nisar, K.S.; Abdelwahab, S.F. A study on canine distemper virus (CDV) and rabies epidemics in the red fox population via fractional derivatives. Results Phys. 2021, 25, 104281. [Google Scholar] [CrossRef]
  19. Asamoah, J.K.K.; Okyere, E.; Yankson, E.; Opoku, A.A.; Adom-Konadu, A.; Acheampong, E.; Arthur, Y.D. Non-fractional and fractional mathematical analysis and simulations for Q fever. Chaos Solitons Fractals 2022, 156, 111821. [Google Scholar] [CrossRef]
  20. Khan, H.; Tunc, C.; Chen, W.; Khan, A. Existence theorems and Hyers-Ulam stability for a class of hybrid fractional differentialequations with p-Laplacial operator. J. Appl. Anal. Comput. 2018, 8, 1211–1226. [Google Scholar]
  21. Omame, A.; Nwajeri, U.K.; Abbas, M.; Onyenegecha, C.P. A fractional order control model for Diabetes and COVID-19 co-dynamics with Mittag-Leffler function. Alex. Eng. J. 2022, 61, 7619–7635. [Google Scholar] [CrossRef]
  22. Ali, W.; Turab, A.; Nieto, J.J. On the novel existence results of solutions for a class of fractional boundary value problems on the cyclohexane graph. J. Inequalities Appl. 2022, 2022, 5. [Google Scholar] [CrossRef]
  23. Baleanu, D.; Etemad, S.; Mohammadi, H.; Rezapour, S. A novel modeling of boundary value problems on the glucose graph. Commun. Nonlinear Sci. Numer. Simul. 2021, 100, 105844. [Google Scholar] [CrossRef]
  24. Rezapour, S.; Tellab, B.; Deressa, C.T.; Etemad, S.; Nonlaopon, K. H-U-type stability and numerical solutions for a nonlinear model of the coupled systems of Navier BVPs via the generalized differential transform method. Fractal Fract. 2021, 5, 166. [Google Scholar] [CrossRef]
  25. Culshaw, R.V.; Ruan, S. A delay-differential equation model of HIV infection of CD4+ T-cells. Math. Biosci. 2000, 165, 27–39. [Google Scholar] [CrossRef]
  26. Cai, L.; Li, X.; Ghosh, M.; Guo, B. Stability analysis of an HIV/AIDS epidemic model with treatment. J. Comput. Appl. Math. 2009, 229, 313–323. [Google Scholar] [CrossRef] [Green Version]
  27. Mohyud-Din, S.T.; Nazir, A.; Almohsin, B.; Ahmed, N.; Khan, U.; Waheed, A.; Hussain, T. On mathematical model of HIV CD4+ T-cells. Alex. Eng. J. 2021, 60, 995–1000. [Google Scholar] [CrossRef]
  28. Perlson, A.S. Modeling the interaction of the immune system with HIV. Lec. Notes Biomath. 1989, 83, 350–370. [Google Scholar]
  29. Perlson, A.S.; Kirschner, D.E.; Boer, R.D. Dynamics of HIV infection of CD4+ T cells. Math. Biosci. 1993, 114, 81–125. [Google Scholar] [CrossRef] [Green Version]
  30. Ding, Y.; Ye, H. A fractional-order differential equation model of HIV infection of CD4+ T-cells. Math. Comput. Model. 2009, 50, 386–392. [Google Scholar] [CrossRef]
  31. Arafa, A.A.M.; Rida, S.Z.; Khalil, M. Fractional modeling dynamics of HIV and CD4+ T-cells during primary infection. Nonlinear Biomed. Phys. 2012, 6, 1. [Google Scholar] [CrossRef] [Green Version]
  32. Bulut, H.; Kumar, D.; Singh, J.; Swroop, R.; Baskonus, H.M. Analytic study for a fractional model of HIV infection of CD4+T lymphocyte cells. Math. Nat. Sci. 2018, 2, 33–43. [Google Scholar] [CrossRef]
  33. Lichae, B.H.; Biazar, J.; Ayati, Z. The fractional differential model of HIV-1 infection of CD4+ T-cells with description of the effect of antiviral drug treatment. Comput. Math. Methods Med. 2019, 2019, 4059549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Nazir, G.; Shah, K.; Debbouche, A.; Khan, R.A. Study of HIV mathematical model under nonsingular kernel type derivative of fractional order. Chaos Solitons Fractals 2020, 139, 110095. [Google Scholar] [CrossRef]
  35. Wang, W.; Wang, X.; Feng, Z. Time periodic reaction-diffusion equations for modeling 2-LTR dynamics in HIV-infected patients. Nonlin. Anal. Real World Appl. 2021, 57, 103184. [Google Scholar] [CrossRef]
  36. Atangana, A. Fractal-fractional differentiation and integration: Connecting fractal calculus and fractional calculus to predict complex system. Chaos Solitons Fractals 2017, 102, 396–406. [Google Scholar] [CrossRef]
  37. Shah, K.; Arfan, M.; Mahariq, I.; Ahmadian, A.; Salahshour, S.; Ferrara, M. Fractal-fractional mathematical model addressing the situation of Corona virus in Pakistan. Results Phys. 2020, 19, 103560. [Google Scholar] [CrossRef] [PubMed]
  38. Gomez-Aguilar, J.F.; Cordova-Fraga, T.; Abdeljawad, T.; Khan, A.; Khan, H. Analysis of fractal-fractional Malaria transmission model. Fractals 2020, 28, 2040041. [Google Scholar] [CrossRef]
  39. Ali, Z.; Rabiei, F.; Shah, K.; Khodadadi, T. Qualitative analysis of fractal-fractional order COVID-19 mathematical model with case study of Wuhan. Alex. Eng. J. 2021, 60, 477–489. [Google Scholar] [CrossRef]
  40. Asamoah, J.K.K. Fractal-fractional model and numerical scheme based on Newton polynomial for Q fever disease under Atangana-Baleanu derivative. Res. Phys. 2022, 34, 105189. [Google Scholar] [CrossRef]
  41. Khan, H.; Alam, K.; Gulzar, H.; Etemad, S.; Rezapour, S. A case study of fractal-fractional tuberculosis model in China: Existence and stability theories along with numerical simulations. Math. Comput. Simul. 2022, 198, 455–473. [Google Scholar] [CrossRef]
  42. Ahmad, S.; Ullah, A.; Akgul, A.; De la Sen, M. Study of HIV disease and its association with immune cells under nonsingular and nonlocal fractal-fractional operator. Complexity 2021, 2021, 1904067. [Google Scholar] [CrossRef]
  43. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2003. [Google Scholar]
  44. Atangana, A.; Araz, S.I. New Numerical Scheme with Newton Polynomial: Theory, Methods, and Applications; Academic Press: Cambridge, MA, USA, 2021. [Google Scholar]
  45. Baleanu, D.; Mohammadi, H.; Rezapour, S. Analysis of the model of HIV-1 infection of CD4+ T-cell with a new approach of fractional derivative. Adv. Differ. Equ. 2020, 2020, 71. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Numerical trajectories by varying both fractal order and fractional order δ = σ = 1.00 , 0.98 , 0.96 , 0.94 , 0.92 , 0.90 via the Adams-Bashforth method. (a) Amount of susceptible CD 4 + T-cells. (b) Amount of infectious CD 4 + T-cells. (c) Free particles of the infection of the HIV-1 virus in the blood.
Figure 1. Numerical trajectories by varying both fractal order and fractional order δ = σ = 1.00 , 0.98 , 0.96 , 0.94 , 0.92 , 0.90 via the Adams-Bashforth method. (a) Amount of susceptible CD 4 + T-cells. (b) Amount of infectious CD 4 + T-cells. (c) Free particles of the infection of the HIV-1 virus in the blood.
Mathematics 10 01366 g001
Figure 2. Numerical trajectories by varying both fractal order and fractional orders δ = σ = 1.00 , 0.99 via the Adams-Bashforth method. (a) A 0.01 -variation in fractal-fractional orders of susceptible CD 4 + T-cells. (b) A 0.01 -variation in fractal-fractional orders of infectious CD 4 + T-cells. (c) A 0.01 -variation in fractal-fractional orders of free particles of the infection of the HIV-1 virus in the blood.
Figure 2. Numerical trajectories by varying both fractal order and fractional orders δ = σ = 1.00 , 0.99 via the Adams-Bashforth method. (a) A 0.01 -variation in fractal-fractional orders of susceptible CD 4 + T-cells. (b) A 0.01 -variation in fractal-fractional orders of infectious CD 4 + T-cells. (c) A 0.01 -variation in fractal-fractional orders of free particles of the infection of the HIV-1 virus in the blood.
Mathematics 10 01366 g002
Figure 3. Numerical trajectories by varying both fractal order and fractional order δ = σ = 1.00 , 0.98 , 0.96 , 0.94 , 0.92 , 0.90 via the Newton polynomials method. (a) Amount of susceptible CD 4 + T-cells. (b) Amount of infectious CD 4 + T-cells. (c) Free particles of the infection of the HIV-1 virus in the blood.
Figure 3. Numerical trajectories by varying both fractal order and fractional order δ = σ = 1.00 , 0.98 , 0.96 , 0.94 , 0.92 , 0.90 via the Newton polynomials method. (a) Amount of susceptible CD 4 + T-cells. (b) Amount of infectious CD 4 + T-cells. (c) Free particles of the infection of the HIV-1 virus in the blood.
Mathematics 10 01366 g003
Figure 4. Numerical trajectories by varying both fractal order and fractional order δ = σ = 1.00 , 0.99 via the Newton polynomials method. (a) A 0.01 -variation in fractal-fractional orders of susceptible CD 4 + T-cells. (b) A 0.01 -variation in fractal-fractional orders of infectious CD 4 + T-cells. (c) A 0.01 -variation in fractal-fractional orders of free particles of the infection of the HIV-1 virus in the blood.
Figure 4. Numerical trajectories by varying both fractal order and fractional order δ = σ = 1.00 , 0.99 via the Newton polynomials method. (a) A 0.01 -variation in fractal-fractional orders of susceptible CD 4 + T-cells. (b) A 0.01 -variation in fractal-fractional orders of infectious CD 4 + T-cells. (c) A 0.01 -variation in fractal-fractional orders of free particles of the infection of the HIV-1 virus in the blood.
Mathematics 10 01366 g004
Figure 5. Numerical trajectories by varying the average number of infected particles N under the fractal order and fractional order δ = σ = 0.95 via the Adams-Bashforth method. (a) The effects of N on T . (b) The effects of N on U . (c) The effects of N on V .
Figure 5. Numerical trajectories by varying the average number of infected particles N under the fractal order and fractional order δ = σ = 0.95 via the Adams-Bashforth method. (a) The effects of N on T . (b) The effects of N on U . (c) The effects of N on V .
Mathematics 10 01366 g005
Figure 6. Numerical trajectories by varying the supply rate for new T-cells θ under the fractal order and fractional order δ = σ = 0.95 via the Adams-Bashforth method. (a) The effects of θ on T . (b) The effects of θ on U . (c) The effects of θ on V .
Figure 6. Numerical trajectories by varying the supply rate for new T-cells θ under the fractal order and fractional order δ = σ = 0.95 via the Adams-Bashforth method. (a) The effects of θ on T . (b) The effects of θ on U . (c) The effects of θ on V .
Mathematics 10 01366 g006
Figure 7. Numerical trajectories in the Comparison of the numerical solutions of the Adams-Bashforth method with Newton polynomials method under the fractal dimension and fractional order δ = σ = 0.90 . (a) Comparison of the graphs w.r.t. two different numerical algorithms on T . (b) Comparison of the graphs w.r.t. two different numerical algorithms on U . (c) Comparison of the graphs w.r.t. two different numerical algorithms on V .
Figure 7. Numerical trajectories in the Comparison of the numerical solutions of the Adams-Bashforth method with Newton polynomials method under the fractal dimension and fractional order δ = σ = 0.90 . (a) Comparison of the graphs w.r.t. two different numerical algorithms on T . (b) Comparison of the graphs w.r.t. two different numerical algorithms on U . (c) Comparison of the graphs w.r.t. two different numerical algorithms on V .
Mathematics 10 01366 g007
Table 1. Results of the two numerical schemes for the susceptible CD 4 + T-cells T (s) under the fractal order and fractional order δ = σ = 0.95 with the step size h = 0.1 .
Table 1. Results of the two numerical schemes for the susceptible CD 4 + T-cells T (s) under the fractal order and fractional order δ = σ = 0.95 with the step size h = 0.1 .
Time: (s)0.0001100200300400
Adams-Bashforth method1000 300.8396 315.4653 318.7946 318.8066
Newton polynomials method1000 300.7881 315.4643 318.7951 318.8066
Table 2. Results of the two numerical schemes for the infectious CD 4 + T-cells U (s) under the fractal order and fractional order δ = σ = 0.95 with the step size h = 0.1 .
Table 2. Results of the two numerical schemes for the infectious CD 4 + T-cells U (s) under the fractal order and fractional order δ = σ = 0.95 with the step size h = 0.1 .
Time: (s)0.0001100200300400
Adams-Bashforth method0 37.7693 46.0612 44.2611 43.7066
Newton polynomials method0 37.7778 46.0645 44.2612 43.7066
Table 3. Results of the two numerical schemes for free particles of the infection of the HIV virus in the blood V (s) under the fractal order and fractional order δ = σ = 0.95 with the step size h = 0.1 .
Table 3. Results of the two numerical schemes for free particles of the infection of the HIV virus in the blood V (s) under the fractal order and fractional order δ = σ = 0.95 with the step size h = 0.1 .
Time: (s)0.0001100200300400
Adams-Bashforth method0 1787.4317 2169.7278 2082.8420 2082.8420
Newton polynomials method0 1787.8611 2169.8866 2082.8597 2082.8472
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Najafi, H.; Etemad, S.; Patanarapeelert, N.; Asamoah, J.K.K.; Rezapour, S.; Sitthiwirattham, T. A Study on Dynamics of CD4+ T-Cells under the Effect of HIV-1 Infection Based on a Mathematical Fractal-Fractional Model via the Adams-Bashforth Scheme and Newton Polynomials. Mathematics 2022, 10, 1366. https://doi.org/10.3390/math10091366

AMA Style

Najafi H, Etemad S, Patanarapeelert N, Asamoah JKK, Rezapour S, Sitthiwirattham T. A Study on Dynamics of CD4+ T-Cells under the Effect of HIV-1 Infection Based on a Mathematical Fractal-Fractional Model via the Adams-Bashforth Scheme and Newton Polynomials. Mathematics. 2022; 10(9):1366. https://doi.org/10.3390/math10091366

Chicago/Turabian Style

Najafi, Hashem, Sina Etemad, Nichaphat Patanarapeelert, Joshua Kiddy K. Asamoah, Shahram Rezapour, and Thanin Sitthiwirattham. 2022. "A Study on Dynamics of CD4+ T-Cells under the Effect of HIV-1 Infection Based on a Mathematical Fractal-Fractional Model via the Adams-Bashforth Scheme and Newton Polynomials" Mathematics 10, no. 9: 1366. https://doi.org/10.3390/math10091366

APA Style

Najafi, H., Etemad, S., Patanarapeelert, N., Asamoah, J. K. K., Rezapour, S., & Sitthiwirattham, T. (2022). A Study on Dynamics of CD4+ T-Cells under the Effect of HIV-1 Infection Based on a Mathematical Fractal-Fractional Model via the Adams-Bashforth Scheme and Newton Polynomials. Mathematics, 10(9), 1366. https://doi.org/10.3390/math10091366

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop