Uncertainty always lives with us. We cannot take exact measurement of initial conditions or parameters values in a mathematical model. As humans, we are remaining alive in an environment where the uncertainties lie in the modelling of physical phenomena. There might be some incomplete information or estimation of the parameter or initial values. To handle uncertainty, we use fuzzy operators rather than classical operators. In this paper, we study a model of HIV-1 infection by taking uncertainty in the initial data under Caputo fractional operator. We explore the existence and uniqueness of the results through fixed-point theory. We study the Ulam–Hyres stability of the considered model. By using the fuzzy Laplace Adomian decomposition method, numerical results are obtained for specific fuzzy initial conditions. To better understand the behaviour of the fuzzy solution, we present the obtained numerical results graphically for various fractional orders where the uncertainty lies in [0, 1].

1. Introduction

AIDS (acquired immune deficiency syndrome) is one of the most dangerous diseases caused by a pathogen virus called HIV (human immunodeficiency virus) and leads a person to death. In 1982, a patient in the USA was infected by HIV. At the end of 2019, about 38 million people, identified by the WHO (World Health Organization), were surviving with HIV. In 2019, 0.69 million people died from diseases associated with AIDS. To investigate how AIDS can spread or control, mathematical models are valuable instruments as they provide a short or a long-term prediction of the frequency of the disease. Many researchers have investigated the dynamics of HIV infection models [13]. A lot of methods have been used by researchers to find a solution to HIV infections [46]. In [7], the authors presented an HIV-1 infection model which contains five compartments, denoted by (the uninfected CD4+ T cells), (the concentration of infected cells), (the concentration of double cells), and and (the densities of pathogen viruses), respectively. The integer-order model of HIV-1 infection is given bywhere the host cells are produced at the rate . Moreover, are the rates of death of host cells, infected cells, double infected cells, and pathogen virus, respectively. is the rate of infection of healthy CD4+ T cells. The rate of production of HIV-1 by a cell is defined by , is infection rate of removal of HIV-1 infected cells, is the rate of production of recombinant by a double-infected cell, and is the removal rate of recombinant. The mathematician has taken an interest in fractional calculus in recent years as it has many applications in engineering and biological sciences [810]. In the last decade, the researcher converted many models from integer order to fractional order which provides best results as compared to integer-order derivatives. As fractional operator has a greater degree of freedom, different researchers studied the existence theory of the solution [1113]. The fractional operators also have been used for modelling of different infectious diseases. Sweilam et al. investigated the optimal control for variable order fractional HIV/AIDS and malaria mathematical models with multitime delay [14]. The Legendre spectral-collocation method was used for solving fractional optimal control of HIV infection of CD4+T cells’ mathematical model in [15]. A complex fractional-order HIV infection model with drug resistance during therapy has been studied by Sweilam et al. [16]. Some other applications of the fractional-order operators in various field of sciences can be found in [1719]. It is worthy to note here that the fractional derivative in the Caputo sense is the most commonly used definition among the definitions of the fractional derivative. The definition of Caputo is mathematically rigorous than the Riemann–Liouville definition. The Caputo derivative exists in the whole interval (0, 1). In addition, Caputo definition is very welcome in applied science and engineering. Furthermore, properties of the Caputo derivative are helpful in translating the higher fractional-order differential systems into lower ones. For a comparison between Caputo and Riemann–Liouville operators, the interested reader is referred to [20, 21]. Therefore, we will consider the fractional operator in Caputo sense.

Confusion and uncertainty still live with us, and that is a reality. Many human beings are exposed to doubting anything surrounding them and questioning why for themselves or others. Because they are not clear, and their reports are incomplete or inaccurate. Now, suppose that we are in a situation where uncertainty occurs with all this inaccurate knowledge. It is a reality that, in many of our real questions where inaccurate data are involved, we do not know how to answer. For scientists, this mindset and attitude of uncertainty is very necessary. Our aim is to sort out how to understand it and operate through it rather than trying to fight ambiguity because the advancement, resources, and life that you want are uncertain. To overcome this situation, Zadeh introduced fuzzy set in 1965 [22]. The concept of fuzzy mapping and control has been addressed by Zadeh and Chang [23]. Furthermore, researchers extended the ordinary and fractional operators to fuzzy operators [24, 25]. These operators have been used by various authors to study different models [2629]. If the information or data is imprecise, then fuzzy operators model physical phenomena easily and accurately. We are interested in investigating model (1) for fractional-order derivative under the uncertain initial conditions. For fractional order , consideralong with fuzzy initial conditions,such that .

2. Preliminaries

This section provides the basic concepts of fuzzy sets and fuzzy fractional calculus. Let FFD, FFI, and LT denote the fuzzy fractional derivative, fuzzy fractional integral, and Laplace transform, respectively.

Definition 1 (see [29, 30]). A fuzzy set is called a fuzzy number if(i) is fuzzy convex.(ii) is normal, i.e., for some .(iii) is upper semicontinuous on .(iv)The closure of is compact.

Definition 2 (see [29, 30]). The parametric form of a fuzzy number is given by , for , if and only if(i) is increasing, left on and right continuous at 0, respectively.(ii) is decreasing, left on and right continuous at 0, respectively.(iii).

Definition 3 (see [31]). Let be a function, and and are any two fuzzy numbers. Then, the H-distance between and is given byIn , satisfies the properties given below:(i).(ii).(iii).(iv) is a complete metric space.

Definition 4 (see [31]). Let be a continuous fuzzy function on , and a FFI is defined asFurthermore, if , where and are the spaces of fuzzy continuous functions and fuzzy Lebesgue integrable functions, respectively, then the FFI can be written aswhere

Definition 5 (see [31]). A fuzzy function is such that , ; then, the FFD in Caputo sense is given bywherewhere .

Definition 6 (see [32]). Let . Then, the fuzzy LT of is defined asWhile, in the parametric form, the fuzzy LT is given byHence,

Theorem 1 (see [31]). The LT of Caputo fractional derivative is defined as

3. Existence and Uniqueness Results

In this section, the existence and uniqueness of the solution of the fuzzy fractional model via fixed-point theory are discussed. Let us write the right-hand side of system (2) as

We can write the given system (2) as

Now, using fractional integral, we obtain

Now, consider a Banach space as , under the norm,

We can write equation (16) aswhere and We take some assumption on nonlinear function as(1)Lipschtiz condition: there exists a constant such that, for each , (2)Growth condition: there exists constants and such that

Theorem 2. If the growth condition is satisfied, then system (3) has at least one solution.

Proof. Consider is closed and fuzzy convex, and a mapping is defined byThen, for any , we haveIt follows that . Thus, ; it implies that is bounded. To show is equi-continuity, consider such that . Then,Hence,This show that the operator is equi-continuous, which is completely continuous using Arzela–Ascoli theorem, as is bounded. Therefore, system (3) has at least one solution using Schauder’s fixed-point theorem.

Theorem 3. Let the Lipschitz condition hold. Ifthen system (3) has a unique solution.

Proof. Let . Then,which shows that is a contraction. Thus, system (3) has a unique solution by Banach contraction theorem.

4. Ulam–Hyres Stability

In this section, we show that the numerical results are stable. Take a small change , where depends on the solution of as(i)For any positive , (ii)

Lemma 1. The above perturb problem has a solution,if it satisfies

Proof. The proof is easy.

Theorem 4. Consider assumption (2) and relation (29) hold. Then, the numerical results of the considered model are Ulam–Hyres stable if , by showing that equation (25) is Ulam–Hyres stable.

Proof. Let be any solution and be at most one solution of equation (25); then,From (29), we haveHence, from (30), it follows that system (6) is “UH stable.” Consequently, the considered model is UH stable.

4.1. Procedure for Solution

In this section, we deduce the solution of the given model by using LADM. Applying LT on (3), we obtain

The infinite series solution iswhere are Adomian polynomials and represent nonlinear terms.

Taking inverse Laplace transform, we obtain

Comparing the first two terms of the series in parametric form,

We can find other terms by the similar way. Hence, the system has the series solution as follows:

Based on the procedure of the method, the following theorem discusses the convergence of the proposed method.

Theorem 5. Let be a Banach space, , , and . Let be an exact solution of the considered model. Then, the series solution converges to the exact solution whenever for all ; i.e., for any , there exists a positive number such that .

Proof. Let us define the sequence of asWe need to show that is a Cauchy sequence in . For this, we considerNow, for any , we haveSince and is bounded, letting , we obtainThus, the sequence is a Cauchy sequence in . This implies that such that . Thus, the required solution converges to the exact solution.

5. Numerical Results and Discussion

In this section, we illustrate the obtained solution of the proposed model through MATLAB-17. We have taken the parameter values from [7]. The values of parameters are , , , , , , , , , and . Consider model (2) under initial conditions aswhere . Using the developed technique on model (2) by applying the initial conditions, we obtain