Corrigendum Corrigendum to ( Global Stability of Malaria Transmission Dynamics Model with Logistic Growth )

In the article titled “Global Stability of Malaria Transmission Dynamics Model with Logistic Growth” [1], there was a sign error on the global stability analysis of the endemic equilibrium point where a preceding sign should be minus (-) instead of plus (+) in equation (34) and inequality (35). e corrected equation and inequality are shown as follows: dV dt = nh (1 − nh) (2 − nh∗ nh − sh∗ sh ) + αsh∗ (2 − nh nh∗ − sh sh∗) + βsh∗iV∗ (3 − ih ih∗ − nhiV∗ nh∗iV − shih∗iV sh∗ihiV∗) + βsh∗iV∗ (1 − iV∗ iV )( iV iV∗ − nh nh∗) + ρih∗ (nh∗ nh − 1)( nh nh∗ + ih ih∗) + nV∗iV∗ ( nVih nV∗ih∗ − iV iV∗)(1 − iV∗ iV ) − ε (nV − nV∗) 2 − ih∗iV∗ (1 − ih ih∗)(1 − iV iV∗) . (34)


Introduction
Malaria is a mosquito-borne disease caused by Plasmodium parasite, which is transmitted through the bites of an infected mosquito.In 2017, the World Health Organization report reveals estimations of 216 million malaria cases and 445 thousand deaths due to malaria were registered worldwide in 2016.However, the most malaria cases and deaths were shared by the WHO Africa region, which account for 90% of cases and 91% deaths.The most predominant malaria parasite in the WHO Africa region is Plasmodium falciparum, accounting for 99% of malaria cases in 2016 [1].
Malaria is entirely preventable and treatable disease if the recommended interventions are properly applied.Individuals should have taken some aggressive measurements to decline malaria burden.Personal protection measures are the first line of defense against mosquito-borne diseases.Mosquito repellent is a method used for personal protection; and these are the substances used for exposed skin to prevent human-mosquito contact.Insecticide Treated Bed Nets (ITNs) are used for individuals against malaria to reduce the morbidity of childhood malaria (below five years of age) by 50% and global child mortality by 20-30% [2,3].When used on a large scale, ITNs are supposed to represent efficient tools for malaria vector control but there is a limitation of resistance to insecticides used for a saturated net.The resistance of the most important African malaria Anopheles gambiae to protrude is already widespread in several West African countries [4,5].
Nowadays, mathematical models become an important and popular tools to understand the transmission dynamics of the disease and give an insight to reduce the impact of malaria burden in the society.This is because mathematical modeling can answer the following questions raised by the public health authorities and policy makers to make the correct decisions: (1) how severe will the epidemics be? (2) How long will it last?(3) How effective will an intervention be? (4) What are the effective measures to control and eliminate an endemic disease?The earliest malaria model study originated from Ross in 1911 [6] and later modification made by Macdonald [7].Some further extensions of Ross-Macdonald models for malaria were described in [8][9][10][11][12][13].Tumwiine et al. [13] define the reproduction number,  0 , and show the existence and stability of the disease-free equilibrium and an endemic equilibrium.Recently, many works on host-vector interaction models have been done in [14][15][16][17][18][19][20][21][22][23][24].In [18,20,22,25], global stability of equilibria has been investigated using suitable Lyapunov functions; and their results show that the disease-free and endemic equilibrium points become globally asymptotically stable if  0 ≤ 1 and  0 > 1, respectively.Application of the optimal control theory becomes an important tool for investigating the efficiency of joint control intervention strategies to minimize the impact of malaria disease and cost-effectiveness of implementing them [19,21,23,24].Their studies suggest that the optimal control strategies can effectively reduce the malaria disease.
Motivated by the above studies, we extend the model presented in [14] by taking into account a logistic model with population dependent birth rates for both human and vector populations that describes self-limiting growth of both the human host and mosquito vector populations.We consider logistic malaria model as no population can grow exponentially at all time, in general.A number of populations initially grow exponentially, but, due to competition and limited resources availability, their population size decline, after some time, to a stable size , called the maximum carrying capacity.The competition for limited resources (including food, territory, light, water, and oxygen) decreases the fertility or survival of individuals.Furthermore, this paper presents application of the model to study abrupt and periodic variations of malaria and sensitivity analysis applied to understand the important parameters in transmission and prevalence of the malaria disease.The purpose of this work is to investigate the global stability of both disease-free and endemic equilibrium points.

Malaria Model
We consider that the total human host population,  ℎ , at a time  * , is divided into three disjoint compartments: susceptible  ℎ , infectious  ℎ , and recovered  ℎ .The total vector population,  V , at a time  * , is divided into two mutually exclusive subpopulations of individuals who are susceptible,  V , and infectious,  V .The susceptible human and vector populations are recruited at the rates  ℎ and  V , respectively, where The susceptible human and vector populations decrease due to natural death at a rate  ℎ for humans and  V for vectors, and those that move to the infected classes at a rate  ℎ  V and  V  ℎ , respectively.The infected human population grows as a result of new infection at a rate  ℎ  ℎ  V and decline due to natural mortality, disease induced death, and recovery at a rate  ℎ ,  ℎ , and  ℎ , respectively.For details, see the schematic diagram of the model in Figure 1.The state variables and parameters for the model are described in "State Variables, Parameters, Descriptions, and Their Dimensions of Malaria Model" section.
The model has made the following assumptions: both the total sizes of human and vector populations not being constant; all variables and parameters involving the model assumed to be nonnegative; all newborns susceptible to infection; mosquitoes not dying because of infection; no recovery compartment for infected mosquitoes; and the recovered human population developing permanent immunity.From the schematics diagram of transmission of malaria between human and mosquito (see Figure 1), we have the governed differential equations which describe the dynamics of malaria, with initial conditions At all times,  ℎ =  ℎ +  ℎ +  ℎ and  V =  V +  V .Moreover, their differential equations are satisfying respectively.We may notice that the vector population equation is completely decoupled from the human equations which is physically reasonable.If we eliminate  ℎ and  V and add the total population equations, then we finally have 2.1.Basic Properties.Since the model system (1) involves human and mosquito populations, all its associated variables and parameters are nonnegative.
Theorem 1. Solutions of the model system (1) with positive initial data will remain nonnegative for all time  ≥ 0.

Invariant
Region.The malaria model (1) will be analyzed in biologically feasible region.Thus, the feasible solutions set for the model written by is positively invariant and then the model is biologically meaningful and mathematically well posed in the domain Ω.
The proof is omitted for simplicity.

Model Analysis
To analyze the malaria model in system (4), we use the normalized quantities instead of the actual populations.Since  ℎ and  V may vary, these scales are not suitable for use in the scaling.However, the typical choice for logistic models is to use the sustainable populations  ℎ and  V for the scales.
In the present case, we shall also consider varying  ℎ and  V .It is, therefore, convenient to write  ℎ =  0 ℎ  ℎ () and , where  0 ℎ and  0 V are typical sizes and where  ℎ and  V take care of the time variations.At the moment, we just assume that  ℎ () =  V () = 1.We shall scale the time  * with the quantity 1/ ℎ by setting  =  ℎ  * .The scaling is then The scaled equations then become subject to suitable initial conditions,
The basic reproduction number,  0 , is the single most important parameter in epidemiological modeling.It measures the average number of the secondary infections caused by a single infective in an entirely susceptible population during its whole infectious period [26].To derive the basic reproduction number  0 of model ( 9), we use the next generation matrix approach described in [27][28][29].The infected compartments of system ( 9) are  ℎ and  V .Following [29], the new infection matrix  and the transition matrix  are given, respectively, by Hence, the basic reproduction number,  0 , is the dominant eigenvalue of the next generation matrix  −1 and becomes From (12), it is noted that the reproduction number depends on the product of the number of humans that one mosquito infects through its infectious lifetime,  0Vℎ , and the number of mosquitoes that one human infects through its infectious lifetime,  0ℎV .
Theorem 2. The disease-free equilibrium point,  0 , is locally asymptotically stable if all eigenvalues of the characteristic equation of the variational matrix lie below zero.
Proof.At the equilibrium point,  0 , the variational matrix is given by The characteristic equation may be written as det[( 0 ) − ] = 0.It implies Clearly, we have Thus, this shows that the solution,  0 , is unstable since  1 and  4 lie above zero.
Proof.The variational matrix at the equilibrium point,  1 , becomes Thus, the characteristic equation of the variational matrix is given by where The characteristic polynomial in (17) has roots  1 =  − 1,  2 =  − 1,  3 = −, with negative real parts since 0 < ,  < 1.By Routh-Hurwitz criterion [28], the other roots  4 and  5 have negative real parts if both  0 and  1 lie above zero.
Theorem 4. The disease-free equilibrium point,  1 , is globally asymptotically stable in Ω if  0 ≤ 1; otherwise it is unstable.
Proof.Consider the Lyapunov function where The time derivative of the function  along the solutions of (9) becomes Thus V ≤ 0 if  0 ≤ 1 and the equality V = 0 holds if and only if  ℎ =  V = 0. Therefore, the largest compact invariant set in {( ℎ ,  V ) ∈ Ω : V = 0} is the singleton { 1 }, where  1 is the disease-free equilibrium.LaSalle's Invariant Principle [30] implies that  1 is globally asymptotically stable in Ω.

Stability of Endemic Equilibrium.
To find the endemic equilibrium  2 , we shall keep the assumptions about  V * and  V * from the singular perturbation equations (see the fourth and fifth equations of system ( 9)) and then focus on the first three equations of system (9).We consider the equilibrium solutions, using  ℎ * as the basic quantity.From the third equation of ( 9), we have an expression for  ℎ * : Addition of the second and third equations of system (9) gives Now, let us consider the first equation of ( 9) which connects  ℎ * and  ℎ * .That is, Equation ( 24) has only real and positive solutions for , then (24) has two solutions: for each value of  ℎ * .
One might wonder what happens when the inequality is violated and  ℎ * > ((1−)/2) 2 .The solution of (24) will then be majorized by the equation which has only one unstable equilibrium point at  = ((1 − )/2) and otherwise tends to 0. Moreover, it is clear that  ℎ * needs to be larger or equal to  ℎ * .Let us, therefore, consider  + ℎ * and  − ℎ * in the light of this restriction.The various situations are best described in terms of the graph shown in Figure 2.
In the graph, the upper solution (upper red dot) is acceptable since  + ℎ * >  ℎ * , whereas the lower solution is outside the region and hence unacceptable.Further inspection of the graph shows the following.(3) For (1 − )/2 <  < 1 − , one or both of the solutions are acceptable.
The main conclusion is that there are acceptable solutions with respect to size for all  ℎ * where 0 ≤  ℎ * ≤ (1 − ) 2 /4.
Assume that  ℎ * < (1 − ) 2 /4, where only  + ℎ * is acceptable with respect to size.Then substituting expressions ( 22) and ( 25) into the equation in (23) and simplifying lead to the following quadratic equation: where From (27), it can easily be seen that  > 0. Further, if  0 > 1, then  < 0. Thus, the number of possible positive real roots of ( 27) can depend on the signs of .This can be analyzed using the Descartes rule of signs on the quadratic polynomial (27).The different possibilities for the roots of ( 27) are tabulated in Table 1.
Thus, the malaria model has a unique endemic equilibrium if  0 > 1 and whenever cases 2 and 3 are satisfied.Hence, the endemic equilibrium then becomes and  ℎ * is the unique positive root of ( 27).
Proof.We shall propose the Lyapunov function where The Lyapunov function  is continuous for all  ℎ ,  ℎ ,  ℎ ,  V ,  V > 0. The time derivative of the function  along the solutions of system (9) becomes From the equilibrium point of the malaria model ( 9), we have the following relations: By adding and subtracting  ℎ *  V * ,  ℎ *  V * ( ℎ  V * / ℎ *  V ) and using ( 31) and ( 33) in (32), after intensive simplification, we have ) . ( Since the arithmetic mean is greater than or equal to the geometric mean, then we have Also, if Furthermore, if Hence, it follows from (34), ( 35), (36), and (38) that / ≤ 0 in Ω.Thus, the equality / = 0 holds only when , where  2 is the endemic equilibrium.From the LaSalle's invariant principle [30], the unique equilibrium  2 of system ( 9) is globally asymptomatically stable for  0 > 1.

Sensitivity Analysis.
To understand the relative importance of parameters which are responsible for transmission and prevalence of malaria disease, described in the model ( 9), we perform a sensitivity analysis.Sensitivity indices help us to measure the relative change in a state variable while a parameter changes.The normalized sensitivity index of a variable to a parameter is the ratio of the relative change in the variable to the relative change in the parameter [27].We calculate the sensitivity indices of  0 to assess which parameter has a great impact on  0 and hence the greatest effect in determining whether the disease dies out or persists with population.
Let  be the generic parameter of model (9).We, now only, derive the normalized sensitivity index of  0 to each of the parameters involved in  0 , defined by the ratio of the relative change in  0 to the relative change in the parameter ; that is, This index shows how sensitive  0 is to a change in the parameter .We notice that This indicates that Π  0  does not depend on any parameter value.Similarly, for the other parameters, we have We evaluate the above sensitivity indices, in Table 3, using the parameter values in Table 2.The basic reproduction number,  0 , is most sensitive to the contact rates of human to vector and vector to human, with   0  = 0.5000 and   0 ] = 0.5000 as it can be seen in Table 3.This shows that any increase (decrease) by 10% in  or ] will increase (decrease) by 5% in  0 .The other parameters with highest sensitivity indices are , with   0  = −0.6266,and , with   0  = −0.3812.Increasing (decreasing)  by 10% will decrease (increase) in  0 by 6.266% and increasing (decreasing)  by 10% will decrease (increase) in  0 by 3.812% or vice versa.The rest of parameters,  and , have less significant effect in  0 .
In conclusion, the vector death rate, the human induced death rate, and the contact rates are important parameters in the model which have a significant impact on prevalence and transmission of the malaria disease; these parameters are able to control so that an intensive effort/work has to be done to eradicate the malaria disease from the population.Furthermore, one can understand from the sensitivity indices that vector control is the most effective control strategy.

Application of the Model
In this section, we present more simulations illustrating the abrupt and periodic variations of the model.We fix a reasonable parameter values of the model for numerical simulations.Steps down the vector population by 50%.
We allow the mosquito sustainable level,  V () =  0 V  V (), noted in Section 3, to vary with respect to time.However, we first keep  ℎ fixed in order to investigate the impact of fast variation in  V on the human populations.Periodic variations in  V are shown in these plots for different periods.Plots of the abrupt changes in the sustainable populations,  ℎ =  0 ℎ  ℎ () (see in Section 3) and  V , are located in Figures 3-7, whereas plots of periodic variations are shown in Figures 8  and 9.
Abrupt changes in the human and mosquito populations may, for example, be due to intensive spraying of the mosquitoes some massive emigration (refugee camps) or immigration for the humans.In Figures 3 and 4, steps down and up in vector sustainable population about 50% are plotted.These plots show that the transition occurs very fast and the system adjusts quickly to the new equilibrium.In Figures 5-7 increases about 50% and 100% and decrease about 50% in  ℎ are shown.Step change in humans needs caution since it may lead to unphysical solutions.In Figure 5, after a transient, the solution converges to the new equilibrium.An increase in  ℎ in Figure 6 may lead to an increase in the fraction of the human susceptible population for small time intervals, but not much dramatic change is shown.However, the slow change in human population after transient is shown in Figure 7.In the plots of periodic variation, one observes that the fast variation is quickly adapted by the vector population.This is because of the fast time scale for the mosquito population.Most of the figures show that the humans follow a   slow variations relative to the vectors.Therefore, it is possible to say that fast variations in  V do not imply large variation in human population.In general, for periods shorter than 1/ ℎ , the human population do not in practice show the variations   in the vector populations, but for long periods they do, but apparently weaker.

Conclusions
In this work, we developed and analyzed a logistic malaria model to study the global stability of both disease-free and endemic equilibrium points.Mathematically, we formulated a five-dimensional system of deterministic ordinary differential equations and defined the domain where the model is epidemiologically feasible and mathematically well-posed.The model used the next generation matrix approach to obtain an explicit formula for a reproduction number,  0 , which is the expected number of secondary cases produced by a single infectious individual during its entire period of infectiousness in a fully susceptible populations.
Qualitative analysis of the model determines stability analysis of the equilibrium points.Accordingly, we obtained two diseases-free equilibrium points  0 and  1 .The equilibrium point,  0 , is unstable and unphysical, while the equilibrium point,  1 , becomes both locally and globally stable whenever  0 < 1 and  0 ≤ 1, respectively.We also have shown that the endemic equilibrium point,  2 , is globally asymptotically stable if  0 > 1.Furthermore, sensitivity analysis of the model shows that the human induced death rate, the contact rates (human to mosquito or vice versa), and mosquito death rate have a significant effect on transmission and prevalence of the malaria disease.Moreover, numerical simulations are carried out in the application of the model to investigate how variations in the sustainable level of the vectors affect the human population.One can see from these simulations that fast variations in  V do not lead to large variations in the human population.

Figure 1 :
Figure 1: The compartmental model for malaria transmission.

Figure 2 :
Figure 2: Graph illustrating the behaviour of the solutions of the equation between  ℎ * and  ℎ * .

Figure 4 :
Figure 4: Numerical simulation for the fractions of human and vector population versus time with constant parameter values the same as in Figure 3 and steps up the vector population by 50%.

Figure 5 :
Figure 5: Numerical simulation for the fractions of human and vector population versus time with constant parameter values the same as in Figure 3 and steps up the human population by 50%.

Figure 6 :
Figure 6: Numerical simulation for the fractions of human and vector population versus time with constant parameter values the same as in Figure 3 and increasing the human population by 100%.

Figure 7 :
Figure 7: Numerical simulation for the fractions of human and vector population versus time with constant parameter values the same as in Figure 3 and  = 0.001, and steps down the human population by 50% in the form of a constant drop over a time interval of length 100.

Figure 9 :
Figure 9: Periodic variations in  V with  V = 1 + 0.3 sin(/6), for  ℎ = 1, and the other parameters and initial values remain fixed as in Figure 8 apart from  = 0.5 over the time interval [0 100].
ℎ : Number of susceptible humans at a time  *  ℎ : Number of infected humans at a time  *  ℎ : Number of recovered humans at a time  *  V : Number of susceptible mosquitoes at a time  *  V : Number of infected mosquitoes at a time  *  ℎ : The total human population at a time  *  V : The total mosquito population at a time  *  ℎ : Sustainable level of human population at a time  *  V : Sustainable level of mosquito population at a time  *  ℎ : Per capita birth rate of human population.Dimension: Time −1  ℎ : Per capita natural death rate for humans.Dimension: Time −1  ℎ : Per capita disease-induced death rate for humans.Dimension: Time −1  ℎ : The human contact rate.Dimension: Mosquitoes −1 × Time −1  ℎ : Per capita recovery rate for humans.Time −1  V : Per capita birth rate of mosquitoes.Dimension: Time −1  V : Per capita natural death rate of mosquitoes.Dimension: Time −1  V : The mosquito contact rate.Dimension: Humans −1 × Time −1 .

Table 2 :
Some parameter values of the malaria model; all parameters are nondimensional.

Table 3 :
Sensitivity indices of  0 to parameters for the malaria model.