Abstract

We derive and analyze a mathematical model of smoking in which the population is divided into four classes: potential smokers, smokers, temporary quitters, and permanent quitters. In this model we study the effect of smokers on temporary quitters. Two equilibria of the model are found: one of them is the smoking-free equilibrium and the other corresponds to the presence of smoking. We examine the local and global stability of both equilibria and we support our results by using numerical simulations.

1. Introduction

Tobacco epidemic is one of the biggest public health threats the world has ever faced. It kills up to half of its users. Nearly, each year six million people die from smoking of whom more than 5 million are users and ex-users and more than 600,000 are nonsmokers exposed to second-hand smoke. Tobacco users who die prematurely deprive their families of income, raise the cost of health care, and hinder economic development. In smoking death statistics, there is a death caused by tobacco every 8 seconds; 10% of the adult population who smoke die of tobacco related diseases. The World Health Organization predicts that, by 2030, 10 million people will die every year due to tobacco related illnesses. This makes smoking the biggest killer globally. In Saudi Arabia, the prevalence of current smoking ranges from 2.4 to 52.3% (median = 17.5%) depending on the age group. The results of a Saudi modern study predicted an increase of smokers number in the country to 10 million smokers by 2020. The current number of smokers in Saudi Arabia is approximately 6 million, and they spend around 21 billion Saudi Riyal on smoking annually. Clearly smoking is a prevalent problem among Saudis that requires intervention for eradication. Persistent education of the health hazards related to smoking is recommended particularly at early ages in order to prevent initiation of smoking [1, 2].

In 1997, Castillo-Garsow et al. [3] proposed a simple mathematical model for giving up smoking. They considered a system described by the simplified model with a total constant population which is divided into three classes: potential smokers, that is, people who do not smoke yet but might become smokers in the future (), smokers (), and people (former smokers) who have quit smoking permanently (). Later, this mathematical model was developed by Sharomi and Gumel (2008) [4]; they introduced a new class of smokers who temporarily quit smoking and they described the dynamics of smoking by the following four nonlinear differential equations: They concluded that the smoking-free equilibrium is globally asymptotically stable whenever a certain threshold, known as the smokers generation number, is less than unity and unstable if this threshold is greater than unity. The public health implication of this result is that the number of smokers in the community will be effectively controlled (or eliminated) at steady state if the threshold is made to be less than unity. Such a control is not feasible if the threshold exceeds unity (a global stability result for the smoking-present equilibrium is provided for a special case). In 2011, Lahrouz et al. [5] proved the global stability of the unique smoking-present equilibrium state of the mathematical model developed by Sharomi and Gumel. Also in 2011, Zaman [6] derived and analyzed a smoking model taking into account the occasional smokers compartment, and later [7] he extended the model to consider the possibility of quitters becoming smokers again. In 2012, Ertürk et al. [8] introduced fractional derivatives into the model and studied it numerically.

In this paper, we adopt the model developed and studied in [4, 5] and examine the effect of peer pressure on temporary quitters . By this we mean the effect of smokers on temporary quitters which may be considered as one of the main causes of their relapse. This is done by taking as a parameter dependent on the pressure of smokers. The aim of this work is to analyze the model by using stability theory of nonlinear differential equations and supporting the results with numerical simulation. The structure of this paper is organized as follows. In Section 2, we calculate the smoking generation number, find the equilibria, and investigate their stability. In Section 3, we present numerical simulations to support our results. Finally, we end the paper with a conclusion in Section 4.

2. Formulation of the Model

Let the total population size at time be denoted by . We divide the population into four subclasses: potential smokers (nonsmoker) , smokers , smokers who temporarily quit smoking , and smokers who permanently quit smoking , such that . We describe the dynamics of smoking by the following four nonlinear differential equations: Here, is the contact rate between potential smokers and smokers, is the rate of natural death, is the contact rate between smokers and temporary quitters who revert back to smoking, is the rate of quitting smoking, is the fraction of smokers who temporarily quit smoking (at a rate ), and is the remaining fraction of smokers who permanently quit smoking (at a rate ).

We assume that the class of potential smokers is increased by the recruitment of individuals at a rate . In model 2, the total population is constant and , , , and are, respectively, the proportions of potential smokers, smokers, temporarily quitters, and permanent quitters where . Since the variable of model 2 does not appear in the first three equations, we will only consider the subsystem: From model 3 we find that Let ; then . The initial value problem , with , has the solution , and . Therefore, , which implies that . Thus, the considered region for model 3 is All solutions of model 3 are bounded and enter the region . Hence, is positively invariant. That is, every solution of model 3, with initial conditions in , remains there for all .

3. Equilibria and Smokers Generation Number

Setting the right-hand side of the equations in model 3 to zero, we find two equilibria of the model: the smoking-free equilibrium and smoking-present equilibrium , where and satisfies the equation The smokers generation number is found by the method of next generation matrix [9, 10]. Let ; then model 3 can be rewritten as such that By calculating the Jacobian matrices at , we find that where Thus, the next generation matrix is and the smokers generation number is the spectral radius .

Next we explore the existence of positive smoking-present equilibrium of model 3. Equation (7) has solutions of the form

where and are If , then and . So we have only one positive solution: If , then and there are no positive solutions. If and , then and and there are no positive solutions. These results are summarized below.

Theorem 1. Model 3 always has the smoking-free equilibrium point . As for the existence of a smoking-present equilibrium point , one has three cases:(i)if and , one has no positive equilibrium point;(ii)if , one has one positive equilibrium point ;(iii)if , one has no positive equilibrium point.

4. Stability of the Equilibria

4.1. Local Stability

First, we investigate the local stability of and state the following theorem.

Theorem 2 (local stability of ). If , the smoking-free equilibrium point of model 3 is locally asymptotically stable. If , is locally stable. If , is unstable.

Proof. Evaluating the Jacobian matrix of model 3 at (using linearization method [11]) gives Thus, the eigenvalues of are , . Clearly is negative if . So all eigenvalues are negative if , and hence is locally asymptotically stable. If , then and is locally stable. If , then which means that there exists a positive eigenvalue. So is unstable.

Next, we investigate the local stability of the positive equilibrium by using the following lemma.

Lemma 3 (see [1214]). Let be a real matrix. If , , and are all negative, then all of the eigenvalues of have negative real part.

in the previous lemma is known as the second additive compound matrix with the following definition.

Definition 4 (see [15] (second additive compound matrix)). Let be an real matrix. The second additive compound of is the matrix defined as follows:

Theorem 5 (local stability of ). The smoking-present equilibrium of model 3 is locally asymptotically stable if .

Proof. Linearizing model 3 at the equilibrium gives The second additive compound matrix is Here Hence, by previous lemma, is locally asymptotically stable if .

4.2. Global Stability

We will examine the global stability of when . First, it should be noted that in for all . Consider the following Lyapunov function [16]: We have   for and only if and . Therefore, by LaSalle's Invariance Principle [16], every solution to the equations of the model with initial conditions in approaches as . We have the following result.

Theorem 6 (global stability of ). If , then is globally asymptotically stable in .

Now we examine the global stability of by using the following theorem [12, 17] to prove that model 3 has no periodic solutions, homoclinic loops, and oriented phase polygons inside the invariant region.

Theorem 7. Let be a vector field which is piecewise smooth on and which satisfies the conditions , in the interior of , where is the normal vector to and is a Lipschitz continuous field in the interior of and . Then the differential equation system , , and has no periodic solutions, homoclinic loops, and oriented phase polygons inside .

Let . One can easily prove that , is positively invariant, and . Thus, we state the following theorem.

Theorem 8. Model 3 has no periodic solutions, homoclinic loops, and oriented phase polygons inside the invariant region .

Proof. Let , , and denote the right-hand side of equations in model 3, respectively. Using to rewrite , , and in the equivalent forms, we have Let be a vector field such that Since the alternate forms of , , and are equivalent in , then So on .
Using the normal vector to , we can see that Hence, model 3 has no periodic solutions, homoclinic loops, and oriented phase polygons inside the invariant region .

Consequently, we have the following result.

Theorem 9 (global stability of ). If , then the smoking-present equilibrium point of system (3) is globally asymptotically stable.

Proof. We know that, if in , then is unstable. Also is a positively invariant subset of and the limit set of each solution of (3) is a single point in since there is no periodic solutions, homoclinic loops, and oriented phase polygons inside . Therefore is globally asymptotically stable.

5. Numerical Simulations

In this section, we illustrate some numerical solutions of model 2 for different values of the parameters, and we show that these solutions are in agreement with the qualitative behavior of the solutions. We use the following parameters: , , , and , with two different values of the contact rate between potential smokers and smokers: for we use and for we use . Model 2 is simulated for the following different initial values such that :(1), , , and ;(2), , , and ;(3), , , and ;(4), , , and .

For , Figure 1(a) shows that the number of potential smokers increases and approaches the total population 1. Figure 1(b) shows that the number of the smokers decreases and approaches zero. In Figures 1(c) and 1(d), temporary quitters and permanent quitters increase at first; after that they decrease and approach zero. We see from Figure 1 that, for any initial value, the solution curves tend to the equilibrium , when . Hence, model 2 is locally asymptotically stable about for the previous set of parameters.

In Figure 2, we use the same parameters and initial values as previously with . Figure 2(a) shows that the number of potential smokers decreases at first; then it increases and approaches . Figure 2(b) shows that the number of smokers increases at first; then it decreases and approaches . In Figures 2(c) and 2(d), temporary quitters and permanent quitters increase at first; after that they decrease and approach and . We see from Figure 2 that, for any initial value, the solution curves tend to the equilibrium , when . Hence, model 2 is locally asymptotically stable about for the above set of parameters.

6. Discussion and Conclusions

In this paper, we presented a mathematical model to analyze the behavior of smoking dynamics in a population with peer pressure effect on temporary quitters . Local asymptotic stability for the smoking-free equilibrium state is obtained when the threshold quantity is less than 1 (i.e., when the contact rate between potential smokers and smokers is less than the sum of natural death rate and the quitting rate ). A Lyapunov function is used to show global stability of the smoking-free equilibrium when the contact rate between potential smokers and smokers is less than or equal to the natural death rate (). This means that the number of smokers may be controlled by reducing the contact rate to be less than the natural death rate . On the other hand, if (i.e., when the contact rate between potential smokers and smokers is greater than the contact rate between smokers and temporary quitters who revert back to smoking), then the smoking-present equilibrium state is locally asymptotically stable. By showing that this model has no periodic solutions, homoclinic loops, and oriented phase polygons inside the invariant region , we proved the global asymptotic stability of . This means that, if , smoking will persist. Some numerical simulations are performed to illustrate the findings of the analytical results. For , the results show that the smokers population reaches a steady state of approximately 6% of the total population.

This model may be extended to include, for example, smokers who after quitting smoking may become potential smokers again. We may also consider the public health impact of smoking related illnesses as well.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.