Ecoepidemiological Model and Optimal Control Analysis of Tomato Yellow Leaf Curl Virus Disease in Tomato Plant

. The purpose of this study is to analyze the impact of control strategies, namely, insecticide spray, roguing of a diseased tomato plant, and protective netting to protect tomato plant from tomato yellow leaf curl virus disease (TYLCVD). For this, we formulate and analyze a simple deterministic model for the transmission dynamics of TYLCVD that incorporates these control strategies. We initially take into account the constant control case, we calculate the basic reproduction number, and we investigate the existence and stability of the disease-free and endemic equilibria. We use the Kamgang-Sallet stability to ensure that the disease-free equilibrium point is globally asymptotically stable when the reproduction number R 0 is less than one. This indicates that TYLCVD dies out independent of the initial size of the tomato population. For R 0 < 1 , the disease-free equilibrium becomes unstable, and the endemic equilibrium is globally asymptotically stable. This indicates that TYLCVD propagates. In the nonconstant control case, we use Pontryagin ’ s maximum principle to derive the necessary conditions for the optimal control of the disease. Our ﬁ ndings show that all the combined e ﬀ orts of two out of three strategies can signi ﬁ cantly reduce the power of infectivity of the disease except the combination of the use of insecticide spray and rouging infected tomato plants. Besides our numerical simulations show, the implementation of the combination of roguing diseased plants and protective netting has a similar e ﬀ ect in minimizing or eliminating TYLCV as the use of all strategies. Hence, as resources are always in scarce, we recommend policymakers to adapt the combination of the use of roguing diseased tomato plants and protective netting to eradicate the disease.


Introduction
Tomato (Solanum lycopersicum L.) is one of the most popular and widely grown vegetables in the world.However, it is highly destructed by tomato yellow leaf curl virus (TYLCV) disease [1,2].This disease is mainly transmitted by an insect vector called whitefly Bemisia tabaci (B.tabaci) of biotype B (Gennadius) (Hemiptera: Aleyrodidae) in a circulative and persistent manner [3].The vector damages tomato plant directly by feeding on phloem, excreting honeydew, and causing phytotoxic disorders [4].With increased populations, they secrete large quantities of honeydew.The vector favors the growth of sooty mold on leaf surfaces and reduces the photosynthetic efficiency of the plants [5].The honeydew also contaminates the marketable part of the plant, reducing its market value.Additionally, in severe infestations, the leaves turn yellow and are dropped off from its leaves [3].
The first report on TYLCV infection in tomato was from Israel and other countries in the Middle East in the 1930s.And since then, the virus has further emerged [5].In Africa, The TYLC disease was first seen in Sudan in 1965; but the causal agent was identified as TYLCV in 1997 [6,7].In Ethiopia, at Melka-Werner, about 90% of tomato plant showed leaf curl virus symptoms with reduced size, suspected to be caused by TYLCV [8].However, only recently occurrence of Begomovirus associated with the TYLC disease was reported for the first time from Melkasa [9].
The whitefly vector feeds on an infected host plant and acquires the virus.Viral transmission can occur within hours and may continue for the life span of the vector.Acquisition and transmission thresholds were found to be between 15 and 30 min, and a single B. tabaci whitefly can accumulate 600 million TYLCV genomes [3].From the site of inoculation, viral DNA is first mainly transported to strongly sink organs, such as root and shoot apices, flowers, and fruit, and then moves to leaves [10].Similarly to other plant viruses, TYLCV moves in the existing host transport routes such as plasmodesmata and phloem, along with carbohydrates [11].About 11-13 days after inoculation, maximum amount of viral DNA and capsid or coat protein (CP) accumulates in the youngest tissues of shoots and roots.Four to seven days later, unique symptoms appear for the first time [10].As the systemic infection proceeds in the growing plant, the virus is accumulated in the strongest sink tissues.The level of viral double-stranded DNA (dsDNA) and newly generated single-stranded DNA (ssDNA), as well as CP, further increases in young organs up to several weeks postinoculation.The infection then gradually spreads to older organs of the host and remains strictly confined to the vascular system [11].
Mathematical modeling has been playing an important role in better understanding the epidemiological patterns.It provides deeper insight into the underlying mechanisms for the spread of emerging and reemerging infectious diseases and the suggested effective control strategies [12].Holt et al. [13] on their paper entitled "An epidemiological model incorporating vector population dynamics applied to African Cassava Mosaic Virus Disease" illustrated the African cassava mosaic virus occurrence in cassava which is transmitted by a cassava-specific whitefly strain, which was later seen in Uganda.The virus also propagates routinely from stem cuttings.The use of uninfected cutting tools and roguing of infected plants are among the control alternatives.Utilization of uninfected cutting tools would be more effective if the infected cuttings transmit disease, whereas roguing would be more important in a largely vector-driven epidemic.
In a later paper, Holt et al. [14] developed an epidemiological model, a susceptible-exposed-infective-(SEI-) type epidemic model for the host plant and a susceptible-infective-(SI-) type for the insect vector population.It represents the incidence of TYLCV in tomato plant, mainly relying on the immigration of vectors from alternative hosts which act as a reservoir of both the virus and vector.This is because, unlike the cassava, the tomato was only an occasional host for this whitefly and is spilledover from other perennials and weedy plants driven by vector and virus dynamics.They considered different strategies to reduce the spread of TYLCV and studied the sensitivity analysis of their results to the parameters to explore different disease management options.In this context, the authors asked, "What the best method to control the disease is?"Because most of the vector lifespan occurs on other hosts.The authors adapted a model framework [15] to explain the transmission process of the disease.Since the tomato crop was a sink for whiteflies and TLCV, interventions that reduce vector immigration and their survival were predicted to be most effective.
Alemneh et al. [16] proposed and examined an ecoepidemiological deterministic model for the transmission dynamics of maize streak virus (MSV) disease in the maize plant.Their model depicted that increased parameters namely the infection and predation rates made an increment of the basic reproduction number that leads to the increment of the number of infected maize population.Hence, the authors suggested that to intervene in MSV disease, endeavors should be exerted to reduce the contact of infected maize and susceptible leafhopper.In addition to this, MSV-infected maize should be treated using insecticide chemicals.This enabled us to bring down the infection rate of leafhoppers, and it should be administered before the reach of the leafhopper or uprooting it.Moreover, infected maize should be burnt from the field.
Nowadays, many sophisticated metaheuristics have been initiated to solve the most complex problems by transforming them into problems of optimization.For instance, Farman et al. [17] used an evolutionary Pade approximation (EPA) scheme for the treatment of nonlinear epidemiological hepatitis-B model, instead of numerical methods.Thus, they showed that as compared to a nonstandard finite difference discretization (NSFD) numerical scheme, EPA scheme is more reliable and significant in approximations of state variables that are highly accurate to the governing equations.Naik et al. [18] proposed and examined a nonlinear fractional order SEIR model for the transmission of HIV epidemics.For such nonlinear fractional-order model, they inspired the application of the homotopy analysis method (HAM) to solve highly nonlinear fractional-order problems, describing biological phenomena.They also gave a research direction to use another approximate solution technique such as fractional-order derivatives with Alangana-Gomez, Atangana-Baleanu, Caputo-Fabrizio, and fractal-fractional, to get a nonlinear fractional-order SEIR epidemic model.
Optimal control theory has found wide-range of the applications in biological and ecological problems [19].Particularly, there have been various studies of epidemiological models where optimal control methods have been applied.Berhe et al. [20] formulated a deterministic model to study the effects of implementing continuous controls on the dysentery epidemic model.It examined the cost-effectiveness of the optimal control measures of the disease.They took three control parameters, namely, treatment, sanitation, and educational campaign as a preventive strategy.As a result, they found that the disease can probably be eradicated by implementing continuous controls in a short period of time.However, utilizing a combination of sanitation of the environment and educational campaign was found to be the most cost-effective control strategy.Okosun and Makinde [21] derived and analyzed a deterministic model for the transmission of childhood disease, incorporating optimal control parameters and investigated the cost-effectiveness of the controls, to identify the most effective strategy.They considered control parameters, such as improvement of hygiene due to health educational campaign, improvement of treatment of the infected children, and reduction in the loss of disease immunity due to the improvement of vaccination and treatment efficacy.Thus, utilization of these control strategies has declined the disease from the community.
Similar results are also obtained if educational campaigns as preventive measure and treatment of infected children were used.However, as resources are scarced, the authors proposed that policymakers likely should focus on optimal provision of prevention and treatment being costeffective.Bokil et al. [22] investigated and analyzed optimal control of a vectored plant disease model for crops with continuous replanting.They considered two plant-vector-virus models which take into account frequency replanting and abundance replanting strategies to study the African cassava mosaic virus.They compared the two models with respect to replanting strategies through a combination of mathematical analysis, parameter sensitivity, and optimal control of the disease dynamism.They used optimal control theory to investigate the effects of roguing and insecticide to maximize the healthy plants to be harvested.The simulation results of their models suggested that various optimal control strategies were suitable for the different replanting practices.Hugo et al. [23] studied optimal control and cost-effectiveness analysis of the TYLCV disease epidemic model.Their model was an extended work of [14].The authors incorporated the time-dependent control to the tomato plants and vector populations in analyzing the cost-effectiveness of the control strategies.Thus, they suggested that the use of the combined protective netting and removal of the infected plant is the cost-effective optimal control strategy.It was sufficient to combat the epidemic of the tomato disease with limited resources.A preprint of our paper has previously been published [24].
The major advantage of these early models was to provide a suitable control strategy through the transmission threshold criterion, which is based on the reproductive capacity of the parasite, R 0 .To wind up, this paper is focused on optimal control strategies analysis of the tomato yellow leaf curl virus disease model, which is adopted from the African cassava mosaic virus disease's model [13] with the inclusion of exposed class into tomato plant population and incorporate three-time dependent controls, representing the interventions.
The organization of the paper is as follows.In Section 2, we presented a model consisting of ordinary differential equations that describe the transmission dynamics of tomato yellow leaf curl virus disease and the underlying assumptions.Section 3 is devoted to the stability analysis of the model.Section 4 is contained by numerical simulations.Our conclusions are given in Section 5.

Mathematical Model
The model that is considered here is a small modification of the model for plant-virus transmission considered in [13].It is a standard model of SEI type for tomato plants and SI for whiteflies Bemisia tabaci B-type insect vector.
The model subdivides the total tomato plant population into the following subclasses: healthy or susceptible (S p ), latently infected (E p ) and infective (I p ), K is the carrying capacity of the tomato farm.Thus, the total population size of tomato plant is N p = S p + E p + I p .The total insect vector (whitefly Bemisia tabaci B-type) population is subdivided in respect to tomato plant into susceptible (virus-free) (S v ) and infective vector (I v ).The latent period in the vector between the acquisition of the TYLCV and the ability to transmit the virus is roughly 30 min [25].Thus, it is assumed to be negligible, i.e., no latent subclass was defined for the whitefly vector in our model.The whiteflies remain infective for their lifetime.Hence, the total population size of the vector (whitefly) is The net replanting rate of tomato plants is rS p ð1 − S p + E p + I p /KÞ where r is the rate replanting healthy tomato.This is because, the model assumed that healthy tomatoes responded proportionately inverse to the extended intraspecific pressure for the healthy, exposed, and infected tomato plants.Replanting tomato is restricted by maximum tomato plant availability and the harvesting of healthy, exposed, and infected tomato plant reflects the continual turnover of the tomato plant population.The tomato fruit is either harvested or removed at a rate of g or move to the exposed stage by inoculating through contact with infective whiteflies at a rate of β p .Moreover, all stages of tomatoes are assumed to be harvested or removed at the same constant rate.Thus, it is assumed that the force of infection at time t is given by β p I v ðtÞS p ðtÞ. Latently infected tomato plants propagate to the infectious stage at a rate of a, corresponding to a mean latent period in a tomato plant population of 1/a.Here, it is assumed that the infective tomatoes remain infected forever.Loss of tomatoes due to nature-and disease-related reduction is rated as b.
The population of whiteflies is assumed to be generated at a rate of cðS v + I v Þð1 − S v + I v /mðS p + I p ÞÞ, where c is the vector birth rate and m is the rate of vector maximum abundance.Besides, we assumed that the vector does not immigrate to other hosts in the case of low tomato plant abundance.The vector dies either from natural causes at a rate of e or propagates to the infective class by acquiring tomato yellow leaf curl virus through contacts with infected tomatoes at a rate of β v I p ðtÞS v ðtÞ, where β v is the rate of virus acquisition by a susceptible vector (whitefly) during a single visit to an infectious tomato plant.The protections target mainly the following: (i) insecticide on tomato plant population: minimizes the inoculation efficiency of the vector, i.e., it reduces to β p ; (ii) roguing diseased tomato plants: good practice for reducing the source of primary infection, i.e., reduced to g; (iii) protective netting, prevent the entry of whitefly vectors, B. tabaci, into tomato plots, i.e., reduced to β p and β v .This is because the only way of controlling TYLCV is by controlling the vector [5].Thus, the efforts made on these three intervention mechanisms enable to control the tomato infections due to TYLCV.
Suppose that the control function like u 1 ðtÞ represents insecticide spray with efficacy q, u 2 ðtÞ represents roguing diseased tomato plants, and u 3 ðtÞ represents protective netting at any time t.Besides, controlling u 1 , u 2 , and u 3 are assumed to be bounded and integrable functions.
It is further assumed that the transmission of the virus by the whitefly vectors is by circulative and persistent mode 3 Journal of Applied Mathematics [5].Moreover, it is also assumed that infective whiteflies stay infective for life.
Based on the above assumptions, the following vectorplant dynamical system is formulated: with initial conditions

Stability Analysis of Free Disease Equilibrium Point
In this section, we study the positivity and boundedness of the solution and the stability of the disease-free equilibrium of model system 1.As the variables in model system 1 represent tomato and whitefly population densities, positivity indicates survival of the population, and boundedness may be interpreted as a natural restriction to growth as a result of limited resources.The stability of the disease-free equilibrium point tells us TYLCV disease dies out from the population.

Positivity and Boundedness of Solution.
For the TYLCV transmission model ( 1) to be biologically meaningful, it is important to prove that all solutions with nonnegative initial data will remain nonnegative for all the time as is presented in [26].
Theorem 1.Let S p ðtÞ ≥ 0, E p ðtÞ ≥ 0, I p ðtÞ ≥ 0, S v ðtÞ ≥ 0 and I v ðtÞ ≥ 0. The solutions S p , E p , I p , S v , I v of the system of differential Equation ( 1) are positive for all t ≥ 0. Besides, the region Ω is positively invariant and all solutions starting in Ω approach, enter, or stay in Ω.
Proof.By adding the first four equations in the system (1), we found that the rate at which the total population of tomato plant changes is given by the following: Since N p = S p + E p + I p and S p ≤ N p , we have Similarly, by adding the last three equations of system (1), we obtain that the rate at which the total population of whitefly vector changes is given by the following: Since N v = S v + I v and K > S p + I p , then Equation ( 6) can be written as If we factorize Equation ( 7), then we can obtain The region This implies that all solutions of tomato plants population are only confined in the feasible region Ω p , and all solutions of the whiteflies population are confined in Ω v .
Therefore, the biological feasibility of system (1) is studied in the following region: Thus, Ω is positively invariant.This means that solutions of the model system with positive initial data remain positive for all time t ≥ 0 and are bound in the region Ω.Therefore, the model is mathematically and epidemiologically well-posed.

Analysis of the Model with Constant Controls.
In this section, it is assumed that the control parameters are constant and determine the basic reproductive number, the steady states, and their stability.

Local and Global Stability.
The disease-free equilibrium of the tomato yellow leaf curl virus disease model (1) exists and is given by The basic reproduction number, R 0 , is calculated by using the next-generation matrix [27].To obtain R 0 for model (1), let the vector of the disease states Then, the model ( 1) can be written as where Calculate the Jacobian matrix (F and V) of F and V by derivating with respect to the infected classes ðE p , I p , I v Þ at the disease-free equilibrium point E 0 .This gives

So that
The basic reproduction number, Theorem 2 below follows from Theorem 2 of [28].
The epidemiological implication of Theorem 2 is that the transmission of TYLCV can be controlled by having R 0 < 1 if the initial total numbers of the subpopulations involved in Equation ( 1) are in the basin of attraction of E 0 .To ensure that eliminating the disease is independent of the initial size of the subpopulation, the disease-free equilibrium must be globally asymptotically stable when R 0 < 1.This is what we are checking next.
Proof.To prove the theorem, we use the Kamgang-Sallet stability theorem in [29].
. Then, the system (1) can be written as where We show that the five sufficient conditions of the Kamgang-Sallet theorem are satisfied as follows: (1) The system (1) is a dynamical system on Ω.This is proved in Theorem 1.
(2) The equilibrium X * 1 is GAS for the subsystem _ . This is obvious from the structure of the involved matrix.
(3) As can be seen from the elements, the matrix A 2 ðXÞ is Metzler (i.e., all the off-diagonal elements are nonnegative) and irreducible for any given X ∈ Ω.
(4) There exists an upper-bound matrix A 2 for the set More precisely, is an upper bound of M. ( Thus, all eigenvalues of A are negative for R 0 < 1 in Equation (19).
Hence, by the Kamgang-Sallet stability theorem, the disease-free equilibrium is globally asymptotically stable for R 0 < 1.
For any initial data, Theorem 3 implies that any solution of the system (1) converges to the DFE when R 0 < 1.In addition, the theorem implies that the model is without backward bifurcation for R 0 < 1.In this case, the classical approach making R 0 < 1 to eliminate TYLCV disease from the farm is sufficient.
Next, we want to check that the system (1) has at least one endemic equilibrium (EE) point for R 0 > 1.Let be an EE of system (1).By setting the right-hand side of (1) equal to zero, we obtain where It can be shown that the equilibria of the model satisfy the following polynomial where Journal of Applied Mathematics Theorem 4. When R 0 > 1, the model ( 1) has at least one endemic equilibrium, which is locally asymptotically stable for R 0 close to one.
The stability of the EE is guaranteed by using the center manifold theorem in [30].

Sensitivity Analysis of Model
Parameters.Sensitivity analysis assists to build confidence in the model by studying the uncertainty associated with parameters in the model.This is because many parameters in the system dynamic models characterize quantities that are very difficult or even impossible to measure accurately in the real world.It helps to comprehend the dynamics of the system under study.In general, sensitivity analysis is carried out to establish which input parameters contribute most to output variability [31].Now let us carry out the sensitivity analysis in order to identify the parameters that have a high impact on the basic reproductive number (R 0 ).Definition 5.The normalized sensitivity index of a variable, R 0 , that depends differentiably on a parameter, p, is defined as

Sensitivity Indices of Basic Reproductive Number.
Here, the sensitivity of R 0 to every parameter of the model is derived.Hence, the sensitivity index of R 0 with respect to K and is equal to 1.It is equal to 0.5 with respect to m, β p , and β v .The rest of the parameters are the following: Since most of the expressions for sensitivity indices are complex, the sensitivity indices are evaluated at the baseline parameter values given in Table 1.The sensitivity index of R 0 with respect to r, for example, is as follows: The detail sensitivity indices of R 0 , resulting from the evaluation of the eight different parameters of the model, are shown below.
The parameters are arranged from most sensitive to least.The most sensitive parameters are the rate of replanting of healthy tomato, the rate at which tomato fruits are either harvested or removed, the carrying capacity of tomato farm, the rate of vectors' maximum abundance, the rate of inoculation of the healthy tomato plant, and the rate of virus acquisition by susceptible vectors.The least sensitive parameter is the loss rate of tomatoes due to infection.The parameters that reduce R 0 can be used as control parameters.

Interpretation of Sensitivity Indices.
The sensitivity indices of the basic reproductive number with respect to the main parameters are found in Table 2.The parameters with the most important that have positive indices are β p , β v , r, K, m, c, and a, and those with negative indices are g, e, and b.The results show that when the inoculation of the healthy tomato plant rate and virus acquisition rate by susceptible vectors increases, the basic reproduction number increases as a result of TYLCV disease propagated since the infected plants and infective vectors are both infectious.When the rate of replanting of a healthy tomato and whitefly vector birth rate is increased, more tomato plants and whitefly vectors are exposed to TYLCV, and these increase their probability of catching the virus and contribute to the spread of the disease.The propagation rate from the exposed class to the infectious class increases the number of infectious tomato plants and vectors.These contribute to the disease spread.On the other hand, tomato fruit harvested or removed and death rate of vectors are reduced to the size of tomato plant population and whitefly vector population.

Analysis of Optimal Control
Let xðtÞ represent the tomato plant population to be protected via insecticides, cultural techniques, and virusresistant cultivars; insecticides reduce the number and movement of the whitefly vector.Cultural techniques (roguing, avoidance, crop residue disposal) reduce the amount of secondary spread within a tomatoes' field.Virus-resistant cultivars reduce the loosed ones, due to infections by 7 Journal of Applied Mathematics TYLCV [32].Thus, the aim of this is to minimize the multiple objective cost functional J, considering the costs of control methods of exposed and infected tomato plants.
If q is the efficacy of insecticide spray, c 1 and c 2 are cost factors due to the size of the infectious tomato plants and whiteflies population; p 1 , p 2 , and p 3 represent the weight attached on the cost control methods, and then the cost rate at which TYLCV disease is controlled at any time t can be given by the following: Since the implementation of any intervention has decreased costs, it is customary to take a nonlinear cost function.Hence, the simplest nonlinear function (the quadratic) in modeling the cost of the interventions is used.Therefore, an optimal control u * = ðu * 1 , u * 2 , u * 3 Þ searches the following: where is the set of admissible controls.
To sum up, the optimal control problem has the following form: 4.1.Pontryagin's Maximum Principle.Since our model has no terminal constraints, it is a normal optimal control problem, and hence, the Hamiltonian takes the following form: where F i ðx, u, tÞ is the right-hand side of the differential equations of the i th state variable.By using Pontryagin's Maximum Principle and the existence of results obtained for optimal control, we obtain: There exists an optimal control u * 1 , u * 2 , u * 3 and corresponding solution, Furthermore, there exist adjoint functions λ i , i = 1, 2, 3, 4, 5 such that with transversality condition By using Pontryagin's maximum principle and the existence result for the optimal control (Makinde et al. [33]), we arrive at the following theorem: 8

Journal of Applied Mathematics
The optimal control, u Proof (Fleming et al. [34]).Provide the existence of optimal control due to the convexity of integrand with respect to (u 1 , u 2 , u 3 ), a priori boundedness of the state solutions, and the Lipschitz property of the state system with respect to the state variables.Employing Pontryagin's maximum principle, we have where λ i , i = 1, 2, 3, 4, 5. Therefore, the adjoint function with each state variable is calculated as follows: Applying Pontryagin's maximum principle, the Hamiltonian equation can be written as follows: Considering the existence of adjoint functions λ i , i = 1, 2, 3, 4, 5 satisfying with transversality condition λ i ðTÞ = 0, i = 1, 2, 3, 4, 5 for the control set u i , hence computed at the optimal control pair and respective corresponding states, which leads to the stated adjoint systems (37) and (38), [35].By taking into account the optimality conditions, and to determine the values for u * 1 , u 2 , u 3 , subject to the constraints, the characterizations (39) can be obtained,

Numerical Simulation
Simulations are carried out to determine the behavior of the system (1).For this purpose, the parameter values listed in Table 1 were used.Most parameter values were taken directly from [14], and the rest were estimated from the data found in [5,36].The estimated parameters were calculated as follows.Since the life span of whiteflies is 20-50 days, the death rate was calculated as per the inverse of the average life span, i.e., e = 1/ðð20 + 50Þ/2Þ = 1/35 = 0:0286day −1 .Depending on cultivars, tomato fruits could be made ready to harvest at about 75 to 90 days after transplanting [36], and tomato fruits' are harvested or removed at the rate of g, which can be calculated as g = 1/ð75 + 90Þ/2 = 0:0121day −1 .Since the exposed period of the tomato plant is 10-14 days [5], the propagation rate from exposed plant to infected plant is calculated as per a = 1/ð10 + 14Þ/2 = 0:0833day −1 .
The main objective of this study is to examine the impact of control strategies such as insecticides, cultural practices, and virus-resistant cultivars on the transmission of TYLCV.In order to support the analytical results, graphical representations of various strategies are visualized to determine their impact whenever the controls are applied to the system [14].

Optimal Control
Effect on the Model.Now it is time to look at the effect of different optimal control strategies on the propagation of TYLCV disease.It is well known that there is no single management option to control the disease.This makes the management of TYLCV challenging and costly [5].A combination of management options is necessary to successfully manage the disease and limit the losses of tomato fruits' production [5].For instance, a combination of cultural and chemical uses is required [32].Therefore, the following optimal control strategies on the propagation of TYLCV disease in the tomato plant population are numerically investigated below.The optimal control is obtained by solving the optimality system (38), (39), and (43).An iterative scheme is used for solving the optimality system.Thus, the state equations are solved with an assumption for the controls over the simulated time using the fourth-order Runge-Kutta scheme.Because of the transversality conditions (43), the adjoint equations are solved by a backward fourth-order Runge-Kutta scheme using the current iteration solutions of the state equation.Then, the controls are updated by using a convex combination of the previous controls and the value taken from the characterizations (46).Generally, it can be written as where k is the current iteration and 0 < α < 1 [19].This process is repeated, and the iterations are also stopped if the values of the unknown variables in the previous iterations are very close in contrast to the iterations [19].
We assume that p 1 > p 3 > p 2 .This assumption is based on the facts that the cost associated with u 1 , u 2 , and u 3 which is the cost of spraying insecticides is applied five times per season; and the use of roguing diseased tomato plants mainly labors cost [32] and the cost associated with the protective netting.Thus, c 1 = 1, c 2 = 0:025, p 1 = 0:006, p 2 = 0:003, and p 3 = 0:005 are chosen and used as parameter values in Table 1.The initial state variables are chosen as S p ð0Þ = 384, E p ð0Þ = 17, I p ð0Þ = 69, S v ð0Þ = 768, and I v ð0Þ = 138.Figure 1 shows that the number of healthy tomato plants is increased gradually while the infected tomato plant population decreased with time in the presence of the control.On the other hand, without control, the number of infected tomatoes escalates while the healthy tomato plant population is reduced.This is logical, because insecticide spray and virus-resistant cultivars help to reduce the incidence of TYLCV disease.
Figure 2 depicted that without control, both susceptible and infected insect vectors increase, but in the presence of control strategies, the insect vector population is decreased.This is because we assume that the tomato plant is the only host for the insect vector and thus suffers from lack of food.

10
Journal of Applied Mathematics This is justifiable because the virus-resistant cultivars resisted not to be infected, and insecticide spray is effective against the TYLCV infection.11 Journal of Applied Mathematics without control and the deescalation of infected tomato plants to some threshold with the use of control.This implies that a combination of the use of insecticide spray and roguing tomato plant control reduced the tomato yellow leaf curl virus disease to some threshold.
Figure 4 depicts that the number of susceptible and infected whitefly vector population is increased in the absence of control.However, according to Figures 4(a) and 4(b), both the susceptible and infective whitefly vectors increased.Even though the infected whitefly vectors increase slightly in the case of control, this may be attributed to roguing infected tomato plants, decreasing the amount of secondary spread of the disease within a farm when incidences of the disease are low.However, if there are higher rates of infection (>10%), roguing may not be successful.In this study, for the simulation purpose, the number    5(a) represented that without the healthy tomato plant, the population decreased with time.Figure 5(b) reveals that the infected tomato plant gradually increased without control.It is decreased with time, providing that the control is utilized.This implies that the use of insecticide spray against the spread of tomato yellow leaf     13 Journal of Applied Mathematics curl virus may be less efficient as it is used to hinder the foundation of the vector B. tabaci within the crop to protect plants from direct devastation.
According to Figure 6, the susceptible and infected whitefly vector populations escalate without control, whereas in the case of control, virus-transmitting whitefly also increased, but at a decreasing rate.This indicates that controlling virus-transmitting whitefly vectors using insecticide spray and roguing diseased tomato plant is hard [37].It is observed in Figure 7(a) that the control strategies resulted in an escalation of a healthy tomato plant population in the presence of control strategies.It dropped off to zero without control.Figure 7(b) portrays a significant decrease in the numbers of infected tomato plants in the case of control.It is dramatically increased in the absence of control.This implies that a combination of insecticide spray, rouging diseased tomato plant, and virus-resistant cultivars declined TYLCV disease.
Figures 8(a) and 8(b) explain that with the application of control, the susceptible and infected whitefly vectors are decreased significantly in the case of control usage, whereas it is increased in the case of no control.

Conclusion
In this paper, a simple deterministic model for the transmission of TYLCV disease that incorporates the strategies of the use of insecticide spray, roguing the diseased tomato plants, and virus-resistance cultivars is used.It performs optimal control analysis of the model.The basic reproduction number is calculated, and the existence of local and global stability of equilibria is analyzed.We initially take into account the constant control case, we calculate the basic reproduction number, and we investigate the existence and stability of the disease-free and endemic equilibria.We use the Kamgang-Sallet stability to ensure that the disease-free equilibrium point is globally asymptotically stable when the reproduction number R 0 is less than one.This indicates that TYLCVD dies out independent of the initial size of the tomato population.For R 0 < 1, the disease-free equilibrium becomes unstable, and the endemic equilibrium is globally asymptotically stable.This indicates that TYLCVD spreads.The sensitivity analysis of the basic reproduction number shows that TYLCV disease has a positive relationship with the rate of harvested or removed tomato fruits, inoculation of healthy tomato plants, and virus acquisition by susceptible vectors.Thus, these parameters are those that have to be targeted mostly by policymakers to combat the TYLCV disease.Hence, the optimal control analysis of the model is made, using Pontryagin's maximum principle.From the simulation results, it can be concluded that all the combined efforts of the two out of the three strategies are listed as insecticide spray, rouging diseased tomato plants, and virus resistance cultivars.These can significantly reduce the disease, except the strategic combination of the use of insecticide spray and rouging infected tomato plants.Since the strategic combination of the use of roguing diseased tomato plants and virus-resistance cultivars has a similar effect as the strategic combination of the use of all strategies, we recommend that policymakers ought to adopt the combination of the use of the combination of roguing diseased tomato plants and virus-resistance cultivars as a strategy.
(i) Strategy I: combination of the use of insecticide spray and virus-resistant cultivars (ii) Strategy II: combination of the use of roguing diseased tomato plants and virus-resistant cultivars (iii) Strategy III: combination of the use of insecticide spray and roguing diseased tomato plants (iv) Strategy IV: combination of the use of insecticide spray, roguing diseased tomato plants, and virusresistant cultivars

5. 1 . 1 .
Strategy I: Combination of the Use of Insecticide Spray and Virus-Resistant Cultivars.The insecticide spray control u 1 and virus-resistant cultivars control u 3 are used to optimize the objective function J while roguing diseased tomato plant control strategy u 2 is set to zero.

Figure 4 :Figure 3 :
Figure 4: The impact of roguing diseased tomato plants and virus-resistant cultivars on (a) susceptible whitefly vector and (b) infective whitefly vector.

Figure 6 :
Figure 6: The impact of insecticide spray and roguing diseased tomato plant on (a) susceptible whitefly vector and (b) infective whitefly vector.

Figure 5 :
Figure 5: The impact of insecticide spray and roguing diseased tomato plant on (a) healthy tomato plant and (b) infected tomato plant.

5. 1 . 4 .
Strategy IV: Combination of the Use of Insecticide Spray, Roguing Diseased Tomato Plants, and Virus-Resistant Cultivars.All the three controls u 1 , u 2 and u 3 are used to optimize the objective function J.

u 1 u 1 Figure 7 :Figure 8 :
Figure 7: The impact of insecticide spray, roguing diseased tomato plant, and virus-resistance cultivars on (a) healthy tomato plant and (b) infected tomato plant.

Table 1 :
Model parameters and values used in the simulation.