Approximating a Giving Up Smoking Dynamic on Adolescent Nicotine Dependence in Fractional Order

In this work, we consider giving up smoking dynamic on adolescent nicotine dependence. First, we use the Caputo derivative to develop the model in fractional order. Then we apply two different numerical methods to compute accurate approximate solutions of this new model in fractional order and compare their results. In order to do this, we consider the generalized Euler method (GEM) and multi-step generalized differential transform method (MSGDTM). We also show the unique positive solution for this model and present numerical results graphically.


Introduction
Mathematical modeling is a tool which have a lot of applications in applied sciences that require a deep study for the different verities of methods used in applied mathematics. Smoking is the prodigious cause of many diseases specially of different type of cancers. As cigarette smoke contains over forty thousands dangerous chemicals, which cause harmful infections to human health. The life of smokers is ten to thirteen years shorter than that of non smokers and at the rat of 1:2 of smokers die from diseases launched by the cigarette smoking. Besides, smokers have 70% more chances of heart attack than non smokers. The incidence rate of lung cancer is 10% high in smokers. To secure the life expectancy every scientist, doctor and mathematician try to control smoking. Mathematicians try to make different smoking models for the best representation of cigarette smoking phenomena. So several authors proposed different smoking models, for example see the first model presented by Castillo-Garsow et al., [1] in which they studied different classes of smokers (potential smokers (P), smokers (S), and quit smokers (Q)). Then Sharami et al., [2] modified the model presented in [1] and introduced a new class named chain smokers. In their work, they presented the development and public health impact of smoking related diseases. Zeb et al., [3] presented the square-root dynamics of a giving up smoking model. In this work they discussed local and global stability of the model and its general solutions in which the interaction of occasional and potential smokers occur. Nowadays, every researcher tries to establish different epidemic models in fractional order. Perhaps, the reality of nature could be better translated by fractional calculus. Fractional calculus is used in many fields of sciences [4][5][6][7][8][9][10][11][12][13][14][15]. In this paper, we present the fractional order derivative and find analytic numeric solution of model presented in [16], and the model is as follow: under the initial conditions: where P, S 1 , S 2 , Q t and Q p all functions of t denote the numbers of potential, occasional and quit smokers. Here Λ is number of incoming population per unit time, b i is incoming rate per unit time i = 1, 2, 3, 4, 5 ðb ¼ P 5 i¼1 b i ¼ 1Þ: Average rate of incoming for a year, d i : Outgoing per unit time i = 1, 2, 3, 4, 5, μ is natural death rate or death related to other diseases not smoking and d represents death rate of all classes caused by smoking diseases.
By establish fractional order derivative into the model [16], we obtain the following fractional order model: Here we consider the Caputo sense fractional derivative and α is the order of the fractional time-derivative, subject to the initial conditions given in Eq (2). In the present paper, we take in account the model presented in [16]. First, we will use the Caputo derivative to develop the model in fractional order. Then we will apply two different numerical methods to compute accurate approximate solutions of this new model in fractional order and compare their results. In order to do this, we consider generalized Euler method (GEM) and the multi-step generalized differential transform method (MSGDTM). We also show the unique positive solution for this model and present numerical results graphically. This paper is organized as follows. In Section 2, we present positivity of the new model. The idea of GEM and MSGDTM for solution of the proposed model are presented shortly in Sections 3. Numerical simulation results graphically of GEM and MSGDTM are presented in Section 4 with comparisons with the results of Runge-Kutta method (RKM). The conclusion is given in Section 5. Appendix is devoted to present some basic definitions and results which is needed in this work are given in the supporting file.

Non-negative Solutions
Let R 5 þ ¼ fX 2 R 5 : X ! 0g and X = (P, S 1 , S 2 , Q t , Q p ) T . We need lemma [Generalized Mean Value Theorem] [17] that help us in the proof of subsequent theorem. (3), are positive for all t > 0.
Proof 1 If the above condition does not satisfied, then at least one of the individuals may be negative. Then the individuals will satisfy one of the following conditions.
1. There exists a first time t 1 which is contradiction to the above supposition.

2.
There exists a first time t 2 , which is contradiction to the above supposition.

Generalized Euler Method (GEM) and Multistep Generalized Differential Transform Method (MSGDTM)
Generally, it is impossible to determine the analytic solution of nonlinear differential equations. However, different authors used different methods for the approximate analytic and numeric solutions of nonlinear differential equations [7][8][9][10][11]. But all these methods are applicable for short time interval. So we have in mind to find the approximate solutions of problems Eqs (2) and (3) via the GEM and MSGDTM. The GEM is derived to get numerical solution of initial value problems with Caputo derivatives by Odibat et al., [11]. This method is a generalization of the classical Euler's method for detail see [11]. Also in this section, we use the MSGDTM that to find numerical solution of the system Eq (3) of fractional order differential equations and the detail analysis of this method is found in [12]. The basics steps of the GDTM can be found in [12][13][14][15]. In this paper, we are not going in detail of these methods but we use them only for the numerical solution of the system Eq (3). Applying the MSGDTM, we obtain the following system: Here, P(k), S 1 (k), S 2 (k), Q t (k) and Q p (k) are the differential transformation of P(t), S 1 (t), S 2 (t), Q t (t) and Q p (t), respectively. Also, PS 1 S 2 N(k) is the differential transformation of function PS 1 S 2 NðkÞ ¼ PðS 1 ðtÞ þ CS 2 ðtÞÞ NðtÞ and defined as follows: In next section, we will present the numerical results.

Numerical and Simulation Results
We applied the GEM and MSGDTM to solve the system Eq (1) for α = 1. In order to demonstrate the effectiveness of these two methods as an approximate tool for solving the nonlinear system of fractional differential Eq (3) for larger time t, we apply these two methods on the interval [0, 30]. It is to be noted that GEM results are obtained for h = 0.00001 while the MSGDTM results are obtained when K = 10 and M = 30000 and RKM for h = 0.001. All the results are calculated by using computer algebra package Mathematica and in this paper, we show only the graphically obtained results.
We assume the parameters of the system Eq (1) shown in Table 1. Figs 1-5 show the approximate solutions for P(t), S 1 (t), S 2 (t), Q t (t) and Q p (t) obtained for different values of α using GEM. From the graphical results given in Figs 1-5, it can be seen that the results obtained using GEM match the results of Runge-Kutta for integer case(α = 1,) very well when α = 1, which implies that GEM can predict the behavior of these variables accurately for the region under consideration.

Discussion
In this paper, we used Generalized Euler method (GEM) and multi-step generalized differential transform method (MSGDTM) as a reasonable basis for studying the dynamics of a new giving up smoking model and modified the integer order model Eq (1) into a fractional-order model Eq (3). The results obtained show that when α ! 1 the solutions of fractional model, {P} α (t), {S 1 } α (t), {S 2 } α (t), {Q t } α (t) and {Q p } α (t), reduce to the standard solutions P(t), S 1 (t), S 2 (t), Q t (t) and Q p (t). Finally, the recent appearance of fractional differential equations as models in some        Theorem 2 (Generalized Taylor's formula) [17] Suppose that D ka a f ðxÞ 2 Cða; b for k = 0, 1, . . ., n + 1, where 0 < α 1. Then we have f ðxÞ ¼ For mathematical properties of fractional derivatives and integrals one can consult the mentioned references in this paper.