# Mathematical Analysis of Malaria-Schistosomiasis Coinfection Model

Authors: E. A. Bakare; C. R. Nwozo
Journal: Epidemiology Research International
 (2016)

---

## Abstract

We formulated and analysed a mathematical model to explore the cointeraction between malaria and schistosomiasis. Qualitative and comprehensive mathematical techniques have been applied to analyse the model. The local stability of the disease-free and endemic equilibrium was analysed, respectively. However, the main theorem shows that ifRMS<1, then the disease-free equilibrium is locally asymptotically stable and the phase will vanish out of the host and if RMS>1, a unique endemic equilibrium is also locally asymptotically stable and the disease persists at the endemic steady state. The impact of schistosomiasis and its treatment on malaria dynamics is also investigated. Numerical simulations using a set of reasonable parameter values show that the two epidemics coexist whenever their reproduction numbers exceed unity. Further, results of the full malaria-schistosomiasis model also suggest that an increase in the number of individuals infected with schistosomiasis in the presence of treatment results in a decrease in malaria cases. Sensitivity analysis was further carried out to investigate the influence of the model parameters on the transmission and spread of malaria-schistosomiasis coinfection. Numerical simulations were carried out to confirm our theoretical findings.


---

## Body

## 1. Introduction

Malaria is highly endemic in various parts of sub-Saharan Africa in which 85% of global malaria cases and 90% of malaria deaths occur [1].Schistosoma mansoni (the causative agent of intestinal schistosomiasis) is also prevalent in many sub-Saharan African countries [2, 3], accounting for approximately one-third of the total cases of schistosomiasis in the region [4]. The disease is a major contributor to disease burden globally and affects low income countries with climates suitable for transmission seriously. It is a life-threatening disease caused by parasites that are transmitted to people by the bites of infected mosquitoes [4, 5]. The bites by mosquitoes have resulted in the death of a child from malaria every 30 secs according to the report by the World Health Organization (WHO) fact sheet (2009) [4, 5].Plasmodium falciparum andPlasmodium vivax are the two common species andPlasmodium falciparum is the most deadly.Plasmodium falciparum malaria remains a major cause of mortality and morbidity in the tropics and subtropics areas of the globe [4, 6]. According to the 2009 world report, half of the world’s population is at risk of malaria, with an estimated 247 million cases that led to about 863,000 death in 2008 mostly among African children, a slight drop from 2006 statistics with the estimation that over 2000 young are lost every day across the globe [4]. This population made malaria the dominant parasitic disease of the tropics and one of the top three killer communicable diseases [4, 5, 7]. Malaria makes development to be very slow in several ways; it affects fertility, population growth, savings and investments, and worker productivity and causes absenteeism and premature mortality [4, 5, 7]. Malaria also affects fetal development during early stage of pregnancy in women due to loss of immunity. However, malaria is preventable and curable when treatment and prevention measures are sought early [4, 5, 7].The disease, schistosomiasis, also known as bilharziasis or snail fever, is a parasitic disease that was first namedbilharzia [8, 9] and it is prevalent in several regions of the developing world, predominantly Africa, South America, and Asia, with about 650 million people living in the endemic areas [10]. It is known that estimated 207 million people are infected, where 85 percent lives in underdeveloped areas of Africa [11], resulting in about 15,000 deaths annually presently [2]. Children below the age of 14 are the major victims of schistosomiasis infection in many parts of the world [10]. The basis of illness in victims is the eggs laid by the parasitic flat worms, that is, blood flukes of the genusSchistosoma [8]. The species of the water-borne flatworm or blood flukes known asschistosomes is the main type that initiates the human schistosomiasis, butSchistosoma mansoni,Schistosoma japonicum, andSchistosoma haematobium are the three major species that are found everywhere [8]. The urinary tract and kidneys as well as the reproductive systems are affected by theSchistosoma haematobium, and they are intense in Africa and the Middle East [8]. The most widely spread species is theSchistosoma mansoni whileSchistosoma japonicum is chiefly found in Asia and these two cause chronic hepatic and intestinal fibrosis [8, 10]. When skin comes in contact with contaminated freshwater in which certain types of snails that carry the parasite are living then the infection can be established [8]. Whenever infected people urinate or defecate in the water, freshwater becomes contaminated bySchistosoma eggs [8]. The eggs hatch, and the parasites infect, mature, and reproduce inside the snails when the appropriate species of snails exist in the water [12]. The parasite eventually leaves the snail and go into the water where it can persist for about 48 hours [8, 12].Schistosoma parasites, when wading, swimming, bathing, or washing, can enter the skin of anyone who comes in contact with contaminated freshwater [8, 12]. The parasites migrate through host tissue and develop into adult worms inside the blood vessels of the body for over numerous weeks [8, 12]. The worms mate and females produce eggs after maturity [8, 12]. Several of these eggs eventually travel to the bladder or intestine and are finally passed into the urine or stool [8, 12]. The schistosomiasis symptoms are caused by the body’s reaction to the eggs but not by the worms themselves [12, 13]. Eggs shed by the adult worms that do not pass out of the body can become lodged in the intestine or bladder, causing inflammation or scarring [12, 13]. Repeatedly infected children can acquire anemia, malnutrition, and learning difficulties [12, 13]. The parasite can as well damage the liver, intestine, spleen, lungs, and bladder even several years after infection [12, 13]. It is known at present that both malaria and intestinal schistosomiasis contribute to common epidemiological distributions and are currently posing a great task to public health and socioeconomic development throughout the tropical region [14]. The interactive pathology between malaria andS. mansoni has received increased investigation in the recent time, as a result of their coendemicities [1, 11, 15, 16]. It has been discovered that considerableS. mansoni infections are linked with a major increase in the incidence of malaria among school-age children [11]. In individuals infected withS. mansoni the technique responsible for the magnification of malaria is not yet fully understood [1, 9]. Thus, it is observed that the interface between the two diseases is perhaps set in motion by contradicting effects; the parasites possess the immunological cytokines; that is, the balance betweenTh1 andTh2 type immune responses which reduces immunological control of malaria may be altered byS. mansoni, while other methods are probable [1, 15, 17–19].It is our view that this study represents the very first modeling work that presents a mathematical analysis of the qualitative dynamics of malaria-schistosomiasis coinfection. There are few studies done on the malaria-schistosomiasis coinfection model so far. In [20], a coepidemic model of malaria andS. mansoni transmission dynamics is established, where the model reports major epidemiological coupling between the two diseases in terms of aggravated malaria incidence among individuals withS. mansoni extreme egg output. Their model was factored forS. mansoni extreme-risk endemic areas, applying epidemiological and clinical data of the relationship betweenS. mansoni and malaria among children in sub-Saharan Africa. They also assessed the potential influence of theS. mansoni malaria interface and mass treatment of schistosomiasis on malaria prevalence in coendemic areas.In this paper, we develop a mathematical model of the interplay between malaria andS. mansoni in which we have modeled the malaria transmission and theS. mansoni together as coendemic deterministic model. Our aim here is to study and analyse a mathematical model of malaria-schistosomiasis transmission model. Additionally, there are some important differences between the model in [16] and the one in this paper. This paper is organized as follows: we present a malaria-schistosomiasis coinfection transmission model formulation in Section 2, where the general mathematical framework, notations, and model equations were analysed with the basic properties of the models and their analysis. In Section 3, we present the existence of steady state solution. In Section 4, the basic reproduction number and stability were derived and carried out. Sensitivity analysis of the model was performed to determine the most important parameters that influence R0 in Section 5. In Section 6, we show our numerical simulation results while, in Section 7, we discussed our conclusions and recommendations.


## 2. Model Formulation

In this model, we denote the total human population byNh and subdivide it into the following subclasses of individuals who are susceptible (SH), individuals with malaria symptoms only (i.e., who are already infected and infective with malaria parasite) (Im), individuals infected with schistosomiasis only (Ihs), individuals infected with both malaria and schistosomiasis (Vms), individuals who recovered from malaria only (Rm), individuals who recovered from schistosomiasis only (Rhs), and individuals who recovered from both malaria and schistosomiasis such that Nh=SH+Im+Ihs+Vms+Rm+Rhs+Rms. The total snail population is denoted by Ns, which comprises susceptible snails (Ss) and infected as well as infectious snails (Is). That is, Ns=Ss+Is. The total mosquito population is denoted by Nv, which comprises susceptible mosquitoes (Sv) and infected as well as infectious mosquitoes (Iv). That is, Nv=Sv+Iv.The population of susceptible humans is generated through birth (at a constant per capita ratebH), by the loss of immunity to the malaria disease only (at a constant per capita rate γ), loss of immunity to the schistosomiasis disease only at a rate k, and loss of immunity to malaria and schistosomiasis disease at a rate φ. It is reduced by natural death (at a rate dH) and through the rate of acquiring malaria through contact with infectious mosquitoes (at a rate β1ϵhσ), where β1 is the transmission probability per bite, ϵh is the per capita biting rate of mosquitoes, and σ is the contact rate of mosquito per human per unit time. It is also reduced by rate of acquiring schistosomiasis through contact with infected snails (at a rate β2). Hence, the rate of change of population of susceptible humans is given by(1)dSHdt=bH+γRm+kRhs+φRms-β1ϵhσSHIv-β2SHIs-dHSH. The rate of change of the population of individuals with malaria only is increased by the rate of acquiring malaria through contact with infectious mosquitoes (at a rate β1ϵhσ) and by the rate of acquiring schistosomiasis through contact with infectious snail (at a rate β2) It is also reduced by human spontaneous recovery (at a rate ω). It is also reduced by the disease induced death rate (at per capita rate θ) and by the natural death rate (at per capita rate dH). Hence it is given by(2)dImdt=β1ϵhσSHIv-β2IsIm-ω+θ+dHIm.The rate of change of the population of individuals infected with schistosomiasis only is increased by the rate of acquiring schistosomiasis through contact with infectious snail (at a rateβ2) and decreased by infected mosquitoes (at a rate β1ϵhσ) and by human spontaneous recovery from schistosomiasis only (at a rate qH). It is also reduced by the disease induced death rate (at per capita rate ρH) and by the natural death rate (at per capita rate dH). Hence it is given by(3)dIhsdt=β2SHIs-β1ϵhσIvIhs-qH+ρH+dHIhs.The rate of change of the population of individuals infected with schistosomiasis and malaria is increased by the rate of acquiring malaria by infected mosquitoes (at a rateβ1ϵhσ) and schistosomiasis through contact with infectious snails (at a rate β2) and reduced by human spontaneous recovery from schistosomiasis only (at a rate α). It is also reduced by the malaria disease induced death rate (at per capita rate m) and schistosomiasis induced death rate (at per capita rate r) and by the natural death rate (at per capita rate dH). Hence it is given by(4)dVmsdt=β1ϵhσIvIhs+β2IsIm-α+r+m+dHVms.The individuals who recovered from malaria only is generated following a human spontaneous recovery (at a rateω) and by the dually infected individuals who recovered from malaria only at a rate (ϕ(1-α)) decreased by loss of immunity (at a rate γ) and by natural death (at a rate dH). Then(5)dRmdt=ωIm-γ+dHRm+ϕ1-αVms.The individuals who recovered from schistosomiasis only are generated following a human spontaneous recovery (at a rateqH) and by the dually infected individuals who recovered from schistosomiasis only at a rate ((1-ϕ)(1-α)) decreased by loss of immunity (at a rate γ) and by natural death (at a rate dH). Then(6)dRhsdt=qHIhs-k+dHRhs+1-ϕ1-αVms.The individuals who recovered from malaria and schistosomiasis are generated following a human spontaneous recovery (at a rateα) decreased by loss of immunity (at a rate φ) and by natural death (at a rate dH). Then(7)dRmsdt=αVms-φ+dHRms.Susceptible snail population is generated by the birth of snails (at a per capita rate ofbs). It is reduced by rate of acquiring schistosomiasis through contacts with infected humans at a rate β3, where η1 is a modification parameter. It is also reduced by natural death (at a rate ds). Thus,(8)dSsdt=bs-β3Ihs+η1VmsSs-dsSs. The population of infected snail is increased by rate of acquiring schistosomiasis through contacts with infected humans at a rate β3 and decreased by the natural death rate (at a rate ds), where η1 is a modification parameter. Hence, it is given by(9)dIsdt=β3Ihs+η1VmsSs-ds+ρsIs.Susceptible mosquito population is generated by the birth of mosquitoes (at a per capita rate ofbv). It is reduced by rate of acquiring malaria through contacts with infected humans at a rate β4ϵvσ, where β4 is probability for a vector (mosquito) to get infected by an infectious human, where η2 is a modification parameter. It is also reduced by natural death (at a rate dv). Thus,(10)dSvdt=bv-β4ϵvσIm+η2VmsSv-dvSv. The population of infected mosquito is increased by rate of acquiring malaria through contacts with infected humans at a rate β2ϵvσ and decreased by the natural death rate (at a rate dv), where η2 is a modification parameter. Hence, it is given by(11)dIvdt=β4ϵvσIm+η2VmsSv-dvIv.


### 2.1. The Full Schistosomiasis-Malaria Coinfection Model

Bringing the above formulation and assumptions together leads to the following set of ordinary differential equations which may be a new malaria-schistosomiasis coinfection model:(12)dSHdt=bH+γRm+kRhs+φRms-β1ϵhσSHIv-β2SHIs-dHSHdImdt=β1ϵhσSHIv-β2IsIm-ω+θ+dHImdIhsdt=β2SHIs-β1ϵhσIvIhs-qH+ρH+dHIhsdVmsdt=β1ϵhσIvIhs+β2IsIm-α+r+m+dHVmsdRmdt=θIh-γ+dHRm+ϕ1-αVmsdRhsdt=qHIhs-g+dHRhs+1-ϕ1-αVmsdRmsdt=αVms-φ+dHRmsdSsdt=bs-β3Ihs+η1VmsSs-dsSsdIsdt=β3Ihs+η1VmsSs-ds+ρsIsdSvdt=bv-β4ϵvσIm+η2VmsSv-dvSvdIvdt=β4ϵvσIm+η2VmsSv-dvIvsubject to the initial conditions SH(0)=SH,0, Im(0)=Im,0, Ihs(0)=Ihs,0, Vms(0)=Vms,0, Rm(0)=Rm,0, Rhs(0)=Rhs,0, Rms(0)=Rms,0, Ss(0)=Ss,0, Is(0)=Is,0, Sv(0)=Sv,0, and Iv(0)=Iv,0.We describe the associated model variables and parameters in the following list and Table1.Table 1
Table showing numerical values of sensitivity Indices.
Parameter Parameter values Sensitivity toRMS θ 0.05 - 0.00000000000186 β 1 0.8333 + 1.0 β 2 0.406 + 1.0 β 3 0.004 + 1.0 β 4 0.09 + 0.99 b v 1000 + 0.74 d H 0.00004 - 2.0 ϵ v 0.2 + 1.0 σ 0.502 + 2.0 ϵ h 0.2 + 1.0 q H 0.56 - 0.99 ρ S 0.0004012 - 0.876 ρ H 0.0039 - 0.0069 d v 0.1429 - 0.23 b s 200 + 1.0 b H 100 + 2.0 d s 0.0000569 - 0.11 ω 0.05 - 0.91Variables in the Model S H: susceptible human I m: human infected with malaria only I h s: human infected with schistosomiasis only V m s: human infected with malaria and schistosomiasis only R m: human recovered from malaria only R h s: human recovered from schistosomiasis only R m s: human recovered from malaria and schistosomiasis only S s: susceptible snail I s: infected snail S v: susceptible mosquito I v: infected mosquito


### 2.2. Qualitative Analysis of the Full Schistosomiasis-Malaria Coinfection Model

The schistosomiasis-malaria coinfection model (12) will be analysed in a biologically feasible region for both humans, snail, and mosquito populations. Hence, for it to be epidemiologically well posed, it is necessary to show that all its state variables are nonnegative for all time t>0.


### 2.3. Positivity of the Solution

Theorem 1.
If the initial value forSH≥0, Im≥0, Ihs≥0, Vms≥0, Rm≥0, Rhs≥0, Rms≥0, Ss≥0, Is≥0, Sv≥0, and Iv≥0 then the solutions (SH(t),Im(t),Ihs(t),Rhs(t),Vms(t),Rm(t),Rhs(t),Rms(t),Ss(t),Is(t),Sv(t),Iv(t)) of the schistosomiasis only model (16) are nonnegative for all t>0.Proof.
We letΓ=sup⁡{t>0:SHt>0,Imt>0,Ihst>0, Vmst>0,Rmt>0,Rhst>0,Rmst>0,Sst>0,Ist>0,Svt>0,Ivt>0}.
Since the variablesSH(0)>0, Im(0)>0, Ihs(0)>0, Vms(0)>0, Rm(0)>0, Rhs(0)>0, Rms(0)>0, Ss(0)>0, Is(0)>0, Sv(0)>0, and Iv(0)>0 then, Γ>0. If Γ<∞, then SH,Im,Ihs,Vms,Rm,Rhs,Rms,Ss,Is,Sv,Iv are equal to zero at Γ. It follows from the first equation of the system (12) that (13)dSHdt=bH+γRm+kRhs+φRms-β1ϵhσSHIv-β2SHIs-dHSH.Therefore,(14)ddtSHtexp⁡β1ϵhσIv+β2Is+dHt=bH+kRhs+γRm+φRmsexp⁡β1ϵhσIv+β2Is+dHt.Hence(15)SHΓexp⁡β1ϵhσIv+β2Is+dHt-SH0=∫0ΓbH+kRhs+γRm+φRmsexp⁡β1ϵhσIv+β2Is+dHzdz,(16)SHΓ=SH0exp⁡-β1ϵhσIv+β2Is+dHΓ+exp⁡-β1ϵhσIv+β2Is+dHΓ∫0ΓbH+kRhs+γRm+φRmsexp⁡β1ϵhσIv+β2Is+dHydy≥0.Also(17)dImdt=β1ϵhσSHIv-β2IsIm-ω+θ+dHIm.Now,(18)ddtImtexp⁡β2Is+ω+θ+dHt=∫0Γβ1ϵhσSHIvexp⁡β2Is+ω+θ+dHydy,ImΓ1=Im0exp⁡-β2Is+ω+θ+dHΓ+exp⁡-β2Is+ω+θ+dHΓ∫0Γβ1ϵhσSHIvexp⁡β2Is+ω+θ+dHydy≥0.
And then(19)dIhsdt=β2SHIs-β1ϵhσIvIhs-qH+ρH+dHIhs.
Now,(20)ddtIhstexp⁡qH+ρH+dHt=∫0Γβ2SHIsexp⁡β1ϵhσIv+qH+ρH+dHydy,IhsΓ=Ihs0exp⁡-β1ϵhσIv+qH+ρH+dHΓ+exp⁡-β1ϵhσIv+qH+ρH+dHΓ∫0Γβ2SHIsexp⁡β1ϵhσIv+qH+ρH+dHydy≥0.
Furthermore,(21)dVmsdt=β1ϵhσIvIhs+β2IsIm-α+r+m+dHVms.
Now,(22)ddtVmstexp⁡α+r+m+dHt=∫0Γβ1ϵhσIvIhs+β2IsImexp⁡α+r+m+dHydy,VmsΓ=Vms0exp⁡-α+r+m+dHΓ+exp⁡-α+r+m+dHΓ∫0Γβ1ϵhσIvIhs+β5IhsImexp⁡α+r+m+dHydy≥0.
Also(23)dRmdt=θIh-γ+dHRm+ϕ1-αVms.
Therefore(24)ddtRmtexp⁡γ+dHt=∫0ΓθIh+ϕ1-αVmsexp⁡γ+dHydy. Hence (25)RmΓexp⁡γ+dHt-Rhs0=∫0Γ1θIh+ϕ1-αVmsexp⁡γ+dHydy so that(26)RmΓ=Rhs0exp⁡-γ+dHΓ+exp⁡-γ+dHΓ∫0ΓθIh+ϕ1-αVmsexp⁡γ+dHydy≥0.
Also,(27)dRhsdt=qHIhs-k+dHRhs+1-ϕ1-αVms. Therefore (28)ddtRhstexp⁡k+dHt=∫0ΓqHIhs+1-ϕ1-αVmsexp⁡k+dHydy. Hence (29)RhsΓexp⁡k+dHt-Rhs0=∫0ΓqHIhs+1-ϕ1-αVmsexp⁡k+dHydyso that(30)RhsΓ=Rhs0exp⁡-k+dHΓ+exp⁡-k+dHΓ1∫0ΓqHIhs+1-ϕ1-αVmsexp⁡k+dHydy≥0.
Also,(31)dRmsdt=αVms-φ+dHRms.
Then,(32)ddtRmstexp⁡φ+dHt=∫0ΓαVmsexp⁡φ+dHydy. Hence (33)RhsΓexp⁡φ+dHt-Rms0=∫0ΓαVmsexp⁡φ+dHydy so that(34)RmsΓ=Rms0exp⁡-φ+dHΓ+exp⁡-φ+dHΓ∫0ΓαVmsexp⁡φ+dHydy≥0.
Also,(35)dSsdt=bs-β3Ihs+η1VmsSs-dsSs. Therefore, (36)Sstexp⁡β3Ihs+η1Vms+dst=∫0Γbvexp⁡β3Ihs+η1Vms-dsydy. Then, we have(37)SsΓ=Ss0exp⁡-β3Ihs+η1Vms+dsΓ+exp⁡-β3Ihs+η1Vms+dsΓ∫0Γbsexp⁡β3Ihs+η1Vms+dsydy≥0. Similarly, (38)dIsdt=β3Ihs+η1VmsSs-ds+ρsIs. Therefore, (39)ddtIstexp⁡ds+ρst=β3Ihs+η1Vmsexp⁡ds+ρsydy so that(40)IsΓ=Is0exp⁡-ds+ρsΓ+exp⁡-ds+ρsΓ∫0Γβ3Ihs+η1Vmsexp⁡ds+ρsydy≥0for all t>0.
Also,(41)dSvdt=bv-β4ϵvσIm+η2VmsSv-dvSv. Therefore, (42)Svtexp⁡β4ϵvσIm+η2Vms+dvt=∫0Γbvexp⁡β4ϵvσIm+η2Vms+dvzdz. Then, we have (43)SvΓ=Sv0exp⁡-β4ϵvσIm+η2Vms+dvΓ+exp⁡-β4ϵvσIm+η2Vms+dvΓ∫0Γbvexp⁡β4ϵvσIm+η2Vms+dvzdz≥0.Similarly, (44)dIvdt=β4ϵvσIm+η2VmsSv-dvIv. Therefore, (45)ddtIvtexp⁡dvt=β4ϵvσIm+η2Vmsexp⁡dvzdz so that(46)IvΓ=Iv0exp⁡-dvΓ+exp⁡-dvΓ∫0Γβ4ϵvσIm+η2Vmsexp⁡dvzdz≥0for all t>0.


### 2.4. Boundedness of the Solution

Theorem 2.
All solutions(SH(t),Im(t),Ihs(t),Vms(t),Rm(t),Rhs(t),Rms(t),Ss(t),Is(t),Sv(t),Iv(t)) of the malaria-schistosomiasis coinfection model (12) are bounded. Therefore, from (12) if (47)limt→∞⁡sup⁡Nht≤bHdH,limt→∞⁡sup⁡Nst≤bsds,limt→∞⁡sup⁡Nvt≤bvdvthen Ns=Ss+Is, Nv=Sv+Iv, and Nh=SH+Im+Ihs+Vms+Rhs+Rm+Rhs+Rms.Proof.
For the proof of boundedness, we note that0<Im(t)≤Nh(t), 0<Ihs(t)≤Nh(t), 0<Vms(t)≤Nh(t), 0<Is(t)≤Ns(t), and 0<Iv(t)≤Nv(t). We add the first eight equations and the last two equations of the malaria-schistosomiasis coinfection model (12) and yield(48)dNhdt=bh-dHNh-ωIm-ρHIhs-r+mVmsdNsdt=bs-dsNs-ρsIsdNvdt=bv-dvNv.All solutions of model (12) are bounded. The feasible region for the human population is given by Γh={SH,Im,Ihs,Vms,Rhs,Rm,Rhs,Rms∣SH + Im+Ihs+Vms+Rhs+Rm+Rhs+Rms≤bH/dH,0 ≤ SH≤SH(t)≤bH/dH,Ihs≥0,Rhs≥0}. And the feasible region for the snail population is given by(49)Γs=Ss,Is∣Ss+Is≤bsds,0≤Ss≤Sst≤bsds,Is≥0.And the feasible region for the mosquito population is given by(50)Γv=Sv,Iv∣Sv+Iv≤bvdv,0≤Sv≤Svt≤bvdv,Iv≥0.Therefore,(51)bH-dHNh-ω+ρH+r+mNht≤dNHtdt≤bH-dHNhtbs-ds+ρsNst≤dNstdt≤bs-dsNstbv-dv+ρvNvt≤dNvtdt≤bv-dvNvt. Hence, (52)bHdH+ρH≤liminft→∞⁡Nht≤limsupt→∞⁡Nht≤bHdH,bsds+ρs≤liminft→∞⁡Nst≤limsupt→∞⁡Nst≤bsdsbvdv≤liminft→∞⁡Nvt≤limsupt→∞⁡Nvt≤bvdv.


### 2.5. Invariant Region

Here, we analysed the malaria-schistosomiasis only model (12) in a biologically feasible region. Hence the system of (12) is split into three parts, namely. We have the human population (Nh: with Nh=SH+Im+Ihs+Vms+Rhs+Rm+Rhs+Rms), the snail (vector) population (Ns: with Ns=Ss+Is), and the mosquito (vector) population (Nv: with Nv=Sv+Iv). Let us consider the feasible region as (53)Γ=Γh∪Γs∪Γv⊂R+7×R+2×R+2 with (54)Γh=SH,Im,Ihs,Vms,Rhs,Rm,Rhs,Rms∈R+7:SH+Im+Ihs+Vms+Rhs+Rm+Rhs+Rms≤bHdHΓs=Ss,Is∈R+2:Ss+Is≤bsdsΓv=Sv,Iv∈R+2:Sv+Iv≤bvdv. We take the following steps to reaffirm the positive invariance of Γ (i.e., solution in Γ remains in Γ for all t>0). The rate of change of humans, snails, and mosquitoes population is given in (12) and it follows that (55)dNhtdt≤bH-dHNhdNstdt≤bs-dsNsdNvtdt≤bv-dvNv. By the standard comparison theorem we obtain (56)Nht≤Nh0e-dHt+bHdH1-e-dHtNst≤Ns0e-dst+bsds1-e-dstNvt≤Nv0e-dvt+bvdv1-e-dvt.In particular, Nh(t)≤bH/dH if Nh(0)≤bH/dH, Ns(t)≤bs/ds if Ns(0)≤bs/ds, and Nv(t)≤bv/dv if Nv(0)≤bv/dv. Therefore, the region Γ is positively invariant. Hence, it is enough to consider the dynamics of the flow generated by (12) in Γ. Thus, the model can be considered to be epidemiologically and mathematically well posed in this region. Therefore every solution of the malaria-schistosomiasis coinfection model (12) with initial data (conditions) in Γ remains in Γ. The summary of this result is given below.Theorem 3.
The regionΓ=Γh∪Γs∪Γv⊂R+7×R+2×R+2 is positively invariant for the malaria-schistosomiasis coinfection model (12) with nonnegative initial conditions in R+11.


## 2.1. The Full Schistosomiasis-Malaria Coinfection Model

Bringing the above formulation and assumptions together leads to the following set of ordinary differential equations which may be a new malaria-schistosomiasis coinfection model:(12)dSHdt=bH+γRm+kRhs+φRms-β1ϵhσSHIv-β2SHIs-dHSHdImdt=β1ϵhσSHIv-β2IsIm-ω+θ+dHImdIhsdt=β2SHIs-β1ϵhσIvIhs-qH+ρH+dHIhsdVmsdt=β1ϵhσIvIhs+β2IsIm-α+r+m+dHVmsdRmdt=θIh-γ+dHRm+ϕ1-αVmsdRhsdt=qHIhs-g+dHRhs+1-ϕ1-αVmsdRmsdt=αVms-φ+dHRmsdSsdt=bs-β3Ihs+η1VmsSs-dsSsdIsdt=β3Ihs+η1VmsSs-ds+ρsIsdSvdt=bv-β4ϵvσIm+η2VmsSv-dvSvdIvdt=β4ϵvσIm+η2VmsSv-dvIvsubject to the initial conditions SH(0)=SH,0, Im(0)=Im,0, Ihs(0)=Ihs,0, Vms(0)=Vms,0, Rm(0)=Rm,0, Rhs(0)=Rhs,0, Rms(0)=Rms,0, Ss(0)=Ss,0, Is(0)=Is,0, Sv(0)=Sv,0, and Iv(0)=Iv,0.We describe the associated model variables and parameters in the following list and Table1.Table 1
Table showing numerical values of sensitivity Indices.
Parameter Parameter values Sensitivity toRMS θ 0.05 - 0.00000000000186 β 1 0.8333 + 1.0 β 2 0.406 + 1.0 β 3 0.004 + 1.0 β 4 0.09 + 0.99 b v 1000 + 0.74 d H 0.00004 - 2.0 ϵ v 0.2 + 1.0 σ 0.502 + 2.0 ϵ h 0.2 + 1.0 q H 0.56 - 0.99 ρ S 0.0004012 - 0.876 ρ H 0.0039 - 0.0069 d v 0.1429 - 0.23 b s 200 + 1.0 b H 100 + 2.0 d s 0.0000569 - 0.11 ω 0.05 - 0.91Variables in the Model S H: susceptible human I m: human infected with malaria only I h s: human infected with schistosomiasis only V m s: human infected with malaria and schistosomiasis only R m: human recovered from malaria only R h s: human recovered from schistosomiasis only R m s: human recovered from malaria and schistosomiasis only S s: susceptible snail I s: infected snail S v: susceptible mosquito I v: infected mosquito


## 2.2. Qualitative Analysis of the Full Schistosomiasis-Malaria Coinfection Model

The schistosomiasis-malaria coinfection model (12) will be analysed in a biologically feasible region for both humans, snail, and mosquito populations. Hence, for it to be epidemiologically well posed, it is necessary to show that all its state variables are nonnegative for all time t>0.


## 2.3. Positivity of the Solution

Theorem 1.
If the initial value forSH≥0, Im≥0, Ihs≥0, Vms≥0, Rm≥0, Rhs≥0, Rms≥0, Ss≥0, Is≥0, Sv≥0, and Iv≥0 then the solutions (SH(t),Im(t),Ihs(t),Rhs(t),Vms(t),Rm(t),Rhs(t),Rms(t),Ss(t),Is(t),Sv(t),Iv(t)) of the schistosomiasis only model (16) are nonnegative for all t>0.Proof.
We letΓ=sup⁡{t>0:SHt>0,Imt>0,Ihst>0, Vmst>0,Rmt>0,Rhst>0,Rmst>0,Sst>0,Ist>0,Svt>0,Ivt>0}.
Since the variablesSH(0)>0, Im(0)>0, Ihs(0)>0, Vms(0)>0, Rm(0)>0, Rhs(0)>0, Rms(0)>0, Ss(0)>0, Is(0)>0, Sv(0)>0, and Iv(0)>0 then, Γ>0. If Γ<∞, then SH,Im,Ihs,Vms,Rm,Rhs,Rms,Ss,Is,Sv,Iv are equal to zero at Γ. It follows from the first equation of the system (12) that (13)dSHdt=bH+γRm+kRhs+φRms-β1ϵhσSHIv-β2SHIs-dHSH.Therefore,(14)ddtSHtexp⁡β1ϵhσIv+β2Is+dHt=bH+kRhs+γRm+φRmsexp⁡β1ϵhσIv+β2Is+dHt.Hence(15)SHΓexp⁡β1ϵhσIv+β2Is+dHt-SH0=∫0ΓbH+kRhs+γRm+φRmsexp⁡β1ϵhσIv+β2Is+dHzdz,(16)SHΓ=SH0exp⁡-β1ϵhσIv+β2Is+dHΓ+exp⁡-β1ϵhσIv+β2Is+dHΓ∫0ΓbH+kRhs+γRm+φRmsexp⁡β1ϵhσIv+β2Is+dHydy≥0.Also(17)dImdt=β1ϵhσSHIv-β2IsIm-ω+θ+dHIm.Now,(18)ddtImtexp⁡β2Is+ω+θ+dHt=∫0Γβ1ϵhσSHIvexp⁡β2Is+ω+θ+dHydy,ImΓ1=Im0exp⁡-β2Is+ω+θ+dHΓ+exp⁡-β2Is+ω+θ+dHΓ∫0Γβ1ϵhσSHIvexp⁡β2Is+ω+θ+dHydy≥0.
And then(19)dIhsdt=β2SHIs-β1ϵhσIvIhs-qH+ρH+dHIhs.
Now,(20)ddtIhstexp⁡qH+ρH+dHt=∫0Γβ2SHIsexp⁡β1ϵhσIv+qH+ρH+dHydy,IhsΓ=Ihs0exp⁡-β1ϵhσIv+qH+ρH+dHΓ+exp⁡-β1ϵhσIv+qH+ρH+dHΓ∫0Γβ2SHIsexp⁡β1ϵhσIv+qH+ρH+dHydy≥0.
Furthermore,(21)dVmsdt=β1ϵhσIvIhs+β2IsIm-α+r+m+dHVms.
Now,(22)ddtVmstexp⁡α+r+m+dHt=∫0Γβ1ϵhσIvIhs+β2IsImexp⁡α+r+m+dHydy,VmsΓ=Vms0exp⁡-α+r+m+dHΓ+exp⁡-α+r+m+dHΓ∫0Γβ1ϵhσIvIhs+β5IhsImexp⁡α+r+m+dHydy≥0.
Also(23)dRmdt=θIh-γ+dHRm+ϕ1-αVms.
Therefore(24)ddtRmtexp⁡γ+dHt=∫0ΓθIh+ϕ1-αVmsexp⁡γ+dHydy. Hence (25)RmΓexp⁡γ+dHt-Rhs0=∫0Γ1θIh+ϕ1-αVmsexp⁡γ+dHydy so that(26)RmΓ=Rhs0exp⁡-γ+dHΓ+exp⁡-γ+dHΓ∫0ΓθIh+ϕ1-αVmsexp⁡γ+dHydy≥0.
Also,(27)dRhsdt=qHIhs-k+dHRhs+1-ϕ1-αVms. Therefore (28)ddtRhstexp⁡k+dHt=∫0ΓqHIhs+1-ϕ1-αVmsexp⁡k+dHydy. Hence (29)RhsΓexp⁡k+dHt-Rhs0=∫0ΓqHIhs+1-ϕ1-αVmsexp⁡k+dHydyso that(30)RhsΓ=Rhs0exp⁡-k+dHΓ+exp⁡-k+dHΓ1∫0ΓqHIhs+1-ϕ1-αVmsexp⁡k+dHydy≥0.
Also,(31)dRmsdt=αVms-φ+dHRms.
Then,(32)ddtRmstexp⁡φ+dHt=∫0ΓαVmsexp⁡φ+dHydy. Hence (33)RhsΓexp⁡φ+dHt-Rms0=∫0ΓαVmsexp⁡φ+dHydy so that(34)RmsΓ=Rms0exp⁡-φ+dHΓ+exp⁡-φ+dHΓ∫0ΓαVmsexp⁡φ+dHydy≥0.
Also,(35)dSsdt=bs-β3Ihs+η1VmsSs-dsSs. Therefore, (36)Sstexp⁡β3Ihs+η1Vms+dst=∫0Γbvexp⁡β3Ihs+η1Vms-dsydy. Then, we have(37)SsΓ=Ss0exp⁡-β3Ihs+η1Vms+dsΓ+exp⁡-β3Ihs+η1Vms+dsΓ∫0Γbsexp⁡β3Ihs+η1Vms+dsydy≥0. Similarly, (38)dIsdt=β3Ihs+η1VmsSs-ds+ρsIs. Therefore, (39)ddtIstexp⁡ds+ρst=β3Ihs+η1Vmsexp⁡ds+ρsydy so that(40)IsΓ=Is0exp⁡-ds+ρsΓ+exp⁡-ds+ρsΓ∫0Γβ3Ihs+η1Vmsexp⁡ds+ρsydy≥0for all t>0.
Also,(41)dSvdt=bv-β4ϵvσIm+η2VmsSv-dvSv. Therefore, (42)Svtexp⁡β4ϵvσIm+η2Vms+dvt=∫0Γbvexp⁡β4ϵvσIm+η2Vms+dvzdz. Then, we have (43)SvΓ=Sv0exp⁡-β4ϵvσIm+η2Vms+dvΓ+exp⁡-β4ϵvσIm+η2Vms+dvΓ∫0Γbvexp⁡β4ϵvσIm+η2Vms+dvzdz≥0.Similarly, (44)dIvdt=β4ϵvσIm+η2VmsSv-dvIv. Therefore, (45)ddtIvtexp⁡dvt=β4ϵvσIm+η2Vmsexp⁡dvzdz so that(46)IvΓ=Iv0exp⁡-dvΓ+exp⁡-dvΓ∫0Γβ4ϵvσIm+η2Vmsexp⁡dvzdz≥0for all t>0.


## 2.4. Boundedness of the Solution

Theorem 2.
All solutions(SH(t),Im(t),Ihs(t),Vms(t),Rm(t),Rhs(t),Rms(t),Ss(t),Is(t),Sv(t),Iv(t)) of the malaria-schistosomiasis coinfection model (12) are bounded. Therefore, from (12) if (47)limt→∞⁡sup⁡Nht≤bHdH,limt→∞⁡sup⁡Nst≤bsds,limt→∞⁡sup⁡Nvt≤bvdvthen Ns=Ss+Is, Nv=Sv+Iv, and Nh=SH+Im+Ihs+Vms+Rhs+Rm+Rhs+Rms.Proof.
For the proof of boundedness, we note that0<Im(t)≤Nh(t), 0<Ihs(t)≤Nh(t), 0<Vms(t)≤Nh(t), 0<Is(t)≤Ns(t), and 0<Iv(t)≤Nv(t). We add the first eight equations and the last two equations of the malaria-schistosomiasis coinfection model (12) and yield(48)dNhdt=bh-dHNh-ωIm-ρHIhs-r+mVmsdNsdt=bs-dsNs-ρsIsdNvdt=bv-dvNv.All solutions of model (12) are bounded. The feasible region for the human population is given by Γh={SH,Im,Ihs,Vms,Rhs,Rm,Rhs,Rms∣SH + Im+Ihs+Vms+Rhs+Rm+Rhs+Rms≤bH/dH,0 ≤ SH≤SH(t)≤bH/dH,Ihs≥0,Rhs≥0}. And the feasible region for the snail population is given by(49)Γs=Ss,Is∣Ss+Is≤bsds,0≤Ss≤Sst≤bsds,Is≥0.And the feasible region for the mosquito population is given by(50)Γv=Sv,Iv∣Sv+Iv≤bvdv,0≤Sv≤Svt≤bvdv,Iv≥0.Therefore,(51)bH-dHNh-ω+ρH+r+mNht≤dNHtdt≤bH-dHNhtbs-ds+ρsNst≤dNstdt≤bs-dsNstbv-dv+ρvNvt≤dNvtdt≤bv-dvNvt. Hence, (52)bHdH+ρH≤liminft→∞⁡Nht≤limsupt→∞⁡Nht≤bHdH,bsds+ρs≤liminft→∞⁡Nst≤limsupt→∞⁡Nst≤bsdsbvdv≤liminft→∞⁡Nvt≤limsupt→∞⁡Nvt≤bvdv.


## 2.5. Invariant Region

Here, we analysed the malaria-schistosomiasis only model (12) in a biologically feasible region. Hence the system of (12) is split into three parts, namely. We have the human population (Nh: with Nh=SH+Im+Ihs+Vms+Rhs+Rm+Rhs+Rms), the snail (vector) population (Ns: with Ns=Ss+Is), and the mosquito (vector) population (Nv: with Nv=Sv+Iv). Let us consider the feasible region as (53)Γ=Γh∪Γs∪Γv⊂R+7×R+2×R+2 with (54)Γh=SH,Im,Ihs,Vms,Rhs,Rm,Rhs,Rms∈R+7:SH+Im+Ihs+Vms+Rhs+Rm+Rhs+Rms≤bHdHΓs=Ss,Is∈R+2:Ss+Is≤bsdsΓv=Sv,Iv∈R+2:Sv+Iv≤bvdv. We take the following steps to reaffirm the positive invariance of Γ (i.e., solution in Γ remains in Γ for all t>0). The rate of change of humans, snails, and mosquitoes population is given in (12) and it follows that (55)dNhtdt≤bH-dHNhdNstdt≤bs-dsNsdNvtdt≤bv-dvNv. By the standard comparison theorem we obtain (56)Nht≤Nh0e-dHt+bHdH1-e-dHtNst≤Ns0e-dst+bsds1-e-dstNvt≤Nv0e-dvt+bvdv1-e-dvt.In particular, Nh(t)≤bH/dH if Nh(0)≤bH/dH, Ns(t)≤bs/ds if Ns(0)≤bs/ds, and Nv(t)≤bv/dv if Nv(0)≤bv/dv. Therefore, the region Γ is positively invariant. Hence, it is enough to consider the dynamics of the flow generated by (12) in Γ. Thus, the model can be considered to be epidemiologically and mathematically well posed in this region. Therefore every solution of the malaria-schistosomiasis coinfection model (12) with initial data (conditions) in Γ remains in Γ. The summary of this result is given below.Theorem 3.
The regionΓ=Γh∪Γs∪Γv⊂R+7×R+2×R+2 is positively invariant for the malaria-schistosomiasis coinfection model (12) with nonnegative initial conditions in R+11.


## 3. Existence of Steady States Solution

We analysed the model equation (12) in this section qualitatively to investigate the condition of existence of equilibrium points. We would like to know what will eventually happen to the disease in the long run. The question that will arise is (i) will the disease (malaria-schistosomiasis disease) die out? Or (ii) will it be established in the population and become endemic? In order to answer these questions we have to investigate the long-term behaviour of the solutions. This behaviour depends largely on the equilibrium points, that is, time-independent solutions of the system. Since these solutions do not depend on time, we set (57)dSHdt=dImdt=dIhsdt=dVmsdt=dRmdt=dRhsdt=dRmsdt=0dSsdt=dIsdt=dSvdt=dIvdt=0. Therefore, we have the equilibrium solution (58)DMS=SH=bHdH,Im=0,Ihs=0,Vms=0,Rhs=0,Rm=0,Rhs=0,Rms=0,Ss=bsds,Is=0,Sv=bvdv,Iv=0.This equilibrium exists for all values of the parameters. Since Ihs=Rhs=Is=0 it implies that the disease will disappear from the population. Thus, this equilibrium is referred to as the disease-free equilibrium (DFE).


### 3.1. Existence of Disease-Free Equilibrium Point (DFE)

We analysed the system qualitatively by studying the system of equation in closed set(59)Ω=SH,Im,Ihs,Vms,Rhs,Rm,Rhs,Rms,Ss,Is,Sv,Iv∈R+11∣0≤SH≤Nh,0≤Ihs≤Nh,0≤Rhs≤Nh,0≤Ss≤Ns,0≤Is≤Ns,0≤Sv≤Nv,0≤Iv≤Nv.A qualitative study of the system exists in two forms, namely,(i)
the disease-free (dies out),(ii)
endemic.When the malaria-schistosomiasis disease dies out naturally, the solution asymptotically tends to a disease-free equilibriumDMS of the form DMS=(SH=bH/dH,Im=0,Ihs=0,Vms=0,Rhs = 0,Rm=0,Rhs=0,Rms=0,Ss=bs/ds,Is=0,Sv=bv/dv,Iv=0). Hence, the threshold that determines the stability of this equilibrium is the reproduction number (which is the average number of secondary infections caused by a single infectious individual during their entire infectious lifetime).


## 3.1. Existence of Disease-Free Equilibrium Point (DFE)

We analysed the system qualitatively by studying the system of equation in closed set(59)Ω=SH,Im,Ihs,Vms,Rhs,Rm,Rhs,Rms,Ss,Is,Sv,Iv∈R+11∣0≤SH≤Nh,0≤Ihs≤Nh,0≤Rhs≤Nh,0≤Ss≤Ns,0≤Is≤Ns,0≤Sv≤Nv,0≤Iv≤Nv.A qualitative study of the system exists in two forms, namely,(i)
the disease-free (dies out),(ii)
endemic.When the malaria-schistosomiasis disease dies out naturally, the solution asymptotically tends to a disease-free equilibriumDMS of the form DMS=(SH=bH/dH,Im=0,Ihs=0,Vms=0,Rhs = 0,Rm=0,Rhs=0,Rms=0,Ss=bs/ds,Is=0,Sv=bv/dv,Iv=0). Hence, the threshold that determines the stability of this equilibrium is the reproduction number (which is the average number of secondary infections caused by a single infectious individual during their entire infectious lifetime).


## 4. The Basic Reproduction NumberRMS

WheneverRMS<1, each individual releases on average less than one infected individual and therefore the disease (malaria-schistosomiasis only) dies out. Whenever RMS<1, each individual release more than one new infected individual and therefore the disease is able to invade the susceptible population. This allows us to determine the effectiveness of control measures. Epidemiologically, the reproductive number tells us how many secondary cases will one infected individual or vector produce in an entirely susceptible population of hosts and vectors, in other words, to determine the number of infected people that are generated by the introduction of a single infected person into a susceptible population. We need the computation of the basic reproduction number RMS to assess the stability of the disease-free equilibrium (DFE) and the endemic equilibrium point (EEP). The basic reproduction number for particular infections is dependent on the biological characteristics of the disease and on the behavioural pattern of the population. When there is a high transmission of the disease per unit time and durations of the infectious period, the basic reproduction number increases.Proposition 4.
The value of the basic reproduction number isRMS=β2β3bHbs/dHds(dH+ρH+qH)(ds+ρs).Proof.
The argument uses the approach of the next generation matrix (see [21]). Suppose(60)F=∂f1DMS∂Im∂f1DMS∂Ihs∂f1DMS∂Vms∂f1DMS∂Is∂f1DMS∂Iv∂f2DMS∂Im∂f2DMS∂Ihs∂f2DMS∂Vms∂f2DMS∂Is∂f2DMS∂Iv∂f3DMS∂Im∂f3DMS∂Ihs∂f3DMS∂Vms∂f3DMS∂Is∂f3DMS∂Iv∂f4DMS∂Im∂f4DMS∂Ihs∂f4DMS∂Vms∂f4DMS∂Is∂f4DMS∂Iv∂f5DMS∂Im∂f5DMS∂Ihs∂f5DMS∂Vms∂f5DMS∂Is∂f5DMS∂Ivwhich results into(61)F=0000β1ϵHσbHdH000β2bHdH0000000β3bsdsβ3η1bsds00β4ϵvσbvdv0β4ϵvσbvdv00,where F is the rate of appearance of new infections in one compartment and V is the transfer of individuals and mosquitoes into one compartment which is a Jacobian matrix evaluated at D0s and the Jacobian matrix of V is evaluated at DMS to give(62)V=ω+θ+dH00000ρH+qH+dH00000α+r+m+dH00000ds+ρs00000dvV-1=1ω+θ+dH000001ρH+qH+dH000001α+r+m+dH000001ds+ρs000001dv.
The product of matrixF and V-1 gives(63)FV-1=0000β1ϵhσbHdHdv000β2bHdHds+ρs0000000β3bsdsρH+qH+dHβ3bsη1α+r+m+dH00β4ϵvσbvdvω+θ+dH0β4ϵvσbvη2dvα+r+m+dH00.
Finding the eigenvalue ofFV-1, we obtain RMS for the malaria-schistosomiasis coinfection model (12) is given by(64)RM12=bH2β4ϵvϵhβ1σ2bvdHω+dH+θdv2RS12=β3bsβ2bHdsdHds+ρsdH+ρH+qHRMS=β3bsβ2bH2β4ϵvϵhβ1σ2bvdsdH2ds+ρsdH+ρH+qHω+dH+θdv2RMS2=β3bsβ2bH2β4ϵvϵhβ1σ2bvdsdH2ds+ρsdH+ρH+qHω+dH+θdv2.Considering the parameters for our model from Table 1, the value of RMS appears to be far greater than 1. This implies that the disease (malaria-schistosomiasis coinfection model) will not die out but remains endemic in the population until some control strategies are applied properly. The disease-free equilibrium is locally stable if RMS<1. We can achieve the global stability for the disease-free equilibrium by following the Lyapunov argument for RMS<1. Hence, when the disease-free equilibrium is unstable, there exists an endemic equilibrium. Next we discuss the local stability of the disease-free equilibrium.


### 4.1. Local Stability of the Disease-Free Equilibrium of Malaria-Schistosomiasis Coinfection Model

The disease-free equilibrium for the malaria-schistosomiasis coinfection model given byDMS=(SH∗,Im∗,Ihs∗,Vms∗,Rm∗,Rhs∗,Rms∗,Ss∗,Is∗,Sv∗,Iv∗) = (bH/dH,0,0,0,0,0,0,0,bs/ds,0,bv/dv,0) exists for all values of dH>0, ds>0, and dv>0. If the threshold value is RMS, then the disease-free equilibrium is locally asymptotically stable and the disease cannot invade or spread in the population or community. The summary of the result is given by Theorem 5.Theorem 5.
The disease-free equilibrium pointDMS is locally asymptotically stable if RMS<1 and unstable if RMS>1.It means that the malaria-schistosomiasis coinfection can be eliminated from the population wheneverRMS<1. RMS<1; then, averagely an infected person produces less than one newly infected person over the entire period of his infectiousness and malaria-schistosomiasis coinfection dies out. For RMS to be less than one, the β1, β2, β3, and β4 (rates of transmission) must decrease without bounds. Otherwise, if RMS>1, then each infected person produces less, averagely more than one new infection, and the malaria-schistosomiasis coinfection can spread or invade the population. Consider the model equation (12). At the equilibrium state DMS, the Jacobian is given by(65)JDMS=-dH000γkφ0-m10-m20-m700000000m200-m800000m100000-m900000000ω0ϕ1-α-γ+dH00000000qH00-m1000000000α00-φ+dH000000-m3-m4000-ds00000m3m40000-ds+ρs000-m50-m600000-dv00m50m6000000-dv,where m1=β2bH/dH,   m2=β1ϵhσbH/dH, m3=β3bs/ds, m4=β3η1bs/ds, m5=β4ϵvσbv/dv, m6=β4ϵvση2bv/dv, m7=ω+dH+θ, m8=dH+ρH+qH, m9=α+r+m+dH, and m10=g+dH+1-ϕ1-α.The local stability ofDMS is determined by the signs of the eigenvalues of the Jacobian matrix (65). The disease-free equilibrium point, DMS, is said to be locally asymptotically stable if the real parts of the eigenvalues of the Jacobian matrix (65) are all negative; otherwise it is said to be unstable. By considering the Jacobian matrix (65), we obtain(66)JDMS=-dH000γkφ0-m10-m20-m700000000m200-m800000m100000-m900000000ω0ϕ1-α-γ+dH00000000qH00-m1000000000α00-φ+dH000000-m3-m4000-ds00000m3m40000-ds+ρs000-m50-m600000-dv00m50m6000000-dv=0,where m1=β2bH/dH, m2=β1ϵhσbH/dH, m3=β3bs/ds, m4=β3η1bs/ds, m5=β4ϵvσbv/dv, m6=β4ϵvση2bv/dv, m7=(ω+dH+θ), m8=(dH+ρH+qH), m9=(α+r+m+dH), and m10=(g+dH)+(1-ϕ)(1-α), such that(67)λ1=-dH<0λ2=-α+dH+r+m<0λ3=-γ+dH<0λ4=-g-dH-1+α+ϕ-ϕα<0λ5=-φ+dH<0λ6=-ds<0λ7=-dv<0.The remaining characteristic polynomial corresponding toJDMS is(68)F1λ=λ2+dv+ω+dH+θλ-dvω+dH+θ-βvϵvσ2bvβ1ϵhbHdvdH+dvω+dH+θ=0,where(69)A1=1A2=dv+ω+dH+θA3=dvω+dH+θ-βvϵvσ2bvβ1ϵhbHdvdH+dvω+dH+θ.The Routh-Hurwitz conditions to establish that all roots of (F1) have negative real parts are A1>0, A2>0, A3>0, and A1A2>A3. Note that all model parameters are positive. For A3 to be positive set (70)dvω+dH+θ-βvϵvσ2bvβ1ϵhbHdvdH+dvω+dH+θ>0 which implies that (71)dvω+dH+θ1-βvϵvσ2bvβ1ϵhbHdv2dH+dvω+dH+θ>0. This leads to (1-RM12)>0⇒RM1<1.(72)dvω+dH+θ1-βvϵvσ2bvβ1ϵhbHdv2dH+dvω+dH+θ>0which is true. Thus, by the Routh-Hurwitz condition, all the eigenvalues have negative real parts whenever RM1<1, such that DMS is locally asymptotically stable. F2(λ)=λ2+(ds+ρs+dH+ρH+qH)λ-β3bsβ2bH/dsdH+(ds+ρs)(dH+ρH+qH), where(73)E1=1E2=ds+ρs+dH+ρH+qHE3=ds+ρsdH+ρH+qH-β3bsβ2bHdsdH. The Routh-Hurwitz conditions to make certain that all roots of (F2) have negative real parts are E1>0, E2>0, E3>0, and E1E2>E3. Note that all model parameters are positive. For E3 to be positive set (74)ds+ρsdH+ρH+qH-β3bsβ2bHdsdH>0 which implies that (75)1-β3bsβ2bHdsdHds+ρsdH+ρH+qH>0. This leads to (1-RS12)>0⇒RS1<1. It can be verified that E1E2-E3>0; that is, (76)ds+ρs+dH+ρH+qH-ds+ρsdH+ρH+qH1-RS12>0; hence it is true. Accordingly, by the Routh-Hurwitz criteria, all the eigenvalues have negative real parts whenever RS1<1, such that DMS is locally asymptotically stable.


### 4.2. Impact of Schistosomiasis on Malaria and Vice Versa

The impact of schistosomiasis on malaria and vice versa can be analysed by expressingR0m in terms of R0s. We solve dH to obtain the following.Suppose(77)P1=dsqH+ρHds+ρsP2=β3bsβ2bHdsds+ρs∴ (78) d H = - P 1 R 0 s + P 1 2 R 0 s 2 + 4 P 2 2 d s d s + ρ s R 0 s .Let(79)P12R0s2+4P2=P3R0s+P4.Therefore(80)dH=-P12R0s+P3R0s+P42dsds+ρsR0s.Hence(81)dH=P3-P1R0s+P42dsds+ρsR0s.We substitutedH in R0m to obtain(82)R0m2=β4ϵvbvβ1ϵhσ2ds2ds+ρs2R0sP3-R0sP1+P42ωds2+2ωdsρs+P3-P1+P4+2θds2+2θdsρsR0sdv2(83)∂R0m2∂Rs=-β4ϵvbvβ1ϵhσ2ds2ds+ρs22ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρsP3-P1+1P3R0s-P1R0s+P42ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρs2(84)∂R0m2∂Rs=-β4ϵvbvβ1ϵhσ2ds2ds+ρs22ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρsP1-P3+1P3R0s-P1R0s+P42ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρs2(85)∂R0m2∂Rs=-β4ϵvbvβ1ϵhσ2ds2ds+ρs2P3-P12ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρs-mP3-mP1+P4dv2P3R0s-P1R0s+P42ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρs2.If(ω+θ)ds2+ρs(ω+θ)ds+1/2P3-1/2P1(P1-P3) ≥ 1/1((ω+θ)ds2+ρs(ω+θ)ds-P1+P3)P4 is strictly positive, then the malaria enhances schistosomiasis infection. This means that whenever (85) is greater than zero, an increase in schistosomiasis will result in an increase of malaria cases in the community. If (85) is equal to zero, it implies that schistosomiasis cases will have no significant impact on the transmission dynamics of malaria. Accordingly, by expressing dH in terms of R0M (impact of malaria on schistosomiasis) (86)dH=-1/2dvω+1/2dvθR0M+1/2dv2R0M2ω2+2dv2Rv2ωθ+dv2R0M2θ2+4β4bvβ1ϵhbHσ22dvR0M, where (87)Q1=dvω+dvθQ2=β4bvβ1ϵhbHσ2dH=-Q1R0M+Q12R0M2+4Q24dvR0M.By assigning(88)Q12R0M2+4Q2=Q3R0M+Q4 therefore (89)dH=-Q1R0M+Q3R0M+Q44dvR0MR0S=4β3bsbHdv2R0M2dsR0MQ1-R0MQ3-Q4ds+ρsR0MQ1-Q4-R0MQ3-2ρHdvR0M-2qHdvR0M.By differentiating R0S partially with respect to R0M we obtain (90)∂R0S2R0M=8β3bsβ2bHdv2R0MQ4R0MρHdv+qHdv-Q1+Q3+Q4dsR0MQ1-R0MQ3-Q42ds+ρsR0MQ1-Q4-R0MQ3-2ρHdvR0M-2qHdvR0M2.If (90) is greater than zero, it means that an increase in the malaria accounts for an increase of schistosomiasis cases in the population. The impact of malaria treatment on schistosomiasis is evaluated by partially differentiating R0M with respect to ω. We obtain ∂R0M/∂ω=-1/θ+ω+dH. Since R0M is a decreasing function of ω, the treatment of schistosomiasis will have a positive impact on the dynamics of malaria.


## 4.1. Local Stability of the Disease-Free Equilibrium of Malaria-Schistosomiasis Coinfection Model

The disease-free equilibrium for the malaria-schistosomiasis coinfection model given byDMS=(SH∗,Im∗,Ihs∗,Vms∗,Rm∗,Rhs∗,Rms∗,Ss∗,Is∗,Sv∗,Iv∗) = (bH/dH,0,0,0,0,0,0,0,bs/ds,0,bv/dv,0) exists for all values of dH>0, ds>0, and dv>0. If the threshold value is RMS, then the disease-free equilibrium is locally asymptotically stable and the disease cannot invade or spread in the population or community. The summary of the result is given by Theorem 5.Theorem 5.
The disease-free equilibrium pointDMS is locally asymptotically stable if RMS<1 and unstable if RMS>1.It means that the malaria-schistosomiasis coinfection can be eliminated from the population wheneverRMS<1. RMS<1; then, averagely an infected person produces less than one newly infected person over the entire period of his infectiousness and malaria-schistosomiasis coinfection dies out. For RMS to be less than one, the β1, β2, β3, and β4 (rates of transmission) must decrease without bounds. Otherwise, if RMS>1, then each infected person produces less, averagely more than one new infection, and the malaria-schistosomiasis coinfection can spread or invade the population. Consider the model equation (12). At the equilibrium state DMS, the Jacobian is given by(65)JDMS=-dH000γkφ0-m10-m20-m700000000m200-m800000m100000-m900000000ω0ϕ1-α-γ+dH00000000qH00-m1000000000α00-φ+dH000000-m3-m4000-ds00000m3m40000-ds+ρs000-m50-m600000-dv00m50m6000000-dv,where m1=β2bH/dH,   m2=β1ϵhσbH/dH, m3=β3bs/ds, m4=β3η1bs/ds, m5=β4ϵvσbv/dv, m6=β4ϵvση2bv/dv, m7=ω+dH+θ, m8=dH+ρH+qH, m9=α+r+m+dH, and m10=g+dH+1-ϕ1-α.The local stability ofDMS is determined by the signs of the eigenvalues of the Jacobian matrix (65). The disease-free equilibrium point, DMS, is said to be locally asymptotically stable if the real parts of the eigenvalues of the Jacobian matrix (65) are all negative; otherwise it is said to be unstable. By considering the Jacobian matrix (65), we obtain(66)JDMS=-dH000γkφ0-m10-m20-m700000000m200-m800000m100000-m900000000ω0ϕ1-α-γ+dH00000000qH00-m1000000000α00-φ+dH000000-m3-m4000-ds00000m3m40000-ds+ρs000-m50-m600000-dv00m50m6000000-dv=0,where m1=β2bH/dH, m2=β1ϵhσbH/dH, m3=β3bs/ds, m4=β3η1bs/ds, m5=β4ϵvσbv/dv, m6=β4ϵvση2bv/dv, m7=(ω+dH+θ), m8=(dH+ρH+qH), m9=(α+r+m+dH), and m10=(g+dH)+(1-ϕ)(1-α), such that(67)λ1=-dH<0λ2=-α+dH+r+m<0λ3=-γ+dH<0λ4=-g-dH-1+α+ϕ-ϕα<0λ5=-φ+dH<0λ6=-ds<0λ7=-dv<0.The remaining characteristic polynomial corresponding toJDMS is(68)F1λ=λ2+dv+ω+dH+θλ-dvω+dH+θ-βvϵvσ2bvβ1ϵhbHdvdH+dvω+dH+θ=0,where(69)A1=1A2=dv+ω+dH+θA3=dvω+dH+θ-βvϵvσ2bvβ1ϵhbHdvdH+dvω+dH+θ.The Routh-Hurwitz conditions to establish that all roots of (F1) have negative real parts are A1>0, A2>0, A3>0, and A1A2>A3. Note that all model parameters are positive. For A3 to be positive set (70)dvω+dH+θ-βvϵvσ2bvβ1ϵhbHdvdH+dvω+dH+θ>0 which implies that (71)dvω+dH+θ1-βvϵvσ2bvβ1ϵhbHdv2dH+dvω+dH+θ>0. This leads to (1-RM12)>0⇒RM1<1.(72)dvω+dH+θ1-βvϵvσ2bvβ1ϵhbHdv2dH+dvω+dH+θ>0which is true. Thus, by the Routh-Hurwitz condition, all the eigenvalues have negative real parts whenever RM1<1, such that DMS is locally asymptotically stable. F2(λ)=λ2+(ds+ρs+dH+ρH+qH)λ-β3bsβ2bH/dsdH+(ds+ρs)(dH+ρH+qH), where(73)E1=1E2=ds+ρs+dH+ρH+qHE3=ds+ρsdH+ρH+qH-β3bsβ2bHdsdH. The Routh-Hurwitz conditions to make certain that all roots of (F2) have negative real parts are E1>0, E2>0, E3>0, and E1E2>E3. Note that all model parameters are positive. For E3 to be positive set (74)ds+ρsdH+ρH+qH-β3bsβ2bHdsdH>0 which implies that (75)1-β3bsβ2bHdsdHds+ρsdH+ρH+qH>0. This leads to (1-RS12)>0⇒RS1<1. It can be verified that E1E2-E3>0; that is, (76)ds+ρs+dH+ρH+qH-ds+ρsdH+ρH+qH1-RS12>0; hence it is true. Accordingly, by the Routh-Hurwitz criteria, all the eigenvalues have negative real parts whenever RS1<1, such that DMS is locally asymptotically stable.


## 4.2. Impact of Schistosomiasis on Malaria and Vice Versa

The impact of schistosomiasis on malaria and vice versa can be analysed by expressingR0m in terms of R0s. We solve dH to obtain the following.Suppose(77)P1=dsqH+ρHds+ρsP2=β3bsβ2bHdsds+ρs∴ (78) d H = - P 1 R 0 s + P 1 2 R 0 s 2 + 4 P 2 2 d s d s + ρ s R 0 s .Let(79)P12R0s2+4P2=P3R0s+P4.Therefore(80)dH=-P12R0s+P3R0s+P42dsds+ρsR0s.Hence(81)dH=P3-P1R0s+P42dsds+ρsR0s.We substitutedH in R0m to obtain(82)R0m2=β4ϵvbvβ1ϵhσ2ds2ds+ρs2R0sP3-R0sP1+P42ωds2+2ωdsρs+P3-P1+P4+2θds2+2θdsρsR0sdv2(83)∂R0m2∂Rs=-β4ϵvbvβ1ϵhσ2ds2ds+ρs22ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρsP3-P1+1P3R0s-P1R0s+P42ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρs2(84)∂R0m2∂Rs=-β4ϵvbvβ1ϵhσ2ds2ds+ρs22ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρsP1-P3+1P3R0s-P1R0s+P42ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρs2(85)∂R0m2∂Rs=-β4ϵvbvβ1ϵhσ2ds2ds+ρs2P3-P12ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρs-mP3-mP1+P4dv2P3R0s-P1R0s+P42ωds2R0s+2ωdsR0sρs+R0sP3-R0sP1+P4+2θds2R0s+2θdsR0sρs2.If(ω+θ)ds2+ρs(ω+θ)ds+1/2P3-1/2P1(P1-P3) ≥ 1/1((ω+θ)ds2+ρs(ω+θ)ds-P1+P3)P4 is strictly positive, then the malaria enhances schistosomiasis infection. This means that whenever (85) is greater than zero, an increase in schistosomiasis will result in an increase of malaria cases in the community. If (85) is equal to zero, it implies that schistosomiasis cases will have no significant impact on the transmission dynamics of malaria. Accordingly, by expressing dH in terms of R0M (impact of malaria on schistosomiasis) (86)dH=-1/2dvω+1/2dvθR0M+1/2dv2R0M2ω2+2dv2Rv2ωθ+dv2R0M2θ2+4β4bvβ1ϵhbHσ22dvR0M, where (87)Q1=dvω+dvθQ2=β4bvβ1ϵhbHσ2dH=-Q1R0M+Q12R0M2+4Q24dvR0M.By assigning(88)Q12R0M2+4Q2=Q3R0M+Q4 therefore (89)dH=-Q1R0M+Q3R0M+Q44dvR0MR0S=4β3bsbHdv2R0M2dsR0MQ1-R0MQ3-Q4ds+ρsR0MQ1-Q4-R0MQ3-2ρHdvR0M-2qHdvR0M.By differentiating R0S partially with respect to R0M we obtain (90)∂R0S2R0M=8β3bsβ2bHdv2R0MQ4R0MρHdv+qHdv-Q1+Q3+Q4dsR0MQ1-R0MQ3-Q42ds+ρsR0MQ1-Q4-R0MQ3-2ρHdvR0M-2qHdvR0M2.If (90) is greater than zero, it means that an increase in the malaria accounts for an increase of schistosomiasis cases in the population. The impact of malaria treatment on schistosomiasis is evaluated by partially differentiating R0M with respect to ω. We obtain ∂R0M/∂ω=-1/θ+ω+dH. Since R0M is a decreasing function of ω, the treatment of schistosomiasis will have a positive impact on the dynamics of malaria.


## 5. Sensitivity Analysis of the Malaria-Schistosomiasis Coinfection Model

Sensitivity analysis is the study of how the uncertainty in the output of a mathematical model or system can be distributed to different sources of uncertainty in its outputs. It is a scientific procedure used in determining how different values of an independent variable will affect a particular dependent variable under a given set of hypotheses. This method is applied within certain limits that will depend on more input variables. The sensitivity index is a measure of the relative change in a variable with respect to the relative change in its parameters. It is known that there are uncertainties which are associated with the estimation of certain parameter values. So, it will be useful to study the effect of these parameters on the basic reproduction numberRMS. The understanding from this study will help us to identify those parameters that cause the most increase/decrease in RMS and by this we will be able to recognize the type of control measure(s) that is most suitable and effective in controlling the disease. Hence, we compute the normalized forward sensitivity index of RMS with respect to the parameters in the model equation (12).The normalized forward sensitivity index of a variable,u, that depends differentiably on index of a parameter, z, is defined as Kzv=∂v/∂z×z/v. As we have an implicit formula for RMS we derive analytically for the sensitivity of RMS as KzRMS=∂RMS/∂z×z/RMS to each parameter involved in RMS.Important indices are given thus:(91)Kβ1RMS=∂RMS∂β1×β1RMS=+1.0Kβ2RMS=∂RMS∂β2×β2RMS=+1.0Kβ3RMS=∂RMS∂β3×β3RMS=+1.0Kβ4RMS=∂RMS∂β4×β4RMS=+0.99KbvRMS=∂RMS∂bv×bvRMS=+0.74KdHRMS=∂RMS∂dH×dHRMS=-2.0KϵvRMS=∂RMS∂ϵv×ϵvRMS=+1.0KϵhRMS=∂RMS∂ϵh×ϵhRMS=+1.0KωRMS=∂RMS∂ω×ωRMS=-0.91KθRMS=∂RMS∂θ×θRMS=-0.00000000000186KσRMS=∂RMS∂γ×σRMS=+2.0KdvRMS=∂RMS∂dv×dvRMS=-0.23KbHRMS=∂RMS∂bH×bHRMS=+2.0KbsRMS=∂RMS∂bs×bsRMS=+1.0KdsRMS=∂RMS∂ds×dsRMS=-0.11KρsRMS=∂RMS∂ρs×ρsRMS=-0.876KρHRMS=∂RMS∂ρH×ρHRMS=-0.0069KqHRMS=∂RMS∂qH×qHRMS=-0.99.


### 5.1. Interpretation of Sensitivity Indices

It is observed that the positive signs of sensitivity indices of the basic reproduction number to the model parameters reveal that an increase or reduction in the value of each of the parameter in this case will point to an increase or decrease in the basic reproduction number of the disease [22]. The case of the KϵvRMS=1.0 shows that increasing or reducing the mosquito biting rate by 10 percent increases or reduces the basic reproduction number, RMS, by 10 percent. However, the negative signs of the sensitivity indices of the basic reproduction number to the model parameters show that an increase or reduction in the value of each of the parameter in this class indicates a corresponding reduction or increase in the basic reproduction number of the disease [22]. Hence, the sensitivity analysis of the malaria-schistosomiasis coinfection model (12) provides an effective means of knowing how the dynamics of the coinfection model respond to changes in parameters of the model [22]. More importantly, it helps as a guide to policy makers and public health authorities in focusing on an appropriate intervention strategy for preventing and controlling the spread of the disease. The sensitivity analysis is applied to study the relative importance of model parameters to malaria-schistosomiasis coinfection. The sensitivity analysis reveals how each parameters influence the malaria-schistosomiasis coinfection output. However, the analysis will be carried out with response to the most crucial and sensitive parameters to the malaria-schistosomiasis coinfection model.


## 5.1. Interpretation of Sensitivity Indices

It is observed that the positive signs of sensitivity indices of the basic reproduction number to the model parameters reveal that an increase or reduction in the value of each of the parameter in this case will point to an increase or decrease in the basic reproduction number of the disease [22]. The case of the KϵvRMS=1.0 shows that increasing or reducing the mosquito biting rate by 10 percent increases or reduces the basic reproduction number, RMS, by 10 percent. However, the negative signs of the sensitivity indices of the basic reproduction number to the model parameters show that an increase or reduction in the value of each of the parameter in this class indicates a corresponding reduction or increase in the basic reproduction number of the disease [22]. Hence, the sensitivity analysis of the malaria-schistosomiasis coinfection model (12) provides an effective means of knowing how the dynamics of the coinfection model respond to changes in parameters of the model [22]. More importantly, it helps as a guide to policy makers and public health authorities in focusing on an appropriate intervention strategy for preventing and controlling the spread of the disease. The sensitivity analysis is applied to study the relative importance of model parameters to malaria-schistosomiasis coinfection. The sensitivity analysis reveals how each parameters influence the malaria-schistosomiasis coinfection output. However, the analysis will be carried out with response to the most crucial and sensitive parameters to the malaria-schistosomiasis coinfection model.


## 6. Numerical Simulation and Graphical Illustration of the Malaria-Schistosomiasis Coinfection Model

Numerical simulation and graphical illustrations are carried out in order to verify some of the analytical results. Different initial starts have been used from Table2 to perform the computer simulations for different cases and displayed graphically in Figures 1–4. Figure 1 shows the variation of susceptible human (Sh), infected human with malaria only (Im), infected human with schistosomiasis only (Ihs), and humans infected with both malaria and schistosomiasis (Vms). It is observed that, in Figure 1(a), the susceptible human population decreases with time due to the malaria-schistosomiasis coinfections and in Figure 1(b), the population of malaria-infected human decreases with time due to treatment of the infected population. In Figure 1(c), the population of the schistosomiasis-infected human population also decreases due to the treatment while the population of the coinfected decreases to a particular level and then reaches its equilibrium point due to treatments. Figure 2 shows the variation of humans who recovered from malaria only (Rm), humans who recovered from schistosomiasis only (Rhs), and humans who recovered from both malaria and schistosomiasis (Rms). It is observed that, in Figure 2(a), the population of malaria-recovered human population decreases with time to a level and in Figure 2(b), the population of the schistosomiasis-only-recovered human population increases with time and then decreases with time until it reaches its equilibrium point while, in Figure 2(c), the population of malaria-schistosomiasis recovered human population increases with time and then reaches its equilibrium point. Figure 3 shows the variation of susceptible mosquito (Sv) and infected mosquito with malaria (Iv). It is observed that, in Figure 3(a), the susceptible snail population decreases with time due to infection with schistosomiasis while, in Figure 3(b), the population of the infected snail maintains a constant equilibrium and then increases with time. Figure 4 shows the variation of susceptible mosquito (Sv) and infected mosquito with malaria (Iv). It is observed that, in Figure 4(a), the population of susceptible mosquito decreases and then increases to maintain a steady state while the infected mosquito population increases with time because of the presence of the infection.Table 2
Table showing numerical values of parameters used in the model.
Parameter Symbol Value Source Disease induced death rate due to malaria only θ 0.05 [23] Probability of human getting infected with malaria β 1 0.8333 [24, 25] Rate of acquiring schistosomiasis through contact with infected snails β 2 0.406 [26] Rate of acquiring schistosomiasis through contacts with infected humans β 3 0.0004 [26] Probability of mosquito getting infected by an infectious human β 4 0.09 [24, 25] Malaria schistosomiasis immunity waning rate φ 0.0005 [24] Schistosomiasis induced death rate r 0.02 [assumed] Malaria disease induced death rate m 0.06 [assumed] Coinfected population who recovered from malaria only ϕ 0.12 [assumed] Disease induced death rate of human due to schistosomiasis only ρ H 0.0039 [27] Disease induced death rate of snails due to schistosomiasis only ρ s 0.0004012 [28] Human spontaneous recovery ω 0.005 [13] Rate of loss of immunity to the schistosomiasis disease only k 0.7 [assumed] Rate of loss of immunity to malaria and schistosomiasis disease φ 0.0005 [24, 25] Human spontaneous recovery from schistosomiasis only q H 0.56 [assumed] Recovery rate of coinfected individual α 0.7 [assumed] Modification parameter η 1 1.3 [assumed] Modification parameter η 2 1.4 [assumed] Per capita birth rate of mosquitoes b v 1000 [24, 29] Per capita birth rate of snails b s 200 [8] Natural death rate of humans d H 0.00004 [10] Per capita biting rate of mosquitoes ϵ v 0.2 [30–32] Contact rate of vector per human per unit time σ 0.502 [33] Per capita biting rate of humans ϵ h 0.2 [30, 34] Natural death rate of mosquitoes d v 0.1429 [24] Natural death rate of snails d s 0.0000569 [26] Rate of loss of immunity to the malaria disease only γ 0.7902 [24, 29, 35] Per capita birth rate of humans b h 100 [36]Figure 1
Simulation showing the variation of susceptible human(Sh) in (a), infected human with malaria only (Im) in (b), infected human with schistosomiasis only (Ihs) in (c), and humans infected with both malaria and schistosomiasis (Vms) in (d).
(a) (b) (c) (d)Figure 2
Simulation showing the variation of humans who recovered from malaria only(Rm) in (a), humans who recovered from schistosomiasis only (Rhs) in (b), and humans who recovered from both malaria and schistosomiasis (Rms) in (c).
(a) (b) (c)Figure 3
Simulation showing the variation of Susceptible Snail(Ss) in (a), and Infected Snail with schistosomiasis (Is) in (b).
(a) (b)Figure 4
Simulation showing the variation of susceptible mosquito(Sv) in (a) and infected mosquito with malaria (Iv) in (b).
(a) (b)


## 7. Concluding Remarks

In this paper, we developed and analysed a deterministic model for the transmission of malaria-schistosomiasis coinfection without intervention strategies. The model was rigorously analysed to gain insights into its qualitative dynamics. We obtained the following results: (i) the malaria-schistosomiasis coinfection model has a locally stable disease-free equilibrium whenever the associated reproduction number is less than unity. (ii) We observed from our analysis of the effect of schistosomiasis on malaria that schistosomiasis infection may be associated with an increased risk of malaria. (iii) We show that the interaction between schistosomiasis and malaria may reduce the effectiveness of malaria treatment for controlling malaria transmission. Sensitivity analysis of the reproduction number to the model parameter was performed to investigate the parameters that possess greater influence on the model. Sensitivity analysis also showed that mosquito biting rateϵv and rate of acquiring schistosomiasis through contacts with infected snails β3 among other parameters contribute most significantly to the transmission and spread of the malaria disease. Hence the intervention strategy that inhibits the human-mosquito contact and human-snail contact should be encouraged in order to achieve a malaria-schistosomiasis-free population.



---

*Source: 3854902-2016-11-23.xml*