Measles, a recurrent virus infection has a short term outbreak but its impact is devastating especially among children under five. Severe measles results in pneumonia, which is the disease with high mortality rate in Ghana. Found in , the graphs show the mortality during an outbreak of measles. In this research paper, the authors looked at a case study of the dynamics of measles disease in the case of interrupted vaccination in Sunyani, a thriving metropolis in Ghana, a country where only 57% deliveries are handled by qualified health personnel . In the Sunyani city and its immediate environs, new deliveries are commonly done in the hospital, clinics and maternity homes but significantly some cases of new births are successfully done without the presence of qualified trained health practitioner, usually in their homes. The health program in Ghana makes it necessary for regular attendant mothers with newborns to vaccinate their children against the six childhood killer diseases including measles in any health facility. The first dose of MMR (mumps, measles, rubella) vaccine is given at 12 - 13 months of age. The second dose is usually taken when the child is 3 years and 4 months old. This can be taken in later years in the management of a disease outbreak. The irregular attendants do not successfully vaccinate their children against measles. New deliveries done outside a health facility make it difficult to track these children for vaccination. There are some children who missed the first and second doses, and there are some still who also missed the second dose because of irregularity in attending a health facility and for immigration, too . There is also among the population a strong believe that MMR vaccines are leaving children with autism. Thereby, there are many who do not have a strong immunity against measles and many have their immunity wane in time. Their development of a strong immunity is interrupted so several children are susceptible to the disease and several still may have secondary infection. We seek to investigate the nature of a possible outbreak of measles in a population with dynamics such as this. Many researchers have studied the global dynamics of infectious diseases in general (see  ) while others also studied the dynamics of specific infectious diseases such as measles, chickenpox, mumps and rubella (see  ). In this paper, we also study the dynamics of measles with a relapse in a varying population where birth, death and immigration dynamics are considered.    studied the global stability of disease models and inculcated immigration and/or births and death dynamics into the population system. Again much work has already been done (see  ) defining what measles is and the symptoms and the effect of chronic infection of measles. It is mentioned that measles is severe among children under five (5) years.
Mathematical models are strong instruments used extensively to study the spread and control of infectious diseases. One important measure that determines the dynamics of disease models is the threshold parameter known as basic reproductive number . This parameter measures the number of infectives generated by a single infectious individual introduced into the susceptible population. In measuring the foregoing, researchers in recent literature use the NGM approach. When (usually when the disease-free equilibrium exists) the introduction of infectious individual can not generate a large outbreak. On the other hand, when , (where the endemic equilibrium exists), the disease will persist in the population.   have mentioned in their works that establishing global properties of dynamical systems using Lyapunov functions is not a trivial problem because there is no systematic procedure in constructing Lyapunov functions for infectious diseases. But one of the most successful procedures used to construct Lyapunov functions by many researchers is the quadratic function of the form (see ) or the nonlinear Goh-Volterra function (see  ).
In this paper, the matrix theoretic method, a type of Lyapunov function, is used to establish the stability in the disease free system. Later in this study, the Goh-Volterra method is used to establish the stability of the endemic equilibrium.   used this matrix theoretic method to study the global dynamics of several specific disease models.
2. The Model
We consider a susceptible-exposed-infectious-recovered-susceptible (SEIRS) model with relapse since there is a vaccination parameter that provides a near-permanent immunity against the disease with non-negative initial conditions. E represents the number of latent(exposed) individuals, represents the rate at which the exposed individuals become infectious (i.e.) represents the average latent period. represents the rate at which the susceptible individual become immune and gain near-permanent immunity and also represents the rate at which the infectious individuals go through the successful treatment process of the disease. Near-permanent immunity because there exist a rate of the recovered individuals falling back into the susceptible compartment. It is considered that this individual have a waning immunity over time so they fall back susceptible to the infectious disease. Already,  studied the “Global Stability Results in an SVIR Epidemic Model with Immunity Loss Rate Depending on the Vaccine-Age” where it is assumed that the waning rate of vaccine-induced immunity is depending on the vaccine-age. represents a vital dynamics-birth and immigration rates. The total number of new infections at a time t is given by . It is also considered that the natural immunity of some exposed individuals moves them into the recovered population. While is the natural death rate for the susceptible and recovered compartments, and represent the death rates for the exposed and infected classes respectively. is a growing population since births and immigration rates are added at any time t. The schematic diagram as shown in Figure 1 is drawn or developed to depict the interactions in the thriving metropolis, Sunyani of Ghana. The schematic diagram above in Figure 1 has the following as its system of equations;
3. Model Transformation & Analysis
The total population size can be determined by
Figure 1. Schematic diagram of the SEIRS model.
or from the differential equation , which is derived by adding the equations in (1). Let , , and denote the fractions of the classes S, E, I and R in the population, respectively. Then for , , and respectively, the following assignments are appropriate;
From the transformation above, the following equations are derived
subject to the restriction that . Let us hereon refer to (3) as the full model. It is also worthy to note that the total population does not appear in Equation (2). We can attribute that to the homogeneity in the system (1). Also, since r appears in the first equation of the Equation (3), we substitute into the first equation to get:
so that we can focus our attention on the subsystem (4).
The transformation has revealed some dynamics and interactions embedded in our model which is not intended to be measured in this study. Let us rewrite system (5) by making these substitutions; ; , ; and to obtain the system of equations:
For biological considerations we study the system (6) in the region
It can be shown that is positively invariant. In the latter part of this study, we show that the transformed model (3) is similar to the new model (6) by numerical simulation.  also did this transformation in their study and used the full model, but for the purpose of this study, with the established equality of the models (3) and model (6) by numerical computation we turn our attention to the latter model.
4. Equilibria & Stability
It can be seen that Equation (6) has a disease-free equilibrium is . The Jacobian matrix evaluated at the DFE is given thus,
The characteristic polynomial of the DFE system is given as,
From the characteristic polynomial above,
By the Routh-Hurwitz’s criterion, for the characteristic equation to have negative eigenvalues, then , and . The Routh-Hurwitz’s criterion is a tool used to establish the negativity or otherwise of the roots of a characteristic equation. Let us first introduce our threshold parameter .
4.1. Basic Reproductive Number, R0
Following the theory made explicit in  on above, we use the features of the model under study to carve out an appropriate measure for the . Let , F, V and be defined as in . From the system (6),
From (11), we proceed to find F and V as defined in  above.
It can be easily shown that
Another construction of the measure mostly considered by researchers long before the introduction of the new generation matrix (NGM) method was introduced is
where is the endemic equilibrium.
We remark that, the NGM’s measure of (as shown in 14) is equivalent to this measure (15) if and only if . But for the purpose of this study, we choose (14) as our measure for .
Theorem 1. Denote . When , then the system (6) has only a on the set . When then besides the DFE, there exist a unique endemic equilibrium EE where the disease is persistent.
4.2. Routh-Hurwitz’s Criterion
Theorem 2 (Routh-Hurwitz’s Criterion) Consider the characteristic equation
determining the n eigenvalues of a real square matrix A, where I is the identity matrix. Then the eigenvalues all have negative real parts if , , , , where
The above theorem (2) is reproduced from  and used in the article .
Theorem 3. The disease-free equilibrium DFE is locally stable as and unstable as .
Proof. It is easy to show from (7) that and whenever . Now for ,
And since ,
therefore . The DF system has negative eigenvalues, hence the DF system is locally stable when . Again, when , then the condition 2 of the Routh-Hurwitz’s criterion is violated; thus may become less than 0 and therefore the system becomes unstable.
4.3. Global Stability of the DFE
The matrix theoretic method is used to prove the global stability of the DFE.
A Matrix-Theoretic Method
The matrix-theoretic method is used to prove the sharp threshold statement in theorem (5). It is a systematic method, and it is presented to guide the construction of a Lyapunov function. Taking the same path as   , let us set
Then the equation for the disease compartment can be written as
Let be the left eigenvector of the non-negative matrix corresponding to the eigenvalue . The following result provides a general method to construct a Lyapunov function for [1.1].   used this Lyapunov function involving the Perron eigenvectors to study the global dynamics of several specific disease models while  used it to consider a general case for infectious diseases. In this paper, this method reproduces to establish the global stability of our system.
Theorem 4 Let F, V and be defined as in (1.2) and (2.1) in  respectively. If , in the , , and , then the function is a Lyapunov function for the system (1.1) (again in ) on .
Proof. The proof as followed in  gives
Since , , and in the region , the last term is non-positive. If , then in and thus is a Lyapunov function for the system [1.1] as defined in .
 has proven that the Lyapunov function used to prove the global stability of the DFE in can also be extended to establish a uniform persistence and thus establish the existence of an EE in . Find the theorem 2.2 in  and the proof thereof.
Theorem 5. The DFE is globally stable.
Proof. Taking the left side expansion of (19), the Lyapunov function as defined in theorem (4) therefore is
it is easy to prove that
It is obvious from this equation that when then and therefore the DFE is globally asymptotically stable and unstable when .
4.4. Existence and Stability of the Endemic Equilibrium
Referring to theorem (1), there exists a unique EE. In this section, we establish the stability of the EE using Routh-Hurwitz’s criterion for the proof of local stability and Lyapunov functions for the global stability.
4.4.1. Local Stability of Endemic Equilibrium
Theorem 6. The EE is locally stable.
Proof. The Jacobian of the EE is given thus
Then the characteristic polynomial of the EE would be
which is of the form where , and are simplified from Equation (23) as
when . It can be seen that , and are all positive. For ,
Now and by extension . The rest of the terms in the left hand side of this inequality are non-negative and therefore obviously greater than . Therefore the EE is locally stable.
4.4.2. Global Stability Analysis of the Endemic Equilibrium
A general form of Lyapunov functions coined from the first integral form of the Goh-Volterra system which is often used in the literature of mathematical biology is used to prove the global stability of the EE. This function takes the form
where x is the variables and are carefully selected constants. This criterion has been used many times in establishing the stability or otherwise of many disease models and also present in [4, 13, 15].
Theorem 7. The EE is globally stable.
Proof. Let and ,
and the equilibrium relation for
, whenever and in the region (30)
In line (26), we define the Goh-Volterra function for the first variable s and differentiate it in the second line (27). In the following line, is substituted for its relation in (6) and evaluated at the equilibrium relation for in line (29). Consider that and , (which is a necessary condition for (26) to hold). Then the inequality in line (29) is justified. Therefore, defined above is a Lyapunov function.
Again, let and
then from the system under study, (6), and
The next Lyapunov function for the e variable is given in line (31), differentiated and evaluated in lines (32) and (33) respectively. The equilibrium relation for is introduced in lines (35) and (36). Again, the inequality is obviously justified whenever . The function defined as is successful Lyapunov function.
Lastly, let us set and ,
The Lyapunov function for the last variable understudy is introduced in line (38) and differentiated in the following line (39). The equilibrium relation for is substituted in lines (42) and (43) thereby justifying the inequality.
Therefore defined as
is a Lyapunov function for the system (6). Arbitrary constants can be chosen from and any linear combination of would be a Lyapunov function for the system. Hence the proof.
5. Numerical Simulations
In this section, we simulate our model to demonstrate the theoretic results found by this study. We only talk about susceptibility and recovery over time since they explain the rest of the compartments. Note that from the schematic diagram (figure ), individuals in the exposed compartments are recruited from the susceptible and the infectious individuals are also recruited from the exposed compartment. We wish you to see that we have shown by this numerical simulation that the first system (1) which was transformed to give the full system (3) is mathematically the same as the reduced model (6). Parameters used for this simulation are shown in Table 1. With all other parameters held constant, we vary vaccination from 0.2 to 0.8 with a step size of 0.2 to see when the disease would die out of the system.
From Figure 2, we can see that the higher the vaccination cover in the population, the lesser the susceptibility to the measles disease in the population. Some significant fraction of the population remains susceptible to measles at a 20% yearly vaccination cover after ten-year span. This means then that in the wake of an outbreak there would still be some that would be impacted by the disease. At a 40% vaccination cover yearly, the susceptible individuals reduce sharply. The same can be said of a vaccination cover of 60% and 80% vaccination cover. It can be drawn from here that when vaccination cover is aimed at 80%, the number of susceptible individuals approaches zero.
In Figure 3, it can be seen that awareness of the disease is already in the population and control steps for treatment and vaccination do not allow for an outbreak to get worse before it becomes better. It can be seen that no matter the
Table 1. Parameter and variable definition & values.
Figure 2. Plot of the susceptible individuals with time.
Figure 3. Plot of the exposed individuals with time.
vaccination cover in the population, the disease will persist in the population for a while and die out. When the time is extended, it would be seen that the disease will completely die out of the system as the exposed individuals approaches zero. It makes sense because, from the schematic diagram (4.1) of the model, the exposed individuals are assumed to go through the infectious stage or attained a natural cure against the disease and move into the recovery stage. Again, we see that the full model and the reduced model are the same.
In Figure 4 also, with time, measles would no longer be seen in the population. At any level of vaccination cover, the disease will die out of the population. Only that, at a higher vaccination rate, more susceptible individuals do not have to get the disease but will gain a strong immunity against the disease.
In Figure 5, the higher the vaccination rate, the more susceptible individuals gain immunity against the disease so that more people fall into the safe haven of being recovered. But it can be seen at 60% vaccination cover is comparable to an 80% vaccination cover. It would therefore be a good practice if vaccination in a population whose system can be compared to that which is studied in this paper vaccinate 80% of its people, newborns and new immigrants so that in any case of an outbreak of measles, the disease will be eradicated in the shortest possible time.
Figure 4. Plot of the infectious individuals with time.
Figure 5. Plot of the recovered individuals with time.
In this paper, the global dynamics of measles has been studied in a varying population system such as could be compared to the thriving metropolis of Sunyani in Ghana. A matrix theoretic method and Goh-Volterra types of Lyapunov function were used to establish the stability of the disease-free and endemic equilibria respectively. The basic reproductive number is defined using the New Generation Matrix method and proved as the threshold parameter. The DFE is globally stable whenever and unstable otherwise. The EE is globally stable when . Also, the transformation of the original system (1) uncovered some unnecessary measures embedded in the model which was removed to give the approximate system (6). The latter system was used to replace the former and towards the end of this paper, we showed by numerical simulation that approximate system (3) is the same as (6). The numerical simulation also showed the theoretic dynamics discussed in this paper. It showed that more people will gain immunity against the measles and by extension, an outbreak of measles would not impact the community if vaccination cover is 80% annually.