Abstract

In this study, a Caputo fractional derivative is employed to develop a model of malaria and HIV transmission dynamics with optimal control. Also, the model’s basic properties are shown, and the basic reproduction number is computed using the next-generation matrix method. Additionally, the order of fractional derivative analysis shows that the infected group decreases at the beginning for the higher-order of fractional derivative. Moreover, the early activation of memory effects through public health education reduces the impact of malaria and HIV infections on further progression and transmission. On the other hand, effective optimal controls reduce the occurrence and prevalence of HIV and malaria infections from the beginning to the end of the investigation. Finally, the numerical simulations are done for the justification of analytical solutions with numerical solutions of the model. Moreover, the MATLAB platform is incorporated for numerical simulation of the solutions.

1. Introduction

Human immunodeficiency virus (HIV) is a human harmful virus that causes a fatal disease called acquired immunodeficiency syndrome (AIDS) [1]. There is no medication that cures AIDS disease, but effective antiretroviral (ARV) drugs could control the transmission and progression of the infection. HIV can be prevented using the prominent rule stated as abstinence-be-faithful-use-condom which is abbreviated as ABC rule. HIV can be transmitted to individuals through contact of fluids with HIV exposed objects. HIV can be transmitted vertically from mother to child and horizontally from person to person in an unsafe exposure to HIV exposed environment. Unsafe sexual practice is one of the major modes of HIV transmission. Malaria infection is a curable infection caused by a parasite called Plasmodium if the infected mosquito bites humans to suck the blood [2]. Globally, the consequence of deforestation is causing an increase in global warming that facilitates the risk of malaria infection. In the most part of sub-Saharan Africa, malaria cases and deaths are largely visible. The malaria and HIV diseases are studied separately for a long period of time as fatal diseases. However, the co-infection of HIV and malaria killed about 2 million each year [1].

Recently, mathematical models have become indispensable tools to describe the transmission and progression dynamics of disease in the human population [311]. The first mathematical model of malaria transmission dynamics is introduced by Ross [7]. Then, different mathematical models are constructed to describe transmission dynamics of malaria infection. Ngwa and Shu developed the malaria model with variable human and mosquito population [9]. Chitins et al. [10, 11] studied the transmission of malaria, considering that human and vector species follow a logistic population. They also studied the bifurcation analysis and the effects of seasonality that facilitate the spreading of the mosquito and feeding on the host. Li [8] studied the transmission of malaria through the stage-structured model and provides a basic analysis that shows backward bifurcation.

Boukhouima et al. [12, 13] studied the dynamics of virus in the human body through a fractional derivative. Hattaf [14] developed the new generalized fractional derivative and applied it to the analysis of HIV dynamics in the body. Wu et al. [15] studied the cellular dynamics of the HIV-1 with uncertainty in the initial data. Dokuyucu and Dutta [16] studied the TB-HIV dynamics with the fractional order derivative. The order of the fractional derivative contributes to the memory effects. Okyere et al. [17] studied the transmission dynamics of malaria using the optimal fractional order derivative. Okyere et al. [18] studied the invasion of disease in the ecological population, whereas the diffusive model of COVID-19 is studied in [19]. As a basis for the memory effect, awareness intervention has contributed to reducing disease transmission [20]. In addition, the work done in [2131] can be used to illustrate how mathematical models have recently played a part in providing information for the public through fitting data and new methods. In addition, in [3134], co-infection models are constructed to describe the dynamics of HIV, cholera, and COVID-19 infections. Although some models are developed to study the dynamics of HIV and malaria, the fractional order derivative is incorporated to study the co-dynamics of malaria and HIV transmission dynamics. In this study, we are motivated by the works done in [1] and extend it to the fractional derivative model.

In addition, the remaining portions of the work are organized as follows. Section 2 discusses the creation of the mathematical model. Section 3 presents a basic analysis of the model without control. Section 4 presents the optimal control problem. Numerical simulation is presented in Section 5. Results and discussions are presented in Section 6. The conclusion is presented in Section 7.

2. Mathematical Model Formulation

In this study, the total population is classified into compartments as follows: (i) susceptible : It consists of all human who are infection-free and have a chance to be infected in the future; (ii) malaria-infected : This consists of only malaria-infected individuals; (iii) HIV-infected: It consists of only HIV-infected individuals; (iv) HIV-malaria-infected: It consists of all individuals infected with both malaria and HIV infections before the AIDS stage; (v) AIDS patients: It consists of individuals infected with HIV at the AIDS stage; (vi) AIDS-malaria-infected: It consists of individuals infected with malaria at the AIDS stage; (vii) susceptible mosquito : It consists of noninfectious mosquito population with chance to become infectious if they suck the blood of malaria infected human individuals; and (viii) infectious mosquito : It consists of infectious mosquito vector.

Moreover, the state variables and parameters used in the model are described in Tables 1 and 2, respectively.

Considering the structural representation given in Figure 1, we have constructed the fractional model for HIV and malaria transmission dynamics without control as given by the following equation:with initial conditions .

3. Mathematical Analysis of the Model

3.1. Invariant Region

Theorem 1. The solution of Caputo fractional derivative model (1) is invariant in the invariant set such thatwhere .

Proof. Considering the human population equation and adding the right and left expressions of the fractional derivative part giveswhich impliesApplying the Laplace transform on both sides and solving the preceding expression as t reaches infinity, we obtainHence, the solution of the developed is bounded for all time t.

3.2. Non-Negativity of Solution

Theorem 2. The solution of Caputo fractional derivative model (1) is non-negative in the invariant region for all time .

Proof. To show the non-negativity of solution, consider model (1) along the axis where state variables vanish as follows:Hence, by Caputo generalized mean value theorem, the solution of constructed model is non-negative.

3.3. Existence and Uniqueness of Solutions

Theorem 3. The solution of Caputo fractional derivative model (1) exists and is unique in the defined invariant region .

Proof. The proof can be shown using the fixed point theory [5].

3.4. Disease-Free Equilibrium

The disease-free equilibrium of Caputo fractional derivative model (1) is a steady state point of disease extinction. It is computed and given by

3.5. Basic Reproduction Number

The basic reproduction number is extensively applied in mathematical epidemiology to describe the status of infections invading the population. The biological significance of basic reproduction number is that it shows the status of the disease in the population, whether the disease in the population is extinct if or persists if [3537]. Recently, a next-generation matrix method is broadly incorporated in the modern research on infections transmission dynamics. Here also, similar to the works done in [3537], we have applied the next-generation matrix to compute the basic reproduction number. Hence, from model (1), we construct a matrix that comprises of newly infected individuals arriving in the compartments, and a matrix encompasses the transition terms in the infected compartments.

First, from only the malaria model, we construct basic reproduction using the subsequent vectors.

At disease-free equilibrium, the Jacobian matrices and obtained from preceding vectors and , respectively, are given by

The next-generation matrix is constructed from the preceding matrices and as follows:

The eigenvalues of the next-generation matrix are computed as follows:

Furthermore, the basic reproduction number is the spectral radius of next-generation matrix . Hence, we have

Second, from only the malaria model, we construct basic reproduction using the subsequent vectors.

At disease-free equilibrium, the Jacobian matrices and obtained from preceding vectors and , respectively, are given by

The next-generation matrix is constructed from the preceding matrices and as follows:

The eigenvalues of the next-generation matrix are computed as follows:

Furthermore, the basic reproductions number is the spectral radius of next-generation matrix . Hence, we have

Moreover, according to the works done in [1], the basic reproduction number of the full model (1) is given by

3.6. Global Stability of Disease-Free Equilibrium

Theorem 4. The disease-free equilibrium of fractional model is globally asymptotically stable if and unstable if .

Proof. To show global stability, we follow the works of Castillo–Chavez and Song done in [2]. Now, we consider the infected compartments of model (1) as given below. First, let be vector of variables for infection-free compartments and be vector of infected compartments such thatwhich implieswhere the subsequent matrix is a Metzler matrix.Hence, the conditions of Castillo–Chavez and song are satisfied. Therefore, the disease-free equilibrium of model (1) is globally asymptotically stable.

4. Extension to Optimal Control Problem

In this section, we study the optimal control problem of HIV and malaria co-infection transmission dynamics for . Since the optimal control problem is very methodical and hypothetical, it is applied in real life by taking feedback from the target community based on numerical solutions of model. Moreover, the control functions are taken in the study can be described as follows: (i) HIV prevention control (): The intervention with this control reduces the chance of getting with HIV infections, (ii) malaria prevention control (): This control function helps in the prevention of malaria infection, (iii) malaria treatment control (): This control function is applied to treat individuals infected with malaria infection, (iv) HIV progression control (): This is HIV treatment control function applied for slowing HIV infection to the advanced stage.

4.1. Objective Function

The objective functional for minimization of cost of using controls and the number of infected compartments is defined by

4.2. Hamiltonian Function

Pontryagin’s maximum principle is applied for the construction of the Hamiltonian function that minimizes the objective function, as given by subsequent equation (38). Let be adjoint variables corresponding to state variables , respectively. Hence, we define Hamiltonian function as follows:which implies

4.3. Adjoint Equations

The adjoint equations are computed from the subsequent partial derivatives, that is,

Therefore, the adjoint equations are computed as follows:with transversal conditions .

4.4. Optimal Control Functions

The optimal controls and are obtained by solving the subsequent equations.

To obtain , we solve, from the subsequent equation, for .

The above implies

To obtain , solve, from the subsequent equation, for .

The above implies,

To obtain , solve, from the subsequent equation, for .

The above implies,

To obtain , solve, from the subsequent equation, for .

The above implies,

In compact form

4.4.1. Optimality System

5. Numerical Simulations

In this section, we have used the compacted fde12 function which is designed to solve Caputo fractional derivative [39]. MATLAB platform is used for numerical simulations of population size to describe effects on memory and applied control measures. The some of the constants used in the simulations are . Moreover, the applied initial value of state variables are and parameters values are given in the subsequent Table 3.

6. Results and Discussion

The fractional derivative model and the optimal control problem are analyzed with the support of numerical simulations.

In Figure 2, the susceptible population size decrease with high order of fractional derivative for long period of time. In Figure 3, the numerical simulations show that the size of only malaria infected human population decrease with high order of fractional derivative at the beginning of simulations. However, after a few years, for long time of interval, the size of malaria infected population increase as the order of fractional derivative increase. Moreover, at the end of numerical simulation, there is an inclination that high order of fractional derivative reduces size of malaria infected population.

In Figure 4, at the beginning of simulations, for high order of fractional derivative the size of only HIV-infected population decreases. Moreover, the small value of fractional derivative reduces a little whereas the HIV-infected population size decrease and increase for high order of fractional derivative. In Figure 5, the high order of fractional derivative reduces HIV-malaria infected population size only at the beginning of infection. In Figure 6, AIDS individuals’ size is low for high order of fractional derivative after some specified time span of simulation. In Figure 7, the high order of fractional derivative reduces the infected population in the beginning of simulation. Moreover, for long period of time the high order of fractional derivative increase the infected individuals.

In Figure 8, the size of susceptible mosquito population decreases with high order of fractional derivative starting from the beginning of simulation time. In Figure 9, at the beginning, the size of infected mosquito population decrease for high value of order of fractional derivative. In Figure 10, it is diagnosed that the size of malaria infected population decrease due to intervention with control measures.

In Figure 11, it is observed that the presence of applied controls reduces the size of HIV-infected population effectively at the beginning and at the end of simulations time. In Figure 12, HIV-malaria infected population size decrease with presence control intervention. In Figure 13, the size of AIDS population decreases due to interventions with control measures for all time of simulation. In Figure 14, the simulation results show that the intervention with applied optimal controls reduces the AIDS-malaria infected population size whereas in Figure 15, the number of infected mosquito population decrease with applied available controls. In general, the study demonstrates that early intervention to manage the dynamics of HIV and malaria transmission by prevention and treatment can lower infection-related mortality and morbidity. Additionally, early treatment intervention can lower the number of deaths brought on by HIV and malaria. As a result, public health education is necessary to raise awareness within the population about HIV and malaria preventive and treatment measures as well as the effects of HIV and malaria coinfections.

7. Conclusion

The high order of fractional derivative reduces the size of infected population at the beginning of the infection. Thus, early activation of population memory effects can reduce the size of infected population. On the other hand, the intervention with optimal control reduces the number of the infected population. Therefore, the results obtained from numerical simulations show that it is advisable to apply public health education for reducing the impact of HIV and malaria infections at the beginning or before the occurrence of infections whereas applying optimal control exhaustingly from the beginning to the entire period of infection reduces the impact of infections. [4146].

Data Availability

All data used in the manuscript are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.