Dynamical Analysis on a Malaria Model with Relapse Preventive Treatment and Saturated Fumigation

Malaria has produced health issues in many parts of the world. One of the reason is due to the recurrence phenomenon, which can happen years after the main infection has appeared in the human body. Furthermore, the fumigation intervention, which has become a major worry in several regions of the world, has yielded unsatisfactory results, as seen by the high number of cases reported each year in several African countries. We present a novel mathematical model that integrates tafenoquine treatments to prevent relapse in the human population and saturation fumigation to control mosquito populations in this study. The endemic threshold, also known as the basic reproduction number, is calculated analytically, as is the existence and local stability of the equilibrium points. Through careful investigation, we discovered that the malaria-free equilibrium is locally asymptotically stable if the basic reproduction number is less than one and unstable if it is greater than one. According to the sensitivity analysis, the utilization of tafenoquine treatment is inversely proportional to the basic reproduction number. Although our model never exhibits a backward bifurcation at the basic reproduction number equal to one, we have demonstrated that it is possible; when the basic reproduction number is greater than one, two stable malaria-endemic equilibrium can exist. As a result, when the basic reproduction number is more than one, the final state will be determined by the initial condition of the population. As a result, enormous temporal fumigation can shift the stability of our malaria model from a big endemic size to a smaller endemic size, which is more advantageous in terms of the malaria prevention strategy. Despite the fact that this is not a case study, the numerical results presented in this article are intended to support any theoretical analysis of current malaria eradication tactics in the field.


Introduction
Malaria is a vector-borne disease caused by the bite of a female mosquito that has been infected with Plasmodium. Of more than 100 species, only five Plasmodium species cause malaria, namely, Plasmodium vivax, Plasmodium malaria, Plasmodium falciparum, Plasmodium knowlesi, and Plasmodium ovale [1]. When this Plasmodium has entered the human bloodstream, it will attack several vital organs in the human body, especially the liver and red blood cells [2]. People who have been infected with malaria will show a variety of symptoms, including chills, fever, and headache, which can even result in death in most cases in the pediatric population.
Until now, there have been many interventions launched by governments in various countries in the world to tackle the spread of malaria. These interventions include the use of vaccines, treatment, use of insecticide-treated bed nets (ITN), and vector control with fumigation and larvicides [3]. Among these mentioned interventions, vector control with fumigation is considered as the most promising and easiest policy to implement [1]. However, several problems in its implementation arise, such as the tendency of mosquitoes to become resistant to fumigants when the intervention is not well controlled [4] or the problem of limited implementation costs. In some cases in the field, the high fumigation intensity needs to be reduced when infected people are too high. This is due to the difficulty of implementation in the field when intervention costs must also be allocated to treatment for infected individuals in the hospitals.
Vaccines for malaria have become one of the main concerns of governments in many parts of the world and the World Health Organization [5]. In 2021, the R21/Matrix-M vaccine has been investigated to be the second malaria vaccine, and it is stated that it has reached the minimum efficacy limit required by WHO, which is a minimum efficacy of 75% [6]. This type of vaccine has an efficacy level of 77% to reduce the chance of successful infection in humans due to an infected female Anopheles mosquito bite. In addition to vaccination, treatment interventions are also needed to prevent the severity or incidence of relapse in malaria patients. Until now, it was stated that primaquine was the primary drug used to avoid relapse in individuals infected with malaria. However, because this drug has to be taken on a regular basis (every 14 days), it has resulted in many treatments not being successful [7]. Therefore, MMV and GlaxoSmithKline (GSK) collaborated to develop a new malaria drug known as tafenoquine, which is a single dose treatment for preventing relapse in malaria-infected individuals [8].
The complexity of malaria has attracted the attention of many researchers to take part in efforts to understand the mechanism of spread and the best intervention for malaria. This is due to the complexity of its infection mechanism, such as recurrence phenomena (relapse, reinfection, and recrudescence), to the problem of the most appropriate intervention. Among these researches, mathematical modeling would play an essential role. Many authors have introduced mathematical models to understand how malaria spreads among human and mosquito populations. The first mathematical model for malaria was introduced by Ross in the early 19 th century [9], where he introduced the mechanism of malaria transmission involving mosquito and human populations in his model. Ross' research was then continued by Macdonald [10] where he introduced the concept of basic reproduction number in his model. Since then, many mathematical models have been introduced by researchers to understand various important factors in the spread of malaria. Authors in [11] proposed a malaria model considering immunological memory which boost protection of reinfection phenomenon. Two host types in malaria transmission are discussed by author in [12]. Furthermore, a two-age class model for malaria transmission is discussed in [13]. A periodic biting rate of malaria mosquitoes is modelled by author in [14]. They used Floquet theory to analyze the stability of their model. Recently, author in [15] proposed a malaria model with optimal control on saturated treatment rate. Another new strategy of transmission blocking drugs for malaria is modelled by Wu and Hu in [16]. They found that increasing the transmission blocking drugs is a more pronounced effect compared to treatment intervention. Another important factors have been discussed through mathematical models such as vector-bias effect [15,17], relapse [18,19], reinfection [20,21], fumigation [15,22], temperature and seasonality [23][24][25], impact of Wolbachia [26], and coinfection [27]. However, the best that we know, there is no mathematical model considering the impact of potential new treatment (tafenoquine) into their model.
In this paper, we introduce and investigate a new mathematical model on malaria transmission. In this model, we divide the human population into five epidemiological classes based on their health status while the mosquito population into two epidemiological classes. Several vital factors were introduced into our model: the effect of a potential new treatment for malaria to prevent relapse; vector-bias phenomena where mosquitoes are more attracted to bite the infected individuals; and fumigation intervention which depends on the number of infected individuals at time t. We perform our mathematical analysis to show the existence of a forward bifurcation and forward hysteresis phenomena on our model, which allows the possibility of existence of three different endemic equilibrium, where two of them is locally stable. Based on this phenomenon, we show from a numerical simulation that a massive fumigation intervention in a limited time window can change the dynamic direction of the system from a large endemic equilibrium to a smaller endemic point. We also show a sensitivity analysis to determine the most influential parameter to our model. This paper is organized as follows. We formulate our model in Section 2. The stability of the malaria-free equilibrium point and the form of the respected basic reproduction number are shown in Section 3. In Section 4, we analyze the existence of the malaria-endemic equilibrium point. In addition, we show the possible forward hysteresis from our model in this section. Some numerical experiments on the proposed model are shown in Section 5. Finally, some relevant conclusions are given in the last section.

Mathematical Model Formulation
The proposed dynamic model for malaria transmission in this article is inspired by our previous work in [15], by taking into account two important factors. The first factor involved in our new model is the involvement of a malaria relapse prevention drug intervention (for example, tafenoquine [28]). The second factor involves fumigation intervention which is not a monotonous function. We assume that when the number of humans infected with malaria is approaching the outbreak, the intervention given can be quite large. However, when the number of infected people continues to grow, efforts for fumigation will be reduced because of the difficulty of intervention during the outbreak.
This model divides the human population based on their health status and whether they received any malaria treatment or not. Hence, let human population be divided into five epidemiological classes: susceptible ðSÞ, latent ðEÞ, infected ðIÞ, exposed treated ðTÞ, and recovered ðRÞ. On the other hand, we classify mosquito population only into two classes: susceptible ðUÞ and infected ðWÞ. The latent individual is an individual who has already been exposed to malaria and has Plasmodium in their lever. If an individual in E gets treated with tafenoquine to prevent relapse, then they will be classified into the class of T. Only infected individual (I) can transmit the Plasmodium into the susceptible mosquito. Based on this assumption, we have the total human population which is given by Computational and Mathematical Methods in Medicine and the total mosquito population is given by The malaria model is governed by the following system of ordinary differential equations: where Π h ðN, WÞ and Π v ðN, UÞ are the infection rate in human and mosquito population, respectively, while ΨðI, u 2 Þ presents the fumigation effectiveness factors. The per capita of birth on humans and mosquitoes is denoted by Λ h and Λ v , respectively. The natural death rate of humans and mosquitoes is given by μ h and μ v , respectively. Furthermore, parameters u 1 and u 2 present medical treatment intervention with tafenoquine and vector control with fumigation, respectively. Letpbe the proportion of exposed individuals who get tafenoquine and succeeded in avoiding relapse afterκ −1 duration of treatment. On the other hand, we assume that the 1 − p proportion of individuals in T failed in treatment. Hence, we have ð1 − pÞδT as the transition from T to I due to treatment failure, where δ −1 is the incubation period of Plasmodium with the effect of tafenoquine. We denote that the recovery rate from malaria is γ, while ξ −1 is the duration of temporal immunity.
We construct the force of infection in human ðΠ h ðN, WÞÞ as follows. Let b be the average bite per mosquito per day. In our model, we take into account the preference of mosquito to be more attracted to bite infected human, rather than noninfected human. This phenomenon is commonly known as "vector-bias" phenomenon [29]. Based on this "vector-biased" assumption, the probability of a mosquito encountering a susceptible human is given by S/ðS + E + T + αI + RÞ, where α > 1 is the vector-bias parameter. Hence, total bite of all mosquito per day is given by bWðS/ðS + E + T + αI + RÞÞ. Assuming ν h as the probability that the bite of infected mosquito succeeded in infecting susceptible human, then bν h WðS/ðS + E + T + αI + RÞÞ present the total of susceptible human who get infected by malaria per time. Since b and ν h are constant parameters with a dimension of bite/day and 1/ðbite × mosquitoÞ, respectively, we assume β h ≔ bν h . Therefore, we have that Using a similar approach, let ν v be the probability of successful infection in mosquitoes; the force of infection on mosquitoes is given by where β v ≔ bν v with a dimension of bite/day and 1/ð bite × humanÞ for b and ν v , respectively. Now, we construct our fumigation term ΨðI, u 2 Þ. We assume that the fumigation intervention depends on the number of infected individuals. Indicators of the endemic of malaria in the field can not be seen from the number of infected mosquitoes, but it can be identified by the high number of infected individuals which is reported in the media. When the number of infected individuals increases, then the intensity of fumigation will increase. However, when the number of infected individuals increases more significantly, then the effectiveness of fumigation will decrease since the policymaker may concentrate more on the number of infected individuals in the hospital, which makes them overwhelmed to control vector population in the field. Hence, we assume that ΨðI, u 2 Þ should have the following properties: (i) When the number of the infected individual is zero, then the fumigation intervention is zero. Hence, we have Ψð0, u 2 Þ = 0 (ii) The fumigation intervention increases at the beginning when the number of infected individual start to increase but will decrease when the number of infected individual is sufficiently large. Hence, we have that ðð∂ΨðI, u 2 ÞÞ/∂IÞ > 0 for I ∈ ð0, I critical Þ and ðð∂ΨðI, u 2 ÞÞ/∂IÞ ≤ 0 for I ∈ ½I critical ,∞Þ. Note that I critical denote the critical number of I when the policymaker is already overwhelmed to conduct an effective fumigation intervention in the field Based on the above assumption, we model our fumigation intervention as where a > 0 is the saturated coefficient.

Computational and Mathematical Methods in Medicine
According to the mentioned assumptions on the infection and fumigation functions, system (3) now is read as with a nonnegative initial conditions given at time t = 0. Figure 1 depicts the flow chart of our malaria transmission model. Biological interpretation and the unity of all parameters in system (7) are given in Table 1.
Let system (7) have an initial condition in the following set: To describe the feasible solution of system (7) and its biological interpretation, we have the following theorem. (8), malaria model in system (7) has a unique solution and remains in D for all time t ≥ 0.

Theorem 1. For initial values in
Proof. Please see Appendix A for the proof.

Malaria-Free Equilibrium and the Basic Reproduction Number
The first equilibrium point of our model is the malaria-free equilibrium point. This equilibrium present a situation where all nonsusceptible population do not exist in the equilibrium condition. For this reason, let E = 0, T = 0, I = 0, R = 0, and W = 0, and then, malaria-free equilibrium ðMFEÞ is obtained by the following subsystem: Taking the right hand side of the above system, it follows that the malaria-free equilibrium of system (7) is given by To conduct further analysis on the qualitative behaviour of our model, it is important to determine the related basic reproduction number of our proposed model. In many epidemiological models, basic reproduction number holds an important role in determining that the diseases die out or exist in the population [34][35][36][37][38]. Basic reproduction number is defined as the expected number of secondary cases caused by one primary case during infection period in a completely susceptible population [39,40]. The basic reproduction number is calculated using the next-generation matrix approach [41]. From system (7), we have that the infected compartments consist of E, T, I, and W. The Jacobian matrix of subsystem of infected compartment on system (7) evaluated in MFE can be written as F + V , where Figure 1: Transmission diagram of malaria model in (7).

Computational and Mathematical Methods in Medicine
where F and V present the transmission and transition terms. Using formula in [41], we have the next-generation matrix ðNGMÞ of system (7) which is given by where E T is the transpose of E, with : ð13Þ Note that each column of F can be spanned by each column of E. Hence, the basic reproduction number of system (7) is given by More example on the method of next-generation matrix method to determine the basic reproduction number in various epidemiological models can be seen in [42][43][44]. The above expression can be rewritten as a multiplication between four important component on malaria transmission on system (7) as follows.
where C 1 = β h /ðu 1 + η + μ h Þ present the number of new latent infected human per infection period of E, C 2 = αβ v /μ v present the number of new infected mosquitoes per infection period of W, C 3 = N v /N h present the ratio of mosquitoes and human, and C 4 = η + u 1 ð1/ð1 + ððpκ + μ h Þ/ðð1 − pÞδÞÞÞÞ present the impact of tafenoquine intervention.
According to Theorem 2 in [45], we have the following theorem regarding the local stability criteria of the malariafree equilibrium of system (7).
The malaria-free equilibrium of system (7) is locally asymptotically stable if R 0 < 1 and unstable if R 0 > 1.

Sensitivity Analysis on the Basic Reproduction Number.
In many mathematical epidemiology models, understanding the impact of key parameters in determining the size of the basic reproduction number is essential to find the best optimal strategy. Therefore, we study the normalized sensitivity analysis of the basic reproduction number using the following formula [30]: where ρ is any key parameter in malaria model in system (7). In our paper, we are only interested in the following parameters: β h , β v , α, u 1 , u 2 , η, δ, p, κ, γ, ξ, and a. Furthermore, we do not show the partial derivative of these parameters since it has a long expressions. Using parameter values as in Table 1, u 1 = 0:2, u 2 = 0, and p = 0:8; the normalized sensitivity of R 0 is given in Table 2 and visualized in Figure 2.
The normalized indices in Table 2 are a nondimensional value, which present the percentage change of R 0 for each increase value of parameter ρ for 1%. For an example, since Γ R 0 p = −0:6477, then increasing probability of individuals in T to succeed in treatment for 10% will reduce R 0 for 6.477%. On the other hand, since Γ R 0 β h = 0:5, then increasing β h for 10% will increase R 0 for 5%. From Figure 2, we can see that β h , β v , α, η, and δ are proportional to R 0 . Increasing these mentioned values will increase R 0 . On the other hand, parameters p, γ, u 1 , and κ are inversely proportional to R 0 . Therefore, increasing the value of p, γ, u 1 , and κ will reduce R 0 . In addition, we can see that fumigation ðu 2 Þ, rate of loss of immunity ðξÞ, and saturated parameter ðaÞ do not affect R 0 . Figure 2 shows the most to the less influential parameter to R 0 in descending order, from left to the right. Figure 3 shows the level set of R 0 with respect to u 1 , α, and p. From Figure 3(a), we can see that increasing the value of p reduces R 0 . It means that more people succeed due to treatment with tafenoquine; then, the possibility to achieve malaria-free equilibrium is bigger. Same interpretation to the rate of treatment u 1 . We can see that more intense intervention of tafenoquine will reduce R 0 . In addition, we can see clearly that better quality of tafenoquine will reduce the burden of intervention in providing tafenoquine treatment to achieve malaria-free conditions. The effect of vector-bias on the success of tafenoquine intervention to reduce R 0 can be seen in Figure 3(b). We can see that more bias the mosquito to be more preferring infected human will increase the R 0 , which makes the intervention of tafenoquine should be given more intense to reduce the value of R 0 .

The Malaria-Endemic Equilibrium
4.1. Existence of Malaria-Endemic Equilibrium. The malariaendemic equilibrium of system (7) is given by where with Note that I † is taken from the positive root of the following polynomial: where Ω is the set of parameter in system (7), and while k 5 , k 4 , k 3 , k 2 , and k 1 have a complex form to be written in this article. It can be seen that whenever I † > 0, then E † , T † , R † , and U † are also positive. On the other hand, S † is always positive since N h ≤ ðΛ h /μ h Þ (see the proof of Theorem 1). On the other hand, since Computational and Mathematical Methods in Medicine and N v ≤ ðΛ v /μ v Þ (see the proof of Theorem 1), then we can guarantee that W † is also positive. From the expression of polynomial in (19), k 6 is always negative since α > 1, k 0 > 0 ⇐ R 0 > 1, while another coefficient is difficult to be determined, whether it was positive or negative. Hence, using the Descartes rules of sign [46], there exists at least one positive root of polynomial (19) whenever R 0 > 1. According to this result and the expression of MEE, we have the following result. Theorem 3. System (7) has at least one malaria-endemic equilibrium point if R 0 > 1. (19) is a six-degree polynomial, it is possible that system (7) have more than one malariaendemic equilibrium point. We use Descartes rules of sign [46] to analyze the maximum possibility of positive root of polynomial in (19). The result is given in Table 3 for the case when R 0 > 1, and Table 4 for the case when R 0 < 1.

Since polynomial in
From Table 2, we can confirm the result in Theorem 3 that we always have at least one malaria-endemic equilibrium when R 0 > 1. If R 0 > 1, then we always have an odd number possibility of the positive root of polynomial (19), i.e., 1, 3, or 5 positive roots. On the other hand, malaria-endemic equilibrium is possible to vanish only when R 0 < 1. However, we still possible to have 2, 4, or 6 positive roots of polynomial (19) when R 0 < 1.

Bifurcation Analysis.
In this section, we perform the bifurcation analysis of our proposed malaria model in system (7). To do this analysis, we use the well-known Castillo-Song bifurcation theorem [47] (please see [48][49][50][51] for more examples on the use of this theorem on epidemiological models). First, for numerical calculation purposes, let us redefine our proposed system (7) as follows: x 6 , where x i for i = 1, 2, ⋯7 present S, E, T, I, R, U, and W, respectively, Next, we determine our bifurcation parameter to replace R 0 . By solving R 0 = 1 with respect to β h , we obtain the bifurcation parameter, namely, β h = β * , as follows: The linearization of MFE of system (22) at β h = β * is given by with
For the bifurcation indicators, we calculate A for system (22) which is expressed by We can confirm that A is always negative (please see the expression of A in Appendix D). Meanwhile, B is given by Since all parameters are positive, and ð1 − pÞ > 0, then we have that B > 0. According to Castillo-Chavez and Song theorem [47], since the quantity of A is negative and B is positive, then system (22)) indicates a forward bifurcation at R 0 = 1. We state the result in the following theorem.

Numerical Experiments on Theorem 4.
In this section, we show the numerical interpretation of Theorem 4. The first numerical experiment is for the bifurcation diagram of system (7), which is given in Figure 4. We use parameter values as mentioned in Table 1, except that it states differently. With this set of parameter values, we have R 0 = 1 when β h = 0:0004079. For the case of a = 400, u 1 = 0:2, and u 2 = 0, the bifurcation diagram is shown in Figure 4(a). It can be seen that the forward bifurcation phenomenon appears, which indicates there always exists a unique endemic equilibrium point when R 0 > 1, and no endemic equilibrium when R 0 < 1. Furthermore, we can see that the malariaendemic equilibrium is always stable (solid red) when R 0 > 1. The autonomous simulation for various initial conditions is shown in Figure 5. We use Runge-Kutta adaptive step size method in MATLAB to run the autonomous simulation in this article [52] (please see [53] for further detail on the method and its algorithm). It can be seen that when R 0 = 0:8 < 1, then the solution from all different initial conditions tends to the malaria-free equilibrium point ( Figure 5). On the other hand, when R 0 > 1, then all trajectories tend to the malaria-endemic equilibrium ( Figure 6).
The autonomous simulation of system (7) when forward hysteresis (Figure 4(b)) appears is given in Figures 7 and 8. The numerical results is using the same parameter values as in Figure 4(b). We only conduct two cases for this scenario, namely, when R 0 > 1 but close to one in which only one stable malaria-endemic appears (Figure 7) and when two stable malaria-endemic equilibrium appears ( Figure 8) when R 0 > 1, but not too far from 1. In the first case, as shown in Figure 7, we can see that all trajectories from all different initial conditions tend to the same malaria-endemic equilibrium. However, when hysteresis starts to appear, which causes two stable malaria-endemic equilibrium, the solutions will tend to two different stable malaria-endemic equilibrium points, depending on their initial conditions. We can see that when the initial condition is close enough to the bigger malaria-endemic equilibrium (blue curve), then the solution tends to the bigger size of malaria-endemic equilibrium. The same thing happens when the initial value of infection is small enough, and then, the solution leads to the smallest stable malaria-endemic  11 Computational and Mathematical Methods in Medicine equilibrium. These simulation results indicate that fumigation may trigger the existence of multiple stable malariaendemic equilibrium for some value when R 0 > 1. Figure 9 confirms the statement. It can be seen that an increase in fumigation rate increases the interval when multiple stable malaria-endemic equilibrium appears.

Autonomous Simulation
From the previous mathematical analysis, we found that our proposed malaria model always exhibits a forward bifurcation at R 0 = 1. These results indicate that the basic reproduction number becomes the only endemic indicator on

12
Computational and Mathematical Methods in Medicine our proposed model. However, our model may show a multiple stable endemic equilibrium when R 0 > 1. This phenomenon is called a forward hysteresis [54]. We found that this phenomenon was affected by the intensity of fumigation ðu 2 Þ and the level of population awareness ðaÞ. Furthermore, our sensitivity analysis indicates how important is the intervention of tafenoquine to prevent the occurrence of relapse and fumigation to control the number of Anopheles mosquitoes in the environment. To visualize our mentioned results, we perform several numerical simulations on our autonomous simulations for several scenarios.

Effect of Vector-Bias.
In malaria transmission, vectorbias has an important role in determining the endemic condition of the population [29]. The larger the vector-bias values, the more mosquito attracted to hunt infected humans for their meal. Figure 10 depicts the dynamic of the solution of our malaria model in (7) for several values of vector-bias parameter. We use the same parameter values as in Table 2, except u 1 = 0:1, u 2 = 0:05, p = 0:8, and varying α from 1 to 5. With these parameters, R 0 is always larger than 1, which makes the solution of system (7) tends to the malaria-endemic equilibrium. We can see that an increased value of the vector-bias parameter at the malariaendemic equilibrium situation will increase the total population in infected humans but reduce the size of the infected mosquito population. This means that the more Anopheles mosquitoes attracted to bite infected humans than healthy humans can negatively impact the human population, where the endemic size can increase. Therefore, efforts to control the mosquito population are essential in this situation.

Effect of Fumigation Saturation
Parameter. The first autonomous simulation was conducted to show the impact of the fumigation saturation parameter a. As we mentioned before, a smaller value a indicates a more prepared community to the increasing number of infected individuals. From the expression of R 0 in (15), it can be seen that a does not 13 Computational and Mathematical Methods in Medicine appear in R 0 . Hence, we conclude that a does not impact the size of R 0 . However, as we have shown in Figure 11, a smaller value of a reduces the size of total infected humans and mosquitoes in the malaria-endemic equilibrium point. Therefore, it can be concluded that although the level of community readiness to carry out fumigation does not affect the final state of population (endemic or not), it is clear that the higher the community readiness (the smaller the value of a), then the smaller the total size of the infected population in malaria-endemic equilibrium.

Effect of Different Fumigation Strategy.
As we have mentioned in sensitivity analysis on R 0 , we find that fumigation does not affect the size of R 0 , but it can reduce the size of malaria-endemic equilibrium when fumigation intervention increases, as shown in Figure 12. Now, we conduct our simulation with three different scenarios, based on the measured fumigation control depending on the implementation time. In the 1 st and the 2 nd scenarios, we choose u 2 to be changed depending on the time interval, using the following step function: while the 3 rd scenario when u 2 = 0:3 for all time t ∈ ½0, 10000. The result is given in Figure 13. We can see from Figure 13 that when there exist two stable malaria-endemic equilibrium points, then proper fumigation intervention may change the direction of stability of the system, which in our numerical experiment is from the large endemic size into small endemic size. When the improvement of fumigation is given too early (2 nd scenario), then after the fumigation intervention loosened back into 0:3, then the dynamic of total infected human goes back to the large endemic size. On the other hand, when the intervention is given several times after the first implementation (1 st scenario), then the dynamic of total infected humans is continuously going to the small endemic equilibrium. Based on this, it is necessary 15 Computational and Mathematical Methods in Medicine to consider the time for implementing an appropriate increase in fumigation intervention so that the solution dynamics can be directed to a smaller endemic point if the bistability phenomenon appears.

Conclusions
Malaria has long been a critical health problem in various parts of the world. Every year, hundreds of millions of people are at risk of becoming infected with malaria, with the majority of cases occurring in Africa. The disease is spread due to the bite of a female Anopheles mosquito and is caused by five different types of Plasmodium. Different types of Plasmodium that infect give different symptoms/serious illness that appears in patients with malaria. Various interventions have been and are being researched, such as vaccination, treatment, vector control with fumigation, and use of insecticide-treated bed net. The high number of cases in various parts of the world until now indicates that our understanding of malaria is still not sufficient to help us optimally control the spread of malaria.
In this research, we introduce a new malaria model that considers two important factors: the use of a new treatment (tafenoquine) to prevent relapse and a saturated fumigation function. The fundamental properties, the existence and stability criteria of the equilibrium points, and how they relate to the basic reproduction number are analyzed in detail. We use Descartes's rule of signs to show a possible number of malaria-endemic equilibrium points when the basic reproduction is less or larger than one. We find that it is possible to have more than one endemic equilibrium when the basic reproduction number is larger than one. Our bifurcation analysis shows how our model consistently exhibits a forward bifurcation at the basic reproduction number equal to one. However, our numerical simulations show forward bifurcation phenomena with hysteresis. This phenomenon results in the emergence of three malaria-endemic equilibrium for a basic reproduction number larger than one.
Our sensitivity analysis shows that tafenoquine has a big potential to control the spread of malaria by preventing the possibility of exposed individuals from relapsing. Furthermore, we also find that although fumigation does not affect the basic reproduction number, it can reduce the number of infected individuals at malaria-endemic equilibrium. Furthermore, a numerical investigation on implementing a high intensity of fumigation in a short time intervention interval may lead to a final switching condition if the forward with hysteresis phenomena appears. We find that when fumigation is implemented in a proper time interval, the direction of endemic equilibrium can be "kicked down" into the smaller size of malaria-endemic equilibrium, which is easier to control with other intervention strategies. We hope that the results of our research in this article can provide another perspective in evaluating the possibility of implementing tafenoquine and fumigation in the field.