Next Article in Journal
New Results for Multivalued Mappings in Hausdorff Neutrosophic Metric Spaces
Previous Article in Journal
Latent Multi-View Semi-Nonnegative Matrix Factorization with Block Diagonal Constraint
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamical Behaviour of a Modified Tuberculosis Model with Impact of Public Health Education and Hospital Treatment

by
Fatima Sulayman
1,2,† and
Farah Aini Abdullah
1,*,†
1
School of Mathematical Sciences, Universiti Sains Malaysia (USM), Penang 11800, Malaysia
2
Department of Mathematical Sciences, Ibrahim Badamasi Babangida University, Lapai 911101, Nigeria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 19 September 2022 / Revised: 1 November 2022 / Accepted: 2 November 2022 / Published: 13 December 2022
(This article belongs to the Topic Mathematical Modeling)

Abstract

:
Tuberculosis (TB), caused by Mycobacterium tuberculosis is one of the treacherous infectious diseases of global concern. In this paper, we consider a deterministic model of TB infection with the public health education and hospital treatment impact. The effective reproductive number, R p h , that measures the potential spread of TB is presented by employing the next generation matrix approach. We investigate local and global stability of the TB-free equilibrium point, endemic equilibrium point, and sensitivity analysis. The analyses of the proposed model show that the model undergoes the phenomenon of backward bifurcation when the effective reproduction number ( R p h ) is less than one, where two stable equilibria, namely, the DFE and an EEP coexist. Further, we compute the sensitivity of the impact of each parameter on the effective reproductive number of the model by employing a normalized sensitivity index formula. Numerical simulation of the proposed model was conducted using Maple 2016 and MatLab R2020b software and compared with the theoretical results for illustration purposes. The investigation results can be useful in providing information to policy makers and public health authorities in mitigating the spread of TB infection by public health education and hospital treatment.

1. Introduction

Tuberculosis (TB) is one of the most hazardous infectious diseases that has become a significant widespread phenomenon, claiming more lives than any other contagious disease every day, according to [1]. Approximately 1/3 of the total population has a TB infection, resulting in millions of deaths and new cases annually (World Health Organization report). The report corroborates that TB is one of the top ten causes of mortality globally of both human and animal populations [2,3,4,5,6]. In 2020, for instance, 10 million individuals developed tuberculosis (TB), and more than 1.5 million died from it, including 214,000 HIV-positive persons (Human Immunodeficiency Virus) [7]. Typically, the signs may not be instantaneous when an individual contracts the disease. Thus, the individual remains asymptomatic for a long time or is latently infected for life [8]. Young adults may become infected by TB when they are most active [5]. Generally, TB-related deaths often happen in middle-income countries, for example, India, which leads the count, followed by Indonesia, China, the Philippines, Pakistan, Nigeria, Bangladesh, and South Africa. Such countries account for over 87 % of the entire TB trouble in the world. As a result, it is vital to implement techniques and methods that make it simple to understand how this disease spreads and predict its progression.
TB is a communicable disease caused by Mycobacterium tuberculosis affecting mostly the lungs [9,10]. However, it can also attack different organs including the brain, kidney, spine, central nervous system, or the lymphatic system [2,5,11,12,13]. It is important to note that the active lung TB disease typically begins with a cough, with sputum or blood on occasion, chest pains, fatigue, unexpected weight reduction, fever, and night sweats, which may last at least three or more weeks at a time. It has been reported by [14] that a pregnant woman who is infected may infect the foetus in some situations. Only individuals who have active TB can spread the disease. The latently infected individuals do not spread the bacteria [5,7]. Transmission starts with one individual, then onto the next, and relies upon the number of infected and expelled drops, the period of contaminated risk exposure, the virulence of the Mycobacterium tuberculosis, and the activity of environmental ventilation [7,10,15].
It was reported that behavioural change played a significant role in the transmission of infections [16,17,18]. Public health educational campaigns on TB disease plays a vital role in TB management and prevention [16,19], and can increase health literacy and awareness of TB among the population. The identification case rate and cure rate of TB patients can be improved by health education on TB disease, as it is an effective choice to minimize the spread of TB. As a result of the emergence of multi-resistance, treatment of TB is very challenging. Drug-resistant TB is a major public health concern across many developing countries, while treatment takes longer and needs more expensive drugs [5]. However, in many developing countries, the treatment of TB is not entirely free. Some TB patients cannot afford the full cost of treatment. As a result, some TB patients are treated at home to save money, as they can not afford to stay in the hospital. Some patients who were not cured choose to be discharged and continue their treatment at home [20].
Several mathematical models have been constructed to study the dynamical behaviour of TB infection. Abimbade et al. [21] formulated and analyzed an optimal control analysis of a TB transmission with incomplete treatment and exogenous re-infection, where they further divided the infected compartment into two groups, namely, uninformed and enlightened. Their findings revealed that both single controls and combinations of three controls have a positive influence on TB burden reduction. Recently, Ojo et al. [22] proposed a mathematical model for TB transmission with control. Their findings demonstrated that minimizing effective contact with infected people and increasing the rate of vaccinating susceptible persons with high vaccine efficacy will bring down the TB burden in a population. Xueyong et al. [16] formulated a TB model with healthy education and treatment. The qualitative analysis of their model exhibited the phenomenon of backward bifurcation, where the stable disease-free equilibrium co-exists with a stable endemic equilibrium, when the reproduction number is less than one. Huo and Zou [20] investigated a TB dynamics model with two types of treatment, namely, treatment at home and treatment in hospital. Their findings revealed that home treatment has a significant detrimental impact on tuberculosis spread. Okuonghae contributed significantly to TB epidemiology. They investigated the dynamics of TB and developed different models with comprehensive observations and results. Okuonghae and Egonmwan [23] presented a deterministic model that explores the impact of diagnosing and treating both latent and active TB infections on the disease’s infection dynamics in a population. Castillo and Song [6] studied a comprehensive review of the literature on TB dynamics. They compiled many dynamical models of TB and established a theoretical framework. Recently, Mustapha et al. [24] introduced a TB model with hospitalization and reinfection. A mathematical analysis of the equilibrium points, the basic reproduction number, and sensitivity analysis was discussed. However, the authors of [24] did not include public health education in their model. This factor continues to be one of the most important aspects in people’s lives. Public health services contributed to reducing the spread of TB infection.
Inspired by the above discussions, this paper seeks to fill a gap in the references cited, focusing on both public health education and hospital treatment. Infectious individuals in this model are divided into two classes: Infected individuals receiving treatment at home and infected individuals receiving treatment at the hospital. In this paper, we consider the impact of public health education and hospital treatment, comparing with the previous paper of [20]. The key differences are, in this paper, both the susceptible and infected compartment are divided into two classes, which are susceptible and educated susceptible, and infected individuals at home and in the hospital. The other parts of the paper are assembled as follows: Section 2 is devoted to the model formulation along with its fundamental properties. Section 3 deals with analysis of the existence of equilibria and how it relates to the effective reproduction number. Section 4, bifurcation analysis was conducted theoretically and supported with some numerical experiments to give a visualization of the results. Sensitivity analysis is presented in Section 5. Numerical results and discussions are given in Section 6. Finally, conclusions and future work are given in Section 7.

2. Materials and Methods

In this section, a TB model with the impact of public health education and hospital treatment was considered by stratifying total population into epidemiological classes represented as N ι ( t ) , into six compartments, namely, susceptible individuals ( S ι ) , Educated Susceptible ( P ι ) , Exposed individuals ( E ι ) , Infectious individuals at home ( I ι 1 ) , Infectious individuals receiving treatment at hospital ( I ι 2 ) and Recovered individuals ( R ι ) , so that:
N ι = S ι + P ι + E ι + I ι 1 + I ι 2 + R ι .

2.1. Clinical Assumptions of the TB Model with Public Health Education and Hospital Treatment

The assumptions used in this model are as follows:
  • There is a constant recruitment rate to the susceptible population, and natural cause of death affects individuals in all compartments, with an extra TB-induced death rate in the infected class.
  • We assume that, at any moment, the educated susceptible group may act as ignorantly and enter the class of susceptible at a constant rate ϕ , i.e., loss of temporary protection [16].
  • We further assumed that some uncured TB patients wanted to be discharged and continue treatment at home at a rate η 2 , [20].
  • Both infected individuals at home and in the hospital experience the infections effect at the rate λ .
  • The infected compartment is divided into two groups, namely, infected individuals who received treatment at home I ι 1 , and infected individuals receiving treatment in hospital I ι 2 .
  • The recovered individual may be again infected by an infectious individual [25].
In this model, it is assumed that the number of susceptible population is generated via recruitment of individuals into the population at a rate Λ . Susceptible individuals are educated at the rate ψ , and then transferred into the educated class. We further assume that education programs provide “temporary protection” at the per-capita rate ϕ [19,26]. We assume that both infected individuals at home and in the hospital experience the infections effect at the rate λ , where
λ = β ( I ι 1 + z I ι 2 ) N ι .
In (2), the parameter β , is the transmission rate, while 0 z 1 is the modification parameter associated with reduced infectiousness of individuals in the hospital ( in the I ι 1 compartment ) compared to individuals at home. Individuals who have received public health education have a high level of awareness. As a result, educated people are a low-risk population. Education reduces the risk of infection by a factor of ν ( 0 ν 1 ) . The scenario ν = 0 indicates that the education is totally effective to prevent the infection, whereas ν = 1 models the case where the education program is totally ineffective [19,26]. The populations of susceptible and educated susceptible individuals are further decreased as a result of natural death rate μ . Therefore, the rates of change of the populations of susceptible and educated susceptible individuals are given, respectively, by
d S ι d t = Λ + ϕ P ι λ S ι ψ S ι μ S ι ,
d P ι d t = ψ S ι ν λ P ι ϕ P ι μ P ι .
The population of exposed individuals is generated via the infection of susceptible ( S ι ) , educated susceptible ( P ι ) and recovered ( R ι ) , at the rate λ , ν λ and σ λ , and decreased as a result of disease symptoms at a rate κ 1 , κ 2 , and natural death at the rate μ , and σ ( 0 σ 1 ) is taken as the reinfection rate so that
d E ι d t = λ S ι + ν λ P ι + σ λ R ι ( κ 1 + κ 2 + μ ) E ι .
Similarly, the population of infected individuals at home is generated as result of disease symptoms at a rate κ 1 , the rate of TB patients who are not cured from the hospital at a rate η 2 and diminished due to progression to hospital at a rate η 1 , recovery by TB patients at home at a rate γ 1 , natural death cause μ , and TB-induced death at the rate δ 1 , respectively. Thus,
d I ι 1 d t = κ 1 E ι + η 2 I ι 2 ( η 1 + γ 1 + μ + δ 1 ) I ι 1 .
In addition, the population of infected individuals receiving treatment at the hospital is increased as a result of the disease symptoms at a rate κ 2 , and hospital admittance of infected individuals receiving treatment at home at a rate η 1 , and diminished due to their returning back home to I ι 1 at a rate η 2 , natural death cause μ , recovery by TB patients at hospital at a rate γ 2 , and TB induced death at the rate δ 2 . Therefore,
d I ι 2 d t = κ 1 E ι + η 1 I ι 1 η 1 + μ + γ 2 + δ 2 I ι 2 .
Finally, the population of recovered individuals is generated via recovery by both infected individuals treated at home γ 1 and at hospital γ 2 at a rate γ 1 , γ 2 , respectively, and diminished as a result of natural death at a rate μ and progression to the exposed class at a rate σ λ .
d R ι d t = γ 1 I ι 1 + γ 2 I ι 2 σ λ R ι μ R ι .
Therefore, based on the above description and assumptions, the model of TB with public health education and hospital treatment lead to following a system of non-linear differential equations; the schematic diagram Figure 1 below, and the parameters indicated in the diagram are explained in Table 1.
d S ι d t = Λ + ϕ P ι ( λ + ψ + μ ) S ι , d P ι d t = ψ S ι ( ν λ + ϕ + μ ) P ι , d E ι d t = λ S ι + ν λ P ι + σ λ R ι ( κ 1 + κ 2 + μ ) E ι , d I ι 1 d t = κ 1 E ι + η 2 I ι 2 ( η 1 + γ 1 + μ + δ 1 ) I ι 1 , d I ι 2 d t = κ 2 E ι + η 1 I ι 1 η 2 + γ 2 + μ + δ 2 I ι 2 , d R ι d t = γ 1 I ι 1 + γ 2 I ι 2 ( σ λ + μ ) R ι ,
where
λ = β ( I ι 1 + z I ι 2 ) N ι ,
with the initial conditions given by Equation (5)
S ι ( 0 ) = S ι 0 > 0 , P ι ( 0 ) = P ι 0 > 0 , E ι ( 0 ) = E ι 0 > 0 , I ι 1 ( 0 ) = I ι 10 > 0 , I ι 2 ( 0 ) = I ι 20 > 0
satisfies the equation
d N ι d t = Λ μ N ι δ 1 I ι 1 δ 2 I ι 2 .

2.2. Basic Properties of the TB Model with Public Health Education and Hospital Treatment

The basic properties of the TB model (23) are explained in this subsection. In an invariant region, the model is proved to be positive and bounded. When studying the dynamical behaviour of an epidemiologically model, this analysis is crucial because it demonstrates whether the model is epidemiologically relevant and mathematically well-posed, that is, whether the model and its predictions are certain [27].

Positivity of Solutions

To established that the TB infection model system (23) is epidemiologically realistic, all the stated variables are positive all the time.
Theorem 1. 
Let initial data be { ( S ι , P ι , E ι , I ι 1 , I ι 2 , R ι ) 0 } Φ . Thus, the solution of { S ι ( t ) , P ι ( t ) , E ι ( t ) , I ι 1 ( t ) , I ι 2 ( t ) , R ι ( t ) } of the model system (23) is non-negative for all t > 0 .
Proof of Theorem 1 
The method described in [28,29], is applied. We use the first equation to consider the non-linear system of (23), which clearly shows that
d S ι d t + ( λ ( t ) + ψ + μ ) S ι > 0 ,
utilizing an integrating factor gives
d d t S ι e x p 0 t ( λ ( ϵ ) + ψ + μ d ϵ > 0 .
Using the initial conditions (5) and integrating (7) results in
S ι ( t ) > S ι 0 e x p 0 t ( λ ( ϵ ) + ψ + μ d ϵ > 0 , t > 0 .
Similarly, it can be shown that P ι ( t ) > 0 , E ι ( t ) > 0 , I ι 1 ( t ) > 0 , I ι 2 ( t ) > 0 , R ι ( t ) > 0 , t > 0. □

2.3. Invariant Region

Theorem 2. 
The solution of the TB model system (23) is enclosed in the region Φ subset R + 6 , given by
Φ = S ι , P ι , E ι , I ι 1 , I ι 2 , R ι R + 6 , N ι Λ μ ,
for the initial conditions (5) in Φ .
Proof of Theorem 2 
The change in the total population is given by
d N ι d t = Λ μ N ι I ι 1 δ 1 I ι 2 δ 2 ,
In the absence of TB infection, there is no death from TB transmission, (that is, δ 1 = δ 2 = 0) [30], hence the rate of change of the total population size in Equation (8) is given as
d N ι d t Λ μ N ι ,
The solution of Equation (9) is solved by using the same approach presented in [31,32], and is given by
N ι ( t ) Λ μ Λ μ N ι 0 e μ t ,
where N ι 0 = N ι ( 0 ) .
Using [33] we note that if N ι 0 < Λ μ , N ι Λ μ asymptotically as t in Equation (10) the total population size N ι Λ μ , which means that 0 N ι Λ μ . Therefore, all the feasible solutions in the model converge in the region Φ [34]. □

3. Existence of Equilibrium

3.1. Tuberculosis Free Equilibrium (TFE)

The T F E state denoted as T 0 , is generally described as a state in which no disease exists in a given population. The infected population can be defined as a disease type. Using the first fourth equation of system model (23) with E ι = I ι 1 = I ι 2 = R ι = 0 into consideration, we arrive at:
T 0 = ( S ι 0 , P ι 0 , E ι 0 , I ι 1 0 , I ι 2 0 , R ι 0 ) = Λ μ μ + ϕ μ + ψ + ϕ , Λ μ ψ μ + ψ + ϕ , 0 , 0 , 0 , 0 .

3.2. Calculation of Effective Reproduction Number ( R p h ) For the System Model (23)

This section studies the effective reproduction number, which is a threshold parameter that controls the spread of a disease. To obtain the effective reproduction number, we apply the next generation approach described in [35]. The associated incidence matrix ( F ) and the transition matrix ( V ) of system model (23), are respectively obtained as:
F = 0 β ( S ι 0 + ν P ι 0 ) N ι 0 β z ( S ι 0 + ν P ι 0 ) N ι 0 0 0 0 0 0 0 , V = A 1 0 0 κ 1 A 2 η 2 κ 2 η 1 A 3 ,
where
A 1 = κ 1 + κ 2 + μ , A 2 = η 1 + μ + γ 1 + δ 1 , A 3 = η 2 + μ + γ 2 + δ 2 .
Thus, the effective reproductive number of system model (23) is calculated from the spectral radius ρ ( F V 1 ) as:
R p h = β ( ϕ + μ + ν ψ ) ( κ 1 A 3 + κ 2 η 2 ) + z ( κ 2 A 2 + η 1 κ 1 ) A 1 ( A 2 A 3 η 1 η 2 ) ( ψ + ϕ + μ ) .
After substituting the values of A 1 , A 2 and A 3 , Equation (14) becomes
R p h = β ν ψ + μ + ϕ η 2 + μ + γ 2 + δ 2 κ 1 + κ 2 η 2 + z η 1 + μ + γ 1 + δ 1 κ 2 + κ 1 η 1 C 0 + C 1 + μ η 1 + μ + γ 1 + δ 1 η 2 + μ + γ 2 + δ 2 η 1 η 2 ψ + μ + ϕ .
where
C 0 = κ 1 η 1 + μ + γ 1 + δ 1 η 2 + μ + γ 2 + δ 2 η 1 η 2 ψ + μ + ϕ ,
C 1 = κ 2 η 1 + μ + γ 1 + δ 1 η 2 + μ + γ 2 + δ 2 η 1 η 2 ψ + μ + ϕ .
The term effective reproduction number, represented by R p h , is defined as the expected number of secondary cases generated by a single infected individual during the period of infectiousness in a population of susceptible individuals where public health education and hospital treatment are incorporated [36].

3.3. Local Stability of TB Free Equilibrium

To prove the local stability of the TFE, the Jacobian of the proposed model system (23) is used. After that, the Jacobian is used to derive the characteristic equation, from which the eigenvalue result is obtained.
Theorem 3. 
The TFE of the proposed model system (23) is locally asymptotically stable when R p h < 1 .
Proof of Theorem 3 
To prove the system’s local stability, the Jacobian of the proposed model system (23) is investigated at TFE, which is then given by:
J ( T 0 ) = ( μ + ψ ) ϕ 0 β S ι 0 N ι 0 β z S ι 0 N ι 0 0 ψ ( ϕ + μ ) 0 ν β P ι 0 N ι 0 ν β z P ι 0 N ι 0 0 0 0 A 1 β ( S 0 + ν β P 0 ) N 0 β z ( S 0 + ν P ι 0 ) N ι 0 0 0 0 κ 1 A 2 η 2 0 0 0 κ 2 η 1 A 3 0 0 0 0 γ 1 γ 2 μ .
The characteristic equation of (16) is given by:
( λ + μ ) ( λ + ψ + μ ) ( λ + ϕ + μ ) ( λ + A 1 ) ( λ + A 2 ) ( λ + A 3 ) = 0 .
It can be seen from (17) that:
λ 1 = μ < 0 , λ 2 = ( ψ + μ ) < 0 , λ 3 = ( ϕ + μ ) < 0 , λ 4 = A 1 < 0 , λ 5 = A 2 < 0 , λ 6 = A 3 < 0 .
The result reveals that all the eigenvalues λ 1 , λ 2 , λ 3 , λ 4 , λ 5 , λ 6 are negative given that all the parameters values are greater than zero. As a result, according to the principle of linearized stability [37], the TFE is asymptotically stable. □

3.4. Global Stability of the TB Free Equilibrium

This section investigates the global stability of the T F E point T 0 of the model system (23). To investigate the global stability of the system model (23), we employ the techniques implemented by [38].
d X d t = F ( X ι , Z ι ) , d Z d t = G ( X ι , Z ι ) ; G ( X ι , 0 ) = 0 .
Here, the components X ι = ( S ι , P ι , R ι ) and Z ι = ( E ι , I ι 1 , I ι 2 ) , where X ι R 3 denotes the uninfected population and Z ι R 3 denotes the infected population.
The TB free equilibrium is defined by T 0 = ( X ι * , 0 ) .
The fixed point T 0 = ( X ι * , 0 ) is a globally asymptotically stable equilibrium for the system model (23) provided that R p h , which is locally asymptotically stable, and the following two conditions must be satisfied:
( H 1 ) : for d X ι d t = F ( X ι , 0 ) , X *
is globally asymptotically stable (GAS),
( H 2 ) : G ( X ι , Z ι ) = A Z ι G ^ ( X ι , Z ι ) , G ^ ( X ι , Z ι ) 0 for ( X ι , Z ι ) Φ ,
If the system model (23) meets the given two criteria, then the following theorem holds.
Theorem 4. 
The TB free equilibrium point T 0 = ( X * , 0 ) of the system model (23) is globally asymptotically stable provided R p h < 1 and the conditions ( H 1 ) and ( H 2 ) are satisfied.
Proof of Theorem 4 
From system model (23) we can get F ( X ι , Z ι ) and G ( X ι , Z ι ) :
F ( X ι , Z ι ) = Λ + ϕ P ι ( λ + ψ + μ ) S ι ψ S ι ( ν λ + ϕ + μ ) P ι γ 1 I ι 1 + γ 2 I ι 2 ( σ λ + μ ) R ι ,
G ( X ι , Z ι ) = λ S ι + ν λ P ι + σ λ R ι ( κ 1 + κ 2 + μ ) E ι κ 1 E ι + η 2 I ι 2 ( η 1 + γ 1 + μ + δ 1 ) I ι 1 κ 2 E ι + η 1 I ι 1 ( η 2 + γ 2 + μ + δ 2 ) I ι 2 .
At T 0 , Equation (20)
d X d t = F ( X ι , 0 ) = Λ + ϕ P ι ( ψ + μ ) S ι ψ S ι ( ϕ + μ ) P ι μ R ι .
From the system (22) above, we see that X ι * = Λ ( μ + ϕ ) μ ( μ + ψ + ϕ ) , Λ ψ μ ( μ + ψ + ϕ ) , 0 is globally asymptotical stable. This can be verified from the solutions, namely,
S ι ( t ) = Λ ( μ + ϕ ) μ ( μ + ψ + ϕ ) + S ι ( 0 ) Λ ( μ + ϕ ) μ ( μ + ψ + ϕ ) e x p ( ψ + μ + ϕ ) t . P ι ( t ) = Λ ψ μ ( μ + ψ + ϕ ) + P ι ( 0 ) Λ ψ μ ( μ + ψ + ϕ ) e x p ( ψ + μ + ϕ ) t . R ι ( t ) = R ι ( 0 ) e x p μ t .
As t , the solution S ι ( ) , P ι ( ) , R ι ( )   Λ ( μ + ϕ ) μ ( μ + ψ + ϕ ) , Λ ψ μ ( μ + ψ + ϕ ) , 0 , which implies the global convergence of (22) in Φ , and this satisfies condition H 1 .
Next, applying the second condition of the theorem H 2 .
From H 2 , we have G ( X ι , Z ι ) = A Z ι G ^ ( X ι , Z ι ) , G ^ ( X ι , Z ι ) 0 for ( X ι , Z ι ) Φ .
Therefore, G ^ ( X ι , Z ι ) = A Z ι G ( X ι , Z ι ) . where A is an n × n matrix, Z ι is a column vector and G ( X ι , Z ι ) is a column vector formed from the infectious compartments. We already know that
G ( X ι , Z ι ) = G 1 ( X ι , Z ι ) G 2 ( X ι , Z ι ) G 3 ( X ι , Z ι ) = λ S ι + ν λ P ι + σ λ R ι ( κ 1 + κ 2 + μ ) E ι κ 1 E ι + η 2 I ι 2 ( η 1 + γ 1 + μ + δ 1 ) I ι 1 κ 2 E ι + η 1 I ι 1 η 2 + γ 2 + μ + δ 2 I ι 2 .
Now let compute A
A = ( κ 1 + κ 2 + μ ) β ( ϕ + μ + ν ψ ) ϕ + μ + ψ β z ( ϕ + μ + ν ψ ) ϕ + μ + ψ κ 1 ( η 1 + γ 1 + μ + δ 1 ) η 2 κ 2 η 1 ( η 2 + γ 2 + μ + δ 2 ) .
The matrix A is a metzler matrix because all its off-diagonal elements are non-negative.
A = ( κ 1 + κ 2 + μ ) β ( ϕ + μ + ν ψ ) ϕ + μ + ψ β z ( ϕ + μ + ν ψ ) ϕ + μ + ψ κ 1 ( η 1 + γ 1 + μ + δ 1 ) η 2 κ 2 η 1 ( η 2 + γ 2 + μ + δ 2 ) E ι I ι 1 I ι 2 .
A Z ι = β ( ϕ + μ + ν ψ ) ( ϕ + μ + ψ ) ( I ι 1 + z I ι 2 ) ( κ 1 + κ 2 + μ ) E ι κ 1 E ι ( η 1 + γ 1 + μ + δ 1 ) I ι 1 + η 2 I ι 2 κ 2 E ι η 1 I ι 1 ( η 2 + γ 2 + μ + δ 2 ) I ι 2 .
Thus, using G ^ ( X ι , Z ι ) = A Z ι G ( X ι , Z ι ) , one obtains the following
G ^ ( X ι , Z ι ) = β ( I ι 1 + z I ι 2 ) ( ϕ + μ + ν ψ ) ( ϕ + μ + ψ ) S ι + ν P ι + σ R ι N ι 0 0 .
Since G ^ ( X ι , Z ι ) 0 , H 2 is not satisfied. This suggests that backward bifurcation may occur at T 0 when R p h < 1 .

3.5. Endemic Equilibrium State ( T * )

The stability of ( T 0 ) drives the disease’s short-term outbreaks. The stability at the endemic equilibrium points ( T * ) characterises its dynamics over a longer duration. We discovered that long-term behaviour has crucial epidemiological implications, such as whether an outbreak of a disease will lead to an endemic scenario or whether the infection will die out. We shall conduct an endemic analysis in this section. T * can now be determined by equating all of the model system in (23) to zero:
Theorem 5. 
The endemic equilibrium state of the model system Equation (23) exists if the effective reproduction number R p h > 1 .
The proof of Theorem 5 is presented in Appendix A.
Theorem 6. 
The TB model with public health education and hospital treatment (23) has:
1. 
One or more endemic equilibria when R p h < 1 .
2. 
A unique endemic equilibrium when R p h > 1 .
3. 
No endemic equilibrium otherwise.

4. Local Stability of Endemic Equilibrium

The local stability of endemic equilibrium ( T * ) of the model system can be investigated by employing the concept of center manifold presented in [6].
Theorem 7. 
The TB model with public health education and hospital treatment (23) undergoes backward bifurcation with bistability of T * at R p h < 1 and R p h > 1 .
The proof of Theorem 7 is presented in Appendix B.
Figure 2 demonstrates a backward bifurcation diagram of force of infection λ against the effective reproduction number R p h of the system model (23). We can see from the figure that as R p h increases to one, the disease also increases and this occurs when R p h < 1 . At R p h = 1 , DFE and EEP coexists; this means that the disease cannot be eradicated from the population as a result of a high-level endemic. We can also discover that when R p h > 1 , the disease persists.

5. Sensitivity Analysis

It is crucial to investigate how sensitive the TB model (23) is to changes in each of its parameters in order to identify control strategies that will assist in the lowering of the infection trajectory. In other words, conducting sensitivity analysis will assist in identifying what should be done or ignored in order to stop the spread of the TB transmission [29,39,40]. When a parameter changes, we can use sensitivity indices to calculate the relative change in a state variable. We employ the normalised forward sensitivity index of a variable to a parameter method provided in [41] in order to perform sensitivity analysis. The ratio of relative change in the variable to relative change in the parameter is known as the sensitivity index. The sensitivity index may also be defined using partial derivatives when the variable is a differentiable function of the parameter. Therefore, as stated in [41], the normalised forward sensitivity index ( S I ) of a threshold, R p h , that is differentiable with respect to a parameter, η , is defined as:
Υ η R p h = R p h η × η R p h .
The sensitivity of R p h to each of the 14 parameters presented in Table 2 is calculated using the effective reproduction number given below
R p h = β ( ϕ + μ + ν ψ ) ( κ 1 A 3 + κ 2 η 2 ) + z ( κ 2 A 2 + η 1 κ 1 ) A 1 ( A 2 A 3 η 1 η 2 ) ( ψ + ϕ + μ ) ,
where A 1 , A 2 , and A 3 are given in Equation (14).
Sensitivity index implied the parameters values obtained in Table 2.
From Table 2, we can see that the parameters have both positive and negative effects on R p h . Positive S I values, such as β , z , μ , ν , ϕ , γ 2 , η 2 and κ 1 , as shown in Table 2 reveal that an increase in these parameters values increases R p h , which brings about the infection attacking the population. While the parameters κ 2 , μ , ψ , γ 1 , η 1 , δ 1 , and δ 2 , have negative S I , that is, an increase in these parameters values decrease R p h , and as a result, the infection gradually fades from the population.
Figure 3 demonstrates the relationship between the effective reproduction number R p h , ψ , η 1 , κ 1 , and κ 2 .  Figure 3a, illustrates a linear relationship between R p h and the awareness rate ψ . We observed from Figure 3a, that R p h decreases with increasing the awareness rate ψ of susceptible individuals. This means that susceptible individuals who are aware of the TB disease and understand how it is transmitted have a great positive impact on the spread of TB infection.
Figure 3b demonstrates a linear relationship between R p h and the rate at which infected individuals progress from I ι 1 class to I ι 2 class. This indicates that R p h decreases with the increasing the rate at which infected individuals I ι 1 at home undergo treatment at the hospital. This means that infected individuals treated at the hospital have a great positive impact on the spread of TB infection too.
Figure 3c,d depicts a relationship between R p h , κ 1 and κ 2 . It is clearly seen from the Figure 3c,d that R p h increases with increasing κ 1 , while R p h decreases with increasing κ 2 . Biologically, this shows that infected individuals receiving treatment at home have a great negative impact on the dynamics of TB disease.

6. Numerical Results and Discussions

The numerical simulation of the model system (23), is performed with the values of the parameters described in Table 3. We used ode45 solver in Matlab, which depends on the Runge-Kutta technique to stimulate model system (23) with the parameter values presented in Table 3, as well as the following initial conditions:
S ι ( 0 ) = 0.6 , P ι ( 0 ) = 0.2 , E ι ( 0 ) = 0.1 , I ι 1 ( 0 ) = 0.1 , I ι 2 ( 0 ) = 0.1 , R ι ( 0 ) = 0 .
Figure 4 is a graphical representation of system model (23) showing the differences between susceptible and educated susceptible individuals. We can see from Figure 4 that, as time increases, the number of educated susceptible individuals’ health continues to increase until it reaches the carrying capacity, whereas, the number of susceptible individuals continues to decrease with time as more and more of them become infected.
Figure 5 depicts the impact of the various values of awareness rate ( ψ ) , on the dynamics of S ι , E ι , I ι 1 and I ι 2 . In general, this figure reflects that when ψ , increases, the number of E ι , I ι 1 a n d I ι 2 decreases rapidly. This demonstrates that public health education of susceptible individuals has the effect of limiting TB disease spread. As a result, public policymakers must focus on increasing the value of awareness rate ψ in order to prevent and control TB infection in a population.
Figure 6 investigates the impact of varying the rate at which the educated susceptible lose the awareness ϕ . This figure demonstrates a worst-case scenario in which an increase in loss of awareness in the susceptible individual leads to an increase in the proportion of infectious individuals.
Figure 7 investigates the effect of varying the reduction of infection by infected individuals receiving treatment at the hospital z . In this figure, z , is varied between 0 to 1. Of course, lowering the value of z, from 1 to 0, as expected, reduces the proportion of infected individuals, although at different rates.
Figure 8 illustrates the impact of varying the rate at which infected individuals at home go out for treatment in the hospital, η 1 between 0 and 1 on the individuals treated at home I ι 1 and infected individuals receiving treatment at the hospital, respectively. Obviously, as expected, an increase in the value of η 1 , from 0 to 1 greatly reduced the number of infected individuals both at home and at hospital. Epidemiologically, this shows that rate at which infected individuals at home go out for treatment in the hospital have a great impact on reducing the transmission of TB disease in the population.

7. Conclusions

This paper presented a new deterministic model of TB infection subject to the use of public health education and hospital treatment. To gain insight into its dynamic features, the model was rigorously analyzed. The analyses of the model, which contains six mutually-exclusive epidemiological partitions, show that the model undergoes the phenomenon of backward bifurcation when the effective reproduction number ( R p h ) is less than one, where two stable equilibria, namely, the DFE and an EEP coexist when the corresponding effective reproductive number is less than one. This backward bifurcation phenomenon of this article is very vital, and this occurs only under education of susceptible individuals and treatment of TB-infected individuals in the hospital. This is mostly telling us that, even if the effective reproduction number is less than unity, while necessary, it is not sufficient for efficiently controlling the spread of a TB epidemic, which is against classical epidemiological theory. In this scenario, TB elimination will depend upon the initial sizes of the sub-populations of the model. The parameters are sensitive to the transmission dynamics of TB diseases, either negatively ( μ , ν , ψ , η 1 , γ 1 , γ 2 , δ 1 , δ 2 , and κ 2 ) or positively ( β 1 , β 2 and κ 1 ) . According to the numerical investigation, increasing the public health education ψ on the S ι class, and the progression rate η 1 from I ι 1 to I ι 2 , have a significant great effect on reducing the prevalence of TB burden see Figure 5 and Figure 8. Given that public health education and hospital treatment can reduce the spread of TB infection, this programme should be continued and improved.
In light of the model study in this paper, some gaps in this paper need to be filled; these proposed gaps will allow for possible extension and improvement of the paper. The proposed model can be extended and improved in the future by:
  • Considering a stochastic model approach. This will result in more realistic TB model dynamics.
  • Since the spread of tuberculosis affects all age groups, it is crucial to consider the dynamics of the TB model by incorporating an age-structured model.
  • Real data will also be considered because collecting data for TB patients is difficult in epidemiological models; as a result, we use data collected or estimated from literature sources. Once we have real-world data for TB patients, we can compare it to theoretical outcomes.
  • Analyzing the dynamics of the TB model using a fractional order differential equation (FODE). It will be extremely interesting to use a FODE to examine the dynamics of TB model.

Author Contributions

Conceptualization, F.S. and F.A.A.; methodology, F.S. and F.A.A.; software, F.S.; validation, F.S. and F.A.A.; formal analysis, F.S.; data curation, F.S. and F.A.A.; writing—original draft preparation, F.S.; writing—review and editing, F.S.; visualization, F.A.A.; supervision, F.A.A.; project administration, F.A.A.; funding acquisition, F.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Ministry of Higher Education (MoHE), Malaysia; Fundamental Research Grant Scheme (FRGS)—Grant code (FRGS/1/2021/STG06/USM/02/9) and Article Processing Charge (APC) Fund from Research Creativity and Management Office (RCMO), Universiti Sains Malaysia.

Data Availability Statement

All the data and models used and/or generated during the study contains in the submitted article.

Acknowledgments

We would like to acknowledge the funder from the Ministry of Higher Education (MoHE), Malaysia; Fundamental Research Grant Scheme (FRGS)—Grant code (FRGS/1/2021/STG06/ USM/02/9) and Article Processing Charge (APC) Fund from Research Creativity and Management Office (RCMO), Universiti Sains Malaysia. We would also like to thank the School of Mathematical Sciences, Universiti Sains Malaysia for research assistance and computing facilities.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Theorem 5 
At the endemic state, the model system Equation (23) has an equilibrium point called TB endemic equilibrium point.
d S ι * d t = d P ι * d t = d E ι * d t = d I ι 1 * d t = d I ι 2 * d t = d R ι * d t = 0 .
By setting the derivatives of the left hand side of Equation (23) to zero and solving simultaneously we have
S ι * = Λ ( ν λ * + ϕ + μ ) ν λ * 2 + A 4 λ * + A 5 , P ι * = ψ λ * ν λ * 2 + A 4 λ * + A 5 , E ι * = Λ ( ν λ * + ϕ + μ + ν ψ ) ( A 11 λ * A 10 λ * ) ( ν λ * 2 + A 4 λ * + A 5 ) ( A 9 λ * + A 10 ) , I ι 1 * = κ 1 Λ A 1 A 3 ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * ( ν λ * 2 + A 4 λ * + A 5 ) ( A 9 λ * + A 10 ) , I ι 2 * = κ 1 Λ A 1 ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * ( ν λ * 2 + A 4 λ * + A 5 ) ( A 9 λ * + A 10 ) , R ι * = κ 1 Λ A 0 ( γ 1 A 3 + γ 2 ) ( ν λ * + ϕ + μ + ν ψ ) λ * ( ν λ * 2 + A 4 λ * + A 5 ) ( A 9 λ * + A 10 ) .
Further
λ * = β ( I ι 1 * + z I ι 2 * ) N ι * ,
First, the nominator
β ( I ι 1 * + z I ι 2 * ) = β κ 1 Λ A 1 A 3 ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * ( ν λ * 2 + A 4 λ * + A 5 ) ( A 9 λ * + A 10 ) + z κ 1 Λ A 1 ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * ( ν λ * 2 + A 4 λ * + A 5 ) ( A 9 λ * + A 10 )
= β ( A 3 + z ) κ 1 Λ A 1 ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * ( ν λ * 2 + A 4 λ * + A 5 ) ( A 9 λ * + A 10 ) ,
and the denominator is N ι = S ι * + P ι * + E ι * + I ι 1 * + I ι 2 * + R ι * . After substituting the values of S ι * + P ι * + E ι * + I ι 1 * + I ι 2 * + R ι * in Equation (A2), we obtain
N ι = Λ ( ν λ * 2 + A 4 λ + A 5 ) ( A 9 λ * + A 10 ) ( ν λ * + ϕ + μ + ψ ) ( A 9 λ * + A 10 ) + A 12 ( λ * ) ,
where
A 12 ( λ * ) = ( ν λ * + ϕ + μ + ν ψ ) ( A 11 λ * + A 10 ) λ * + κ 1 A 1 A 3 ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * + κ 1 A 1 ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * κ 1 A 1 ( γ 1 A 3 + γ 2 ) ( ν λ * + ϕ + μ + ν ψ )
= ( ν λ * + ϕ + μ + ν ψ ) λ * ( A 13 λ * + A 14 ) ,
where
A 13 = A 11 + κ 1 A 1 A 3 + κ 1 A 1 σ , A 14 = A 10 + μ κ 1 A 1 A 3 + μ κ 1 A 1 + κ 1 A 1 ( γ 1 A 3 + σ 2 ) .
Thus
N ι = Λ [ ( ν λ * + ϕ + μ + ψ ) ( A 9 λ * + A 10 ) + ( ν λ * + ϕ + μ + ψ ) ( A 13 λ * + A 14 ) ] ( ν λ * 2 + A 4 λ * + A 5 ) ( A 9 λ * + A 10 ) .
Hence, from (A3) using (A5) and (A8), to get
λ * = κ 1 A 1 β ( A 3 + z ) ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * [ ( ν λ * + ϕ + μ + ψ ) ( A 9 λ * + A 10 ) + ( ν λ * + ϕ + μ + ν ψ ) ( A 13 λ * + A 14 ) λ * ] ,
λ * κ 1 A 1 β ( A 3 + z ) ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) λ * [ ( ν λ * + ϕ + μ + ψ ) ( A 9 λ * + A 10 ) + ( ν λ * + ϕ + μ + ν ψ ) ( A 13 λ * + A 14 ) λ * ] = 0 ,
λ * = 0 or 1 κ 1 A 1 β ( A 3 + z ) ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) [ ( ν λ * + ϕ + μ + ψ ) ( A 9 λ * + A 10 ) + ( ν λ * + ϕ + μ + ν ψ ) ( A 13 λ * + A 14 ) λ * ] = 0 ,
( ν λ * + ϕ + μ + ψ ) ( A 9 λ * + A 10 ) + ( ν λ * + ϕ + μ + ν ψ ) ( A 13 λ * + A 14 ) λ * ] κ 1 A 1 β ( A 3 + z ) ( σ λ * + μ ) ( ν λ * + ϕ + μ + ν ψ ) .
ν λ * ( A 9 λ * + A 10 ) + ( ϕ + μ + ψ ) ( A 9 λ * + A 10 ) + ν λ * 2 ( A 13 λ * + A 14 ) + ( ϕ + μ + ν ψ ) ( A 13 λ * + A 14 ) λ * κ 1 A 1 β ( A 3 + z ) ( σ λ * + μ ) ν λ * κ 1 A 1 β ( A 3 + z ) ( σ λ * + μ ) ( ϕ + μ + ν ψ ) = 0 .
ν A 9 λ * 2 + ν A 10 λ * + A 9 ( ϕ + μ + ψ ) λ * + A 10 ( ϕ + μ + ψ ) + ν A 13 λ * 3 + ν A 14 λ * 2 + A 13 ( ϕ + μ + ν ψ ) λ * 2 + A 14 ( ϕ + μ + ν ψ ) λ * κ 1 A 1 β ( A 3 + z ) σ λ * 2 μ κ 1 A 1 β ( A 3 + z ) ν λ * ( ϕ + μ + ν ψ ) κ 1 A 1 β ( A 3 + z ) σ λ * μ ( ϕ + μ + ν ψ ) κ 1 A 1 β ( A 3 + z ) = 0 ,
ν A 13 λ * 3 + A 15 λ * 2 + A 16 λ * + A 17 = 0 .
where
A 15 = ν A 9 + ν A 14 + A 13 ( ϕ + μ + ν ψ ) σ ν κ 1 A 1 β ( A 3 + z ) A 16 = ν A 10 + A 9 ( ϕ + μ + ψ ) + A 14 ( ϕ + μ + ν ψ ) κ 1 A 1 β ( A 3 + z ) ( μ ν + σ ( ϕ + μ + ν ψ ) ) A 17 = A 10 ( ϕ + μ + ψ ) ( 1 R p h ) .
Here, R p h is the effective reproduction number given in (15).
Clearly, it is evident that ( I ι ) is given by the positive real roots of the polynomial (A13). The number of possible positive real roots of the cubic polynomial (A13) depends on the signs of A 15 , A 16 , and A 17 . The following theorem can be established. □

Appendix B

Proof of Theorem 7 
To apply the method in presented in Theorem 4.1 of [6], the following simplification and change of variables are made on the model system (23). Let
S ι = x 1 , P ι = x 2 , E ι = x 3 , I ι 1 = x 4 , I ι 2 = x 5 , R ι = x 6 ,
so that
N ι = x 1 + x 2 + x 3 + x 4 + x 5 + x 6 ,
Also, by utilizing the vector notation X = ( x 1 , x 2 , x 3 , x 4 , x 5 , x 6 ) T , then, the model system (23) becomes d X d t = ( f 1 , f 2 , f 3 , f 4 , f 5 , f 6 ) T as follows
d x 1 d t = f 1 = Λ + ϕ x 2 ( λ + ψ + μ ) x 1 , d x 2 d t = f 2 = ψ x 1 ( ν λ + ϕ + μ ) x 2 , d x 3 d t = f 3 = λ x 1 + ν λ x 2 + σ λ x 6 ( κ 1 + κ 2 + μ ) x 3 , d x 4 d t = f 4 = κ 1 x 3 + η 2 x 5 ( η 1 + γ 1 + μ + δ 1 ) x 4 , d x 5 d t = f 5 = κ 2 x 3 + η 1 x 4 η 2 + γ 2 + μ + δ 2 x 5 , d x 6 d t = f 6 = γ 1 x 4 + γ 2 x 5 ( σ λ + μ ) x 6 ,
where
λ = β ( x 4 + z x 5 ) N ι ,
The Jacobian of the system (A17), evaluated at the T F E , T 0 (denoted by J ( T 0 ) ), is given by
J ( T 0 ) = ( μ + ψ ) ϕ 0 β x 1 0 N ι 0 β z x 1 0 N ι 0 0 ψ ( ϕ + μ ) 0 ν β x 2 0 N ι 0 ν β z x 2 0 N ι 0 0 0 0 A 1 β ( x 1 0 + ν x 2 0 ) N ι 0 β z ( x 1 0 + ν x 2 0 ) N ι 0 0 0 0 κ 1 A 2 η 2 0 0 0 κ 2 η 1 A 3 0 0 0 0 γ 1 γ 2 μ ,
where A 1 , A 2 , and A 3 are as in (13), from which it has been shown in (14) that the effective reproduction number, R p h is given by
R p h = β ( ϕ + μ + ν ψ ) ( κ 1 A 3 + κ 2 η 2 ) + z ( κ 2 A 2 + η 1 κ 1 ) A 1 ( A 2 A 3 η 1 η 2 ) ( ψ + ϕ + μ )
Consider the case when R p h = 1 . Suppose, further, that β = β 0 is chosen as a bifurcation parameter, since R p h is often inconvenient to use directly as bifurcation parameter. Solving for β gives R p h = 1 when
β = β 0 = A 1 ( A 2 A 3 η 1 η 2 ) ( ψ + ϕ + μ ) β ( ϕ + μ + ν ψ ) ( κ 1 A 3 + κ 2 η 2 ) + z ( κ 2 A 2 + η 1 κ 1 )
The linearized system of the transformed model system (23) with β = β 0 chosen as a bifurcation parameter has a simple zero eigenvalue. We then calculate the right eigenvector W and the left eigenvector V which are associated with the zero eigenvalue of the Jacobian of (A19) at (denoted by J β 0 ) chosen such that J ( T 0 ) W = 0 and V J ( T 0 ) = 0 with V W = 1 , where
W = [ w 1 , w 2 , w 3 , w 4 , w 5 , w 6 ] , V = [ v 1 , v 2 , v 3 , v 4 , v 5 , v 6 ] .
Then
J ( T 0 ) W i = ( μ + ψ ) ϕ 0 β x 0 N ι 0 z β x 0 N ι 0 0 ψ ( ϕ + μ ) 0 ν β x 0 N ι 0 z ν β x 0 N ι 0 0 0 0 A 1 β x 1 0 + ν β x 2 0 N ι 0 z β x 1 0 + ν z β x 2 0 N ι 0 0 0 0 κ 1 A 2 η 2 0 0 0 κ 2 η 1 A 3 0 0 0 0 γ 1 γ 2 μ w 1 w 2 w 3 w 4 w 5 w 6 = 0 0 0 0 0 0 .
i.e.,
( μ + ψ ) w 1 + ϕ w 2 β x 1 0 N ι 0 w 4 z β x 1 0 N ι 0 w 5 = 0 , ψ w 1 ( μ + ϕ ) w 2 ν β x 2 0 N ι 0 w 4 z ν β x 2 0 N ι 0 w 5 = 0 , A 1 w 3 + β x 1 0 + ν β x 2 0 N ι 0 w 4 + z β x 1 0 + ν z β x 2 0 N ι 0 w 5 = 0 , κ 1 w 3 + A 2 w 4 + η 2 w 5 = 0 , κ 2 w 3 + η 1 w 4 A 3 w 5 = 0 , γ 1 w 4 + γ 2 w 5 μ w 6 = 0 .
Solving (A23), gives
w 1 = A 1 ( η 2 + z A 2 ) { ϕ ( β x 1 0 + ν β x 2 0 ) + β x 1 0 ( ( ψ + μ ) ( ϕ + μ ) ψ ϕ ) } w 5 ( ( ψ + μ ) ( ϕ + μ ) ψ ϕ ) ( A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) κ 1 ) ( ψ + μ ) < 0 , w 2 = A 1 ( η 2 + z A 2 ) x 1 0 + ν x 2 0 ) w 5 ( ( ψ + μ ) ( ϕ + μ ) ψ ϕ ) ( A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) κ 1 ) < 0 , w 3 = β ( x 1 0 + ν x 2 0 ) ( η 2 + z A 2 ) w 5 A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) κ 1 > 0 , w 4 = { η 2 A 1 N ι 0 + κ 1 β ( x 1 0 + ν x 2 0 ) z } w 5 A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) κ 1 > 0 , w 5 > 0 ( can take any value ) , w 6 = { A 1 N ι 0 γ 1 η 2 + ( A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) κ 1 ) γ 2 + β ( x 1 0 + ν x 2 0 ) κ 1 γ 1 z } w 5 μ ( A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) κ 1 ) > 0 .
Similarly, calculating the left eigenvector V = ( v 1 , v 2 , v 3 , v 4 , v 5 , v 6 ) T with V J ( T 0 ) = 0 , gives
J ( T 0 ) V i = ( v 1 , v 2 , v 3 , v 4 , v 5 , v 6 ) ( μ + ψ ) ϕ 0 β x 0 N ι 0 z β x 0 N ι 0 0 ψ ( ϕ + μ ) 0 ν β x 0 N ι 0 z ν β x 0 N ι 0 0 0 0 A 1 β x 0 + ν β x 0 N ι 0 z β x 0 + ν z β x 0 N ι 0 0 0 0 κ 1 A 2 η 2 0 0 0 κ 2 η 1 A 3 0 0 0 0 γ 1 γ 2 μ = 0 0 0 0 0 0 .
i.e.,
( μ + ψ ) v 1 + ψ v 2 = 0 , ϕ v 1 ( μ + ϕ ) v 2 = 0 , A 1 v 3 + κ 1 v 4 + κ 2 v 5 = 0 , β x 1 0 N ι 0 v 1 ν β x 2 0 N ι 0 v 2 + β x 1 0 + ν β x 2 0 N ι 0 v 3 A 2 v 4 + η 1 v 5 + γ 1 v 6 = 0 , z β x 1 0 N ι 0 v 1 ν z β x 2 0 N ι 0 v 2 + z β x 1 0 + ν z β x 2 0 N ι 0 v 3 η 2 v 4 A 3 v 5 + γ 2 v 6 = 0 , μ 1 v 6 = 0 .
Solving (A26), gives
v 1 = v 2 = v 6 = 0 , v 3 = κ 2 { κ 1 β ( x 1 0 + ν x 2 0 ) + ( A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) ) } v 5 A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) > 0 , v 4 = β ( x 1 0 + ν x 2 0 ) ) κ 2 v 5 A 1 A 2 N ι 0 β ( x 1 0 + ν x 2 0 ) κ 1 > 0 , v 5 > 0 ( can take any value ) .
The local dynamics of model system (23) around TB free equilibrium T 0 are totally calculated by a and b represented as
a = k , i , j = 1 n v k w i w j 2 f k ( 0 , 0 ) x i x j ,
b = k , i = 1 4 v k w i 2 f k ( 0 , 0 ) x i φ ( 0 , 0 ) .
Now, to determine the coefficients a and b described in Theorem 4.1 of Castillo Chavez and Song [6] as follow:
a = 2 v 3 β N ι 0 ( w 1 + w 2 ν + w 6 σ ) ( w 4 + w 5 z ) .
and
b = 2 f 3 x 4 0 β 0 ( 0 , 0 ) + 2 f 3 x 5 0 β 0 ( 0 , 0 ) = x 1 0 + ν x 2 0 + σ x 6 0 N ι 0 + z x 1 0 + ν z x 2 0 + z σ x 6 0 N ι 0 > 0 .
Hence, since a > 0 and b > 0 , in this situation, it follows from Theorem 4.1 of Castillo Chavez and Song that the system model (23) or the transformed model (A17), will undergo a phenomenon of backward bifurcation when R p h = 1 . Therefore, establishing that TB endemic equilibrium is locally asymptotically stable if R p h < 1 .

References

  1. World Health Organization. Global Tuberculosis Report; WHO: Geneva, Switzerland, 2019; Volume 214. [Google Scholar]
  2. Ullah, S.; Khan, M.A.; Farooq, M.; Gul, T. Modeling and analysis of tuberculosis (TB) in Khyber Pakhtunkhwa, Pakistan. Math. Comput. Simul. 2019, 165, 181–199. [Google Scholar] [CrossRef]
  3. Sudre, P.; Ten Dam, G.; Kochi, A. Tuberculosis: A global overview of the situation today. Bull. World Health Organ. 1992, 70, 149. [Google Scholar] [PubMed]
  4. Dolin, P.J.; Raviglione, M.C.; Kochi, A. Global tuberculosis incidence and mortality during. Bull. World Health Organ. 1994, 72, 213. [Google Scholar] [PubMed]
  5. World Health Organization. Global Tuberculosis Report 2020; WHO: Geneva, Switzerland, 2020. [Google Scholar]
  6. Castillo-Chavez, C.; Song, B. Dynamical models of tuberculosis and their applications. Math. Biosci. Eng. 2004, 1, 361. [Google Scholar] [CrossRef]
  7. Kasereka Kabunga, S.; Doungmo Goufo, E.F.; Ho Tuong, V. Analysis and simulation of a mathematical model of tuberculosis transmission in Democratic Republic of the Congo. Adv. Differ. Equ. 2020, 2020, 642. [Google Scholar] [CrossRef]
  8. Adebiyi, A.O. Mathematical Modeling of the Population Dynamics of Tuberculosis. Master’s Thesis, University of the Western Cape, Cape Town, South Africa, 2016. [Google Scholar]
  9. Jaramillo, E. Tuberculosis and stigma: Predictors of prejudice against people with tuberculosis. J. Health Psychol. 1999, 4, 71–79. [Google Scholar] [CrossRef]
  10. Daniel, T.M.; Bates, J.H.; Downes, K.A. History of tuberculosis. In Tuberculosis: Pathogenesis, Protection, and Control; Wiley Online Library: Hoboken, NJ, USA, 1994; pp. 13–24. [Google Scholar]
  11. Khan, M.A.; Ahmad, M.; Ullah, S.; Farooq, M.; Gul, T. Modeling the transmission dynamics of tuberculosis in Khyber Pakhtunkhwa Pakistan. Adv. Mech. Eng. 2019, 11, 1687814019854835. [Google Scholar] [CrossRef] [Green Version]
  12. Khajanchi, S.; Das, D.K.; Kar, T.K. Dynamics of tuberculosis transmission with exogenous reinfections and endogenous reactivation. Phys. A Stat. Mech. Appl. 2018, 497, 52–71. [Google Scholar] [CrossRef]
  13. Zhang, J.; Li, Y.; Zhang, X. Mathematical modeling of tuberculosis data of China. J. Theor. Biol. 2015, 365, 159–163. [Google Scholar] [CrossRef]
  14. Bar, B. Tuberculose et grossesse. Bull. Acad. Méd. 1922, 219. [Google Scholar]
  15. World Health Organization. WHO Consolidated Guidelines on Tuberculosis: Tuberculosis Preventive Treatment; WHO: Geneva, Switzerland, 2020. [Google Scholar]
  16. Zhou, X.; Shi, X.; Cheng, H. Modelling and stability analysis for a tuberculosis model with healthy education and treatment. Comput. Appl. 2015, 32, 245–260. [Google Scholar] [CrossRef]
  17. Xiang, H.; Song, N.N.; Huo, H.F. Modelling effects of public health educational campaigns on drinking dynamics. J. Biol. Dyn. 2016, 10, 164–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Del Valle, S.; Hethcote, H.; Hyman, J.M.; Castillo-Chavez, C. Effects of behavioural changes in a smallpox attack model. Math. Biosci. 2005, 195, 228–251. [Google Scholar] [CrossRef] [PubMed]
  19. Xiang, H.; Zou, M.X.; Huo, H.F. Modeling the effects of health education and early therapy on tuberculosis transmission dynamics. Int. J. Nonlinear Sci. Numer. 2019, 20, 243–255. [Google Scholar] [CrossRef]
  20. Huo, H.F.; Zou, M.X. Modelling effects of treatment at home on tuberculosis transmission dynamics. Appl. Math. Model. 2016, 40, 9474–9484. [Google Scholar] [CrossRef]
  21. Abimbade, S.F.; Olaniyi, S.; Ajala, O.A.; Ibrahim, M.O. Optimal control analysis of a tuberculosis model with exogenous reinfection and incomplete treatment. Optim. Control Appl. Methods 2020, 41, 2349–2368. [Google Scholar] [CrossRef]
  22. Ojo, M.M.; Peter, O.J.; Goufo, E.F.D.; Panigoro, H.S.; Oguntolu, F.A. Mathematical model for control of tuberculosis epidemiology. J. Appl. Math. Comput. 2022, 1–19. [Google Scholar] [CrossRef]
  23. Egonmwan, A.O.; Okuonghae, D. Analysis of a mathematical model for tuberculosis with diagnosis. J. Appl. Math Comput. 2019, 59, 129–162. [Google Scholar] [CrossRef]
  24. Mustapha, U.T.; Idris, B.; Musa, S.S.; Yusuf, A. Mathematical modeling and analysis of mycobacterium tuberculosis transmission in humans with hospitalization and reinfection. J. Appl. Math. Comput. Mech. 2022, 21, 55–66. [Google Scholar] [CrossRef]
  25. Kar, T.K.; Mondal, P.K. Global dynamics of a tuberculosis epidemic model and the influence of backward bifurcation. J. Math. Model. Algorithms 2012, 11, 433–459. [Google Scholar] [CrossRef]
  26. Okuonghae, D.; Omosigho, S.E. Analysis of a mathematical model for tuberculosis: What could be done to increase case detection. J. Theor. Biol. 2011, 269, 31–45. [Google Scholar] [CrossRef]
  27. Tang, H.; Li, M.; Yan, X.; Lu, Z.; Jia, Z. Modeling the dynamics of drug spreading in China. Int. J. Environ. Res. Public Health 2021, 18, 288. [Google Scholar] [CrossRef] [PubMed]
  28. Olaniyi, S.; Falowo, O.D.; Okosun, K.O.; Mukamuri, M.; Obabiyi, O.S.; Adepoju, O.A. Effect of saturated treatment on malaria spread with optimal intervention. Alex. Eng. J. in press. 2022. [CrossRef]
  29. Olaniyi, S.; Obabiyi, O.S.; Okosun, K.O.; Oladipo, A.T.; Adewale, S.O. Mathematical modelling and optimal cost-effective control of COVID-19 transmission dynamics. Eur. Phys. J. Plus. 2020, 135, 938. [Google Scholar] [CrossRef] [PubMed]
  30. Cooke, K.L.; Van Den Driessche, P. Analysis of an SEIRS epidemic model with two delays. J. Math. Biol. 1996, 35, 240–260. [Google Scholar] [CrossRef]
  31. Obasi, C.; Mbah, G.C.E. On the stability analysis of a mathematical model of Lassa fever disease dynamics. J. Nig. Soc. Math. Biol. 2019, 2, 135–144. [Google Scholar]
  32. Srivastava, A.; Srivastava, P.K. Nonlinear dynamics of a SIRI model incorporating the impact of information and saturated treatment with optimal control. Eur. Phys. J. Plus. 2022, 137, 1–25. [Google Scholar] [CrossRef]
  33. Rota, G.C.; Birkhoff, G. Ordinary Differential Equations; Blaisdell Publishing Company: Waltham, MA, USA, 1969. [Google Scholar]
  34. Hethcote, H.W. The mathematics of infectious diseases. SIAM Rev. 2000, 42, 599–653. [Google Scholar] [CrossRef] [Green Version]
  35. Van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef]
  36. Falowo, O.D.; Olaniyi, S.; Oladipo, A.T. Optimal control assessment of Rift Valley fever model with vaccination and environmental sanitation in the presence of treatment delay. Model. Earth Syst. Environ. 2022, 1–15. [Google Scholar] [CrossRef]
  37. Keeling, M.J.; Rohani, P. Temporally forced models. In Modeling Infectious Diseases in Humans and Animals; Princeton University Press: Princeton, NJ, USA, 2011; pp. 155–189. [Google Scholar]
  38. Castillo-Chavez, C.; Feng, Z.; Huang, W. On the computation of R0 and its role on. In Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An Introduction; Springer: Berlin/Heidelberg, Germany, 2002; Volume 1, p. 229. [Google Scholar]
  39. Islam, M.S.; Ira, J.I.; Kabir, K.A.; Kamrujjaman, M. COVID-19 Epidemic compartments model and Bangladesh. Preprint 2020, 10. [Google Scholar] [CrossRef]
  40. Olaniyi, S. Dynamics of Zika virus model with nonlinear incidence and optimal control strategies. Appl. Math. Inf. Sci. 2018, 12, 969–982. [Google Scholar] [CrossRef]
  41. Chitnis, N.; Hyman, J.M.; Cushing, J.M. Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model. Bull. Math. Biol. 2008, 70, 1272–1296. [Google Scholar] [CrossRef] [PubMed]
  42. Song, B.; Castillo-Chavez, C.; Aparicio, J.P. Tuberculosis models with fast and slow dynamics: The role of close and casual contacts. Math. Biosci. 2002, 180, 187–205. [Google Scholar] [CrossRef] [PubMed]
  43. Gomes, M.G.M.; Aguas, R.; Lopes, J.S.; Nunes, M.C.; Rebelo, C.; Rodrigues, P.; Struchiner, C.J. How host heterogeneity governs tuberculosis reinfection? Proc. R. Soc. B Biol. Sci. 2012, 279, 2473–2478. [Google Scholar] [CrossRef] [Green Version]
  44. Rodrigues, P.; Gomes, M.G.M.; Rebelo, C. Drug resistance in tuberculosis—A reinfection model. Theor. Popul. Biol. 2007, 71, 196–212. [Google Scholar] [CrossRef] [Green Version]
  45. Bowong, S.; Tewa, J.J. Global analysis of a dynamical model for transmission of tuberculosis with a general contact rate. Commun. Nonlinear. Sci. Numer. Simul. Commun Nonlinear Sci. 2010, 15, 3621–3631. [Google Scholar] [CrossRef]
Figure 1. The flow diagram of TB model with public health education and hospital treatment.
Figure 1. The flow diagram of TB model with public health education and hospital treatment.
Axioms 11 00723 g001
Figure 2. The bifurcation diagram of force of infection against R p h which illustrates a backward bifurcation for the system model (23).
Figure 2. The bifurcation diagram of force of infection against R p h which illustrates a backward bifurcation for the system model (23).
Axioms 11 00723 g002
Figure 3. The relationship between ψ (a), η 1 (b), κ 1 (c), and κ 2 (d) and the effective reproduction number R p h .
Figure 3. The relationship between ψ (a), η 1 (b), κ 1 (c), and κ 2 (d) and the effective reproduction number R p h .
Axioms 11 00723 g003
Figure 4. Simulations of system (23) showing the behaviour of susceptible and educated susceptible. Parameters used are Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Figure 4. Simulations of system (23) showing the behaviour of susceptible and educated susceptible. Parameters used are Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Axioms 11 00723 g004
Figure 5. Impact of awareness rate ( ψ ) on S ι , E ι , I ι 1 , and I ι 2 , varying ψ . Other parameters are given as Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Figure 5. Impact of awareness rate ( ψ ) on S ι , E ι , I ι 1 , and I ι 2 , varying ψ . Other parameters are given as Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Axioms 11 00723 g005
Figure 6. Impact at which educated susceptible lose awareness ( ϕ ) on I ι 1 , and I ι 2 , varying ϕ . Other parameters are given as Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Figure 6. Impact at which educated susceptible lose awareness ( ϕ ) on I ι 1 , and I ι 2 , varying ϕ . Other parameters are given as Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Axioms 11 00723 g006
Figure 7. Impact of reduction of infection rate as a result of I 2 ( z ) on I ι 1 , and I ι 2 , varying z . Other parameters are given as Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Figure 7. Impact of reduction of infection rate as a result of I 2 ( z ) on I ι 1 , and I ι 2 , varying z . Other parameters are given as Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Axioms 11 00723 g007
Figure 8. Impact of progression rate ( η 1 ) from I ι 1 to I ι 2 on I ι 1 , and I ι 2 , varying η 1 . Other parameters are given as Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Figure 8. Impact of progression rate ( η 1 ) from I ι 1 to I ι 2 on I ι 1 , and I ι 2 , varying η 1 . Other parameters are given as Λ = 450,862 , β = 0.86 , ν = 0.25 , ϕ = 0.2 , σ = 0.65 , z = 0.11 , μ = 0.02041 , ψ = 0.5 , η 1 = 0.02 , η 2 = 0.02 , γ 1 = 0.5 , γ 2 = 0.6 , κ 1 = 0.03 , κ 2 = 0.3 , δ 1 = 0.2 , δ 2 = 0.12 .
Axioms 11 00723 g008
Table 1. Interpretation of the Variables and Parameters of the TB Model with Public Education and hospital Treatment (23).
Table 1. Interpretation of the Variables and Parameters of the TB Model with Public Education and hospital Treatment (23).
State VariablesExplanation
S ι The number of individuals who are susceptible
P ι The number of individuals who are educated susceptible
E ι The number of individuals who are exposed to TB
I ι 1 The number of infected individuals at home
I ι 2 The number of infected individuals at hospital
R ι The number of individuals who have  recovered
N ι Human population size
Parameters
Λ Inflow of recruitment into susceptible class
μ Per capita natural death rate of humans
ψ Information dissemination (awareness rate)
ϕ Rate at which the educated susceptible become susceptible
ν Reduction of infection rate as a result of awareness
β Transmission rate
zReduction of infection rate as a result of individuals in the hospital
κ 1 Progression rate from E ι to I ι 1
κ 2 Progression rate from E ι to I ι 2
γ 1 recovery rate for I ι 1
γ 2 recovery rate for I ι 2
σ Mod. parameter for re-infection among the recovered individuals
η 1 Progression rate from I ι 1 to I ι 2
η 2 Progression rate from I ι 2 to I ι 1
δ 1 TB induced mortality rate for I ι 1
δ 2 TB induced mortality rate for I ι 2
Table 2. Sensitivity Indices of the Model Parameters in Relation to R p h .
Table 2. Sensitivity Indices of the Model Parameters in Relation to R p h .
ParametersDescriptionSensitivity Index
μ Per capita natural death rate of humans−0.6970
ψ Information dissemination (awareness rate)−0.2554
ϕ The rate at which susceptible individuals lose awareness+0.2454
ν Reduction in risk of infection due to awareness+0.1103
β Transmission rate for contact with I 1 +1.0000
zReduction of infection rate as a result of I 2 +0.4476
κ 1 Progression rate from E ι to I 1 +0.4228
κ 2 Progression rate from E ι to I 2 −0.2147
γ 1 Treatment rate for I ι 1 −0.2573
γ 2 Treatment rate for I ι 2 +0.5570
η 1 Progression rate from I ι 1 to I ι 2 −0.6832
η 2 Progression rate from I ι 2 to I ι 1 +0.1845
δ 1 TB induced mortality rate for I ι 1 −0.2058
δ 2 TB induced mortality rate for I ι 2 −0.3714
Table 3. The Parameters and Baseline Values of the Model with Public Health Education and Hospital Treatment (23).
Table 3. The Parameters and Baseline Values of the Model with Public Health Education and Hospital Treatment (23).
ParametersBaseline ValuesRangesReferences
Λ 3,768,410 year 1 [3,000,000, 4,000,000][42]
μ 0.02041 year 1 [0.0143, 0.03][23]
ψ Variable[0–1]Assumed
ϕ Variable[0–1]Assumed
ν Variable[0–1]Assumed
β VariableVariedAssumed
σ 0.25[0–1][25,43]
z0.11[0–0.9][16]
κ 1 Variable[0–1]Assumed
κ 2 Variable[0–1]Assumed
γ 1 0.09 year 1 [0–1][20]
γ 2 0.72 year 1 [0–1][20]
η 1 Variable[0–1]assumed
η 2 Variable[0–1]assumed
δ 1 0.2 year 1 -[44]
δ 2 0.02 year 1 -[45]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sulayman, F.; Abdullah, F.A. Dynamical Behaviour of a Modified Tuberculosis Model with Impact of Public Health Education and Hospital Treatment. Axioms 2022, 11, 723. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11120723

AMA Style

Sulayman F, Abdullah FA. Dynamical Behaviour of a Modified Tuberculosis Model with Impact of Public Health Education and Hospital Treatment. Axioms. 2022; 11(12):723. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11120723

Chicago/Turabian Style

Sulayman, Fatima, and Farah Aini Abdullah. 2022. "Dynamical Behaviour of a Modified Tuberculosis Model with Impact of Public Health Education and Hospital Treatment" Axioms 11, no. 12: 723. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11120723

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop