Sensitivity and mathematical model analysis on secondhand smoking tobacco

In this paper, we are concerned with a mathematical model of secondhand smoker. The model is biologically meaningful and mathematically well posed. The reproductive number R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{0}$$\end{document} is determined from the model, and it measures the average number of secondary cases generated by a single primary case in a fully susceptible population. If R0<1,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{0}<1,$$\end{document} the smoking-free equilibrium point is stable, and if R0>1,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{0}>1,$$\end{document} endemic equilibrium point is unstable. We also provide numerical simulation to show stability of equilibrium points. In addition, sensitivity analysis of parameters involving in the dynamic system of the proposed model has been included. The parameters involving in reproductive number measure the relative change in R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{0}$$\end{document} when the value of the parameter changes.

smoke. The highest proportions exposed were estimated in Europe, the Western Pacific, and Southeast Asia, with more than 50% of population exposed. Proportion of people exposed was lowest in Africa [6]. The prevalence of secondhand smoke exposure among adolescents in Ethiopia is highest [11].
Different scholars discussed about secondhand smoke in relation to health problem. For instance, [12] discussed the global health burden of secondhand smoke, and [13] emphasized a new and alarming consequence of smoking in indoor environments.
All studies mentioned above showed secondhand smoke is a big problem in the world as the result of leading peoples to diseases and death. There are a number of studies on mathematical model of smoking by taking different assumption. We direct refer the reader to the papers [14][15][16][17].
In this article, we will see the effect of secondhand smoke in the society. We will prepare a mathematical model and interpret the model graphically. The model, we will develop extends the work of Eduardo L. L. M and Maribel L. B in paper [18]. We consider death (νP) of secondhand smoker by this infection in addition to natural death [19,20], and we discuss sensitivity analysis of R 0 . Different parameters will be used for the infection rate PS from P to S and infection rate δQS from Q to S. The infection rate δQS from Q to S means the infection rate when relapse smoking, that is, restarting smoking after a period of abstinence. In [21], relapse rates within the first year of abstinence ranged from 60 to 90%, while 2 years of continuous cessation indicated a likelihood of 80% to maintain long-term abstinence. A research conducted in the USA has the following results [22]. 1 The prevalence of relapse is 6.8%. 2 Prevalence and odds of relapse were higher among young people compared with elders. 3 Former smokers living in smoke-free homes had 60% lower odds of relapse compared with those living in homes that allowed smoking inside. 4 Odds of relapse were higher among never married, widowed, divorced and separated individuals, compared with the married group. 5 Continuous smoking cessation for 6 months or more significant decreased odds of relapse.
We organized this paper as follows. In Section 2, we prepare the proposed model. In Sections 3, 4 and 5, we discuss free equilibrium point, reproduction number, sensitivity analysis and endemic equilibrium point, respectively. In Section 6, we discuss stability of equilibrium solution. In Section 7, we see the numerical simulation for the equilibrium solution. Finally, conclusion, abbreviations and declaration are included.

Model assumption
We have the following assumptions in preparing the model.
1 There are deaths as a result of smoking and secondhand smoking. 2 The number of healthy and completely recovered people is not directly considered in the system. 3 People can be in one of the three groups: (a) Secondhand smokers or those at risk of others smoking, represented by P; (b) People who have stopped smoking but are at risk due to their smoking habit before, represented by Q. (c) People who are addicted to tobacco and now smoke it, represented by S. 4 There is a constant flow, α , of healthy people who become at risk of active smokers to the state P. 5 People leave from the states P or S or Q under study due to factors such as living in a completely nonsmoking population or the death of the individual.

Description of variables and parameters
In the following table, we describe the variables and parameters to create the mathematical model that represents the dynamics of transmission of the habit of smoking ( Table 1).

The dynamic system
As we see in Fig 1 P increases by α per unit time 2 P decreases due to natural death ( µP people per unit time); the influence of factors that cause a person to move from the population at risk of being an active smoker to the population of healthy people ( σ P people per unit time); death as secondhand smoker ( νP ); and the impact of smoking ( PS people per unit time). 3 S increases by the impact of ( PS and δQS people per unit time). 4 S decreases due to natural death ( µS ), death as a result of smoking tobacco ( κS ), the influence of factors that cause a person to leave the population of active smokers and join the population of people who have stopped smoking. 5 Q increases due to the influence of factors that cause a person to leave the population of active smokers and join the population of people who have stopped smoking 6 Q decreases due to natural death ( µQ ), deaths as a result of smoker some days before ( ηQ ) and the influence of factors that cause a person to move from the population to other areas and also medical treatment to the population of healthy people ( ξ Q). We now describe the dynamic system as follows. (1)

) can be expressed as an inequality
Integrating both sides from 0 to t, the solution is obtained as In the same manner from Eqs. (2) and (3), we obtain and , respectively. The latter one is justified by ζ S(t) ≥ 0 . (1), (2) and (3) are bounded.

The population size N(t) is given by
Differentiating both sides with respect to t, we obtain Substituting (1), (2) and (3) in (4), we have The solution of Eq. (5) is Therefore, the solution of the dynamic system (1), (2) and (3) is bounded. So, let the total population size at time t be denoted by N(t). So, We assume that N(t) is constant, and for the purpose of computation, we take P(t), Q(t) and S(t) as proportions of N(t) where P(t) + Q(t) + S(t) = 1.
Thus, the region that we consider for our model is In this paper, P(t), Q(t) and S(t) are proportions of N(t).

Free equilibrium point
Smoking-free equilibrium points are steady-state solutions of a mathematical model indicating that there is no smoking (only presented secondhand smoker). Thus, setting dP dt = 0, dS dt = 0 and dQ dt = 0, we obtain From (7), we get S = 0. Substituting this value in (6) and (8) and calculating P and Q, we have respectively P = α ν + µ + σ and Q = 0. Thus, the free equilibrium point of the

Reproduction number
The basic reproduction number, R 0 , of an infectious disease is the average number of secondary cases generated by a single primary case in a fully susceptible population [23]. R 0 is the most widely used epidemiological measurement of the transmission potential in a given population. Statistical estimation of R 0 has been performed for various infectious diseases [24,25], aiming toward understanding the dynamics of transmission and evolution and designing effective public health intervention strategies. In particular, R 0 has been used for determining the minimum coverage of immunization, because the threshold condition to prevent a major epidemic in a randomly-mixing. It is primarily used as a threshold parameter: If R 0 < 1, the disease will fade out of the population, but if R 0 > 1 , the disease will persist and become endemic to the population. Furthermore, the larger the magnitude of R 0 , the faster the disease will spread and presumably the more difficult it will be to control.
Even though there are different methods in which R 0 can be calculated, we use the nextgeneration method. This is the most common method of calculating R 0 . In this method, we place appropriate terms from the infected class equations into the vectors F and V. Terms that describe appearances of new infections in each compartment belong in F, and other terms belong in V. The Jacobian matrices obtained by differentiating F and V with respect to the relevant subset of variables are computed and evaluated at a nontrivial disease-free � = {(P, Q, S) ∈ R 3 : P + S + Q = 1, and P, S, Q ≥ 0}.

Sensitivity analysis
Sensitivity analysis tells us how important each parameter is to interpret the model. Sensitivity allows us to measure the relative change in a variable when a parameter changes. The sensitivity of R 0 with respect to a parameter is the ratio of the relative change in the variable ( R 0 ) to the relative change in the parameter. When the variable is a differentiable function of the parameter, the sensitivity index may be alternatively defined using partial derivatives [26]. Thus, for a parameter p, we denote sensitivity of R 0 with respect to p by S R 0 p and as in [27] defined by We calculate the sensitivity analysis of parameters involved in R 0 as follows. We can see that S R 0 and S R 0 α are positive. This indicates that and α are directly proportional to R 0 . This means an increase (or decrease) in and α will cause an increase (or decrease) in R 0 with the same proportion. So, the average number of secondary cases of infection increases (or decreases) in the community. It can also be seen that the remaining five parameters ν, µ, σ , κ and ζ are inversely proportional to R 0 . So, an increase (or decrease) in these parameters leads to minimizing (or maximizing) the endemic nature of smoking in the community. The relationship between R 0 and the above parameters can be described graphically as shown below (Fig. 2).

Endemic equilibrium point
Endemic equilibrium is steady-state solutions when the smoking persists in the population. Assume that the smoking will persist in the population. Thus, in (7), if S = 0, then From (6) and (8), we get Substituting (18) and (19) in (17), we obtain the quadratic equation (in a variable S) (16)  where R 0 is given as in (9), L = ν + µ + σ , M = µ + κ + ζ , and N = η + µ + ξ . From Eq. (20), we can calculate In this case, P and Q can be found by substituting Eq. (22) into Eqs. (18) and (19), respectively. If R 0 < 1, then −V > 0 and W < 0 , and thus either S is negative real number or complex number. In this case, S does not exist and hence P and Q. If R 0 = 1, then V = δζ L − NM and W = 0. Consequently, We write the above discussion as a theorem as follows.

Theorem 3
The endemic equilibrium point of the dynamic system (1), (2) and (3) exists if R 0 > 1 or R 0 = 1 and δζ L > NM and does not exist if R 0 < 1. For R 0 > 1 or R 0 = 1 and δζ L > NM, the endemic equilibrium point is the ordered triple (P, S, Q), where S, P and Q are given as in (22), (18) and (19), respectively.

Stability
The physical stability of the equilibrium solution of the dynamic system (6), (7) and (8) is related to its eigenvalue.

Definition 4
The equilibrium solution of the dynamic system is stable if all the eigenvalues of the Jacobian matrix of the dynamic system evaluated at the equilibrium solution have negative real part, otherwise the solution is unstable.
The Jacobian matrix of the dynamic system is

Free equilibrium point
At the free equilibrium point, the Jacobian matrix (23) becomes The eigenvalues of this matrix are We see that 1 and 3 are negative. If R 0 < 1, then 2 < 0 and hence the free equilibrium point is stable. If R 0 > 1, then 2 > 0 and hence the free equilibrium point is unstable. From this discussion, we now have the following theorem.
Theorem 5 If R 0 < 1, then the free equilibrium point is stable. If R 0 > 1, then the free equilibrium point is unstable.

Endemic equilibrium point
The Jacobian matrix (23) of the dynamic system (6), (7) and (8) at the endemic equilibrium point is where (P * , S * , Q * ) = (P, S, Q), here P, S and Q are given as in (18), (19) and (22). Its characteristic polynomial is where Here, our objective is to show all the roots of the characteristic polynomial (24) have negative real part. However, it is difficult to show this due to the complexity of a 1 , a 2 and a 3 . Later, we will illustrate the stability of the endemic equilibrium solution using numerical simulations, by imposing the value of parameters or variables. It is easy to see that if a 3 > 0 and a 1 a 2 − a 3 > 0, then all the roots of (24) have negative real part. In this case, the endemic equilibrium point is stable by Routh-Hurwitz criteria. By this criteria, if a 3 < 0 or a 1 a 2 − a 3 < 0, then the endemic equilibrium point is unstable.

Numerical simulation
The following initial values (IV) for P, S and Q have been taken at t = 0 . For numerical simulation, we use Runge-Kutta 4-5 methods and MATLAB 2018 software. It will be seen separately for free and endemic equilibrium point.
Since the free equilibrium point does not depend on , in both cases, we have We will construct the graphs in each case for t ∈ [0, 400] as follows.
(i) = 0.035 As it can be seen in Fig. 3, for R 0 < 1 we conclude that the graphs of P, S and Q for different initial values approach P = 2.1413, S = 0 and Q = 0 , respectively, as t → ∞ respectively. We recall that the point (P, S, Q) = (2.1413, 0, 0) is smoking-free equilibrium point of the model. Hence, the smoking-free equilibrium point is stable. (ii) = 0.165 As it is described in Fig. 4, for R 0 > 1 we conclude that the graphs of P, S and Q for different initial values do not approach P = 2.1413, S = 0 and Q = 0 , respectively, as t → ∞ . We again recall that the point (P, S, Q) = (2.1413, 0, 0) is smoking-free equilibrium point of the model. Thus, smoking-free equilibrium point is unstable.

Endemic equilibrium point
We take α = 0.3, ν = 0.01, µ = 0.0035, σ = 0.2, = 0.2, δ = 0.05, κ = 0.01, ζ = 0.2, η = 0.02 and ξ = 0.1. In this case, the endemic equilibrium point is We will construct the graphs for t ∈ [0, 400] as follows. In Fig. 5, we see that for R 0 > 1 we conclude that the graphs of P, S and Q for different initial values approach P = 0.8846, S = 0.7085 and Q = 0.8916 as t → ∞ , respectively. Therefore, the endemic equilibrium point is stable. We will construct the graphs for t ∈ [0, 600] as follows. In Fig. 6, we notice that for R 0 = 1 and δζ L > NM we conclude that the graphs of P, S and Q for different initial values approach P = 0.6125, S = 1.2232 and Q = 0.6121 as t → ∞ , respectively. So, here the endemic equilibrium point is stable.

Conclusion
1 The average number of new smokers from a single addicted smoker is given by Secondhand smoker persists in the population if R 0 > 1 or R 0 = 1 and δζ L > NM. Secondhand smoker will be disappeared over time if R 0 < 1. 2 An increase (or decrease) in leads an increase (or decrease) secondhand smoker in the community.
α (ν + µ + σ )(µ + κ + ζ ) . 3 An increase (or decrease) in the parameters ν, µ, σ , κ and ζ leads to minimizing (or maximizing) the endemic nature of smoking in the community. So secondhand smoker will be minimized or maximized. 4 If the average number of smokers that a single addicted smoker produces is less than one, the population of secondhand smokers disappears over time. 5 If the average number of smokers that a single addicted smoker produces is greater than one, the population of secondhand smokers persists.