Oalib Open Access Library Journal 2333-9721 2333-9705 Scientific Research Publishing 10.4236/oalib.1115019 Oalib-150342 Article Biomedical Life Sciences Business Economics Chemistry Materials Science Computer Science Communications Earth Environmental Sciences Engineering Medicine Healthcare Physics Mathematics Social Sciences Humanities Mathematical Modelling and Analysis of Measles Control Using Vaccination Konlan Musah 1 Gunu Rauf Ibrahim Mohamed 1 Iddrisu Abdul-Karim 1 2 Department of Mathematics and Statistics, University of Energy and Natural Resources, Sunyani, Ghana Department of Statistics, University of Botswana, Gaborone, Botswana

The authors declare no conflicts of interest.

 28 02 2026 02 2026 13 03 1 18 11 02 2026 20 03 2026 23 03 2026 © 2026 by the authors and Scientific Research Publishing Inc. 2026 This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license ( https://creativecommons.org/licenses/by/4.0/ ). https://doi.org/10.4236/oalib.1115019

In this paper, the analysis of a mathematical model for the control of measles using vaccination is presented. The disease reproductive number (

R

0

) is computed and used to derive an expression of a reduction factor in the disease spread resulting from the implementation of immunization program in the community. The Routh-Hurwitz stability conditions are adopted to analyze the model equilibria in terms of their local stability. Lyapunov function candidates are used to determine the global stability status of both the measles-free and measles-persistent equilibrium states. Elasticity index of each parameter embedded into

R

0

is computed to help ascertain the significance of the model parameter contribution to measles epidemics in the society. Furthermore, population simulations are performed to demonstrate how changes in some parameter values influence the dynamical evolution of our model population sub-classes.

 Measles Vaccination Campaign Reproductive Number with/without Vaccination Local and Global Stability Sensitivity Analysis 1. Introduction

Measles is globally known as the most fatal and contagious human disease [1][2]. It is one of the diseases with devastating impact among children under the age of five. Measles is caused by a virus called morbillivirus and is primarily transmitted among humans through aerosolized respiratory droplets [3]-[5]. The initial symptoms of measles include but are not limited to: runny nose, cough, rash all over the body, fever, malaise and conjunctivitis which usually manifest between ten to fourteen days after one’s exposure to the disease [3][6]. If left unattended, severe measles can lead to complications such as: encephalitis, blindness, ear infections, severe diarrhea, pneumonia and death [2][7]. Even though an effective and safe vaccine for the control of measles exists, this disease continues to reemerge in many countries including the developed nations [4][6][8][9]. Thus, measles resurgence still remains a global public health issue that requires attention [2][10]. Regular measles vaccination campaign for children and young adults in measles endemic countries is currently the recommended control strategy for reducing measles disease burden [11]. In the last twenty-three years, about 60.3 million measles deaths have been averted by vaccination. However, between 2022 and 2023, there was an estimated 20% increase in measles cases globally with the number of countries/territories that reported large scale outbreaks increased from 36 to 57 [12]. Key contributing factors to these outbreaks especially in Africa include: inadequate measles surveillance systems, low measles immunization coverage and poor health care systems. In Ghana for example, a total of 2282 measles cases were reported in the year 2023 [12]. Several studies have reported that measles surveillance remains suboptimal in many countries [13]. Zumah et al. [14] evaluated and rated the measles disease surveillance system’s performance in the Bono region of Ghana as suboptimal. Some data from their study is shown in Table 1. Mathematical modelling has been used as an effective tool for providing meaningful insights into infectious disease dynamics and control [2][9][15][16]. Thus, some researchers have developed dynamical models for measles epidemic. Opoku et al. [4] developed a dynamical model for measles control that took into account maternal antibodies protecting new born babies from measles infections and double-dose vaccination program. Their study results suggested that measles first dose vaccination is a potential means of combating measles epidemics. Garba et al. [5] used a nine-compartment model in investigating the significance of simultaneously combining vaccination and treatment in the control of measles epidemics. Wireko et al. [17] proposed a fractional order dynamical system for studying measles. Alemneh and Belay [18] analyzed measles transmission model that took into account the presence of measles virus in the environment. James et al. [6] formulated a deterministic mathematical model to provide insight into measles dynamics in Nigeria. In another study conducted by James et al. [2], it was demonstrated that reducing the contact rate between a susceptible and an

Table 1. Measles surveillance data from the bono region of Ghana.

Year 2019 2018 2017 2016 2015 2014 2013 2012 2011 2010
SMC 184 222 220 166 127 165 77 27 30 19
LCMC 4 0 2 5 0 1 5 0 2 0
% MCV1 - - 89.4 92.3 93.3 94.0 88.2 99.6 96.1 95.6

SMC = Suspected Measles Cases, LCMC = Laboratory Confirmed Measles Cases and %MCV1 = percentage of children vaccinated against the first measles dose vaccine. Source [14].

infectious person and increasing the vaccination rate of an effective measles vaccine were the most effective control measures against measles epidemics. A mathematical model analysis presented by Kuddus et al. [19] in Bangladesh revealed that contact rate is the most significant parameter to the spread of measles. Authors in [20] investigated the infection dynamics and control of measles in Pakistan using a dynamical system based on mass action principle. Their numerical simulation results suggested that improving vaccine efficiency and coverage rate will translate into a drastic reduction in the spread of measles in the country. Researchers in [9] examined the dynamics of measles by extending an SEIR compartmental model to include vaccination and treatment. Their new model was simulated using reported measles data from Indonesia. In this current study, we employ mathematical modeling approach to gain more insights into measles infection dynamics in the presence of vaccination campaign.

 2. Construction of Measles Vaccination Model

To formulate a dynamical model for measles control using vaccination, the population under study is stratified into five sub-population classes. Namely: susceptibles (S), vaccinated (V), exposed (E), infected (I) and recovered (R) classes. Following this description, N ( t ) , the population size at any given time t is made up of S ( t ) , V ( t ) , E ( t ) , I ( t ) and R ( t ) . The susceptible population under consideration is generated at a constant rate π . Some of these susceptible individuals get vaccinated against measles at rate ρ . The model assumes that a portion of those vaccinated attain immunity against measles and move to recovery class at rate τ . Following the waning of the vaccine, some vaccinated individuals return to susceptible sub-class at rate η . Other susceptible individuals get exposed to measles at a force of infection β I N , where β is the effective contact rate between a susceptible and an infected individual. The exposed people progress to infected class at rate α . The infected individuals recover at rate γ . The model assumes that recovered people attain permanent immunity against further measles attacks [8], hence, there is no movement from the recovery class back to susceptible class. The parameters denoting the transfer and removal rates from the model compartments are described in Table 2.

Expressing Figure 1Figure 1 in the form of equations gives:

{ d S d t = π + η V − ( β I N + ρ + μ ) S d V d t = ρ S − ( η + τ + μ ) V d E d t = β I S N − ( μ + α ) E d I d t = α E − ( μ + δ + γ ) I d R d t = γ I + τ V − μ R

Table 2. Measles model parameter description and initial conditions.

Parameter Description Value Source
π Human recruitment rate 123.7 Computed
μ Natural human mortality rate 0.015 [ 17 ]
ρ Vaccination rate 0.842 Estimated
η Vaccine waning rate 0.6 [ 4 ]
τ Immunity rate of vaccinated individuals 0.5 Assumed
β Effective contact rate 0.000402 Estimated
α Progression rate of exposed individuals to infected class 0.25 [ 4 ]
γ Immunity rate of infected individuals 0.6 [ 4 ]
δ Measles induced mortality rate 0.125 [ 4 ]
State Variable Initial Value
S ( 0 ) Initial susceptible number 1000 Assumed
V ( 0 ) Initial vaccinated number 250 Assumed
E ( 0 ) Initial exposed number 500 Assumed
I ( 0 ) Initial infected number 100 Assumed
R ( 0 ) Initial recovered number 0 Assumed

Figure 1. Illustrative diagram for measles transmission dynamics.

To ease our analysis, we use: g 0 = ( ρ + μ ) , g 1 = ( η + τ + μ ) , g 2 = ( α + μ ) , and g 3 = ( γ + μ + δ ) in the rest of the study.

2.1. Positivity of the Model Solutions

Here, we verify the basic properties of the system of differential equations representing model (1).

Lemma 1.The solutions S ( t ) , V ( t ) , E ( t ) , I ( t ) and R ( t ) are non-negative and bounded ∀   t ≥ 0 whenever the initial value set: { S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) , R ( 0 ) } remains positive.

Proof. First, we consider:

d S d t = π + η V − ( Λ + g 0 ) S ,     Λ = β I N ⇒ d S d t ≥ − ( Λ + g 0 ) S ⇒ ∫ 1 S d S ≥ − ∫ ( Λ + g 0 ) d t ⇒ S ( t ) ≥ S ( 0 ) e − ( Λ + g o t ) ≥ 0

Similarly, the following results can be derived:

V ( t ) ≥ V ( 0 ) e − g 1 t ≥ 0 ,     E ( t ) ≥ E ( 0 ) e − g 2 t ≥ 0 ,     I ( t ) ≥ I ( 0 ) e − g 3 t ≥ 0 and     R ( t ) ≥ R ( 0 ) e − μ t ≥ 0.

□

Lemma 2.The positive set D = { ( S , V , E , I , R ) ∈ ℝ + 5 : S + V + E + I + R ≤ π μ } represents the invariant region of the model system of equations.

Proof. For any given time t ≥ 0 , N is:

N = S + V + E + I + R ⇒ d N d t = d d t ( S + V + E + I + R ) ⇒ d N d t = π − μ N − δ I

Thus, in the absence of measles induced mortality,

d N d t ≤ π − μ N ⇒ ∫ d N N − π μ ≤ − ∫ μ d t ⇒ N ≤ π μ     as     t → + ∞ ⇒ N = S + V + E + I + R ≤ π μ

□

 2.2. The Measles-Free Equilibrium Point (DFE)

The DFE point of system (1) denoted by ξ * is the non-trivial solution of system (4) below with the condition E * = I * = 0

{ π − ( β I * N * + g 0 ) S * = 0 ρ S * − g 1 V * = 0 β I * S * N * − g 2 E * = 0 α E * − g 3 I * = 0 γ I * + τ V * − μ R * = 0 ξ * = ( S * , V * , E * , I * , R * ) = ( g 1 π g 0 g 1 − η ρ , ρ π g 0 g 1 − η ρ , 0 , 0 , τ ρ π μ ( g 0 g 1 − η ρ ) )

The Measles Basic Reproductive Number ( R 0 )

The method of next generating matrix is used to derive R 0 . To achieve this, we first expressed the infectious and infected sub-system of system (1) as, d X d t = ( F − G ) X T . Here, X T denotes the transpose of X = ( E , I ) while F and G represent the rates of generation of new infections and transfers respectively. That is,

F = ( β I S N 0 )     and     G = ( g 2 E − α E + g 3 I )

Evaluating the Jacobian matrices F and G of F and G at the DFE gives respectively:

F = ( 0 β S * N * 0 0 ) = ( 0 β μ g 1 g 0 g 1 − η ρ 0 0 )     and     G = ( g 2 0 − α g 3 )

Using F and G from (6) with S * = g 1 π g 0 g 1 − η ρ and N * = π μ , we obtain F G − 1 given by:

F G − 1 = ( 0 β μ g 1 g 0 g 1 − η ρ 0 0 ) ( 1 g 2 0 α g 2 g 3 1 g 3 ) = ( α β μ g 1 g 2 g 3 ( g 0 g 1 − η ρ ) β μ g 1 g 3 ( g 0 g 1 − η ρ ) 0 0 )

Thus, using the equation | F G − 1 − λ I | = 0 where I is a unit matrix and λ an eigenvalue of F G − 1 , we obtain R 0 as the spectral radius of F G − 1 given by:

R 0 = α β μ g 1 g 2 g 3 ( g 0 g 1 − η ρ )

In epidemiology, the magnitude of R 0 provides an extent of the severity of the epidemics in the community [16]. In other words R 0 < 1 , gives hope that the disease will die out with time. On the other hand, if R 0 > 1 , the outbreak will continue to unfold in the population. Using the parameter values presented in Table 2, we obtain R 0 = 1.903 × 10 − 5 .

 2.3. Stability Analysis

2.3.1. Local Stability of Measles Disease-Free Equilibrium (DFE) State

Theorem 1.The point, ξ * = ( g 1 π g 0 g 1 − η ρ , ρ π g 0 g 1 − η ρ , 0 , 0 , τ ρ π μ ( g 0 g 1 − η ρ ) ) admits a local asymptotic stability(LAS)if R 0 < 1 but becomes unstable when R 0 > 1

Proof. We consider the matrix of partial derivatives of the model (1), J ( x ) , that is,

J ( x ) = ( − ( β I * N * + g 0 ) η 0 − β S * N * 0 ρ − g 1 0 0 0 β I * N * 0 − g 2 β S * N * 0 0 0 α − g 3 0 0 τ 0 γ − μ )

Now, evaluating J ( x ) at the DFE ( ξ * ) gives:

J ( ξ * ) = ( − g 0 η 0 − β g 1 μ g 0 g 1 − η ρ 0 ρ − g 1 0 0 0 0 0 − g 2 β g 1 μ g 0 g 1 − η ρ 0 0 0 α − g 3 0 0 τ 0 γ − μ )

It can be easily observed that the matrix in (10) admits three negative eigenvalues: λ 5 = − μ , λ 2 = − g 1 , λ 1 = − g 0 , while the other eigenvalues λ 3 and λ 4 can be obtained from the matrix in (11):

J 1 = ( − g 2 β g 1 μ g 0 g 1 − η ρ α − g 3 )

whose characteristic equation is:

λ 2 + Δ 1 λ + Δ 0 = 0

where:

Δ 1 = g 2 + g 3 ,     Δ 0 = g 2 g 3 ( 1 − R 0 )

Since Δ 1 > 0 and Δ 0 > 0 if R 0 < 1 , matrix J 1 will have negative eigenvalues. Hence, according to the Routh-Hurwitz stability conditions, our proposed model disease-free state ( ξ * ) is a local asymptotic stable equilibrium whenever R 0 < 1 . □

2.3.2. Global Asymptotic Stability

For the purpose of investigating the asymptotic stability of ξ * within its global neighborhood, we define a positive bi-variate function as:

L = α E + g 2 I

Differentiating L gives

d d t ( L ) = d d t ( α E + g 2 I ) = α d E d t + g 2 d I d t = α ( β S I N − g 2 E ) + g 2 ( α E − g 3 I ) = g 2 g 3 ( R 0 − 1 ) I = − g 2 g 3 ( 1 − R 0 ) I

It follows from (15) that d L d t ≤ 0 whenever R 0 ≤ 1 and d L d t = 0 only when I = 0 . This therefore implies that L is a Lyapunov function. Thus, it follows from the Lyapunov version of the LaSale’s Invariant Principle [21] that the measles DFE ( ξ * ) has a GAS if R 0 ≤ 1 but unstable otherwise.

 2.4. Existence of a Unique Measles Endemic Equilibrium Point (EE)

This section is dedicated to the derivation of a unique measles persistence equilibrium state. Let ξ * * = ( S * * , V * * , E * * , I * * , R * * ) be the measles persistent equilibrium (EE) point for model (1) that is ξ * * satisfies the system:

{ π + η V * * − ( β N * * I * * + g 0 ) S * * = 0 ρ S * * − g 1 V * * = 0 β N * * I * * S * * − g 2 E * * = 0 α E * * − g 3 I * * = 0 γ I * * + τ V * * − μ R * * = 0

Solving for S * * , V * * , E * * , I * * and R * * yields, the system:

{ S * * = g 1 π ( g 0 g 1 − η ρ ) R 0 V * * = π ρ ( g 0 g 1 − η ρ ) R 0 E * * = π g 2 ( 1 − 1 R 0 ) I * * = π ( g 0 g 1 − η ρ ) β g 1 μ ( R 0 − 1 ) R * * = β g 1 π μ τ ρ + γ π R 0 ( g 0 g 1 − η ρ ) ( R 0 − 1 ) β g 1 μ 2 R 0 ( g 0 g 1 − η ρ )

It is clear from (17), that ξ * * only exists if R 0 > 1 .

2.4.1. Local Stability Analysis of the Measles Persistent Equilibrium State (EE)

To establish the local stability of the model EE, we consider again the jacobian matrix:

J ( x ) = ( − ( β I * * N * * + g 0 ) η 0 − β S * * N * * 0 ρ − g 1 0 0 0 β I * * N * * 0 − g 2 β S * * N * * 0 0 0 α − g 3 0 0 τ 0 γ − μ )

Now, evaluating J ( x ) at ( ξ * * ) gives:

J ( ξ * * ) = ( − ( m 0 + g 0 ) η 0 − m 1 0 ρ − g 1 0 0 0 m 0 0 − g 2 m 1 0 0 0 α − g 3 0 0 τ 0 γ − μ ) m 0 = 1 g 1 ( g 0 g 1 − η ρ ) ( R 0 − 1 )

The equilibrium state ξ * * will be stable if the matrix in (19) has negative eigenvalues. Expanding the matrix in (19) along the last column gives one negative eigenvalue − μ . Now, the properties of the remaining eigenvalues can be obtained from the reduced matrix

J 2 = ( − ( m 0 + g 0 ) η 0 − m 1 ρ − g 1 0 0 m 0 0 − g 2 m 1 0 0 α − g 3 )

As in [22]-[24], all the eigenvalues of matrix J 2 will be negative or have negative real parts if the trace of ( J 2 ) < 0 and the determinant of ( J 2 > 0 ). It is not hard to see from (20) that

{ trace ( J 2 ) = − ( g 0 + g 1 + g 2 + g 3 ) + 1 g 1 ( g 0 g 1 − η ρ ) ( 1 − R 0 ) determinant ( J 2 ) = ( R 0 − 1 )

It is clear from (21) that trace ( J 2 ) < 0 and determinant ( J 2 ) > 0 if and only if R 0 > 1 . Thus, matrix J 2 admits negative eigenvalues if R 0 > 1 . Consequently the measles persistent equilibrium state is locally asymptotically stable whenever it exists ( R 0 > 1 ).

2.4.2. Global Stability Analysis of the Measles-Persistent Equilibrium State

Theorem 2.If R 0 > 1 , there exists a unique measles-persistent equilibrium state( ξ * * )that is globally asymptotically stable.

Proof. We first defined a positive definite function L as:

L ( ξ * * ) = [ ( S − S * * ) − S * * ln ( S S * * ) ] + [ ( V − V * * ) − V * * ln ( V V * * ) ]     + [ ( E − E * * ) − E * * ln ( E E * * ) ] + [ ( I − I * * ) − I * * ln ( I I * * ) ]     + [ ( R − R * * ) − R * * ln ( R R * * ) ]

Taking the time derivative of L ( ξ * * ) gives:

d L ( ξ ∗ ∗ ) d t = ( S − S ∗ ∗ S ) d S d t + ( V − V ∗ ∗ V ) d V d t + ( E − E ∗ ∗ E ) d E d t + ( I − I ∗ ∗ I ) d I d t + ( R − R ∗ ∗ R ) d R d t                           = ( S − S ∗ ∗ S ) [ π + η V − ( Λ + g 0 ) S ] + ( V − V ∗ ∗ V ) ( ρ S − g 1 V )                                   + ( E − E ∗ ∗ E ) ( Λ S − g 2 E ) + ( I − I ∗ ∗ I ) ( α I − g 3 I ) + ( R − R ∗ ∗ R ) ( γ I + τ V − μ R )                           = π + η V − ( Λ + g 0 ) S − ( π + η V ) S ∗ ∗ S + ( Λ + g 0 ) S ∗ ∗ + ρ S − g 1 V − ρ S V ∗ ∗ V                                   + g 1 V ∗ ∗ + Λ S − g 2 E − Λ S E ∗ ∗ E + g 2 E ∗ ∗ + α E − g 3 I − α E I ∗ ∗ I + g 3 I ∗ ∗ + γ I                                   + τ V − μ R − ( γ I + τ V ) R ∗ ∗ R + μ R ∗ ∗                           = L + − L − where L + = π + η V + ( Λ + g 0 ) S ∗ ∗ + ρ S + g 1 V ∗ ∗ + Λ S + g 2 E ∗ ∗ + α E + g 3 I ∗ ∗ + γ I + τ V + μ R ∗ ∗ L − = ( Λ + g 0 ) S + ( π + η V ) S ∗ ∗ S + g 1 V + ρ S V ∗ ∗ V + g 2 E + λ S E ∗ ∗ E + g 3 I + α E I ∗ ∗ I + μ R               + ( γ I + τ V ) R ∗ ∗ R   where   Λ = β I N

If we now assume L + ≤ L −

then, owing to the positivity of the model parameters, it follows from (22) that d L ( ξ * * ) d t ≤ 0 with the equality holding only when there lations S = S * * , V = V * * , E = E * * , I = I * * and R = R * * are satisfied Thus, following [21][25][24], the model solution set { S ( t ) , V ( t ) , E ( t ) , I ( t ) , R ( t ) } → { S * * , V * * , E * * , I * * , R * * } as t approaches infinity. This therefore suggests that { ξ * * } will be a global asymptotic stable equilibrium whenever it exists ( R 0 > 1 ).□

 3. Local Sensitivity Analysis

In order to help identify the model parameters with high impact on measles dynamics and if possibly target them during any intervention aimed at curtailing measles outbreaks, we carry out sensitivity analysis on R 0 parameters. Using the normalized forward sensitivity index relation:

Γ x R 0 = ∂ R 0 ∂ x × x R 0

Using the formular in (23), we generate the indices of ( R 0 ) parameters as indicated in Table 3 below:

Table 3. The values of the elasticity/sensitivity indices.

Parameter Elasticity/Sensitivity Index
α +0.0566
β +1.0000
η +0.5181
μ +0.8709
δ −0.1689
γ −0.8108
ρ −0.9629
τ −0.5030
4. Evaluating the Impact of Vaccination on Measles Transmission

To evaluate the significance of the measles first dose vaccination campaign, we compare the reproductive number of the disease for the model with and without vaccination. Now, system (1) in the absence of vaccination becomes:

{ d S d t = π − ( β I N + μ ) S d E d t = β I S N − g 2 E d I d t = α E − g 3 I d R d t = γ I − μ R The Measles-Free Equilibrium and <inline-formula> <mml:math> <mml:mrow> <mml:msub> <mml:mi>R</mml:mi> <mml:mn>0</mml:mn> </mml:msub> </mml:mrow> </mml:math> </inline-formula> of the Model without Vaccination

The disease-free equilibrium of system (24) denoted by ξ 0 = ( S 0 , E 0 , I 0 , R 0 ) = ( π μ , 0 , 0 , 0 ) . Substituting S * into (6) by S 0 = π μ and computing the reproductive ratio of the model (24), we obtain the uncontrolled reproductive number denoted by:

R 0 w v = α β g 2 g 3 = α β ( α + μ ) ( γ + μ + δ )

Now, comparing the expression of R 0 from (8) to that of R 0 w v from (25), we see that:

R 0 = α β μ g 1 g 2 g 3 ( g 0 g 1 − η ρ ) = μ g 1 ( g 0 g 1 − η ρ ) α β g 2 g 3 = μ g 1 ( g 0 g 1 − η ρ ) R 0 w v = μ g 1 μ g 1 + ρ ( τ + μ ) R 0 w v = μ ( η + τ + μ ) μ ( η + τ + μ ) + ρ ( τ + μ ) R 0 w v

Since μ ( η + τ + μ ) μ ( η + τ + μ ) + ρ ( τ + μ ) < 1 , it means that R 0 is a decreasing function of R 0 w v . Hence, an effective first dose vaccination campaign for measles reduces the spread of the disease by μ ( η + τ + μ ) μ ( η + τ + μ ) + ρ ( τ + μ ) [24].

 5. Numerical Simulations Results

To investigate the dynamical evolution of the model population sub-classes, we performed some numerical simulations. To achieve this, the following initial population sizes are considered: S ( 0 ) = 1000 , V ( 0 ) = 250 , E ( 0 ) = 500 , I ( 0 ) = 100 and R ( 0 ) = 0 with the parameter values provided in Table 2. Some of these parameters are estimated from the Bono region of Ghana measles surveillance data. The graphical results are shown from Figures 2-11. Figures 2-6 show the behavior of the model sub-classes at different rates of measles vaccination. These results suggest that high vaccination rate is needed to keep measles susceptibility rate down. Also, Figure 7Figures 7-11 give the effect of varying contact rate on the disease dynamics in the community. It is clear from these graphs that high contact rates

Figure 2. Effect of perturbing ρ on susceptible sub-population.

Figure 3. Effect of perturbing ρ on vaccinated sub-population.

Figure 4. Effect of perturbing ρ on exposed sub-population.

Figure 5. Effect of perturbing ρ on infected sub-population.

Figure 6. Effect of perturbing ρ on recovered sub-population.

Figure 7. Effect of perturbing β on susceptible sub-population.

Figure 8. Effect of perturbing ρ on vaccinated sub-population.

Figure 9. Effect of perturbing β on exposed sub-population.

Figure 10. Effect of perturbing ρ on infected sub-population.

Figure 11. Effect of perturbing β on recovered sub-population.

keep the number of exposed and infected individuals high.

 6. Conclusions and Future Research

In this work, a dynamical model for measles control using immunization in a community is formulated and analysed. Concerning the model analysis, we first verified the fundamental qualitative properties of an epidemiological model. We computed the measles infection-free for the model and determined its local and global stability conditions with reference to the magnitude of the model reproductive ratio R 0 . We further computed a unique measles-persistent equilibrium in terms of R 0 , and proved that this endemic equilibrium state always remained locally and globally stable whenever it exists ( R 0 > 1 ). Additionally, we ascertained the significance of using effective vaccination as a control measure for measles by deriving the expression of a reduction factor in the disease spread resulting from the implementation of measles immunization program in a community. The analysis of the sensitivity indices suggested that the contact rate between a susceptible and an infected person followed by measles vaccination rate is most influential in the transmission and control of measles respectively. The local and global asymptotic stability of both the measles-free and persistent equilibria states indicates the feasibility of eradicating measles in the society. Hence, there is a need for combined support from both governments/organizations and the public to strengthen the existing measles surveillance systems and immunization campaign for a measles-free world to be achieved.

For future research, any national data could be used to fit the model. Furthermore, the model could be extended to include the measles vaccine booster, maternally protected babies, and measles treatment.

 Author Contribution(s)

MK: Conceptualization, Data Curation, Formal Analysis, Methodology, Resources, Software, Validation; RIMG: Conceptualization, Data Curation, Formal Analysis, Methodology, Resources, Software, Visualization, Writing—Original Draft Preparation, Writing—Review & Editing; AI: Conceptualization, Supervision, Writing—Original Draft Preparation, Writing—Review & Editing

 Acknowledgments

Authors thankfully acknowledge the contributions from other faculty members and reviewers in shaping this paper.

 Data

Data used for our analysis is obtained from measles published works and have been duly referenced.

 References Nigusie Mitku, S. and Koya, P.R. (2017) Mathematical Modeling and Simulation Study for the Control and Transmission Dynamics of Measles. American Journal of Applied Mathematics, 5, 99-107. https://doi.org/10.11648/j.ajam.20170504.11 10.11648/j.ajam.20170504.11https://doi.org/10.11648/j.ajam.20170504.11 Mitku, S. Koya, P.R. 2017 Mathematical Modeling and Simulation Study for the Control and Transmission Dynamics of Measles American Journal of Applied Mathematics 5 10.11648/j.ajam.20170504.11 Peter, O.J., Panigoro, H.S., Ibrahim, M.A., Otunuga, O.M., Ayoola, T.A. and Oladapo, A.O. (2023) Analysis and Dynamics of Measles with Control Strategies: A Mathematical Modeling Approach. International Journal of Dynamics and Control, 11, 2538-2552. https://doi.org/10.1007/s40435-022-01105-1 10.1007/s40435-022-01105-1https://doi.org/10.1007/s40435-022-01105-1 Peter, O.J. Panigoro, H.S. Ibrahim, M.A. Otunuga, O.M. Ayoola, T.A. Oladapo, A.O. 2023 Analysis and Dynamics of Measles with Control Strategies: A Mathematical Modeling Approach International Journal of Dynamics and Control 11 10.1007/s40435-022-01105-1 World Health Organization (2022) Measles Outbreak Guide. World Health Organization. 2022 Measles Outbreak Guide Opoku, S., Seidu, B. and Akuka, P.N.A. (2023) A Mathematical Analysis of the Impact of Maternally Derived Immunity and Double-Dose Vaccination on the Spread and Control of Measles. Computational and Mathematical Biophysics, 11, Article ID: 20230106. https://doi.org/10.1515/cmb-2023-0106 10.1515/cmb-2023-0106https://doi.org/10.1515/cmb-2023-0106 Opoku, S. Seidu, B. Akuka, P.N.A. 2023 A Mathematical Analysis of the Impact of Maternally Derived Immunity and Double-Dose Vaccination on the Spread and Control of Measles Computational and Mathematical Biophysics 11 202301 ID 10.1515/cmb-2023-0106 Garba, S.M., Safi, M.A. and Usaini, S. (2017) Mathematical Model for Assessing the Impact of Vaccination and Treatment on Measles Transmission Dynamics. Mathematical Methods in the Applied Sciences, 40, 6371-6388. https://doi.org/10.1002/mma.4462 10.1002/mma.4462https://doi.org/10.1002/mma.4462 Garba, S.M. Safi, M.A. Usaini, S. 2017 Mathematical Model for Assessing the Impact of Vaccination and Treatment on Measles Transmission Dynamics Mathematical Methods in the Applied Sciences 40 10.1002/mma.4462 James Peter, O., Ojo, M.M., Viriyapong, R. and Abiodun Oguntolu, F. (2022) Mathematical Model of Measles Transmission Dynamics Using Real Data from Nigeria. Journal of Difference Equations and Applications, 28, 753-770. https://doi.org/10.1080/10236198.2022.2079411 10.1080/10236198.2022.2079411https://doi.org/10.1080/10236198.2022.2079411 Peter, O. Ojo, M.M. Viriyapong, R. Oguntolu, F. 2022 Mathematical Model of Measles Transmission Dynamics Using Real Data from Nigeria Journal of Difference Equations and Applications 28 10.1080/10236198.2022.2079411 Arsal, S.R., Aldila, D. and Handari, B.D. (2020) Short Review of Mathematical Model of Measles. AIP Conference Proceedings, 2264, Article 020003. https://doi.org/10.1063/5.0023439 10.1063/5.0023439https://doi.org/10.1063/5.0023439 Arsal, S.R. Aldila, D. Handari, B.D. 2020 Short Review of Mathematical Model of Measles AIP Conference Proceedings 2264 020003 10.1063/5.0023439 Edward, S., Raymond, K.E., Gabriel, K.T., Nestory, F., Godfrey, M.G. and Arbogast, M.P. (2015) A Mathematical Model for Control and Elimination of the Transmission Dynamics of Measles. Applied and Computational Mathematics, 4, 396-408. https://doi.org/10.11648/j.acm.20150406.12 10.11648/j.acm.20150406.12https://doi.org/10.11648/j.acm.20150406.12 Edward, S. Raymond, K.E. Gabriel, K.T. Nestory, F. Godfrey, M.G. Arbogast, M.P. 2015 A Mathematical Model for Control and Elimination of the Transmission Dynamics of Measles Applied and Computational Mathematics 4 10.11648/j.acm.20150406.12 Jaharuddin and Bakhtiar, T. (2020) Control Policy Mix in Measles Transmission Dynamics Using Vaccination, Therapy, and Treatment. International Journal of Mathematics and Mathematical Sciences, 2020, Article ID: 1561569. https://doi.org/10.1155/2020/1561569 10.1155/2020/1561569https://doi.org/10.1155/2020/1561569 Bakhtiar, T. Vaccination, T 2020 Control Policy Mix in Measles Transmission Dynamics Using Vaccination, Therapy, and Treatment International Journal of Mathematics and Mathematical Sciences 2020 156156 ID 10.1155/2020/1561569 Wang, R., Jing, W., Liu, M. and Liu, J. (2022) Trends of the Global, Regional, and National Incidence of Measles, Vaccine Coverage, and Risk Factors in 204 Countries from 1990 to 2019. Frontiers in Medicine, 8, Article ID: 798031. https://doi.org/10.3389/fmed.2021.798031 10.3389/fmed.2021.79803135127753https://doi.org/10.3389/fmed.2021.798031 Wang, R. Jing, W. Liu, M. Liu, J. Global, R Measles, V 2022 Trends of the Global, Regional, and National Incidence of Measles, Vaccine Coverage, and Risk Factors in 204 Countries from 1990 to 2019 Frontiers in Medicine 8 798031 ID 10.3389/fmed.2021.798031 35127753 Rahmayani, S.A., Aldila, D. and Handari, B.D. (2021) Cost-Effectiveness Analysis on Measles Transmission with Vaccination and Treatment Intervention. AIMS Mathematics, 6, 12491-12527. https://doi.org/10.3934/math.2021721 10.3934/math.2021721https://doi.org/10.3934/math.2021721 Rahmayani, S.A. Aldila, D. Handari, B.D. 2021 Cost-Effectiveness Analysis on Measles Transmission with Vaccination and Treatment Intervention AIMS Mathematics 6 10.3934/math.2021721 Minta, A.A., Ferrari, M., Antoni, S., Lambert, B., Sayi, T.S., Hsu, C.H., et al. (2024) Progress toward Measles Elimination—Worldwide, 2000-2023. MMWR. Morbidity and Mortality Weekly Report, 73, 1036-1042. https://doi.org/10.15585/mmwr.mm7345a4 10.15585/mmwr.mm7345a439541251https://doi.org/10.15585/mmwr.mm7345a4 Minta, A.A. Ferrari, M. Antoni, S. Lambert, B. Sayi, T.S. Hsu, C.H. 2024 Progress toward Measles Elimination—Worldwide, 2000-2023 MMWR. Morbidity and Mortality Weekly Report 73 10.15585/mmwr.mm7345a4 39541251 Minta, A.A., Ferrari, M., Antoni, S., Portnoy, A., Sbarra, A., Lambert, B., et al. (2023) Progress toward Measles Elimination—Worldwide, 2000-2022. MMWR. Morbidity and Mortality Weekly Report, 72, 1262-1268. https://doi.org/10.15585/mmwr.mm7246a3 10.15585/mmwr.mm7246a337971951https://doi.org/10.15585/mmwr.mm7246a3 Minta, A.A. Ferrari, M. Antoni, S. Portnoy, A. Sbarra, A. Lambert, B. 2023 Progress toward Measles Elimination—Worldwide, 2000-2022 MMWR. Morbidity and Mortality Weekly Report 72 10.15585/mmwr.mm7246a3 37971951 Zumah, F., Asem, L., Der, A.D. and Sackey, S. (2023) Evaluation of Measles Surveillance System, Bono Region, Ghana. medRxiv. Zumah, F. Asem, L. Der, A.D. Sackey, S. System, B Region, G 2023 Evaluation of Measles Surveillance System, Bono Region, Ghana Konlan, M., Abassawah Danquah, B., Okyere, E., Osman, S., Amenyo Kessie, J. and Kobina Donkoh, E. (2024) Global Stability Analysis and Modelling Onchocerciasis Transmission Dynamics with Control Measures. Infection Ecology & Epidemiology, 14, Article 2347941. https://doi.org/10.1080/20008686.2024.2347941 10.1080/20008686.2024.234794138736969https://doi.org/10.1080/20008686.2024.2347941 Konlan, M. Danquah, B. Okyere, E. Osman, S. Kessie, J. Donkoh, E. 2024 Global Stability Analysis and Modelling Onchocerciasis Transmission Dynamics with Control Measures Infection Ecology & Epidemiology 14 2347941 10.1080/20008686.2024.2347941 38736969 Konlan, M. (2024) Modeling the Inflow of Exposed and Infected Migrants on the Dynamics of Malaria. European Journal of Mathematical Analysis, 4, Article 7. https://doi.org/10.28924/ada/ma.4.7 10.28924/ada/ma.4.7https://doi.org/10.28924/ada/ma.4.7 Konlan, M. 2024 Modeling the Inflow of Exposed and Infected Migrants on the Dynamics of Malaria European Journal of Mathematical Analysis 4 7 10.28924/ada/ma.4.7 Wireko, F.A., Asamoah, J.K.K., Adu, I.K. and Ndogum, S. (2024) Non-Optimal and Optimal Fractional Control Analysis of Measles Using Real Data. Informatics in Medicine Unlocked, 49, Article 101548. https://doi.org/10.1016/j.imu.2024.101548 10.1016/j.imu.2024.101548https://doi.org/10.1016/j.imu.2024.101548 Wireko, F.A. Asamoah, J.K.K. Adu, I.K. Ndogum, S. 2024 Non-Optimal and Optimal Fractional Control Analysis of Measles Using Real Data Informatics in Medicine Unlocked 49 101548 10.1016/j.imu.2024.101548 Alemneh, H.T. and Belay, A.M. (2023) Modelling, Analysis, and Simulation of Measles Disease Transmission Dynamics. Discrete Dynamics in Nature and Society, 2023, Article ID: 9353540. https://doi.org/10.1155/2023/9353540 10.1155/2023/9353540https://doi.org/10.1155/2023/9353540 Alemneh, H.T. Belay, A.M. Modelling, A 2023 Modelling, Analysis, and Simulation of Measles Disease Transmission Dynamics Discrete Dynamics in Nature and Society 2023 935354 ID 10.1155/2023/9353540 Kuddus, M.A., Mohiuddin, M. and Rahman, A. (2021) Mathematical Analysis of a Measles Transmission Dynamics Model in Bangladesh with Double Dose Vaccination. Scientific Reports, 11, Article No. 16571. https://doi.org/10.1038/s41598-021-95913-8 10.1038/s41598-021-95913-834400667https://doi.org/10.1038/s41598-021-95913-8 Kuddus, M.A. Mohiuddin, M. Rahman, A. 2021 Mathematical Analysis of a Measles Transmission Dynamics Model in Bangladesh with Double Dose Vaccination Scientific Reports 11 No 10.1038/s41598-021-95913-8 34400667 Memon, Z., Qureshi, S. and Memon, B.R. (2020) Mathematical Analysis for a New Nonlinear Measles Epidemiological System Using Real Incidence Data from Pakistan. The European Physical Journal Plus, 135, Article No. 378. https://doi.org/10.1140/epjp/s13360-020-00392-x 10.1140/epjp/s13360-020-00392-x32435550https://doi.org/10.1140/epjp/s13360-020-00392-x Memon, Z. Qureshi, S. Memon, B.R. 2020 Mathematical Analysis for a New Nonlinear Measles Epidemiological System Using Real Incidence Data from Pakistan The European Physical Journal Plus 135 No 10.1140/epjp/s13360-020-00392-x 32435550 LaSalle, J.P. (1976) Stability Theory and Invariance Principles. In: DynamicalSystems, Elsevier, 211-222. LaSalle, J.P. Systems, E 1976 Stability Theory and Invariance Principles In: Dynamical Systems 211 Konlan, M. and Chuaya, R.G. (2025) Stability Analysis of a Mathematical Model for Examination Malpractice Dynamics. European Journal of Mathematical Analysis, 5, 16. https://doi.org/10.28924/ada/ma.5.16 10.28924/ada/ma.5.16https://doi.org/10.28924/ada/ma.5.16 Konlan, M. Chuaya, R.G. 2025 Stability Analysis of a Mathematical Model for Examination Malpractice Dynamics European Journal of Mathematical Analysis 5 10.28924/ada/ma.5.16 Ojo, M.M. and Goufo, E.F.D. (2021) Assessing the Impact of Control Interventions and Awareness on Malaria: A Mathematical Modeling Approach. Communications in Mathematical Biology and Neuroscience, 2021, Article ID: 93. Ojo, M.M. Goufo, E.F.D. 2021 Assessing the Impact of Control Interventions and Awareness on Malaria: A Mathematical Modeling Approach Communications in Mathematical Biology and Neuroscience 2021 93 ID Konlan, M., Danquah, B.A. and Iddrisu, A. (2025) Mathematical Modelling to Assess the Impact of a Hypothetical Vaccine on the Dynamics of Human Onchocerciasis. Communications in Mathematical Biology and Neuroscience, 2025, Article ID: 117. https://doi.org/10.28919/cmbn/9539 10.28919/cmbn/9539https://doi.org/10.28919/cmbn/9539 Konlan, M. Danquah, B.A. Iddrisu, A. 2025 Mathematical Modelling to Assess the Impact of a Hypothetical Vaccine on the Dynamics of Human Onchocerciasis Communications in Mathematical Biology and Neuroscience 2025 117 ID 10.28919/cmbn/9539 Chuaya, R.G. and Konlan, M. (2025) Analysis of a Mathematical Model of Corruption as an Epidemic. OALib, 12, 1-18. https://doi.org/10.4236/oalib.1114184 10.4236/oalib.1114184https://doi.org/10.4236/oalib.1114184 Chuaya, R.G. Konlan, M. 2025 Analysis of a Mathematical Model of Corruption as an Epidemic OALib 12 10.4236/oalib.1114184