Preventions and Controls on Congenital Transmissions of Zika: Mathematical Analysis

Affiliation(s)

^{1}
Department of Mathematics, Gujarat University, Ahmedabad, Gujarat, India.

^{2}
L. D. College of Engineering, Ahmedabad, Gujarat, India.

ABSTRACT

Vector-borne diseases threat lives of millions of people in many countries of the world. Zika is one of the vector-borne diseases which initially spread by the bite of an infected Aedes species mosquito (*Ae. aegypti* and *Ae. albopictus*) and then it transmits vertically from a pregnant woman to her fetus or from an infected human to their sexual partners. The congenital transmission of Zika virus (ZIKV) results in new born with microcephaly and other neurological abnormalities. The control of infected mosquitos is the best efficient way to control spread of ZIKV. Spraying insecticide is the safest and easiest way to control mosquitos, but sometimes it is cost worthy for long period of spraying. Controlled prevention from the vector bites can also help to control disease spread. To control congenital transmission and sexual transmission of ZIKV, preventions should be taken to reduce/stop pregnancy rate and safe heterosexual transmission among adults. Also, there is no specific treatment available for Zika disease. Treatment is aimed at relieving symptoms with rest, fluids and medications. Controlled combinations of rest, fluids and medications will help to recover early. As costs are incorporated with spraying, preventions and treatment, our aim is to minimise the total cost associated by controlling spraying, preventions and treatment. To fulfil this purpose a mathematical model is developed with disease dynamics in nine compartments namely Susceptible human child, Susceptible human male, Susceptible human female, Infected human child, Infected human male, Infected human female, Recovered human, Susceptible vector and Infected vector including vertical transmission of Zika disease. Numerical simulations have been carried out to optimise controls, and basic reproduction number and stability are calculated.

Vector-borne diseases threat lives of millions of people in many countries of the world. Zika is one of the vector-borne diseases which initially spread by the bite of an infected Aedes species mosquito (

1. Introduction

Zika virus (ZIKV) is a Flavivirus. It is initially transmitted to humans by the bites of infected female mosquitoes from the Aedes genus. The pathogen responsible for spread of Zika disease is known as Zika virus (ZIKV).In past two years remarkable changes has been seen in the epidemiology of Zika virus (ZIKV). The transmission of ZIKV has been first reported from continental America and the Caribbean. Also, recent reports indicate an increase in detected cases of congenital malformations and neurological complications associated with ZIKV infection. To treat, prevent, or diagnose ZIKV infection there is no specific treatment, vaccine, or fast diagnostic test available at this time.

Dick et al. [1] , Rodrigez et al. [2] and Macnamara et al. [3] observed that in the Zika forest of Uganda the Zika virus was initially isolated from a rhesus monkey in 1947. In 1954 same virus was isolated from humans in Nigeria. There after only sporadic infected human cases were reported from Africa and Southeast Asia. Duffy et al. [4] observed that first largest outbreak of ZIKV infection was reported in 2007 on Yap Island of the North Pacific. Musso et al. [5] observed that, in French Polynesia during October 2013, 28000 ZIKV infection cases were reported. Pan American Health Organization [6] noted that the recent outbreak began in April 2015 from Brazil, has covered many southern, central American countries and the Caribbean with spread of disease, and more than 140,000 suspected and confirmed cases are reported by the end of February 2016.

Foy et al. [7] reported that infection of ZIKV spread from an infected male to a female during their sexual intercourse. Thereafter, in February 2016 cases of sexually transmitted ZIKV were reported from Dallas [8] County, U.S. and France by CDC [9] , Hills et al. [10] and Mansuy et al. [11] respectively. According to the Toronto Star [12] spread of Zika infection because of sexual activities has been reported in Argentina, Canada, Chile, France, Italy, New Zealand, Peru, Portugal, and the USA from 2015. Hills et al. [10] recorded that the disease has minor impact on sexual activity as ZIKV infections have mild symptoms for two to seven days. As ZIKV is sexually transmissible, CDC [9] issued related guidelines for preventions should be taken for safe sex during a Zika outbreak. Didier et al. [13] noted that ZIKV has been found in semen samples which point out the possibility of transmission of ZIKV through sexual activity.

Gatherer et al. [14] found that like other flaviviruses ZIKV could also be transmitted by blood transfusion.

Mlakar et al. [15] and Cauchemez et al. [16] indicated that ZIKV increases the chances of microcephaly in new-born babies of infected mothers. Cao-Lormeau et al. [17] derived that it causes Guillain-Barre syndrome (GBS) as well. According to Pan American Health Organization, from October 2015 to February 2016 approximately 6000 suspected cases with 139 deaths of microcephaly in new- borns due to ZIKV infection were reported in Brazil, and from December 2015 to February 2016, more than 200 cases and 118 cases with 5 deaths of GBS due to suspected ZIKV infection were recorded in Colombia and in El Salvador respectively. As per estimation of WHO [18] , based on the reported dengue data from 2015, about four million people in the Americas may be infected by ZIKV in 2016.Oliveira et al. [19] observed that in February 2015, abnormal brain development in the foetuses of pregnant and ZIKV infected women were found.

Dalia et al. [20] suggested for prevent pregnancy and effective strategies against the vector will control disease spread. Ebenezer et al. [21] suggested that to control disease the best strategy is to combine all preventive, treatment and Insecticides controls simultaneously. Daozhou et al. [22] noted that to control spread of ZIKV preventive steps are necessary during sexual activities in ZIKV outbreak area.

Together all these facts lead to the increase of the people susceptible to the disease. As there is no vaccination available for the ZIKV disease, it is consider as a severe problem. Mathematical models are essential tools to study the dynamics of the spread of infectious diseases like ZIKV. Basic reproduction number provides information about how infection will be sustained. In this study, a new model with optimal control on spraying insecticides, preventions and treatment is examined. Pontryagin’s maximum principle, established by Pontrayagin et al. [23] is used to determine the optimal control. Result proves that optimal control gives significant reduction in ZIKV spread.

The paper is organised as follow. In Section 2, mathematical model by a system of ordinary differential equations with notations, assumptions and the flow of populations between compartments are described. For autonomous model having fixed rates of controls basic reproduction number of the whole system (human-mosquito combined) is calculated at disease free equilibrium and endemic equilibrium points. In Section 3, Stability of model has been discussed. In Section 4, cost control function for controls on spraying, preventions and treatment is formed and validated for the system of equations obtained in Section 2. In Section 5, numerical simulations are carried out, for both autonomous and control model. In Section 6, conclusions suggest that how to control disease spread with minimal cost.

2. Mathematical Model

The mathematical model is developed with following notations.

A system of non-linear differential equations is formulated to investigate spread of congenital transmission of ZIKV with control spraying, preventions and treatment. Dynamics of human population ${N}_{H}\left(t\right)$ and Vector population ${N}_{V}\left(t\right)$ are developed. Human population is classified amongst eight compartments viz. number of susceptible human child ${S}_{HC}\left(t\right)$ , number of susceptible human female ${S}_{HF}\left(t\right)$ , number of susceptible human male ${S}_{HM}\left(t\right)$ , number of infected human child ${I}_{HC}\left(t\right)$ , number of infected human female ${I}_{HF}\left(t\right)$ , number of infected human male ${I}_{HM}\left(t\right)$ , number of recovered human adult ${R}_{H}\left(t\right)$ , while, vector population ${N}_{V}\left(t\right)$ is divided in to two compartments viz. number of susceptible vectors ${S}_{V}\left(t\right)$ , number of infected vectors ${I}_{V}\left(t\right)$ . The population dynamics of above compartments is shown in Figure 1.

Figure 1. Dynamics of individuals in different compartments.

To investigate effects of congenital transmission of disease and its control by spraying on mosquitoes, preventions and treatment, human population is distributed in human child and human male (adult) and human female (adult) classes.

For model formulation increase in total population of vectors (mosquitoes) is considered, as subpopulation of vectors (mosquitoes) survive from the spraying, get matures and reproduces.

To prepare model, following possibilities of disease spread are considered:

1) Vertical Transmission of infection to new-borns from infected mothers;

2) Horizontal transmission to human child/adult (male and female) from infected vectors (mosquitoes) when they bite to human child/adult (male and female);

3) Heterosexual transmissions amongst human adults (male and female);

4) Horizontal transmission to vectors (mosquitoes) from infected human child/adult (male and female) when vectors (mosquitoes) bite to human child/ adult (male and female).

2.1. Disease Dynamics amongst Human Population

After birth from an infected human female, infected human children will join infected human child class at rate $\left(1-{u}_{3}\right)\delta {I}_{HF}$ . Susceptible human children will either join infected human child class at rate $\left(1-{u}_{1}\right)b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}$ due to bite of an

infected vector (mosquito) or susceptible human male (adult) or female (adult) classes at growth rates ${g}_{M}$ and ${g}_{F}$ respectively. Infected human child will join recovered human class at recovery/treatment rate ${u}_{4}{k}_{HC}$ . Disease induced death rate ${\alpha}_{HC}$ will make very high effect on infected child population as possibility of microcephaly and Guillain-Barre syndrome due to ZIKV transmission occurred during pregnancy.

Susceptible human male (adult) and female (adult) will either join infected human male (adult) and female (adult) class at rates

$\left(1-{u}_{1}\right)b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HF}$ and $\left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HM}$

respectively due to vector bite and sexual intercourse with infected human female (adult) and male (adult) respectively. Infected human male (adult) and female (adult) will join recovered human class at recovery/treatment rate ${u}_{4}{k}_{HM}$ and ${u}_{4}{k}_{HF}$ respectively. Disease induced death rate ${\alpha}_{H}$ will effect on infected human male and female (adult) population. At each stage, the natural death rate in human child/adult classes ${\mu}_{H}$ is taken into account.

2.2. Disease Dynamics amongst Vector (Mosquito) Population

Vector (mosquito) population is divided in to two compartments susceptible vectors and infected vectors. Susceptible vector gets infection from infected human child/adult, when vectors (mosquitoes) bite them, and joins infected vector

class at rate $b\frac{{\beta}_{V}}{{N}_{V}}\left({I}_{HC}+{I}_{HF}+{I}_{HM}\right)$ . A portion of susceptible and infected vec-

tors (mosquitoes) will be eliminated due to either natural death rate in vectors (mosquitoes) ${\mu}_{V}$ or spraying at rate $u{}_{5}r$ . Disease induced death rate ${\alpha}_{V}$ will also affect the infected vector (mosquito) population.

With above discussion the dynamics of Zika disease can be represented by the system of non-linear differential equations as

$\frac{\text{d}{S}_{HC}\left(t\right)}{\text{d}t}={B}_{HC}-\left(1-{u}_{1}\right)b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}{S}_{HC}-{\mu}_{H}{S}_{HC}-\left({g}_{F}+{g}_{M}\right){S}_{HC}$ (1)

$\frac{\text{d}{S}_{HF}\left(t\right)}{\text{d}t}={B}_{HF}-\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HM}\right){S}_{HF}+{g}_{F}{S}_{HC}-{\mu}_{H}{S}_{HF}$ (2)

#Math_33# (3)

#Math_34# (4)

$\frac{\text{d}{I}_{HF}\left(t\right)}{\text{d}t}=\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HM}\right){S}_{HF}-{u}_{4}{k}_{HF}{I}_{HF}-\left({\mu}_{H}+{\alpha}_{H}\right){I}_{HF}$ (5)

$\frac{\text{d}{I}_{HM}\left(t\right)}{\text{d}t}=\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HF}\right){S}_{HM}-{u}_{4}{k}_{HM}{I}_{HM}-\left({\mu}_{H}+{\alpha}_{H}\right){I}_{HM}$ (6)

$\frac{\text{d}{R}_{H}\left(t\right)}{\text{d}t}={u}_{4}\left({k}_{HC}{I}_{HC}+{k}_{HF}{I}_{HF}+{k}_{HM}{I}_{HM}\right)-{\mu}_{H}{R}_{H}$ (7)

$\frac{\text{d}{S}_{V}\left(t\right)}{\text{d}t}={B}_{V}-{\mu}_{V}{S}_{V}-b\frac{{\beta}_{V}}{{N}_{V}}\left({I}_{HC}+{I}_{HF}+{I}_{HM}\right){S}_{V}-{u}_{5}r{S}_{V}$ (8)

$\frac{\text{d}{I}_{V}\left(t\right)}{\text{d}t}=b\frac{{\beta}_{V}}{{N}_{V}}\left({I}_{HC}+{I}_{HF}+{I}_{HM}\right){S}_{V}-\left({\mu}_{V}+{\alpha}_{V}\right){I}_{V}-{u}_{5}r{I}_{V}$ (9)

Adding Equations (1) to (7),

$\begin{array}{l}\frac{\text{d}{N}_{H}}{\text{d}t}={B}_{HC}+{B}_{HF}+{B}_{HM}-{\mu}_{H}\left({S}_{HC}+{S}_{HF}+{S}_{HM}+{I}_{HC}+{I}_{HF}+{I}_{HM}+{R}_{H}\right)\\ \text{}+\left(1-{u}_{3}\right)\delta {I}_{HF}-{\alpha}_{HC}{I}_{HC}-{\alpha}_{H}\left({I}_{HF}+{I}_{HM}\right)\end{array}$

And, adding Equations (8) and (9),

$\frac{\text{d}{N}_{V}}{\text{d}t}={B}_{V}-{\mu}_{V}\left({S}_{V}+{I}_{V}\right)-{u}_{5}r\left({S}_{V}+{I}_{V}\right)-{\alpha}_{V}{I}_{V}$

Hence,

$\frac{\text{d}{N}_{H}}{\text{d}t}\le {B}_{HC}+{B}_{HF}+{B}_{HM}-{\mu}_{H}{N}_{H}\Rightarrow \underset{t\to \infty}{\mathrm{lim}}\text{Sup}{N}_{H}\le \frac{{B}_{HC}+{B}_{HF}+{B}_{HM}}{{\mu}_{H}}$

and,

$\frac{\text{d}{N}_{V}}{\text{d}t}\le {B}_{V}-{\mu}_{V}{N}_{V}\Rightarrow \underset{t\to \infty}{\mathrm{lim}}\text{Sup}{N}_{V}\le \frac{{B}_{V}}{{\mu}_{V}}$ .

So, the feasible region of the system is

$\Omega =\left\{\begin{array}{l}\left({S}_{HC},{S}_{HF},{S}_{HM},{I}_{HC},{I}_{HF},{I}_{HM},{R}_{H},{S}_{V},{I}_{V}\right):\\ {S}_{THC}+{S}_{UHC}+{S}_{THA}+{S}_{UHA}+{I}_{THC}+{I}_{UHC}+{I}_{THA}+{I}_{UHA}\le \frac{{B}_{HC}+{B}_{HF}+{B}_{HM}}{{\mu}_{H}},\\ \text{}{S}_{V}+{I}_{V}\le \frac{{B}_{V}}{{\mu}_{V}},\\ {S}_{HC}>0,{S}_{HF}>0,{S}_{HM}>0,{I}_{HC}>0,{I}_{HF}\ge 0,{I}_{HM}\ge 0,{R}_{H}\ge 0,{S}_{V}>0,{I}_{V}\ge 0\end{array}\right\}$

2.3. Computation of Basic Reproduction Number ${R}_{0}$

As at disease free equilibrium, ${I}_{HC}=0,{I}_{HF}=0,{I}_{HM}=0,{I}_{V}=0$ and ${R}_{H}=0$ .

Hence, let ${X}_{0}=\left(\frac{{B}_{HC}}{{\mu}_{H}+{g}_{F}+{g}_{M}},\frac{{B}_{HF}}{{\mu}_{H}-{g}_{F}},\frac{{B}_{HM}}{{\mu}_{H}-{g}_{M}},0,0,0,0,\frac{{B}_{V}}{{\mu}_{V}+{u}_{5}r},0\right)$ be disease free equilibrium point of system.

Using next generation method [24] , let

${X}^{\prime}={\left({I}_{HC},{I}_{HF},{I}_{HM},{I}_{V},{S}_{HC},{S}_{HF},{S}_{HM},{S}_{V},{R}_{H}\right)}^{\prime}$

where dash denotes derivative

$\therefore {X}^{\prime}=\frac{\text{d}X}{\text{d}t}=\Im \left(x\right)-v\left(x\right)$

where,

$\Im \left(X\right)=\left(\begin{array}{c}\left(1-{u}_{1}\right)b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}{S}_{HC}\\ \left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}{S}_{HF}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HM}{S}_{HF}\\ \left(1-{u}_{1}\right)b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}{S}_{HM}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HF}{S}_{HM}\\ b\frac{{\beta}_{V}}{{N}_{V}}\left({I}_{HC}+{I}_{HF}+{I}_{HM}\right){S}_{V}\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right)$

and

$\upsilon \left(X\right)=\left(\begin{array}{c}-\left(1-{u}_{3}\right)\delta {I}_{HF}+{u}_{4}{k}_{HC}{I}_{HC}+{\mu}_{H}{I}_{HC}+{\alpha}_{HC}{I}_{HC}\\ {u}_{4}{k}_{HF}{I}_{HF}+{\mu}_{H}{I}_{HF}+{\alpha}_{H}{I}_{HF}\\ {u}_{4}{k}_{HM}{I}_{HM}+{\mu}_{H}{I}_{HM}+{\alpha}_{H}{I}_{HM}\\ {\mu}_{V}{I}_{V}+{\alpha}_{V}{I}_{V}+{u}_{5}r{I}_{V}\\ -{B}_{HC}+\left(1-{u}_{1}\right)b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}{S}_{HC}+{\mu}_{H}{S}_{HC}+\left({g}_{F}+{g}_{M}\right){S}_{HC}\\ -{B}_{HF}+\left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}{S}_{HF}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HM}{S}_{HF}+{\mu}_{H}{S}_{HF}-{g}_{F}{S}_{HC}\\ -{B}_{HM}+\left(1-{u}_{1}\right)b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}{S}_{HM}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HF}{S}_{HF}+{\mu}_{H}{S}_{HM}-{g}_{M}{S}_{HC}\\ -{B}_{V}+b\frac{{\beta}_{V}}{{N}_{V}}\left({I}_{HC}+{I}_{HF}+{I}_{HM}\right){S}_{V}+{u}_{5}r{S}_{V}+{\mu}_{V}{S}_{V}\\ {\alpha}_{V}{S}_{V}+{\mu}_{V}{I}_{V}+{u}_{5}r{I}_{V}\end{array}\right)$

Using, $F=\left[\frac{\partial {\Im}_{i}\left({X}_{0}\right)}{\partial {X}_{j}}\right]$ and $V=\left[\frac{\partial {\upsilon}_{i}\left({X}_{0}\right)}{\partial {X}_{j}}\right]\text{for}i,j=1,2,3,\cdots ,9$

Therefore,

$F=\left(\begin{array}{ccccccccc}0& 0& 0& {b}_{14}& 0& 0& 0& 0& 0\\ 0& 0& {b}_{23}& {b}_{24}& 0& 0& 0& 0& 0\\ 0& {b}_{32}& 0& {b}_{34}& 0& 0& 0& 0& 0\\ {b}_{41}& {b}_{42}& {b}_{43}& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& 0& 0& 0\end{array}\right)$

where,

$\begin{array}{l}{b}_{14}=\frac{\left(1-{u}_{1}\right)b{\beta}_{HC}{B}_{HC}}{\left({\mu}_{H}+{g}_{F}+{g}_{M}\right){N}_{H}},\text{}{b}_{23}=\frac{\left(1-{u}_{2}\right)c{\beta}_{S}{B}_{HF}}{\left({\mu}_{H}-{g}_{F}\right){N}_{H}},\\ {b}_{24}=\frac{\left(1-{u}_{1}\right)b{\beta}_{HF}{B}_{HF}}{\left({\mu}_{H}-{g}_{F}\right){N}_{H}},\text{}{b}_{32}=\frac{\left(1-{u}_{2}\right)c{\beta}_{S}{B}_{HM}}{\left({\mu}_{H}-{g}_{M}\right){N}_{H}},\\ {b}_{34}=\frac{\left(1-{u}_{1}\right)b{\beta}_{HM}{B}_{HM}}{\left({\mu}_{H}-{g}_{M}\right){N}_{H}},\text{}{b}_{41}=\frac{b{\beta}_{V}{B}_{V}}{\left({\mu}_{V}+{u}_{5}r\right){N}_{V}},\\ {b}_{42}=\frac{b{\beta}_{V}{B}_{V}}{\left({\mu}_{V}+{u}_{5}r\right){N}_{V}},\text{}{b}_{43}=\frac{b{\beta}_{V}{B}_{V}}{\left({\mu}_{V}+{u}_{5}r\right){N}_{V}}\end{array}$

and

$V=\left(\begin{array}{ccccccccc}{a}_{11}& {a}_{12}& 0& 0& 0& 0& 0& 0& 0\\ 0& {a}_{22}& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& {a}_{33}& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& {a}_{44}& 0& 0& 0& 0& 0\\ 0& 0& 0& {a}_{54}& {a}_{55}& 0& 0& 0& 0\\ 0& 0& {a}_{63}& {a}_{64}& 0& {a}_{66}& 0& 0& 0\\ 0& {a}_{72}& 0& {a}_{74}& 0& 0& {a}_{77}& 0& 0\\ {a}_{81}& {a}_{82}& {a}_{83}& 0& 0& 0& 0& {a}_{88}& 0\\ {a}_{91}& {a}_{92}& {a}_{93}& 0& 0& 0& 0& 0& {a}_{99}\end{array}\right)$

where,

$\begin{array}{l}{a}_{11}={u}_{4}{k}_{HC}+{\mu}_{H}+{\alpha}_{HC},\text{}{a}_{12}=-\left(1-{u}_{3}\right)\delta ,\\ {a}_{22}={u}_{4}{k}_{HF}+{\mu}_{H}+{\alpha}_{H},\text{}{a}_{33}={u}_{4}{k}_{HM}+{\mu}_{H}+{\alpha}_{H}\\ {a}_{33}={u}_{4}{k}_{HM}+{\mu}_{H}+{\alpha}_{H},\text{}{a}_{44}={\mu}_{V}+{\alpha}_{V}+{u}_{5}r,\\ {a}_{54}=\frac{\left(1-{u}_{1}\right)b{\beta}_{HC}{B}_{HC}}{\left({\mu}_{H}+{g}_{F}+{g}_{M}\right){N}_{H}},\text{}{a}_{55}={\mu}_{H}+{g}_{F}+{g}_{M},\\ {a}_{63}=\frac{\left(1-{u}_{2}\right)c{\beta}_{S}{B}_{HF}}{\left({\mu}_{H}-{g}_{F}\right){N}_{H}},\text{}{a}_{64}=\frac{\left(1-{u}_{1}\right)b{\beta}_{HF}{B}_{HF}}{\left({\mu}_{H}-{g}_{F}\right){N}_{H}},\\ {a}_{66}={\mu}_{H}-{g}_{F},\text{}{a}_{72}=\frac{\left(1-{u}_{2}\right)c{\beta}_{S}{B}_{HM}}{\left({\mu}_{H}-{g}_{M}\right){N}_{H}},\\ {a}_{74}=\frac{\left(1-{u}_{1}\right)b{\beta}_{HM}{B}_{HM}}{\left({\mu}_{H}-{g}_{M}\right){N}_{H}},\text{}{a}_{77}={\mu}_{H}-{g}_{M},\\ {a}_{81}=\frac{b{\beta}_{V}{B}_{V}}{\left({\mu}_{V}+{u}_{5}r\right){N}_{V}},\text{}{a}_{82}=\frac{b{\beta}_{V}{B}_{V}}{\left({\mu}_{V}+{u}_{5}r\right){N}_{V}},\\ {a}_{83}=\frac{b{\beta}_{V}{B}_{V}}{\left({\mu}_{V}+{u}_{5}r\right){N}_{V}},\text{}{a}_{88}={\mu}_{V}+{u}_{5}r,\\ {a}_{91}=-{u}_{4}{k}_{HC},\text{}{a}_{92}=-{u}_{4}{k}_{HF},\\ {a}_{93}=-{u}_{4}{k}_{HM},\text{}{a}_{99}={\mu}_{H}\end{array}$

The basic reproduction number ${R}_{0}$ is spectral radius of matrix $F{V}^{-1}$ .

With the parametric values given in Table 1, ${R}_{0}=0.4757$ .

3. Stability

3.1. Stability at Disease Free Equilibrium

If all eigenvalues of Jacobian matrix of the system of differential Equations (1) to (10) have negative real parts at

${X}_{0}=\left(\frac{{B}_{HC}}{{\mu}_{H}+{g}_{F}+{g}_{M}},\frac{{B}_{HF}}{{\mu}_{H}-{g}_{F}},\frac{{B}_{HM}}{{\mu}_{H}-{g}_{M}},0,0,0,0,\frac{{B}_{V}}{{\mu}_{V}+{u}_{5}r},0\right)$

then disease free equilibrium point becomes stable. Jacobian matrix of the system at DFE is,

$J=\left(\begin{array}{ccccccccc}{c}_{11}& 0& 0& 0& 0& 0& 0& 0& {c}_{19}\\ 0& {c}_{22}& 0& 0& 0& {c}_{26}& 0& 0& {c}_{29}\\ 0& 0& {c}_{33}& 0& {c}_{35}& 0& 0& 0& {c}_{39}\\ 0& 0& 0& {c}_{44}& {c}_{45}& 0& 0& 0& {c}_{49}\\ 0& 0& 0& 0& {c}_{55}& {c}_{56}& 0& 0& {c}_{59}\\ 0& 0& 0& 0& {c}_{65}& {c}_{66}& 0& 0& {c}_{69}\\ 0& 0& 0& {c}_{74}& {c}_{75}& {c}_{76}& {c}_{77}& 0& 0\\ 0& 0& 0& {c}_{84}& {c}_{85}& {c}_{86}& 0& {c}_{88}& 0\\ 0& 0& 0& {c}_{94}& {c}_{95}& {c}_{96}& 0& 0& {c}_{99}\end{array}\right)$

where,

$\begin{array}{l}{c}_{11}=-{\mu}_{H}-{g}_{F}-{g}_{M},\text{}{c}_{19}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HC}{B}_{HC}}{\left({\mu}_{H}+{g}_{F}+{g}_{M}\right){N}_{H}},\\ {c}_{21}={g}_{F},\text{}{c}_{22}=-{\mu}_{H},\text{}{c}_{26}=-\frac{\left(1-{u}_{2}\right)c{\beta}_{S}{B}_{HF}}{\left({\mu}_{H}-{g}_{F}\right){N}_{H}},\\ {c}_{29}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HF}{B}_{HF}}{\left({\mu}_{H}-{g}_{F}\right){N}_{H}},\text{}{c}_{31}={g}_{M},\text{}{c}_{33}=-{\mu}_{H},\\ {c}_{35}=-\frac{\left(1-{u}_{2}\right)c{\beta}_{S}{B}_{HM}}{\left({\mu}_{H}-{g}_{M}\right){N}_{H}},\text{}{c}_{39}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HM}{B}_{HM}}{\left({\mu}_{H}-{g}_{M}\right){N}_{H}},\\ {c}_{44}=-{u}_{4}{k}_{HC}-{\mu}_{H}-{\alpha}_{HC},\text{}{c}_{45}=\left(1-{u}_{3}\right)\delta ,\\ {c}_{49}=-{c}_{19},\text{}{c}_{55}=-{u}_{4}{k}_{HF}-{\mu}_{H}-{\alpha}_{H},\text{}{c}_{56}=-{c}_{26},\text{}{c}_{59}=-{c}_{29},\\ {c}_{65}=-{c}_{35},\text{}{c}_{66}=-{u}_{4}{k}_{HM}-{\mu}_{H}-{\alpha}_{H},\text{}{c}_{69}=-{c}_{39},\\ {c}_{74}={u}_{4}{k}_{HC},\text{}{c}_{75}={u}_{4}{k}_{HF},\text{}{c}_{76}={u}_{4}{k}_{HM},\text{}{c}_{77}=-{\mu}_{H},\\ {c}_{84}=-\frac{b{\beta}_{V}{B}_{V}}{\left({\mu}_{V}+{u}_{5}r\right){N}_{V}},\text{}{c}_{85}={c}_{84},\text{}{c}_{86}={c}_{84},\text{}{c}_{88}=-{\mu}_{V}-{u}_{5}r,\\ {c}_{94}=-{c}_{84},\text{}{c}_{95}=-{c}_{84},\text{}{c}_{96}=-{c}_{84},\text{}{c}_{99}={c}_{88}-{\alpha}_{V}.\end{array}$

With the parametric values given in Table 1, $\text{trace}\left(J\right)<0$ and also $\mathrm{det}\left(J\right)>0$ , if ${R}_{0}<1$ .

This implies that disease free equilibrium is locally asymptotically stable if ${R}_{0}<1$ and unstable, if ${R}_{0}>1$ .

3.2. Stability at Endemic Equilibrium Point

For the endemic equilibrium, using

$\begin{array}{l}\frac{\text{d}{S}_{THC}}{\text{d}t}=0,\text{}\frac{\text{d}{S}_{UHC}}{\text{d}t}=0,\text{}\frac{\text{d}{S}_{THA}}{\text{d}t}=0,\text{}\\ \frac{\text{d}{S}_{UHA}}{\text{d}t}=0,\text{}\frac{\text{d}{I}_{THC}}{\text{d}t}=0,\\ \frac{\text{d}{I}_{UHC}}{\text{d}t}=0,\text{}\frac{\text{d}{I}_{THA}}{\text{d}t}=0,\text{}\\ \frac{\text{d}{I}_{UHA}}{\text{d}t}=0,\text{}\frac{\text{d}{S}_{V}}{\text{d}t}=0,\text{}\frac{\text{d}{I}_{V}}{\text{d}t}=0\end{array}$

These give, ${X}^{\ast}=\left({S}_{THC}^{\ast},{S}_{UHC}^{\ast},{S}_{THA}^{\ast},{S}_{UHA}^{\ast},{I}_{THC}^{\ast},{I}_{UHC}^{\ast},{I}_{THA}^{\ast},{I}_{UHA}^{\ast},{S}_{V}^{\ast},{I}_{V}^{\ast}\right)$ as an endemic equilibrium point.

${J}^{\ast}=\left(\begin{array}{ccccccccc}{d}_{11}& 0& 0& 0& 0& 0& 0& 0& {d}_{19}\\ {d}_{21}& {d}_{22}& 0& 0& 0& {d}_{26}& 0& 0& {d}_{29}\\ {d}_{31}& 0& {d}_{33}& 0& {d}_{35}& 0& 0& 0& {d}_{39}\\ {d}_{41}& 0& 0& {d}_{44}& {d}_{45}& 0& 0& 0& {d}_{49}\\ 0& {d}_{52}& 0& 0& {d}_{55}& {d}_{56}& 0& 0& {d}_{59}\\ 0& 0& {d}_{63}& 0& {d}_{65}& {d}_{66}& 0& 0& {d}_{69}\\ 0& 0& 0& {d}_{74}& {d}_{75}& {d}_{76}& {d}_{77}& 0& 0\\ 0& 0& 0& {d}_{84}& {d}_{85}& {d}_{86}& 0& {d}_{88}& 0\\ 0& 0& 0& {d}_{94}& {d}_{95}& {d}_{96}& 0& {d}_{98}& {d}_{99}\end{array}\right)$

where,

Table 1. Description of notation and parametric values.

$\begin{array}{l}{d}_{11}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HC}}{{N}_{H}}{I}_{V}^{\ast}-{\mu}_{H}-{g}_{F}-{g}_{M},\\ {d}_{19}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HC}}{{N}_{H}}{S}_{HC}^{\ast},\text{}{d}_{21}={g}_{F},\\ {d}_{22}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HF}}{{N}_{H}}{I}_{V}^{\ast}-\frac{\left(1-{u}_{2}\right)c{\beta}_{S}}{{N}_{H}}{I}_{HM}^{\ast}-{\mu}_{H},\\ {d}_{26}=-\frac{\left(1-{u}_{2}\right)c{\beta}_{S}}{{N}_{H}}{S}_{HF}^{\ast},\text{}{d}_{29}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HF}}{{N}_{H}}{S}_{HF}^{\ast},\\ {d}_{31}={g}_{M},\text{}{d}_{33}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HM}}{{N}_{H}}{I}_{V}^{\ast}-\frac{\left(1-{u}_{2}\right)c{\beta}_{S}}{{N}_{H}}{I}_{HF}^{\ast}-{\mu}_{H},\\ {d}_{35}=-\frac{\left(1-{u}_{2}\right)c{\beta}_{S}}{{N}_{H}}{S}_{HM}^{\ast},\text{}{d}_{39}=-\frac{\left(1-{u}_{1}\right)b{\beta}_{HM}}{{N}_{H}}{I}_{V}^{\ast},\\ {d}_{41}=-{d}_{39},\text{}{d}_{44}=-{u}_{4}{k}_{HC}-{\mu}_{H}-{\alpha}_{HC},\text{}{d}_{45}=\left(1-{u}_{3}\right)\delta ,\\ {d}_{49}=-{d}_{19},\text{}{d}_{52}=\frac{\left(1-{u}_{1}\right)b{\beta}_{HF}}{{N}_{H}}{I}_{V}^{\ast}+\frac{\left(1-{u}_{2}\right)c{\beta}_{S}}{{N}_{H}}{I}_{HM}^{\ast},\\ {d}_{55}=-{u}_{4}{k}_{HF}-{\mu}_{H}-{\alpha}_{H},\text{}{d}_{56}=-{d}_{26},\text{}{d}_{59}=-{d}_{29},\\ {d}_{63}=\frac{\left(1-{u}_{1}\right)b{\beta}_{HM}}{{N}_{H}}{I}_{V}^{\ast}+\frac{\left(1-{u}_{2}\right)c{\beta}_{S}}{{N}_{H}}{I}_{HF}^{\ast},\text{}{d}_{65}=-{d}_{35},\\ {d}_{66}=-{u}_{4}{k}_{HM}-{\mu}_{H}-{\alpha}_{H},\text{}{d}_{69}=-{d}_{39},\text{}{d}_{74}={u}_{4}{k}_{HC},\text{}\\ {d}_{75}={u}_{4}{k}_{HF},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{d}_{76}={u}_{4}{k}_{HM},\text{}{d}_{77}=-{\mu}_{H},\text{}{d}_{84}=-\frac{b{\beta}_{V}}{{N}_{V}}{S}_{V}^{\ast},\\ {d}_{85}={d}_{84},\text{}{d}_{86}={d}_{84},\text{}{d}_{88}=-{\mu}_{V}-\frac{b{\beta}_{V}}{{N}_{V}}\uff08{I}_{Hc}{}^{*}+{I}_{HF}{}^{*}+{I}_{HM}{}^{*}\uff09-{u}_{5}r,\\ {d}_{94}=-{d}_{84},\text{}{d}_{95}=-{d}_{84},\text{}{d}_{96}=-{d}_{84},\\ {d}_{98}=\frac{b{\beta}_{V}}{{N}_{V}}\uff08{I}_{Hc}{}^{*}+{I}_{HF}{}^{*}+{I}_{HM}{}^{*}\uff09,\text{}{d}_{99}=-{\mu}_{V}-{\alpha}_{V}-{u}_{5}r.\end{array}$

With the parametric values given in Table 1, $\text{trace}\left({J}^{\ast}\right)<0$ and also $\mathrm{det}\left({J}^{\ast}\right)>0$ , if ${R}_{0}<1$ . This implies that endemic equilibrium is locally asymptotically stable if ${R}_{0}<1$ and unstable, if ${R}_{0}>1$ .

4. Optimal Control

For Zika disease an optimal control model is formulated, to derive optimal prevention from mosquito bite ${u}_{1}$ , optimal prevention to stop pregnancy ${u}_{2}$ and sexual transmission ${u}_{3}$ , optimal treatment ${u}_{4}$ and optimal spraying ${u}_{5}$ with minimal implementation cost, in order to minimise the number of infected individuals for model described by Equations (1) to (9) in the time interval $\left[0,T\right]$ with the feasible region same as given by $\Omega $ in Section 2.

Considering the cost-functional as,

$\begin{array}{l}J\left({u}_{1},{u}_{2},{u}_{3},{u}_{4},{u}_{5}\right)\\ ={\displaystyle \underset{0}{\overset{T}{\int}}({w}_{21}{\left({S}_{HC}\right)}^{2}+{w}_{22}{\left({S}_{HF}\right)}^{2}+{w}_{23}{\left({S}_{HM}\right)}^{2}+{w}_{24}{\left({I}_{HC}\right)}^{2}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{w}_{25}{\left({I}_{HF}\right)}^{2}+{w}_{26}{\left({I}_{HM}\right)}^{2}+{w}_{27}{\left({R}_{H}\right)}^{2}+{w}_{28}{\left({S}_{V}\right)}^{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{w}_{29}{\left({I}_{V}\right)}^{2}+{w}_{11}{u}_{1}^{2}+{w}_{12}{u}_{2}^{2}+{w}_{13}{u}_{3}^{2}+{w}_{14}{u}_{4}^{2}+{w}_{15}{u}_{5}^{2})\text{d}t\end{array}$

where, $x=\left({S}_{HC},{S}_{HF},{S}_{HM},{I}_{HC},{I}_{HF},{I}_{HM},{R}_{H},{S}_{V},{I}_{V}\right)$ , $u=\left({u}_{1},{u}_{2},{u}_{3},{u}_{4},{u}_{5}\right)$ and ${w}_{ij}$ are weights to regularise the optimal control. In order minimise the disease spread using, ${w}_{ij}>0$ , for $i=2,j=1,2,\cdots ,9$ . Also to minimise the cost associated with strategies applied for preventions, treatment and spraying in a way that spread of infection can be controlled by choosing weights ${w}_{km}>0$ , for $k=1,m=1,2,\cdots ,5$ . So, the optimal control task reads as $\underset{u}{\mathrm{min}}J\left({u}_{1},{u}_{2},{u}_{3},{u}_{4},{u}_{5}\right)$ such that $P\left(x,u\right)=0$ . Where $P\left(x,u\right)=0$ denotes the system of equations de-

fined in (1) to (9), i.e. $\stackrel{\dot{}}{x}=M\left(t,x,u\right),x\left(0\right)={x}_{0}$ . The optimal control ${u}^{\ast}$ can be obtained from $J\left({u}_{1}^{\ast},{u}_{2}^{\ast},{u}_{3}^{\ast},{u}_{4}^{\ast},{u}_{5}^{\ast}\right)=\underset{u}{\mathrm{min}}J\left({u}_{1},{u}_{2},{u}_{3},{u}_{4},{u}_{5}\right)$ such that $P\left(x,u\right)=0$ .

The control set is

$\Gamma =\left\{\begin{array}{l}\left({u}_{1},{u}_{2},{u}_{3},{u}_{4},{u}_{5}\right)/{u}_{i}\left(t\right)\text{ispiecewisecontinouson}\left[0,T\right],\\ {a}_{i}\le {u}_{i}\left(t\right)\le {b}_{i},i=1,2,3,4,5,{a}_{i},{b}_{i}\text{areconstantsin}\left[0,1\right]\end{array}\right\}$

Using Lagrangian techniques for a problem along with Hamiltonian, the adjoint variable is needed to construct for the optimal control problem given by (1) to (9).

Introducing the Lagrangian to derive the optimality conditions,

$\begin{array}{l}L\left(x,u,\lambda \right)=({w}_{21}{\left({S}_{HC}\right)}^{2}+{w}_{22}{\left({S}_{HF}\right)}^{2}+{w}_{23}{\left({S}_{HM}\right)}^{2}+{w}_{24}{\left({I}_{HC}\right)}^{2}\\ \text{}+{w}_{25}{\left({I}_{HF}\right)}^{2}+{w}_{26}{\left({I}_{HM}\right)}^{2}+{w}_{27}{\left({R}_{H}\right)}^{2}+{w}_{28}{\left({S}_{V}\right)}^{2}+{w}_{29}{\left({I}_{V}\right)}^{2}\\ \text{}+{w}_{11}{u}_{1}^{2}+{w}_{12}{u}_{2}^{2}+{w}_{13}{u}_{3}^{2}+{w}_{14}{u}_{4}^{2}+{w}_{15}{u}_{5}^{2})\end{array}$

To obtain the minimal value of the Lagrangian, defining the Hamiltonian $H$ for the control problem as,

$\begin{array}{l}H\left(x,u,\lambda \right)={w}_{21}{\left({S}_{HC}\right)}^{2}+{w}_{22}{\left({S}_{HF}\right)}^{2}+{w}_{23}{\left({S}_{HM}\right)}^{2}+{w}_{24}{\left({I}_{HC}\right)}^{2}\\ +{w}_{25}{\left({I}_{HF}\right)}^{2}+{w}_{26}{\left({I}_{HM}\right)}^{2}+{w}_{27}{\left({R}_{H}\right)}^{2}+{w}_{28}{\left({S}_{V}\right)}^{2}\\ +{w}_{29}{\left({I}_{V}\right)}^{2}+{w}_{11}{u}_{1}^{2}+{w}_{12}{u}_{2}^{2}+{w}_{13}{u}_{3}^{2}+{w}_{14}{u}_{4}^{2}+{w}_{15}{u}_{5}^{2}\\ +{\lambda}_{1}\left({B}_{HC}-\left(1-{u}_{1}\right)b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}{S}_{HC}-{\mu}_{H}{S}_{HC}-\left({g}_{F}+{g}_{M}\right){S}_{HC}\right)\\ +{\lambda}_{2}\left({B}_{HF}-\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HM}\right){S}_{HF}+{g}_{F}{S}_{HC}-{\mu}_{H}{S}_{HF}\right)\\ +{\lambda}_{3}\left({B}_{HM}-\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HF}\right){S}_{HM}+{g}_{M}{S}_{HC}-{\mu}_{H}{S}_{HM}\right)\\ +{\lambda}_{4}\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}{S}_{HC}+\left(1-{u}_{3}\right)\delta {I}_{HF}-{u}_{4}{k}_{HC}{I}_{HC}-\left({\mu}_{H}+{\alpha}_{HC}\right){I}_{HC}\right)\\ +{\lambda}_{5}\left(\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HM}\right){S}_{HF}-{u}_{4}{k}_{HF}{I}_{HF}-\left({\mu}_{H}+{\alpha}_{H}\right){I}_{HF}\right)\\ +{\lambda}_{6}\left(\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HF}\right){S}_{HM}-{u}_{4}{k}_{HM}{I}_{HM}-\left({\mu}_{H}+{\alpha}_{H}\right){I}_{HM}\right)\\ +{\lambda}_{7}\left({u}_{4}\left({k}_{HC}{I}_{HC}+{k}_{HF}{I}_{HF}+{k}_{HM}{I}_{HM}\right)-{\mu}_{H}{R}_{H}\right)\\ +{\lambda}_{8}\left({B}_{V}-{\mu}_{V}{S}_{V}-b\frac{{\beta}_{V}}{{N}_{V}}\left({I}_{HC}+{I}_{HF}+{I}_{HM}\right){S}_{V}-{u}_{5}r{S}_{V}\right)\\ +{\lambda}_{9}\left(b\frac{{\beta}_{V}}{{N}_{V}}\left({I}_{HC}+{I}_{HF}+{I}_{HM}\right){S}_{V}-\left({\mu}_{V}+{\alpha}_{V}\right){I}_{V}-{u}_{5}r{I}_{V}\right)\end{array}$

To get optimality Pontrayagin’s maximum (minimum) principle, for the model is defined as follow.

If $\left({u}_{1}^{\ast},{u}_{2}^{\ast},{u}_{3}^{\ast},{u}_{4}^{\ast},{u}_{5}^{\ast}\right)$ is optimal solution of an optimal control problem then there exists a nontrivial vector function $\lambda \left(t\right)=\left({\lambda}_{1}\left(t\right),{\lambda}_{2}\left(t\right),{\lambda}_{3}\left(t\right),{\lambda}_{4}\left(t\right),{\lambda}_{5}\left(t\right),{\lambda}_{6}\left(t\right),{\lambda}_{7}\left(t\right),{\lambda}_{8}\left(t\right),{\lambda}_{9}\left(t\right)\right)$ satisfying following equations:

1) The state equation $\frac{\text{d}x}{\text{d}t}=\frac{\partial H\left(t,{u}_{1}^{\ast},{u}_{2}^{\ast},{u}_{3}^{\ast},{u}_{4}^{\ast},{u}_{5}^{\ast},\lambda \left(t\right)\right)}{\partial \lambda}$

2) The optimality condition $0=\frac{\partial H\left(t,{u}_{1}^{\ast},{u}_{2}^{\ast},{u}_{3}^{\ast},{u}_{4}^{\ast},{u}_{5}^{\ast},\lambda \left(t\right)\right)}{\partial u}$

3) The adjoint equation $\frac{\text{d}\lambda}{\text{d}t}=\frac{\partial H\left(t,{u}_{1}^{\ast},{u}_{2}^{\ast},{u}_{3}^{\ast},{u}_{4}^{\ast},{u}_{5}^{\ast},\lambda \left(t\right)\right)}{\partial x}$ .

Using equations stated above for the Hamiltonian defined by (13), using state equations as given by (1) to (9), adjoint equations as

${\stackrel{\dot{}}{\lambda}}_{1}=-2{w}_{21}{S}_{HC}+\left({\lambda}_{1}-{\lambda}_{4}\right)\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}\right)+{\lambda}_{1}\left({\mu}_{H}+{g}_{F}+{g}_{M}\right)$

${\stackrel{\dot{}}{\lambda}}_{2}=-2{w}_{22}{S}_{HF}+\left({\lambda}_{2}-{\lambda}_{5}\right)\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HM}\right)+{\lambda}_{2}\left({\mu}_{H}-{g}_{F}\right)$

${\stackrel{\dot{}}{\lambda}}_{3}=-2{w}_{23}{S}_{HM}+\left({\lambda}_{3}-{\lambda}_{6}\right)\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}+\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{I}_{HF}\right)+{\lambda}_{3}\left({\mu}_{H}-{g}_{M}\right)$

${\stackrel{\dot{}}{\lambda}}_{4}=-2{w}_{24}{I}_{HC}+\left({\lambda}_{8}-{\lambda}_{9}\right)\left(\frac{b{\beta}_{V}}{{N}_{V}}{S}_{V}\right)+\left({\lambda}_{4}-{\lambda}_{7}\right){u}_{4}{k}_{HC}+{\lambda}_{4}\left({\mu}_{H}+{\alpha}_{HC}\right)$

$\begin{array}{l}{\stackrel{\dot{}}{\lambda}}_{5}=-2{w}_{25}{I}_{HF}+\left({\lambda}_{3}-{\lambda}_{6}\right)\left(\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{S}_{HM}\right)+\left({\lambda}_{8}-{\lambda}_{9}\right)\left(\frac{b{\beta}_{V}}{{N}_{V}}{S}_{V}\right)\\ \text{}+\left({\lambda}_{5}-{\lambda}_{7}\right){u}_{4}{k}_{HF}+{\lambda}_{5}\left({\mu}_{H}+{\alpha}_{H}\right)-{\lambda}_{5}\left(1-{u}_{3}\right)\delta \end{array}$

$\begin{array}{l}{\stackrel{\dot{}}{\lambda}}_{6}=-2{w}_{26}{I}_{HF}+\left({\lambda}_{2}-{\lambda}_{5}\right)\left(\left(1-{u}_{2}\right)c\frac{{\beta}_{S}}{{N}_{H}}{S}_{HF}\right)+\left({\lambda}_{8}-{\lambda}_{9}\right)\left(\frac{b{\beta}_{V}}{{N}_{V}}{S}_{V}\right)\\ \text{}+\left({\lambda}_{6}-{\lambda}_{7}\right){u}_{4}{k}_{HM}+{\lambda}_{6}\left({\mu}_{H}+{\alpha}_{H}\right)\end{array}$

${\stackrel{\dot{}}{\lambda}}_{7}=-2{w}_{27}{R}_{H}+{\mu}_{H}{\lambda}_{7}$

${\stackrel{\dot{}}{\lambda}}_{8}=-2{w}_{28}{S}_{V}+\left({\lambda}_{8}-{\lambda}_{9}\right)\left(\frac{b{\beta}_{V}}{{N}_{V}}\left({I}_{HC}+{I}_{HF}+{I}_{HM}\right)\right)+{\lambda}_{8}\left({\mu}_{V}+{u}_{5}r\right)$

$\begin{array}{l}{\stackrel{\dot{}}{\lambda}}_{9}=-2{w}_{29}{I}_{V}+\left({\lambda}_{1}-{\lambda}_{4}\right)\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HC}}{{N}_{H}}{S}_{HC}\right)+\left({\lambda}_{2}-{\lambda}_{5}\right)\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HF}}{{N}_{H}}{S}_{HF}\right)\\ \text{}+\left({\lambda}_{3}-{\lambda}_{6}\right)\left(\left(1-{u}_{1}\right)b\frac{{\beta}_{HM}}{{N}_{H}}{S}_{HM}\right)+{\lambda}_{9}\left({\mu}_{V}+{\alpha}_{V}+{u}_{5}r\right)\end{array}$

andoptimality conditions are as

$\begin{array}{l}2{w}_{11}{u}_{1}+\left({\lambda}_{1}-{\lambda}_{4}\right)\left(b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}{S}_{HC}\right)+\left({\lambda}_{2}-{\lambda}_{5}\right)\left(b1\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}{S}_{HF}\right)\\ +\left({\lambda}_{3}-{\lambda}_{6}\right)\left(b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}{S}_{HM}\right)=0\end{array}$

$2{w}_{12}{u}_{2}+\left({\lambda}_{2}-{\lambda}_{5}\right)\frac{c{\beta}_{S}}{{N}_{H}}{S}_{HF}{I}_{HM}+\left({\lambda}_{3}-{\lambda}_{6}\right)\frac{c{\beta}_{S}}{{N}_{H}}{S}_{HM}{I}_{HF}=0$

$2{w}_{13}{u}_{3}-{\lambda}_{4}\delta {I}_{HF}=0$

$2{w}_{14}{u}_{4}+\left({\lambda}_{7}-{\lambda}_{4}\right){k}_{HC}{I}_{HC}+\left({\lambda}_{7}-{\lambda}_{5}\right){k}_{HF}{I}_{HF}+\left({\lambda}_{7}-{\lambda}_{6}\right){k}_{HM}{I}_{HM}=0$

$2{w}_{15}{u}_{5}-{\lambda}_{8}r{S}_{V}-{\lambda}_{9}r{I}_{V}=0$

Solving optimality conditions for optimal control and the property of control space $u$ give,

$\begin{array}{c}{u}_{1}^{\ast}\left(t\right)=\mathrm{max}\{\mathrm{min}\{1,\frac{1}{2{w}_{11}}[\left({\lambda}_{4}-{\lambda}_{1}\right)\left(b\frac{{\beta}_{HC}}{{N}_{H}}{I}_{V}{S}_{HC}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left({\lambda}_{5}-{\lambda}_{2}\right)\left(b1\frac{{\beta}_{HF}}{{N}_{H}}{I}_{V}{S}_{HF}\right)\text{\hspace{0.17em}}+\left({\lambda}_{6}-{\lambda}_{3}\right)\left(b\frac{{\beta}_{HM}}{{N}_{H}}{I}_{V}{S}_{HM}\right)]\},0\}\end{array}$ (10)

${u}_{2}^{\ast}\left(t\right)=\mathrm{max}\left\{\mathrm{min}\left\{1,\frac{1}{2{w}_{12}}\left({\lambda}_{5}-{\lambda}_{2}\right)\frac{c{\beta}_{S}}{{N}_{H}}{S}_{HF}{I}_{HM}+\left({\lambda}_{6}-{\lambda}_{3}\right)\frac{c{\beta}_{S}}{{N}_{H}}{S}_{HM}{I}_{HF}\right\},0\right\}$ (11)

${u}_{3}^{\ast}\left(t\right)=\mathrm{max}\left\{\mathrm{min}\left\{1,\frac{1}{2{w}_{13}}{\lambda}_{4}\delta {I}_{HF}\right\},0\right\}$ (12)

#Math_177# (13)

${u}_{5}^{\ast}\left(t\right)=\mathrm{max}\left\{\mathrm{min}\left\{1,\frac{1}{2{w}_{15}}\left({\lambda}_{8}r{S}_{V}+{\lambda}_{9}r{I}_{V}\right)\right\},0\right\}$ (14)

5. Numerical Simulations and Observations

For parametric values given in Table 1, the basic reproduction number ${R}_{0}=0.4757<1$ , $\text{trace}\left(J\right)=-4.88<0$ , $\mathrm{det}\left(J\right)=0.0014>0$ , $\text{trace}\left({J}^{\ast}\right)=-4.88<0$ and $\mathrm{det}\left({J}^{\ast}\right)=0.0041>0$ Above values verifies that non-linear system given by Equations (1) to (9) is asymptotically stable at disease free equilibrium point as well as endemic equilibrium point.

In Figure 2, global stability of human population for ${R}_{0}<1$ is shown. It can be seen that in the beginning susceptible human population increases but after 15 days it moves toward stability. Also, infected human population decreases initially but no prevention from vector bite, treatment, sexual activity, pregnancy, and no spraying increase infected population after 10 days. High disease induced death rate controls the infected human population.

In Figure 3, the effects of increase in prevention from mosquito bites on various infected human classes are analysed. It can be observed that prevention from vector bites in human population is an essential tool to control disease spread.

From Figure 4, it can be analysed that how the infected human population can be controlled by giving proper treatment on time. As there is no specific vaccination or treatment available for Zika disease cure, the treatment is considered as cure of symptoms of disease.

Figure 5 indicates that proper preventive steps to stop sexual transmission of disease taken at the time of disease out breaks will help to control disease spread.

From Figure 6, it can be observed that the preventions to stop pregnancy amongst infected female and effected area by ZIKV, will reduce the new child births having microcephaly and hence infected child population.

Figure 7 indicates that how the effective spraying can control infected mosquito population and hence finally spraying on mosquito will control disease spread indirectly in human population.

Thus to investigate effects of congenital transmission of ZIKV and its control

Figure 2. Global stability of human population for ${R}_{0}<1$ .

Figure 3. Effects of increase in prevention from mosquito bites on infected population.

Figure 4. Effects of increase in treatment on infected population.

by spraying on mosquitoes, preventions and treatment on human population, it is required to make the total cost associated with above controls minimise using optimal policy for all controls. (Figure 8)

Figure 5. Effects of increase in prevention from sexual contacts on infected population.

Figure 6. Effects of increase in prevention from pregnancy on infected population.

Figure 7. Effects of increase in spraying (insecticide) on infected mosquito population.

To minimise the total effective cost to control disease spray, the policy is to be designed in a way that during first week of disease outbreaks, 26% preventions on treatment and 40% prevention to stop pregnancy, 15% prevention on sexual activities amongst human, 8% spraying is required to have for 1% of prevention from vector bite.

Figure 8. Controls.

6. Conclusions

In this paper, the spread of ZIKV considered initially as vector-borne infection and then after its spread amongst sexual partners and from mother to child is analysed. The control over spraying on vectors, prevention from vector bite, prevention during sexual activity, prevention to stop pregnancy and prevention for treatment with time incorporated. It is observed that all controls disease effectively in each compartment.

In future, the model can be studied with addition of transmission of ZIKV by blood transmission in human population.

Acknowledgements

The authors thank DST-FIST file # MS1-097 for support to the department of Mathematics.

Cite this paper

Shah, N. , Patel, Z. and Yeolekar, B. (2017) Preventions and Controls on Congenital Transmissions of Zika: Mathematical Analysis.*Applied Mathematics*, **8**, 500-519. doi: 10.4236/am.2017.84040.

Shah, N. , Patel, Z. and Yeolekar, B. (2017) Preventions and Controls on Congenital Transmissions of Zika: Mathematical Analysis.

References

[1] Dick, G.W., Kitchen, S.F. and Haddow, A.J. (1952) Zika Virus (I). Isolations and Serological Specificity. Transaction of the Royal Society of Tropical Medicine and Hygiene, 46, 509-520.

https://doi.org/10.1016/0035-9203(52)90042-4

[2] Rodrigez, A.J. (2015) Zika: The New Arbovirus Threat for Latin America. The Journal of Infection in Developing Countries, 9, 684-685.

[3] Macnamara, F.N. (1954) Zika Virus: A Report on Three Cases of Human Infection during an Epidemic of Jaundice in Nigeria. Transaction of the Royal Society of Tropical Medicine and Hygiene, 48, 139-145.

https://doi.org/10.1016/0035-9203(54)90006-1

[4] Duffy, M.R., et al. (2009) Zika Virus Outbreak on Yap Island, Federated States of Micronesia. The New England Journal of Medicine, 360, 2536-2543.

https://doi.org/10.1056/NEJMoa0805715

[5] Musso, D., Nilles, E.J. and Cao-Lormeau, V.M. (2014) Rapid Spread of Emerging Zika Virus in the Pacific Area. Clinical Microbiology and Infection, 20, 595-596.

https://doi.org/10.1111/1469-0691.12707

[6] Pan American Health Organization (PAHO), Zika Virus Infection.

http://www.paho.org/hq/index.php?option=com

[7] Foy, B.D., et al. (2011) Probable Non-Vector-Borne Transmission of Zika Virus, Colorado, USA. Emerging Infectious Diseases, 17, 880-882.

https://doi.org/10.3201/eid1705.101939

[8] Dallas County Health and Human Services (DCHHS), DCHHS Reports First Zika Virus Case in Dallas County Acquired Through Sexual Transmission. 2 February 2016.

https://www.dallascounty.org/department/hhs/documents/DCHHS_Zika

[9] Centers for Disease Control and Prevention (CDC). Update: Interim Guidelines for Prevention of Sexual Transmission of Zika Virus—United States, 2016. 23 February 2016.

http://emergency.cdc.gov/han/han00388.asp

[10] Hills, S.L., et al. (2016) Transmission of Zika Virus through Sexual Contact with Travellers to Areas of Ongoing Transmission—Continental United States. Morbidity and Mortality Weekly Report, 65, 215-216.

https://doi.org/10.15585/mmwr.mm6508e2

[11] Mansuy, J.M., et al. (2016) Zika Virus: High Infectious Viral Load in Semen, a New Sexually Transmitted Pathogen? The Lancet Infectious Diseases, 16, 405.

https://doi.org/10.1016/S1473-3099(16)00138-9

[12] Toronto Star, Canada’s First Case of Sexually-Transmitted Zika Virus Confirmed in Ontario. 25 April 2016.

https://www.thestar.com/news/gta/2016/04/25/first-canadian-case-of-sexually-transmitted-zika-virus-detected-in-ontario.html

[13] Musso, D., Roche, C., Robin, E., Nhan, T., Teissier, A., Cao-Lormeau, V., et al. (2015) Potential Sexual Transmission of Zika Virus. Emerging Infectious Diseases, 21, 359-361. www.cdc.gov/eid

[14] Gatherer, D. and Kohl, A. (2015) Zika Virus: A Previously Slow Pandemic Spreads Rapidly through the Americas. Journal of General Virology, 97, 269-273.

https://doi.org/10.1099/jgv.0.000381

[15] Mlakar, J., et al. (2016) Zika Virus Associated with Microcephaly. The New England Journal of Medicine, 374, 951-958.

https://doi.org/10.1056/NEJMoa1600651

[16] Cauchemez, S., Besnard, M., Bompard, P., Dub, T., Guillemette-Artur, P., Eyrolle-Guignot, D., Salje, H., Van Kerkhove, M.D., Abadie, V., Garel, C., Fontanet, A. and Mallet, H.P. (2016) Association between Zika Virus and Microcephaly in French Polynesia, 2013-15: A Retrospective Study. The Lancet, 387, 2125-2132.

https://doi.org/10.1016/S0140-6736(16)00651-6

[17] Cao-Lormeau, V.M., Blake, A., Mons, S., Lastère, S., Roche, C., Vanhomwegen, J., Dub, T., Baudouin, L., Teissier, A., Larre, P., Vial, A.L., Decam, C., Choumet, V., Halstead, S.K., Willison, H.J., Musset, L., Manuguerra, J.C., Despres, P., Fournier, E., Mallet, H.P., Musso, D., Fontanet, A., Neil, J. and Ghawché, F. (2016) Guillain-Barré Syndrome Outbreak Associated with Zika Virus Infection in French Polynesia: A Case Control Study. The Lancet, 387, 1531-1539.

https://doi.org/10.1016/S0140-6736(16)00562-6

[18] WHO Report.

http://www.who.int/mediacentre/news/statements/2016/1st-emergency-committee-zika/en/

[19] Oliveira Melo, A.S., Malinger, G., Ximenes, R., Szejnfeld, P.O., Alves, S.S. and Bispo de Filippis, A.M. (2016) Zika Virus Intrauterine Infection Causes Fetal Brain Abnormality and Microcephaly: Tip of the Iceberg? Ultrasound in Obstetrics & Gynecology, 47, 6-7.

https://doi.org/10.1002/uog.15831

[20] Pizza, D.M.M., Loaiza, A.M., Montoya, J.F.A., Manrique, O.A., Sossa, V.A., Munoz, C.A.A., Raigosa, S., Contreras, H.M., Contreras, I.P., Perea, M.E.C. and Enríquez, M.A. (2016) A Model for the Risk of Microcephaly Induced by the Zika Virus (ZIKV). Open Journal of Modelling and Simulation, 4, 109-117.

https://doi.org/10.4236/ojmsi.2016.43010

[21] Bonyah, E. and Okosun, K.O. (2016) Mathematical Modeling of Zika Virus. Asian Pacific Journal of Tropical Disease, 6, 673-679.

[22] Gao, D., Lou, Y., He, D., Porco, T.C., Kuang, Y., Chowell, G. and Ruan, S. (2016) Prevention and Control of Zika as a Mosquito-Borne and Sexually Transmitted Disease: A Mathematical Modeling Analysis. Scientific Reports, 6, Article No. 28070.

[23] Pontryagin, L.S., Boltyanskii, V.G., Gamkrelidze, R.V. and Mishchenko, E.F. (1986) The Mathematical Theory of Optimal Process. Gordon and Breach Science Publishers, New York, 4.

[24] Van den Driessche, P. and Watmough, J. (2002) Reproduction Numbers and Sub-Threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Mathematical Biosciences, 180, 29-48.

[1] Dick, G.W., Kitchen, S.F. and Haddow, A.J. (1952) Zika Virus (I). Isolations and Serological Specificity. Transaction of the Royal Society of Tropical Medicine and Hygiene, 46, 509-520.

https://doi.org/10.1016/0035-9203(52)90042-4

[2] Rodrigez, A.J. (2015) Zika: The New Arbovirus Threat for Latin America. The Journal of Infection in Developing Countries, 9, 684-685.

[3] Macnamara, F.N. (1954) Zika Virus: A Report on Three Cases of Human Infection during an Epidemic of Jaundice in Nigeria. Transaction of the Royal Society of Tropical Medicine and Hygiene, 48, 139-145.

https://doi.org/10.1016/0035-9203(54)90006-1

[4] Duffy, M.R., et al. (2009) Zika Virus Outbreak on Yap Island, Federated States of Micronesia. The New England Journal of Medicine, 360, 2536-2543.

https://doi.org/10.1056/NEJMoa0805715

[5] Musso, D., Nilles, E.J. and Cao-Lormeau, V.M. (2014) Rapid Spread of Emerging Zika Virus in the Pacific Area. Clinical Microbiology and Infection, 20, 595-596.

https://doi.org/10.1111/1469-0691.12707

[6] Pan American Health Organization (PAHO), Zika Virus Infection.

http://www.paho.org/hq/index.php?option=com

[7] Foy, B.D., et al. (2011) Probable Non-Vector-Borne Transmission of Zika Virus, Colorado, USA. Emerging Infectious Diseases, 17, 880-882.

https://doi.org/10.3201/eid1705.101939

[8] Dallas County Health and Human Services (DCHHS), DCHHS Reports First Zika Virus Case in Dallas County Acquired Through Sexual Transmission. 2 February 2016.

https://www.dallascounty.org/department/hhs/documents/DCHHS_Zika

[9] Centers for Disease Control and Prevention (CDC). Update: Interim Guidelines for Prevention of Sexual Transmission of Zika Virus—United States, 2016. 23 February 2016.

http://emergency.cdc.gov/han/han00388.asp

[10] Hills, S.L., et al. (2016) Transmission of Zika Virus through Sexual Contact with Travellers to Areas of Ongoing Transmission—Continental United States. Morbidity and Mortality Weekly Report, 65, 215-216.

https://doi.org/10.15585/mmwr.mm6508e2

[11] Mansuy, J.M., et al. (2016) Zika Virus: High Infectious Viral Load in Semen, a New Sexually Transmitted Pathogen? The Lancet Infectious Diseases, 16, 405.

https://doi.org/10.1016/S1473-3099(16)00138-9

[12] Toronto Star, Canada’s First Case of Sexually-Transmitted Zika Virus Confirmed in Ontario. 25 April 2016.

https://www.thestar.com/news/gta/2016/04/25/first-canadian-case-of-sexually-transmitted-zika-virus-detected-in-ontario.html

[13] Musso, D., Roche, C., Robin, E., Nhan, T., Teissier, A., Cao-Lormeau, V., et al. (2015) Potential Sexual Transmission of Zika Virus. Emerging Infectious Diseases, 21, 359-361. www.cdc.gov/eid

[14] Gatherer, D. and Kohl, A. (2015) Zika Virus: A Previously Slow Pandemic Spreads Rapidly through the Americas. Journal of General Virology, 97, 269-273.

https://doi.org/10.1099/jgv.0.000381

[15] Mlakar, J., et al. (2016) Zika Virus Associated with Microcephaly. The New England Journal of Medicine, 374, 951-958.

https://doi.org/10.1056/NEJMoa1600651

[16] Cauchemez, S., Besnard, M., Bompard, P., Dub, T., Guillemette-Artur, P., Eyrolle-Guignot, D., Salje, H., Van Kerkhove, M.D., Abadie, V., Garel, C., Fontanet, A. and Mallet, H.P. (2016) Association between Zika Virus and Microcephaly in French Polynesia, 2013-15: A Retrospective Study. The Lancet, 387, 2125-2132.

https://doi.org/10.1016/S0140-6736(16)00651-6

[17] Cao-Lormeau, V.M., Blake, A., Mons, S., Lastère, S., Roche, C., Vanhomwegen, J., Dub, T., Baudouin, L., Teissier, A., Larre, P., Vial, A.L., Decam, C., Choumet, V., Halstead, S.K., Willison, H.J., Musset, L., Manuguerra, J.C., Despres, P., Fournier, E., Mallet, H.P., Musso, D., Fontanet, A., Neil, J. and Ghawché, F. (2016) Guillain-Barré Syndrome Outbreak Associated with Zika Virus Infection in French Polynesia: A Case Control Study. The Lancet, 387, 1531-1539.

https://doi.org/10.1016/S0140-6736(16)00562-6

[18] WHO Report.

http://www.who.int/mediacentre/news/statements/2016/1st-emergency-committee-zika/en/

[19] Oliveira Melo, A.S., Malinger, G., Ximenes, R., Szejnfeld, P.O., Alves, S.S. and Bispo de Filippis, A.M. (2016) Zika Virus Intrauterine Infection Causes Fetal Brain Abnormality and Microcephaly: Tip of the Iceberg? Ultrasound in Obstetrics & Gynecology, 47, 6-7.

https://doi.org/10.1002/uog.15831

[20] Pizza, D.M.M., Loaiza, A.M., Montoya, J.F.A., Manrique, O.A., Sossa, V.A., Munoz, C.A.A., Raigosa, S., Contreras, H.M., Contreras, I.P., Perea, M.E.C. and Enríquez, M.A. (2016) A Model for the Risk of Microcephaly Induced by the Zika Virus (ZIKV). Open Journal of Modelling and Simulation, 4, 109-117.

https://doi.org/10.4236/ojmsi.2016.43010

[21] Bonyah, E. and Okosun, K.O. (2016) Mathematical Modeling of Zika Virus. Asian Pacific Journal of Tropical Disease, 6, 673-679.

[22] Gao, D., Lou, Y., He, D., Porco, T.C., Kuang, Y., Chowell, G. and Ruan, S. (2016) Prevention and Control of Zika as a Mosquito-Borne and Sexually Transmitted Disease: A Mathematical Modeling Analysis. Scientific Reports, 6, Article No. 28070.

[23] Pontryagin, L.S., Boltyanskii, V.G., Gamkrelidze, R.V. and Mishchenko, E.F. (1986) The Mathematical Theory of Optimal Process. Gordon and Breach Science Publishers, New York, 4.

[24] Van den Driessche, P. and Watmough, J. (2002) Reproduction Numbers and Sub-Threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Mathematical Biosciences, 180, 29-48.