Analysis and Simulation of Fractional Order Smoking Epidemic Model

In recent years, there are many new definitions that were proposed related to fractional derivatives, and with the help of these definitions, mathematical models were established to overcome the various real-life problems. The true purpose of the current work is to develop and analyze Atangana-Baleanu (AB) with Mittag-Leffler kernel and Atangana-Toufik method (ATM) of fractional derivative model for the Smoking epidemic. Qualitative analysis has been made to `verify the steady state. Stability analysis has been made using self-mapping and Banach space as well as fractional system is analyzed locally and globally by using first derivative of Lyapunov. Also derive a unique solution for fractional-order model which is a new approach for such type of biological models. A few numerical simulations are done by using the given method of fractional order to explain and support the theoretical results.


Introduction
Mathematics was firstly used in biology in the twelfth century when Fibonacci used his popular Fibonacci series to explain a growing population. Daniel Bernoulli used mathematics to describe the effect of small pox. The term biological mathematics was primarily used by Johannes Reinke in 1901. It is aimed at the mathematical image and modeling of biological processes. It is also used to recognize phenomena in the living organism. Bio math has made major progress during the last few decades, and this progress will continue in upcoming decades. Math has played a large role in natural science but now it also will be more useful in biology. We should teach basic concepts of bio math at early stages. The basic steps in mathematical biology are few. The initial step is to explain the biological process and raise a question about the basics. The second step is to build up a mathematical model that represents the underlying biological process. The 3rd step is to apply the method and concepts of math to obtain predictions about the model. The last step is to check whether this prediction answers the raised questions. After this, anyone can further explore the biological question by using mathematical models [1]. Now in the modern world, tobacco smoke is the most inhaled substance. Tobacco is made by mixing its agricultural form with many substances. The smoke is inhaled through the lungs. The most dangerous epidemic in the world is the smoking epidemic. Due to smoking 50% of its users died. Every year, about 60 million people die due to smoking. During the last few decades, there has been a huge boost in deaths. The death rate will rise thrice annually in 2030, almost seventy percent is in developing countries.
According to WHO, 10 million people will die due to smoking. As compared to other diseases, the death ratio is higher than all. The person who uses tobacco dies 14 years earlier than someone who does not smoke [2]. Tobacco smoking is the major reason for cancer and another disease. About 70% of people died due to tobacco-related diseases in developing countries [3]. Three million people died due to smoking yearly.
Nowadays, smoking is the most dangerous habit. The heart attack ratio is 70 percent more than a nonsmoker. There are 900 million men smokers and 200 women smokers in the world. After every 6 seconds, there is a death due to smoking. Smoking is a major cause of lung and heart attack in the world. Due to smoking, the chances of other diseases like heart attack, stroke, especially lung cancer, throat, mouth, esophagus, and pancreas are increased. Tobacco causes many tissue-related diseases. Tobacco smoke is a mixture of several toxic gases. It includes 98 of which are linked with an increased risk of cardiovascular disease, 69 of which are known to be carcinogenic. Daily, a smoker takes 1 or 2 milligrams of nicotine per cigarette. Thus, if a person smokes 5 cigarettes, it takes at least 5-10 milligrams of nicotine. The effect of smoking is not limited to the person who but also adverse effects for other people. It causes 22% of death annually.
The chemical composition of tobacco varies according to the environment. These leaves are mixed with many chemicals. Tobacco smoke contains a large number of different chemicals such as benzopyrene, NNK, aldehydes, carbon monoxide, hydrogen cyanide, phenol, nicotine, and harmala alkoids. The radioactive element polonium 210 is also occurring in tobacco. There are almost 4 thousand noxious substances in smoke which is the main reason for cancer [4,5]. The chemical composition of smoke depends on puff frequency and other materials. Nicotine is the main issue for disturbing the nervous system, rise in heartbeat, raising blood pressure, and shrinking the small blood vessels which are the main basis of wrinkles. The amount of oxygen decreased in the lungs due to carbon monoxide (CO). The natural lungs cleaner that is minuscule hairs is destroyed by hydrogen cyanide. Lead, nickel arsenic, and cadmium are also present in smoke. Some pesticides like DDT are also found in smoke. The major reason for skin and lung cancer is a toxic chemical that is present in smoke. 10 million deaths will occur in the 20 th century and 1 billion in the 21st century due to smoking [6]. Cigarette smoking affects human fertility badly [7].
The generalization of classical calculus is called fractional calculus which is concerned with the operation of integration and differentiation of fractional order. In the 19th century, fractional calculus mathematicians introduced fractional differential equations, fractional dynamics, and fractional geometry. Fractional calculus is used in almost every field of science. It is used to model physical as well as engineering processes. In many cases, standard mathematical models of integer order do not work properly. Due to this reason, fractional calculus made a major contribution to the field of mechanics, chemistry, biology, and image processing. By using fractional calculus, several physical prob-lems are solved. By using integer-order derivatives, the system shows many problems such as history and nonlocal effects. Primarily, all the studies were dependent on Caputo fractional-order and Reimann Liouville fractional (RLF) derivatives. Nowadays, it has been highlighted that these derivatives have the issue, and the issue is they have a singular kernel. That is the reason so many new definitions were presented in the studies [8][9][10][11][12][13][14][15][16]. These new definitions were very impactful because they have nonsingular kernels which are according to their needs. Caputo fractional derivatives [17], the Caputo-Fabrizio derivative [11], and AB [18] fractional derivative have differed from each other only because Caputo is defined by a power law, Fabrizio defined by using exponential decay law, and AB defined by ML law. Tateishi et al. describe the role of fractional time operator derivative in a study of anomalous diffusion [12]. With the help of analytical techniques, Bulut et al. deliberate the role of differential equations of arbitrary order [19]. The key concepts of fractional differential equations and their application are explained by Kilbas et al. [20]. Atangana and Koca examined the Keller-Segel model about a fractional derivative having a nonsingular kernel [21]. Fractional logistic maps are newly introduced by Huang et al. [22]. Zaman studied the qualitative response of the dynamics of giving up smoking [23]. The giving up smoking model linked with Caputo fractional derivative is a probe by Singh et al. [24].
Numerous studies identified sociodemographic, environmental, and behavioral risk factors such as age, sex, occupations, indoor air pollutions, smoking, and alcohol consumption [25][26][27] as being associated with the development of TB in humans. In [28], a simple model for the effect of tobacco smoking in the in-host dynamics of HIV is formulated with the aim of studying how tobacco smoking affects HIV in-host dynamics. Zoonotic tuberculosis (zTB) knowledge, prevention, and control practices via a survey in Bangladesh considering impact of smoking are also in [29].
In this work, we get the approximate solutions of the fractional smoking model by using the Atangana-Toufik method.

Basic Concepts of Fractional Operators
Definition 1. For a function gðtÞ ∈ W 1 2 ð0, 1Þ, b > a and σ ∈ ½0, 1, the definition of AB derivative in the Caputo sense is given by where Computational and Mathematical Methods in Medicine By using ST for (1), we obtain Definition 2. The Laplace transform (LT) of the Caputo fractional derivative of a function gðtÞ of order σ > 0 is defined as where E σ,σ 1 is the two-parameter ML function with σ, σ 1 > 0: Further, the ML function satisfies the following equation [17].
Definition 4. Suppose that gðtÞ is continuous on an open interval ða, bÞ, then the fractal-fractional integral of gðtÞ of order σ having ML type kernel and given by

Model Formulation
We will study the giving up smoking model for overall population at time t. We separate the population into 5 groups, potential smokers PðtÞ, occasional smokers LðtÞ, heavy smokers SðtÞ, temporary quitters QðtÞ, and smokers who quit permanently RðtÞ specified by TðtÞ = PðtÞ + LðtÞ + SðtÞ + QðtÞ + RðtÞ. Due to smoking, the chances of other diseases like heart attack, stroke, especially lung cancer, throat, mouth, esophagus, and pancreas are increased [24]. The model is developed as follows With the initial conditions The rate of change between potential smoker and occasional smokers is represented by b, a represents the rate of natural death, the rate of occasional smoker and temporary smokers by c, the rate of change between quitters and smoker is shown by f , the rate of giving up smoking is shown by d, fraction of temporary giving up smoker is represented by ð1 − eÞ (at the rate of d), e shows the remaining fraction of smokers who give up smoking forever (at a rate d).

Qualitative Analysis
By substituting the values of parameters in given system of differential equations and the rate of change with respect to time is zero, we get By simplifying the above equations, we get disease-free equilibrium, denoted by E 0 , i.e., E 0 = ð1, 0, 0, 0, 0Þ.

Stability Analysis and Reproductive Number.
It is important to find the verge conditions to check the status of population, whether the disease persist or dies out. In case of disease free equilibrium point, R 0 < 1, which shows Computational and Mathematical Methods in Medicine that the disease will die out. In case of endemic equilibrium, R 0 > 1. Consider the Jacobian matrix (JM) as Since the JM is J = F − V where We know that K = FV −1 and using the relation jK − λIj = 0 solving on mathematica for the Eigenvalue λ, which represents the reproductive number R 0 , i.e., Hence, R 0 = 0:431034 < 1, according to the given parameter values.
Proof. E 0 of the given system is locally asymptotically stable if Re ðλÞ < 0 where λ can be evaluated from the relation jJ 0 − λIj = 0.
By using the relation jJ 0 − λIj = 0, we get. Re ðλÞ as All the Eigenvalues are negative real parts which represent that the given system is locally asymptotically stable. Theorem 6. When the reproductive number R 0 > 1, the endemic equilibrium points E 1 of the PLSQR model is globally asymptotically stable.
Proof. The Lyapunov function can be written as Therefore, applying the derivative respect to t on both sides yields Now, we can write their values for derivatives as follows Putting We can organize the above as follows Computational and Mathematical Methods in Medicine To avoid the complexity, the above can be written as where It is concluded that if Σ < Ω, this yields, dM/dt < 0, however when, P = P * , L = L * , S = S * , We can see that the largest compact invariant set for the suggested model in is the point fE 1 g the endemic equilibrium of the considered model. By the help of the Lasalles invariance concept, it follows that E 1 is globally asymptotically stable in Γ if Σ < Ω.

Atangana-Baleanu Caputo Sense with Mittag-Leffler Kernel
By applying AB fractional derivative of order σ and σ ∈ ð0, 1, into ML kernel, then, the system (8) becomes The initial conditions associated with the system (37) are We will discuss the numerical value of solution of system for different values of ρ. With the help of iterative method and the Padè approximation results are obtained.
We Taking ST on both sides of (37), we get

Computational and Mathematical Methods in Medicine
Rearranging, we get Now taking inverse ST on both sides of equation (52), we get ð61Þ We next attain the following recursive formula. ð66Þ And the solution of (62) is for all x, r ∈ X, and 0 ≤ θ < 1: Suppose that H is Picard Hstable. Suppose that system (62), we have where 1 − σ/qðσÞσΓðσ + 1ÞN σ ð−1/1 − σV σ Þ is the fractional Lagrange multiplier.

Theorem 8.
Define K be a self-map is given by Proof. In first step, we will show that K is fixed point 7 Computational and Mathematical Methods in Medicine ð91Þ Applying the properties of the norm and also using the triangular inequality, we obtain Computational and Mathematical Methods in Medicine K fulfills the conditions associated with Theorem 6, when And we add that K is Picard K-stable.
Theorem 9. Prove that system (27) has special solution is unique.
Proof. Let H be the Hilbert space defined as H = L 2 ððp, qÞ × ð0, TÞÞ where So, the following operators are considered −aR + edS: We establish the inner product of where ðS 11 − S 12 , I 21 − I 22 , A 31 − A 32 , T 41 − T 42 , R 51 − R 52 Þ are the special solutions of the system. By using the inner function and the norm, we have

Conclusion
The advanced numerical scheme of fractional differential equation has been investigated in this article for smoking model by using ATM. With the help of fixed point theory uniqueness and stability of the smoking, the model has been examined. System is analyzed qualitatively to verify the steady state position of the dynamic. Proposed system is analyzed locally; also, global stability has been made using first derivative of Lyapunov. The arbitrary derivative of fractional order has been taken in ATM with no singular kernel. Effective results are obtained for the proposed model. Also discussed some theoretical results and proved the efficiency of the proposed techniques. Numerical simulations are carried out to check the actual behavior of the dynamic using the advanced ATM. These results will be helpful to understand further analysis and to control different outbreak caused by smoking.

Data Availability
No data were used to support this study.