Next Article in Journal
Determination of a Nonlinear Coefficient in a Time-Fractional Diffusion Equation
Next Article in Special Issue
Understanding the Role of Intraspecific Disease Transmission and Quarantine on the Dynamics of Eco-Epidemiological Fractional Order Model
Previous Article in Journal
Exact Solutions and Cosmological Constraints in Fractional Cosmology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of Age-Structure Smoking Models with Government Intervention Coverage under Fractal-Fractional Order Derivatives

1
College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan 030024, China
2
Department of Mathematics, Taiyuan University of Technology, Taiyuan 030024, China
3
Department of Mathematics and Statistics, Tshwane University of Technology, Pretoria 0183, South Africa
4
Department of Epidemiology and Biostatistics, School of Public Health, University of Medical Sciences, Ondo City 351101, Nigeria
5
Department of Mathematical and Computer Sciences, University of Medical Sciences, Ondo City 353301, Nigeria
6
Department of Mathematical Science, Montana Technological University, Butte, MT 59701, USA
7
School of Informatics, Computing and Cyber Systems, Northern Arizona University, Flagstaff, AZ 86011, USA
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(5), 370; https://doi.org/10.3390/fractalfract7050370
Submission received: 5 March 2023 / Revised: 27 April 2023 / Accepted: 27 April 2023 / Published: 29 April 2023

Abstract

:
The rising tide of smoking-related diseases has irreparably damaged the health of both young and old people, according to the World Health Organization. This study explores the dynamics of the age-structure smoking model under fractal-fractional (F-F) derivatives with government intervention coverage. We present a new fractal-fractional model for two-age structure smokers in the Caputo–Fabrizio framework to emphasize the potential of this operator. For the existence-uniqueness criterion of the given model, successive iterative sequences are defined with limit points that are the solutions of our proposed age-structure smoking model. We also use the functional technique to demonstrate the proposed model stability under the Ulam–Hyers condition. The two age-structure smoking models are numerically characterized using the Newton polynomial. We observe that in Groups 1 and 2, a change in the fractal-fractional orders has a direct effect on the dynamics of the smoking epidemic. Moreover, testing the inherent effectiveness of government interventions shows a considerable impact on potential, occasional, and temporary smokers when the fractal-fractional order is 0.95. It is the view that this study will contribute to the applicability of the schemes, the rich dynamics of the fractal, and the fractional perspective of future predictions.

1. Introduction

The tobacco epidemic is one of the world’s most serious public health hazards, claiming the lives of more than 8 million people each year. The primary carcinogen that has been identified globally is tobacco. The direct use of tobacco contributes to more than 7 million of these fatalities, whereas second-hand smoke exposure by non-smokers accounts for 1.2 million deaths. According to the World Health Organization’s (WHO) projection, a small decrease in the smoking rate indicates that the WHO will likely miss its 2030 goal of reducing the prevalence of smoking to 20% [1]. In addition, scientists and researchers predict that between 2020 and 2040, the smoking-related cancer mortality rate will increase by 44% in men and 53% in women, resulting in an excess of 8.6 million fatalities. According to the study, 1.4 million fatalities could be avoided if the objective set for 2030 is achieved and the trend persists into 2040 [2]. The WHO projection of smoking has shown a minor increase in the number of young female smokers, a massive loss of life in the working-age population, and a high level of smoke-related sickness among the elderly. This has sparked a large discussion among all stakeholders regarding the future cancer landscape because the world population has been rapidly ageing over the past few decades, due to some countries’ policies to lower fertility rates and raise the life expectancy.
In the last century, mathematical modeling has been shown to be a very effective technique for tracking and managing many diseases by supplying crucial qualitative data. Mathematical models offer a substitute for real data, which are sometimes hard to come by and are costly. The data have a significant impact on the models and play a key role in understanding the dynamics of infectious diseases and developing better methods to prevent their spread in the future [3,4,5,6,7,8]. Numerous researchers have used fractional-order models to study the dynamics of different diseases (see [9] and references therein). The memory effect of fractional calculus facilitates accurate physical system predictions. There are many fractional order derivatives, such as the Riemann–Liouville–Caputo, Atangana–Baleanu–Caputo (ABC), and Caputo–Fabrizio (CF) that have been proposed. Interest in fractional calculus theory and its demonstrated applications are currently growing. According to the literature, the reproduction of fresh tasks that simulate various kinds of real-world circumstances is responsible for the appearance of innovative operators. In this regard, a variety of practical problems have been solved using different fractional operators with different numerical solutions. In [10], the authors examined a recent study on the human liver using the CF fractional derivative, and concluded that the infection rate decreased when the operator value was reduced from 1. The authors in [11] used both the CF and ABC fractional derivatives to study the dynamics of Q fever. In their research, they found that the paths of some fractional and integer orders converge to the same endemic equilibrium point. With regards to smoking, in [12], the authors introduced a mathematical model for giving up smoking for the first time, and in [13] the author expanded on the study in [12], by adding a compartment of occasional smokers in the model, and presented the qualitative behavior of the model. The optimal control model was used in the smoking dynamics described in [14]. They explored two different control variables in the form of education and treatment campaigns aimed at changing people’s attitudes toward smoking, and they showed for the first time, that there is an optimal control for the problem. In [15], the authors extended the smoking model associated with the Caputo fractional derivatives. In [16], the authors also considered the existence of theoretical and numerical solutions for the smoking model under the Caputo–Fabrizio fractional derivative. Researchers have developed new complex fractional operators with multiple operators to mimic real-world problems. These new operators are sufficient to describe extremely complex phenomena that cannot be described by a single operator. The fractal-fractional order models are commonly used in epidemiology to understand the dynamics of infectious diseases. In the realm of mathematical biology, in [17], a listeriosis disease model that was incorporated into the Caputo and ABC fractal-fractional derivatives was investigated. The study in [18] used a fractal-fractional model to study the sensitivity analysis of COVID-19. In [19], the authors studied the dynamics and sensitivity analysis of pine wilt disease with asymptomatic carriers using the ABC fractal-fractional operators, and found that eliminating beetles might significantly lower the number of infections. A new fractal-fractional age-structure model for the omicron SARS-CoV-2 variant under the Caputo–Fabrizio fractional order derivative was investigated by [20]. For many applications of fractal-fractional operators, see [21] and references therein.
In this paper, we use fractal-fractional operators to create a model to capture the repeated patterns and memory of the exponential decay for the two age-structure smoking models. Fractals are naturally found in most biological items. As a result, modeling epidemiology in two dimensions is a true reflection of the projection that results from the study. Memory and heredity traits, which are complex behavioral patterns of biological systems, are the goal of dealing with fractional order systems in our newly designed two age-structure smoking models, which allows us to create a more realistic approach to biological systems. The memory function allows for fractional order models to incorporate more knowledge from the past, thereby allowing for a more accurate prediction and translation. In addition, the hereditary property specifies the genetic profile as well as the age and status of the immune system. From these, we see no single mathematical formulation of the two age-structure smoking models in fractal-fractional dynamics. Comparing our result to [15,16,22], this study reveals some essential qualitative information that is important for the government in the course of developing control interventions for cancer deaths from smoking.
The structure of the paper is as follows: In Section 2, we briefly discuss the essential definitions and their relationship to fractal-fractional operators. In Section 3, we formulate the model in the classical and fractal-fractional cases. The positivity, boundedness, and invariant region for our proposed model are examined in Section 4. In Section 5, we discuss the analysis of the existence and uniqueness using successive iterative sequences, and limit the point principle. We cover the stability findings for the HU stability category in Section 6. In Section 7 and Section 8, we provide numerical and simulation results to verify the effectiveness of the two fractal-fractional age-structure smoking models. Finally, we offer some conclusions in Section 9 with some future recommendations.

2. Preliminaries

The fundamental definitions that will be used throughout the work are contained in this section. The F-F derivative, is covered generally by these definitions.
Definition 1([23,24]).
The fractal-fractional derivatives of function A ( ε ) in the Riemann–Liouville and Caputo senses of the power-type kernel is expressed as
0 F F P D ε α , β [ A ( ε ) ] = 1 Γ ( n ε ) 0 ε d A ( s ) d s β ( ε s ) n α 1 d s ,
where n 1 < α , β n N , with α is the fractional order index, β is the fractal dimension, and
d A ( ε ) d ε β = lim ε s A ( ε ) A ( s ) ε β s β
defines the fractal derivative of function A ( ε ) .
Definition 2 ([23,24]).
Suppose that there is a continuous and fractal differentiable in ( m , w ) , then the F-F derivative of A ( ε ) of order α in the CF, involving exponential decay, is given by
0 F F E D ε α , β [ A ( ε ) ] = M ( α ) 1 α d d ε β 0 ε exp α ( ε s ) ( 1 α ) A ( s ) d s ,
with α > 0 , β n N , and M ( α ) is the normalization function, such that M ( 0 ) = M ( 1 ) = 1 .
Definition 3 ([23,24]).
The fractal-fractional integral associated with the CF-type fractal-fractional derivative is defined as
0 F F E I ε α , β [ A ( ε ) ] = β ( 1 α ) ε β 1 M ( α ) A ( ε ) + α β M ( α ) 0 ε s β 1 A ( s ) d s ,

3. Model Formulation

We consider the entire population N i ( ε ) = N 1 ( ε ) + N 2 ( ε ) and N i = P i + Q i + S i + R i + U i , where i = 1 , 2 , as two age groups, represent Group-1 with individuals under 70 years of age ( N 1 = P 1 + Q 1 + S 1 + R 1 + U 1 ) and Group-2 with individuals aged 70 years and over ( N 2 = P 2 + Q 2 + S 2 + R 2 + U 2 ). We consider the age limit from the real data in [25]. This compartmental model captures the government intervention coverage x ( ε ) that depends on time ε and the proportion that varies between zero and one. We assume that there is no government intervention coverage when x ( ε ) = 0 , while x ( ε ) = 1 represents the perfect or maximum government intervention coverage. In the model, x ( ε ) is replaced by ( 1 x ( ε ) ) , since in real life, government intervention coverage may not be 100% effective over time. Table 1 give detailed description of the variables and parameters in the model, respectively. Considering the inter-relationship, we formulate the following deterministic system of non-linear differential equations;
d P 1 ( ε ) d ε = Λ 1 ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 1 μ 1 P 1 Λ 2 P 1 , d Q 1 ( ε ) d ε = ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 1 λ 1 Q 1 μ 1 Q 1 Λ 2 Q 1 , d S 1 ( ε ) d ε = λ 1 Q 1 + ν 1 Q 1 S 1 ( μ 1 + η 1 ) S 1 Λ 2 S 1 , d R 1 ( ε ) d ε = η 1 ( 1 δ 1 ) S 1 ν 1 Q 1 S 1 μ 1 R 1 Λ 2 R 1 , d U 1 ( ε ) d ε = η 1 δ 1 S 1 μ 1 U 1 Λ 2 U 1 , d P 2 ( ε ) d ε = Λ 2 P 1 ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 2 μ 2 P 2 , d Q 2 ( ε ) d ε = Λ 2 Q 1 + ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 2 λ 2 Q 2 μ 2 Q 2 , d S 2 ( ε ) d ε = Λ 2 S 1 + λ 2 Q 2 + ν 2 Q 2 S 2 ( μ 2 + η 2 ) S 2 , d R 2 ( ε ) d ε = Λ 2 R 1 + η 2 ( 1 δ 2 ) S 2 ν 2 Q 2 S 2 μ 2 R 2 , d U 2 ( ε ) d ε = Λ 2 U 1 + η 2 δ 2 S 2 μ 2 U 2 .
A further change that is made to the aforementioned ordinary differential model (1) is the creation of a fractional-order system of order α . The fractional-order situation is taken into consideration, since it is significantly distinct from other fractional-order systems because it has non-local characteristics and hereditary traits that are not present in the extensively used integer-order differential operators in biology. Non-integer modeling makes predictions for the future by using the past and present. The fractional order is represented by the α value, while the fractal dimension is represented by β . The fractal perspective would improve the comprehension of the phenomena because in reality, items are made up of fractals. Thus, the following is the age-structure smoking model with government intervention under the CF F-F order derivative:
F F E D ε α , β P 1 ( ε ) = Λ 1 ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 1 μ 1 P 1 Λ 2 P 1 , F F E D ε α , β Q 1 ( ε ) = ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 1 λ 1 Q 1 μ 1 Q 1 Λ 2 Q 1 , F F E D ε α , β S 1 ( ε ) = λ 1 Q 1 + ν 1 Q 1 S 1 ( μ 1 + η 1 ) S 1 Λ 2 S 1 , F F E D ε α , β R 1 ( ε ) = η 1 ( 1 δ 1 ) S 1 ν 1 Q 1 S 1 μ 1 R 1 Λ 2 R 1 , F F E D ε α , β U 1 ( ε ) = η 1 δ 1 S 1 μ 1 U 1 Λ 2 U 1 , F F E D ε α , β P 2 ( ε ) = Λ 2 P 1 ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 2 μ 2 P 2 , F F E D ε α , β Q 2 ( ε ) = Λ 2 Q 1 + ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 2 λ 2 Q 2 μ 2 Q 2 , F F E D ε α , β S 2 ( ε ) = Λ 2 S 1 + λ 2 Q 2 + ν 2 Q 2 S 2 ( μ 2 + η 2 ) S 2 , F F E D ε α , β R 2 ( ε ) = Λ 2 R 1 + η 2 ( 1 δ 2 ) S 2 ν 2 Q 2 S 2 μ 2 R 2 , F F E D ε α , β U 2 ( ε ) = Λ 2 U 1 + η 2 δ 2 S 2 μ 2 U 2 ,
where time ε > 0 and F F E D ε α , β is the CF F-F derivative for α and β ( 0 , 1 ] with initial conditions
P 1 ( 0 ) = P 1 ( 0 ) 0 ,   Q 1 ( 0 ) = Q 1 ( 0 ) 0 ,   S 1 ( 0 ) = S 1 ( 0 ) 0 ,   R 1 ( 0 ) = R 1 ( 0 ) 0 ,   U 1 ( 0 ) = U 1 ( 0 ) 0 ,
P 2 ( 0 ) = P 2 ( 0 ) 0 ,   Q 2 ( 0 ) = Q 2 ( 0 ) 0 ,   S 2 ( 0 ) = S 2 ( 0 ) 0 ,   R 2 ( 0 ) = R 2 ( 0 ) 0 ,   U 2 ( 0 ) = U 2 ( 0 ) 0 .

4. Dynamics of the Model

In order to show that the age-structure model given by (1) and (2) are epidemiologically meaningful, we have to prove that the associated state variables of the model stay positive. It can also be explained that the solution of the age-structure model with positive initial conditions will remain positive and bounded every time ε 0 .
Lemma 1.
Let P 1 ( 0 ) 0 ,   Q 1 ( 0 ) 0 ,   S 1 ( 0 ) 0 ,   R 1 ( 0 ) 0 ,   U 1 ( 0 ) 0 ,   P 2 ( 0 ) 0 ,   Q 2 ( 0 ) 0 ,   S 2 ( 0 ) 0 ,   R 2 ( 0 ) 0 ,   U 2 ( 0 ) 0 . Such that the solution set
Ψ = { P 1 ( ε ) , Q 1 ( ε ) , S 1 ( ε ) ,   R 1 ( ε ) , U 1 ( ε ) , P 2 ( ε ) , Q 2 ( ε ) , S 2 ( ε ) ,   R 2 ( ε ) , U 2 ( ε ) } R + 10
of the smoking epidemic model (1) is positive for all ε > 0 . Furthermore,
lim sup t N 1 ( ε ) Λ 1 μ 1 ,   lim sup t N 2 ( ε ) Λ 2 μ 2 .
Proof. 
Let ε 0 = sup { P 1 ( 0 ) 0 ,   Q 1 ( 0 ) 0 ,   S 1 ( 0 ) 0 ,   R 1 ( 0 ) 0 ,   U 1 ( 0 ) 0 ,   P 2 ( 0 ) 0 ,   Q 2 ( 0 ) 0 ,   S 2 ( 0 ) 0 ,   R 2 ( 0 ) 0 ,   U 2 ( 0 ) 0 . } [ 0 , ε ] . Taking all initial values that are greater than zero, then, ε 0 > 0 . Let us now consider the model initial dynamical Equation (1).
d P 1 d ε = Λ 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 μ 1 P 1 Λ 2 P 1 .
For simplicity, we let M 1 * = ( 1 + x ) + β 1 S 1 + β 2 S 2 + Λ 2 + μ 1 . Thus,
d P 1 d ε + M 1 * P 1 ( ε ) = Λ 1 .
This leads to
d d ε ( P 1 ( ε ) e x p ( 0 ε 0 M 1 * ( s ) d s ) ) = ( Λ 1 ) e x p ( 0 ε 0 M 1 * ( s ) d s ) .
So, we obtain
P 1 ( ε ) = P 1 ( 0 ) e x p ( 0 ϵ 0 M 1 * ( s ) d s ) [ 0 ε 0 ( Λ 1 ) e x p ( 0 ε 0 M 1 * ( s ) d s ) ] + e x p 0 ϵ 0 M 1 * ( s ) d s 0 ϵ 0 Λ 1 e x p 0 s M 1 * ( s ) d s
Hence, we proved that P 1 ( 0 ) > 0 for all ε > 0 . In a similar manner, we can prove that Q 1 ( 0 ) 0 ,   S 1 ( 0 ) 0 ,   R 1 ( 0 ) 0 ,   U 1 ( 0 ) 0 ,   P 2 ( 0 ) 0 ,   Q 2 ( 0 ) 0 ,   S 2 ( 0 ) 0 ,   R 2 ( 0 ) 0 ,   U 2 ( 0 ) 0 for all ε > 0 .
Now, for boundedness, we have the following expression d N i d ε = Λ i μ i N i , i = 1 , 2 , such that Λ i μ i N i d N i d ε Λ i μ i N i , and hence, lim inf ε N i ( ε ) Λ i μ i lim sup ε N i ( ε ) Λ i μ i .
Considered invariant, we let Φ = Ω i R + × R + , where
Ψ = { ( P 1 ( ε ) , Q 1 ( ε ) , S 1 ( ε ) ,   R 1 ( ε ) , U 1 ( ε ) , P 2 ( ε ) , Q 2 ( ε ) , S 2 ( ε ) ,   R 2 ( ε ) , U 2 ( ε ) ) R + 10 : N i ( ε ) Λ i μ i ) } .
Now, for the positive invariant of R + 10 ,
d N i d ε = Λ i μ i N i , i = 1 , 2 ,
solving this inequality, we obtain the following
N i ( ε ) N i ( 0 ) e μ i ( ε ) + Λ i μ i ( 1 e μ i ( ε ) ) .
Therefore,
lim sup ε N i ( ε ) Λ i μ i .
Furthermore, from (2) we have
F F E D ε α , β P 1 ( ε ) | P 1 = 0 = Λ 1 0 , F F E D ε α , β Q 1 ( ε ) | Q 1 = 0 = ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 1 0 , F F E D ε α , β S 1 ( ε ) | S 1 = 0 = λ 1 Q 1 + ν 1 Q 1 S 1 0 , F F E D ε α , β R 1 ( ε ) | R 1 = 0 = η 1 ( 1 δ 1 ) S 1 0 , F F E D ε α , β U 1 ( ε ) | U 1 = 0 = η 1 δ 1 S 1 0 , F F E D ε α , β P 2 ( ε ) | P 2 = 0 = Λ 2 P 1 0 , F F E D ε α , β Q 2 ( ε ) | Q 2 = 0 = Λ 2 Q 1 + ( 1 x ( ε ) ) ( β 1 S 1 + β 2 S 2 ) P 2 0 , F F E D ε α , β S 2 ( ε ) | S 2 = 0 = Λ 2 S 1 + λ 2 Q 2 + ν 2 Q 2 S 2 0 , F F E D ε α , β R 2 ( ε ) | R 2 = 0 = Λ 2 R 1 + η 2 ( 1 δ 2 ) S 2 0 , F F E D ε α , β U 2 ( ε ) | U 2 = 0 = Λ 2 U 1 + η 2 δ 2 S 2 0 ,
and from the results in [24], the result of model (2) cannot escape from the hyperplanes, and this concludes the prologue to Lemma 1. Hence, the two age-structure smoking transmission models are deemed to be mathematically well posed Φ . □

5. Existence and Uniqueness Results

In this section, we analyze the existing and unique findings for our fractal-fractional models (2) using a successive iterative technique, and limit the points. In that case, we reformulate model (2) with the help of Definition 2, as follows:
P 1 ( ε ) = P 1 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) Λ 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 μ 1 P 1 Λ 2 P 1 + α β M ( α ) 0 ε s β 1 Λ 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 μ 1 P 1 Λ 2 P 1 d s , Q 1 ( ε ) = Q 1 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 λ 1 Q 1 μ 1 Q 1 Λ 2 Q 1 + α β M ( α ) 0 ε s β 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 λ 1 Q 1 μ 1 Q 1 Λ 2 Q 1 d s , S 1 ( ε ) = S 1 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) λ 1 Q 1 + ν 1 Q 1 S 1 ( μ 1 + η 1 ) S 1 Λ 2 S 1 + α β M ( α ) 0 ε s β 1 λ 1 Q 1 + ν 1 Q 1 S 1 ( μ 1 + η 1 ) S 1 Λ 2 S 1 d s , R 1 ( ε ) = R 1 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) η 1 ( 1 δ 1 ) S 1 ν 1 Q 1 S 1 μ 1 R 1 Λ 2 R 1 + α β M ( α ) 0 ε s β 1 η 1 ( 1 δ 1 ) S 1 ν 1 Q 1 S 1 μ 1 R 1 Λ 2 R 1 d s , U 1 ( ε ) = U 1 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) η 1 δ 1 S 1 μ 1 U 1 Λ 2 U 1 + α β M ( α ) 0 ε s β 1 η 1 δ 1 S 1 μ 1 U 1 Λ 2 U 1 d s , P 2 ( ε ) = P 2 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) Λ 2 P 1 ( β 1 S 1 + β 2 S 2 ) P 2 μ 2 P 2 + α β M ( α ) 0 ε s β 1 Λ 2 P 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 2 μ 2 P 2 d s , Q 2 ( ε ) = Q 2 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) Λ 2 Q 1 + ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 2 λ 2 Q 2 μ 2 Q 2 + α β M ( α ) 0 ε s β 1 Λ 2 Q 1 + ( β 1 S 1 + β 2 S 2 ) P 2 λ 2 Q 2 μ 2 Q 2 d s , S 2 ( ε ) = S 2 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) Λ 2 S 1 + λ 2 Q 2 + ν 2 Q 2 S 2 ( μ 2 + η 2 ) S 2 + α β M ( α ) 0 ε s β 1 Λ 2 S 1 + λ 2 Q 2 + ν 2 Q 2 S 2 ( μ 2 + η 2 ) S 2 d s , R 2 ( ε ) = R 2 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) Λ 2 R 1 + η 2 ( 1 δ 2 ) S 2 ν 2 Q 2 S 2 μ 2 R 2 + α β M ( α ) 0 ε s β 1 Λ 2 R 1 + η 2 ( 1 δ 2 ) S 2 ν 2 Q 2 S 2 μ 2 R 2 d s , U 2 ( ε ) = U 2 ( 0 ) + β ( 1 α ) ε β 1 M ( α ) Λ 2 U 1 + η 2 δ 2 S 2 μ 2 U 2 + α β M ( α ) 0 ε s β 1 Λ 2 U 1 + η 2 δ 2 S 2 μ 2 U 2 d s .
We assume functions A i ( ε ) , where i = 1 , 2 , 3 , , 10 , as given below
A 1 ( ε , P 1 ( ε ) ) = Λ 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 μ 1 P 1 Λ 2 P 1 , A 2 ( ε , Q 1 ( ε ) ) = ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 λ 1 Q 1 μ 1 Q 1 Λ 2 Q 1 , A 3 ( ε , S 1 ( ε ) ) = λ 1 Q 1 + ν 1 Q 1 S 1 ( μ 1 + η 1 ) S 1 Λ 2 S 1 , A 4 ( ε , R 1 ( ε ) ) = η 1 ( 1 δ 1 ) S 1 ν 1 Q 1 S 1 μ 1 R 1 Λ 2 R 1 , A 5 ( ε , U 1 ( ε ) ) = η 1 δ 1 S 1 μ 1 U 1 Λ 2 U 1 , A 6 ( ε , P 2 ( ε ) ) = Λ 2 P 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 2 μ 2 P 2 , A 7 ( ε , Q 2 ( ε ) ) = Λ 2 Q 1 + ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 2 λ 2 Q 2 μ 2 Q 2 , A 8 ( ε , S 2 ( ε ) ) = Λ 2 S 1 + λ 2 Q 2 + ν 2 Q 2 S 2 ( μ 2 + η 2 ) S 2 , A 9 ( ε , R 2 ( ε ) ) = Λ 2 R 1 + η 2 ( 1 δ 2 ) S 2 ν 2 Q 2 S 2 μ 2 R 2 , A 10 ( ε , U 2 ( ε ) ) = Λ 2 U 1 + η 2 δ 2 S 2 μ 2 U 2 .
For proving the desired theorem on the existence property, we have the following assumption: ( ) : For all | { P 1 ( ε ) , Q 1 ( ε ) , S 1 ( ε ) , R 1 ( ε ) , U 1 ( ε ) , P 2 * ( ε ) , Q 2 * ( ε ) , S 2 * ( ε ) , R 2 * ( ε ) , U 2 * ( ε ) L q [ 0 , 1 ] } , there exists positive constants K i > 0 , such that P 1 ( ε ) K 1 , Q 1 ( ε ) K 2 , S 1 ( ε ) K 3 , R 1 ( ε ) K 4 , U 1 ( ε ) K 5 , P 2 ( ε ) K 6 , Q 2 ( ε ) K 7 , S 2 ( ε ) K 8 , R 2 ( ε ) K 9 , U 2 ( ε ) K 10 .
Theorem 1.
Suppose that functions A i , where i = 1 , 2 , 3 , , 10 fulfill the Lipschitz conditions, provided that ( ) holds and satisfies
W = max { ψ 1 , ψ 2 , ψ 3 , ψ 4 , ψ 5 , ψ 6 , ψ 7 , ψ 8 , ψ 9 , ψ 10 } < 1 .
Proof. 
We consider 1 ( ε , P 1 ) . For P 1 ( ε ) , P 1 * ( ε ) , we estimate
| | A 1 ( ε , P 1 ( ε ) ) A 1 ( ε , P 1 * ( ε ) ) | | = | | Λ 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 μ 1 P 1 Λ 2 P 1 Λ 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 * μ 1 P 1 * Λ 2 P 1 * | | , = | | ( 1 x ) ( β 1 S 1 + β 2 S 2 ) + μ 1 + Λ 2 P 1 P 1 * | | , ( 1 x ) β 1 | | S 1 | | + ( 1 x ) β 2 | | S 2 | | + μ 1 + Λ 2 P 1 P 1 * | | , ( 1 x ) β 1 ψ 3 + ( 1 x ) β 2 ψ 8 + μ 1 + Λ 2 P 1 P 1 * | | , ψ 1 | | P 1 P 1 * | | .
Hence, A 1 ( ε , P 1 ( ε ) ) fulfills the Lipschitz property and ψ 1 < 1 .
For Q 1 ( ε ) , Q 1 * ( ε ) , we have
| | A 2 ( ε , Q 1 ( ε ) ) A 2 ( ε , Q 1 * ( ε ) ) | | = | | ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 λ 1 Q 1 μ 1 Q 1 Λ 2 Q 1 ( 1 x ) ( β 1 S 1 + β 2 S 2 ) P 1 λ 1 Q 1 * μ 1 Q 1 * Λ 2 Q 1 * | | , = | | λ 1 + μ 1 + Λ 2 Q 1 Q 1 * | | , λ 1 + μ 1 + Λ 2 | | Q 1 Q 1 * | | , ψ 2 | | Q 1 Q 1 * | | .
Thus, 2 ( ε , Q 1 ( ε ) ) satisfies the Lipschitz condition and ψ 2 < 1 .
Similarly,
| | A 3 ( ε , S 1 ( ε ) ) A 3 ( ε , S 1 * ( ε ) ) | | ψ 3 | | S 1 ( ε ) S 1 * ( ε ) | | , | | A 4 ( ε , R 1 ( ε ) ) A 4 ( ε , R 1 * ( ε ) ) | | ψ 4 | | R 1 ( ε ) R 1 * ( ε ) | | , | | A 5 ( ε , U 1 ( ε ) ) A 5 ( ε , U 1 * ( ε ) ) | | ψ 5 | | U 1 ( ε ) U 1 * ( ε ) | | , | | A 6 ( ε , P 2 ( ε ) ) A 6 ( ε , P 2 * ( ε ) ) | | ψ 6 | | P 2 ( ε ) P 2 * ( ε ) | | , | | A 7 ( ε , Q 2 ( ε ) ) A 7 ( ε , Q 2 * ( ε ) ) | | ψ 7 | | Q 2 ( ε ) Q 2 * ( ε ) | | , | | A 8 ( ε , S 2 ( ε ) ) A 8 ( ε , S 2 * ( ε ) ) | | ψ 8 | | S 2 ( t ) S 2 * ( t ) | | , | | A 9 ( ε , R 2 ( ε ) ) A 9 ( ε , R 2 * ( ε ) ) | | ψ 9 | | R 2 ( ε ) R 2 * ( ε ) | | , | | A 10 ( ε , U 2 ( ε ) ) A 10 ( ε , U 2 * ( ε ) ) | | ψ 10 | | U 2 ( ε ) U 2 * ( ε ) | | .
Hence A i satisfy the Lipschitz condition and Lipschitz constants ψ i < 1 for all i = 1 , 2 , 3 , , 10 , thus the proof is completed.
Now, we consider model (4) and reformulate in the form of kernels A i and, taking the initial conditions as
P 1 ( 0 ) = Q 1 ( 0 ) = S 1 ( 0 ) = R 1 ( 0 ) = U 1 ( 0 ) = P 2 ( 0 ) = Q 2 ( 0 ) = S 2 ( 0 ) = R 2 ( 0 ) = U 2 ( 0 ) = 0 ,
we have
P 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 1 ε , P 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 1 ε , P 1 ( ε ) d s , Q 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 2 ε , Q 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 2 ε , Q 1 ( ε ) d s , S 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 3 ε , S 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 3 ε , S 1 ( ε ) d s , R 1 ( ε ) = β ( 1 α ) ε β 1 A ( α ) A 4 ε , R 1 ( ε ) + α β A ( α ) 0 ε s β 1 A 4 ε , R 1 ( ε ) d s , U 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 5 ε , U 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 5 ε , U 1 ( ε ) d s , P 2 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 6 ε , P 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 6 ε , P 2 ( ε ) d s , Q 2 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 7 ε , Q 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 7 ε , Q 2 ( ε ) d s , S 2 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 8 ε , S 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 8 ε , S 2 ( ε ) d s , R 2 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 9 ε , R 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 9 ε , R 2 ( ε ) d s , U 2 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 10 ε , U 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 10 ε , U 2 ( ε ) d s .
For the iterative scheme of model (9), we have
P 1 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 1 ε , P 1 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 1 ε , P 1 n 1 ( ε ) d s , Q 1 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 2 ε , Q 1 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 2 ε , Q 1 n 1 ( ε ) d s , S 1 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 3 ε , S 1 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 3 ε , S 1 n 1 ( ε ) d s , R 1 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 4 ε , R 1 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 4 ε , R 1 n 1 ( ε ) d s , U 1 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 5 ε , U 1 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 5 ε , U 1 n 1 ( ε ) d s , P 2 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 6 ε , P 2 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 6 ε , P 2 n 1 ( ε ) d s , Q 2 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 7 ε , Q 2 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 7 ε , Q 2 n 1 ( ε ) d s , S 2 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 8 ε , S 2 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 8 ε , S 2 n 1 ( ε ) d s , R 2 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 9 ε , R 2 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 9 ε , R 2 n 1 ( ε ) d s , U 2 n ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 6 ε , P 2 n 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 6 ε , P 2 n 1 ( ε ) d s , P 1 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 1 ε , P 1 n ( ε ) + α β M ( α ) 0 ε s β 1 A 1 ε , P 1 n ( ε ) d s , Q 1 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 2 ε , Q 1 n ( ε ) + α β M ( α ) 0 ε s β 1 A 2 ε , Q 1 n ( ε ) d s , S 1 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 3 ε , S 1 n ( ε ) + α β M ( α ) 0 ε s β 1 A 3 ε , S 1 n ( ε ) d s , R 1 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 4 ε , R 1 n ( ε ) + α β M ( α ) 0 ε s β 1 A 4 ε , R 1 n ( ε ) d s , U 1 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 5 ε , U 1 n ( ε ) + α β M ( α ) 0 ε s β 1 A 5 ε , U 1 n ( ε ) d s , P 2 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 6 ε , P 2 n ( ε ) + α β M ( α ) 0 ε s β 1 A 6 ε , P 2 n ( ε ) d s , Q 2 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 7 ε , Q 2 n ( ε ) + α β M ( α ) 0 ε s β 1 A 7 ε , Q 2 n ( ε ) d s , S 2 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 8 ε , S 2 n ( ε ) + α β M ( α ) 0 ε s β 1 A 8 ε , S 2 n ( ε ) d s , R 2 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 9 ε , R 2 n ( ε ) + α β M ( α ) 0 ε s β 1 A 9 ε , R 2 n ( ε ) d s , U 2 n + 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 6 ε , P 2 n ( ε ) + α β M ( α ) 0 ε s β 1 A 6 ε , P 2 n ( ε ) d s .
Next, we formulate the recursive inequalities for the differences, ;
A 1 ( P 1 n + 1 ( ε ) ) = P 1 n + 1 P 1 n = ( β ( 1 α ) ε β 1 M ( α ) A 1 ( ε , P 1 n ( ε ) ) + α β M ( α ) 0 ε s β 1 A 1 ( s , P 1 n ( s ) ) d s ) ( β ( 1 α ) ε β 1 M ( α ) A 1 ( ε , P 1 n 1 ( ε ) ) + α β M ( α ) 0 t s β 1 A 1 ( s , P 1 n 1 ( s ) ) d s ) .
A 2 ( Q 1 n + 1 ( ε ) ) = Q 1 n + 1 Q 1 n = ( β ( 1 α ) ε β 1 M ( α ) A 2 ( ε , Q 1 n ( ε ) ) + α β M ( α ) 0 ε s β 1 A 2 ( s , Q 1 n ( s ) ) d s ) ( β ( 1 α ) ε β 1 M ( α ) A 2 ( ε , Q 1 n 1 ( ε ) ) + α β M ( α ) 0 t s β 1 A 2 ( s , Q 1 n 1 ( s ) ) d s ) .
A 3 ( S 1 n + 1 ( ε ) ) = S 1 n + 1 S 1 n = ( β ( 1 α ) ε β 1 M ( α ) A 3 ( ε , S 1 n ( ε ) ) + α β M ( α ) 0 ε s β 1 A 3 ( s , S 1 n ( s ) ) d s ) ( β ( 1 α ) ε β 1 M ( α ) A 3 ( ε , S 1 n 1 ( ε ) ) + α β M ( α ) 0 t s β 1 A 3 ( s , S 1 n 1 ( s ) ) d s ) .
Similarly,
A 4 ( R 1 n + 1 ( ε ) ) = R 1 n + 1 R 1 n ,   A 5 ( U 1 n + 1 ( ε ) ) = U 1 n + 1 U 1 n , A 6 ( P 2 n + 1 ( ε ) ) = P 2 n + 1 P 2 n ,   A 7 ( Q 2 n + 1 ( ε ) ) = Q 2 n + 1 Q 2 n , A 8 ( S 2 n + 1 ( ε ) ) = S 2 n + 1 S 2 n ,   A 9 ( R 2 n + 1 ( ε ) ) = R 2 n + 1 R 2 n , A 10 ( U 2 n + 1 ( ε ) ) = U 2 n + 1 U 2 n .
Now, we apply the norm on A i , we obtain
| | A 1 ( P 1 n + 1 ( ε ) ) | | = | | P 1 n + 1 P 1 n | | = | | ( β ( 1 α ) ε β 1 M ( α ) A 1 ( ε , P 1 n ( ε ) ) + α β M ( α ) 0 ε s β 1 A 1 ( s , P 1 n ( s ) ) d s ) ( β ( 1 α ) ε β 1 M ( α ) A 1 ( ε , P 1 n 1 ( ε ) ) + α β M ( α ) 0 ε s β 1 A 1 ( s , P 1 n 1 ( s ) ) d s ) | | = β ( 1 α ) ε β 1 M ( α ) | | A 1 ( ε , P 1 n ( ε ) ) A 1 ( ε , P 1 n 1 ( ε ) ) | | + α β M ( α ) 0 ε s β 1 | | A 1 ( s , P 1 n ( s ) ) A 1 ( s , P 1 n 1 ( s ) ) | | d s .
| | A 2 ( Q 1 n + 1 ( ε ) ) | | = | | Q 1 n + 1 Q 1 n | | = | | ( β ( 1 α ) ε β 1 M ( α ) A 2 ( ε , Q 1 n ( ε ) ) + α β M ( α ) 0 ε s β 1 A 2 ( s , Q 1 n ( s ) ) d s ) ( β ( 1 α ) ε β 1 M ( α ) A 2 ( ε , Q 1 n 1 ( ε ) ) + α β M ( α ) 0 ε s β 1 A 2 ( s , Q 1 n 1 ( s ) ) d s ) | | = β ( 1 α ) ε β 1 M ( α ) | | A 2 ( ε , Q 1 n ( ε ) ) A 2 ( ε , Q 1 n 1 ( ε ) ) | | + α β M ( α ) 0 ε s β 1 | | A 2 ( s , Q 1 n ( s ) ) A 2 ( s , Q 1 n 1 ( s ) ) | | d s .
Similarly,
| | A 3 ( S 1 n + 1 ( ε ) ) | | = | | S 1 n + 1 S 1 n | | ,   | | A 4 ( R 1 n + 1 ( ε ) ) | | = | | R 1 n + 1 R 1 n | | , | | A 5 ( U 1 n + 1 ( ε ) ) | | = | | U 1 n + 1 U 1 n | | ,   | | A 6 ( P 2 n + 1 ( ε ) ) | | = | | P 2 n + 1 P 2 n | | , | | A 7 ( Q 2 n + 1 ( ε ) ) | | = | | Q 2 n + 1 Q 2 n | | ,   | | A 8 ( S 2 n + 1 ( ε ) ) | | = | | S 2 n + 1 S 2 n | | , | | A 9 ( R 2 n + 1 ( ε ) ) | | = | | R 2 n + 1 R 2 n | | ,   | | A 10 ( U 2 n + 1 ( ε ) ) | | = | | U 2 n + 1 U 2 n | | .
Theorem 2.
The fractal-fractional of age-structure smoking epidemic model (3) has at least one solution if we have
W = max { ψ 1 , ψ 2 , ψ 3 , ψ 4 , ψ 5 , ψ 6 , ψ 7 , ψ 8 } < 1 .
Proof. 
For simplicity, we let the function
Ω 1 n ( ε ) = P 1 n + 1 ( ε ) P 1 ( ε ) ,   Ω 2 n ( ε ) = Q 1 n + 1 ( ε ) Q 1 ( ε ) , Ω 3 n ( ε ) = S 1 n + 1 ( ε ) S 1 ( ε ) ,   Ω 4 n ( ε ) = R 1 n + 1 ( ε ) R 1 ( ε ) , Ω 5 n ( ε ) = U 1 n + 1 ( ε ) U 1 n ( ε ) ,   Ω 6 n ( ε ) = P 2 n + 1 P 2 ( ε ) , Ω 7 n ( ε ) = Q 2 n + 1 ( ε ) Q 2 ( ε ) ,   Ω 8 n ( ε ) = S 2 n + 1 S 2 ( ε ) , Ω 9 n ( ε ) = Q 2 n + 1 ( ε ) Q 2 ( ε ) ,   Ω 10 n ( ε ) = S 2 n + 1 S 2 ( ε ) .
Now, we apply the norm on (18), and we obtain
| | Ω 1 n ( ε ) | | = | | P 1 n + 1 ( ε ) P 1 ( ε ) | | = | | ( β ( 1 β ) ε β 1 M ( α ) A 1 ( t , P 1 n ( ε ) ) + α β M ( α ) 0 ε s β 1 A 1 ( s , P 1 n ( s ) ) d s ) ( β ( 1 α ) ε α 1 M ( α ) A 1 ( ε , P 1 ( ε ) ) + α β M ( α ) 0 ε s β 1 A 1 ( s , P 1 ( s ) ) d s ) | | β ( 1 α ) ε β 1 M ( α ) | | A 1 ( ε , P 1 n ( ε ) ) A 1 ( ε , P 1 ( ε ) ) | | + α β M ( α ) 0 ε s β 1 | | A 1 ( s , P 1 n ( s ) ) A 1 ( s , P 1 ( s ) ) | | d s ( β ( 1 α ) ε β 1 M ( α ) + ( α β M ( α ) 0 ε s β 1 ) d s ) | | A 1 ( s , P 1 n ( s ) ) A 1 ( s , P 1 ( s ) ) | | β ( 1 α ) + α M ( α ) | | A 1 ( s , P 1 n ( s ) ) A 1 ( s , P 1 ( s ) ) | | ( β ( 1 α ) + α M ( α ) ) n W n | | A 1 ( s , P 1 n ( s ) ) 1 ( s , P 1 ( s ) ) | | ,
where W < 1 and n . Similarly,
Ω 2 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 2 ( s , Q 1 n ( s ) ) A 2 ( s , Q 1 ( s ) ) | | , Ω 3 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 3 ( s , S 1 n ( s ) ) A 3 ( s , S 1 ( s ) ) | | , Ω 4 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 4 ( s , R 1 n ( s ) ) A 4 ( s , R 1 ( s ) ) | | , Ω 5 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 5 ( s , U 1 n ( s ) ) A 5 ( s , U 1 ( s ) ) | | , Ω 6 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 6 ( s , P 2 n ( s ) ) A 6 ( s , P 2 ( s ) ) | | , Ω 7 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 7 ( s , Q 2 n ( s ) ) A 7 ( s , Q 2 ( s ) ) | | , Ω 8 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 8 ( s , S 2 n ( s ) ) A 8 ( s , S 2 ( s ) ) | | , Ω 9 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 9 ( s , R 2 n ( s ) ) A 9 ( s , R 2 ( s ) ) | | , Ω 10 n ( ε ) = ( β ( 1 α ) + α M ( α ) ) n W n | | A 10 ( s , U 2 n ( s ) ) A 10 ( s , U 2 ( s ) ) | | .
Hence, we say Ω i n ( ε ) 0 as n for i N 1 10 and W < 1 , which completes the proof. □
Theorem 3.
The fractal-fractional age-structure smoking epidemic model (2) has a unique solution if the following holds;
( β ( 1 α ) + α M ( α ) ) ψ i 1 , i N 1 10 .
Proof. 
Let us consider the contradiction that there exists another solution of fractal fractional age-structure smoking epidemic model (2), such that
( P ¨ 1 ( ε ) , Q ¨ 1 ( ε ) , S ¨ 1 ( ε ) ,   R ¨ 1 ( ε ) , U ¨ 1 ( ε ) , P ¨ 2 ( ε ) , Q ¨ 2 ( ε ) , S ¨ 2 ( ε ) ,   R ¨ 2 ( ε ) , U ¨ 2 ( ε ) )
such that
P ¨ 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 1 ε , P ¨ 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 1 ε , P ¨ 1 ( ε ) d s , Q ¨ 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 2 ε , Q ¨ 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 2 ε , Q ¨ 1 ( ε ) d s , S ¨ 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 3 ε , S ¨ 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 3 ε , S ¨ 1 ( ε ) d s , R ¨ 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 4 ε , R ¨ 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 4 ε , R ¨ 1 ( ε ) d s , U ¨ 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 5 ε , U ¨ 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 5 ε , U ¨ 1 ( ε ) d s , P ¨ 2 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 6 ε , P ¨ 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 6 ε , P ¨ 2 ( ε ) d s , Q ¨ 2 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 7 ε , Q ¨ 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 7 ε , Q ¨ 2 ( ε ) d s , S ¨ 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 8 ε , S ¨ 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 8 ε , S ¨ 2 ( ε ) d s , R ¨ 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 9 ε , R ¨ 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 9 ε , R ¨ 2 ( ε ) d s , U ¨ 1 ( ε ) = β ( 1 α ) ε β 1 M ( α ) A 10 ε , U ¨ 2 ( ε ) + α β M ( α ) 0 ε s β 1 A 10 ε , U ¨ 2 ( ε ) d s .
Now, we take the difference of P 1 ( ε ) , P ¨ 1 ( ε ) and also the norm, and we obtain
| | P 1 ( ε ) P ¨ 1 ( ε ) | | = | | ( β ( 1 α ) ε β 1 M ( α ) A 1 ε , P 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 1 ε , P 1 ( ε ) d s ) ( β ( 1 α ) ε β 1 M ( α ) A 1 ε , P ¨ 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 1 ε , P ¨ 1 ( ε ) d s ) | | β ( 1 α ) ε β 1 M ( α ) | | A 1 ε , P 1 ( ε ) A 1 ε , P ¨ 1 ( ε ) | | + α β M ( α ) 0 ε s β 1 | | A 1 ( ε , P 1 ( ε ) A 1 ε , P ¨ 1 ( ε ) | | d s β ( 1 α ) ε β 1 M ( α ) + α β M ( α ) 0 ε s β 1 ψ 1 | | P 1 P ¨ 1 | | ( β ( 1 α ) + α M ( α ) ) ψ 1 | | P 1 P ¨ 1 | | | | P 1 P ¨ 1 | | ( β ( 1 α ) + α M ( α ) ψ 1 ) | | P 1 P ¨ 1 | | 0 ( 1 β ( 1 α ) + α M ( α ) ψ 1 ) | | P 1 P ¨ 1 | | 0 .
The above inequality is true if | | P 1 P ¨ 1 | | = 0 implies P 1 = P ¨ 1 . Similarly,
| | Q 1 ( ε ) Q ¨ 1 ( ε ) | | = | | ( β ( 1 α ) ε β 1 M ( α ) A 2 ε , Q 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 2 ε , Q 1 ( ε ) d s ) ( β ( 1 α ) ε β 1 M ( α ) A 2 ε , Q ¨ 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 2 ε , Q ¨ 1 ( ε ) d s ) | | β ( 1 α ) ε β 1 M ( α ) | | A 2 ε , Q 1 ( ε ) A 2 ε , Q ¨ 1 ( ε ) | | + α β M ( α ) 0 ε s β 1 | | A 2 ( ε , Q 1 ( ε ) A 2 ε , Q ¨ 1 ( ε ) | | d s β ( 1 α ) ε β 1 A ( α ) + α β M ( α ) 0 ε s β 1 ψ 1 | | Q 1 Q ¨ 1 | | ( β ( 1 α ) + α M ( α ) ) ψ 1 | | Q 1 Q ¨ 1 | | | | Q 1 Q ¨ 1 | | ( β ( 1 α ) + α M ( α ) ψ 1 ) | | Q 1 Q ¨ 1 | | 0 ( 1 β ( 1 α ) + α M ( α ) ψ 1 ) | | Q 1 Q ¨ 1 | | 0 .
The above inequality is true if | | Q 1 Q ¨ 1 | | = 0 implies Q 1 = Q ¨ 1 . Similarly, we can conclude that S 1 = S ¨ 1 , R 1 = R ¨ 1 , U 1 = U ¨ 1 , P 2 = P ¨ 2 , Q 2 = Q ¨ 2 S 2 = S ¨ 2 , R 2 = R ¨ 2 , U 2 = U ¨ 2 . Hence, fractal-fractional age-structure smoking epidemic model (2) has a unique solution. □

6. Hyers–Ulam (HU) Stability Results

Definition 4.
The fractal-fractional integral system (9) is HU stable if δ i > 0 and ϕ i > 0 for all i N 1 10 , such that
| P 1 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 1 ε , P 1 ( ε ) + α β M ( α ) 0 t s β 1 A 1 ε , P 1 ( ε ) d s | δ 1 , | Q 1 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 2 ε , Q 1 ( ε ) + α β M ( α ) 0 t s β 1 A 2 ε , Q 1 ( ε ) d s | δ 2 , | S 1 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 3 ε , S 1 ( ε ) + α β M ( α ) 0 t s β 1 A 3 ε , S 1 ( ε ) d s | δ 3 , | R 1 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 4 ε , R 1 ( ε ) + α β M ( α ) 0 t s β 1 A 4 ε , R 1 ( ε ) d s | δ 4 , | U 1 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 5 ε , U 1 ( ε ) + α β M ( α ) 0 t s β 1 A 5 ε , U 1 ( ε ) d s | δ 5 , | P 2 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 6 ε , P 2 ( ε ) + α β M ( α ) 0 t s β 1 A 6 ε , P 2 ( ε ) d s | δ 6 , | Q 2 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 7 ε , Q 2 ( ε ) + α β M ( α ) 0 t s β 1 A 7 ε , Q 2 ( ε ) d s | δ 7 , | S 2 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 8 ε , S 2 ( ε ) + α β M ( α ) 0 t s β 1 A 8 ε , S 2 ( ε ) d s | δ 8 , | R 2 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 9 ε , R 2 ( ε ) + α β M ( α ) 0 t s β 1 A 9 ε , R 2 ( ε ) d s | δ 9 , | U 2 ( ε ) β ( 1 α ) ε β 1 M ( α ) A 10 ε , U 2 ( ε ) + α β M ( α ) 0 t s β 1 A 10 ε , U 2 ( ε ) d s | δ 10 .
There exists an approximate solution of fractal-fractional age-structure smoking epidemic model (3), such that
( P ¨ 1 ( ε ) , Q ¨ 1 ( ε ) , S ¨ 1 ( ε ) ,   R ¨ 1 ( ε ) , U ¨ 1 ( ε ) , P ¨ 2 ( ε ) , Q ¨ 2 ( ε ) , S ¨ 2 ( ε ) ,   R ¨ 2 ( ε ε ) , U ¨ 2 ( ε ) )
such that it satisfies the given model, such that
| P 1 ( ε ) P ¨ 1 ( ε ) | = | ( β ( 1 α ) ε β 1 M ( α ) A 1 ε , P 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 1 ε , P 1 ( ε ) d s ) ( β ( 1 α ) ε β 1 M ( α ) A 1 ε , P ¨ 1 ( ε ) + α β M ( α ) 0 ε s β 1 A 1 ε , P ¨ 1 ( ε ) d s ) | = β ( 1 α ) ε β 1 M ( α ) A 1 ε , P 1 ( ε ) A 1 ε , P ¨ 1 ( ε ) | + α β M ( α ) 0 ε s β 1 | A 1 ( ε , P 1 ( ε ) A 1 ε , P ¨ 1 ( ε ) | d s β ( 1 α ) ε β 1 M ( α ) + α β M ( α ) 0 ε s β 1 ψ 1 | | P 1 P ¨ 1 | | ( β ( 1 β ) + α M ( α ) ) ψ 1 | | P 1 P ¨ 1 | | | | P 1 P ¨ 1 | | ( β ( 1 α ) + α M ( α ) ψ 1 ) | | P 1 P ¨ 1 | | 0 ( 1 β ( 1 α ) + α M ( α ) ψ 1 ) | | P 1 P ¨ 1 | | 0 .
Letting λ 1 = ( 1 β ( 1 α ) + α A ( α ) ) | | P 1 P ¨ 1 | | , ψ 1 = ζ 1 , the above inequalities become | P 1 ( ε ) P ¨ 1 ( ε ) | λ 1 ζ 1 . Similarly,
| Q 1 ( ε ) Q ¨ 1 ( ε ) | λ 2 ζ 2 , | S 1 ( ε ) S ¨ 1 ( ε ) | λ 3 ζ 3 , | R 1 ( ε ) R ¨ 1 ( ε ) | λ 4 ζ 4 , | U 1 ( ε ) U ¨ 1 ( ε ) | λ 5 ζ 5 , | P 2 ( ε ) P ¨ 2 ( ε ) | λ 6 ζ 6 , | Q 2 ( ε ) Q ¨ 2 ( ε ) | λ 7 ζ 7 , | S 2 ( ε ) S ¨ 2 ( ε ) | λ 8 ζ 8 , | R 2 ( ε ) R ¨ 2 ( ε ) | λ 9 ζ 9 , | U 2 ( ε ) U ¨ 2 ( ε ) | λ 10 ζ 10 .
Consequently, by definition fractal-fractional age-structure smoking epidemic model (2) is HU stable. Hence, we complete the proof.

7. Numerical Scheme

In this section, we present the numerical results for the fractal-fractional age-structure smoking epidemic transmission model based on the Newton polynomial. The Cauchy problem of the CF fractal-fractional derivative can be given as
C F F D ε α β A ( α ) = F ( ε , A ( ε ) ) , A ( 0 ) = A 0 .
By employing the initial condition together with operator C F I 0 α β , we convert the Cauchy problem (27) to CF fractal-fractional integral equations, as
A ( ε ) = A ( 0 ) + β ( 1 α ) t β 1 M ( α ) F ( ε , A ( ε ) ) + α β M ( α ) 0 ε s β 1 F ( s , A ( ε ) ) d s .
Taking points ε ( n + 1 ) = ( n + 1 ) h and ε n = n h ,   n = 0 , 1 , 2 , , with h being the time step, then we simply obtain a general Newton numerical scheme for our proposed model (2). For more details about the numerical analysis, see [23,26]
A n + 1 = A n + ( 1 α ) M ( α ) β ε n β 1 F ( ε ( n ) , A ( ε ( n 1 ) ) ) β ε n 1 β 1 F ( ε ( n 1 ) , A ( ε ( n 1 ) ) ) + α M ( α ) m = 2 n [ ( 4 3 β ε m 1 β 1 A ( ε m 1 , A m 1 ) h + 5 12 β ε m 2 β 1 A ( ε m 2 , A m 2 ) h + 23 12 β ε m β 1 A ( ε m , A m ) h ) ] .
For simplicity, we rewrite fractal-fractional age-structure smoking epidemic model (3) as
C F F D ε α , β P 1 ( ε ) = F 1 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β Q 1 ( ε ) = F 2 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β S 1 ( ε ) = F 3 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β U 1 ( ε ) = F 4 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β R 1 ( ε ) = F 5 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β P 2 ( ε ) = F 6 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β Q 2 ( ε ) = F 7 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β S 2 ( ε ) = F 8 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β U 2 ( ε ) = F 9 ( ε , P i , Q i , S i , R i , U i ) , C F F D ε α , β R 2 ( ε ) = F 10 ( ε , P i , Q i , S i , R i , U i ) ,
where i = 1 , 2 and A i are all of the variables, thus i = N 1 10 . Therefore, the fractal-fractional age-structure smoking epidemic model’s numerical scheme based on the Newton polynomial is as follows
A i n + 1 = A i m + ( 1 α ) M ( α ) [ β ε n β 1 F i ( ε n , P i n , Q i n , S i n , R i n , U i n ) β ε m 1 β 1 F i ( ε n 1 , P i n 1 , Q i n 1 , S i n 1 , R i n 1 , U i n 1 ) ] + α M ( α ) m = 2 n [ ( 4 3 β ε m β * 1 F i ( ε m 1 , P i m 1 , Q i m 1 , S i m 1 , R i m 1 , U i r * 1 ) h + 5 12 β ε m 2 β 1 F i ( ε m 2 , P i m 2 , Q i m 2 , S i m 2 , R i m 2 , U i m 2 ) h + 23 12 β ε m β 1 F i ( ε m , P i m , Q i m , S i m , R i m , U i m ) h ) ] .

8. Numerical Results and Discussion

The numerical results and some discussion about an approximate solution of the fractal-fractional age-structure smoking transmission model are presented in this section of the paper. We compute the model associated with the fractal dimension β ( 0 , 1 ] and the CF fractional operator α ( 0 , 1 ] using the fractal-fractional Newton polynomial for the numerical simulation. We consider some suitable initial conditions together with the parameter values, as follows: P 1 ( 0 ) = 30 ;   Q 1 ( 0 ) = 14 ;   S 1 ( 0 ) = 16 ;   R 1 ( 0 ) = 13 ;   U 1 ( 0 ) = 14 ;   P 2 ( 0 ) = 20 ;   Q 2 ( 0 ) = 10 ;   S 2 ( 0 ) = 12 ;   R 2 ( 0 ) = 14 ;   U 2 ( 0 ) = 15 ; Λ 1 = 0.088 ;   Λ 2 = 0.012 ;   μ 1 = 0.0002 ;   μ 2 = 0.040 ;   λ 1 = 0.065 ;   λ 2 = 0.035 ;   β 1 = 0.067 ;   β 2 = 0.050 ;   ν 1 = 0.00078 ;   ν 2 = 0.00022 ;   δ 1 = 0.082 ;   δ 2 = 0.071 ;   η 1 = 0.80 ;   η 2 = 0.50 .
For the F-F age-structure smoking epidemic model, the numerical results are shown in Figure 1 and Figure 2, which show how the model works. In Figure 1 and Figure 2, when the fractional order α and fractal dimension β are reduced from 1, the rate of smokers who have temporarily and permanently quit smoking increases significantly, and the asymptotic behaviors in potential smokers and smoker compartments change in both Group-1 and Group-2. This indicates that the combination of F-F operators directly affects the dynamics of the age-structure smoking epidemic model, and it also shows the hidden behaviors of smoking-related diseases. In Figure 3 and Figure 4, when the fractal-fractional order is α = β = 0.95 when the government intervention parameter is varied, we noticed that the rate at which people join the group of smokers who have potentially, occasionally, and temporarily quit significantly decreases in both Group-1 and Group-2. Moreover, this intervention has a minimal effect on smokers and those who have permanently quit in both Group-1 and Group-2.

9. Conclusions

In this study, we developed fractal-fractional two age-structure smoking epidemic models. The models are initially formulated via a classical integer order differential system. Then, in the rest phase, we extended to the CF fractal-fractional derivatives. The basic mathematical properties, including boundedness, positivity, and model invariance, were explored. By means of the fixed point theorem, we proved the existence and uniqueness sets of the solutions of the fractal-fractional two age-structure smoking epidemic models, and due to a small change in the initial conditions, we proved the stability solutions using the HU stability type. In the second part of the paper, the two new fractal-fractional age-structure smoking epidemic models were solved numerically using the newly built Newton polynomial, followed by numerical simulations. The graphical representation indicated the dynamics of the fractal-fractional operators and the unique dynamism of a single operator when the other operator is kept constant. To be able to control the smoking epidemic, governments and all stakeholders should implement feasible campaigns, taking into account all age groups. In the future, this work could be used to study different fractional operators and fractal-fractional optimal control.

Author Contributions

E.A.: Conceptualization, Methodology, Formal analysis, Writing—original draft, Writing—review & editing. A.A.: Methodology, Formal analysis, Writing—review & editing. O.J.P.: Formal analysis, Writing—original draft, Writing—review & editing. J.O.A.: Software, Validation, Data curation, Visualization, Writing- review & editing. K.O.: Project administration, Investigation, Funding acquisition, Supervision, Methodology, Formal analysis, Writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization Tobacco Report. Fact-Sheets about Tobacco. 2021. Available online: https://www.who.int/news-room/fact-sheets/detail/tobacco (accessed on 11 December 2022).
  2. World Health Organization Tobacco Report. Tobacco in the Western Pacific. 2021. Available online: https://www.who.int/westernpacific/health-topics/tobacco (accessed on 11 December 2022).
  3. Aslam, M.; Murtaza, R.; Abdeljawad, T.; Rahman, G.U.; Khan, A.; Khan, H.; Gulzar, H. A fractional order HIV/AIDS epidemic model with Mittag-Leffler kernel. Adv. Differ. Equ. 2021, 2021, 107. [Google Scholar] [CrossRef]
  4. Sher, M.; Shah, K.; Khan, Z.A.; Khan, H.; Khan, A. Computational and theoretical modeling of the transmission dynamics of novel COVID-19 under Mittag-Leffler Power Law. Alex. Eng. J. 2020, 59, 3133–3147. [Google Scholar] [CrossRef]
  5. Chukwu, W.; Nyabadza, F. A mathematical model and optimal control for Listeriosis disease from ready-to-eat food products. medRxiv 2020. [Google Scholar] [CrossRef]
  6. Osman, S.; Makinde, O.D.A. Mathematical model for co-infection of listeriosis and anthrax diseases. Int. J. Math. Math. Sci. 2018, 2018, 1725671. [Google Scholar] [CrossRef]
  7. Sardar, T.; Rana, S.; Bhattacharya, S.; Al-Khaled, K.; Chattopadhyay, J. A generic model for a single strain mosquito-transmitted disease with memory on the host and the vector. Math. Biosci. 2015, 263, 18–38. [Google Scholar] [CrossRef]
  8. Mensah, J.A.; Dontwi, I.K.; Bonyah, E. Stability analysis of Zika-malaria co-infection model for malaria endemic region. J. Adv. Math. Comput. Sci. 2018, 26, 1–22. [Google Scholar] [CrossRef]
  9. Peter, O.J.; Yusuf, A.; Ojo, M.M. A Mathematical Model Analysis of Meningitis with Treatment and Vaccination in Fractional Derivatives. Int. J. Appl. Comput. Math. 2022, 8, 117. [Google Scholar] [CrossRef]
  10. Baleanu, D.; Jajarmi, A.; Mohammad, H.; Rezapour, S. A new study on the mathematical modelling of human liver with Caputo-Fabrizio fractional derivative. Chaos Solit. Fract. 2020, 134, 109705. [Google Scholar] [CrossRef]
  11. Asamoah, J.K.K.; Okyere, E.; Yankson, E.; Opoku, A.A.; Adom-Konadu, A.; Acheampong, E.; Arthur, Y.D. Non-fractional and fractional mathematical analysis and simulations for Q fever. Chaos Solitons Fractals 2022, 156, 111821. [Google Scholar] [CrossRef]
  12. Castillo-Garsow, C.; Jordan-Salivia, G.; Rodriguez-Herrera, A. Mathematical Models for the Dynamics of Tobacco Use, Recovery, and Replase; Technical Report BU-1505-M; Cornell University: Ithaca, NY, USA, 2000. [Google Scholar]
  13. Zaman, G. Qualitative Behavior of Giving up Smoking Models. Bull. Malays. Math. Sci. Soc. 2011, 34, 403–415. [Google Scholar]
  14. Zaman, G. Optimal Campaign in the Smoking Dynamics. Comput. Math. Methods Med. 2011, 2011, 163834. [Google Scholar] [CrossRef]
  15. Abdullah, M.; Ahmad, A.; Raza, N.; Farman, M.; Ahmad, M. Approximate Solution and Analysis of Smoking Epidemic Model with Caputo Fractional Derivatives. Int. J. Appl. Comput. Math. 2018, 4, 112. [Google Scholar] [CrossRef]
  16. Khan, S.A.; Shah, K.; Zaman, G.; Jarad, F. Existence theory and numerical solutions to smoking model under Caputo-Fabrizio fractional derivative. Chaos 2019, 29, 013128. [Google Scholar] [CrossRef]
  17. Bonyah, E.; Yavuz, M.; Baleanu, D.; Kumar, S. A robust study on the listeriosis disease by adopting fractal-fractional operators. Alex. Eng. J. 2022, 61, 2016–2028. [Google Scholar] [CrossRef]
  18. Ahmad, Z.; Bonanomi, G.; di Serafino, D.; Giannino, F. Transmission dynamics and sensitivity analysis of pine wilt disease with asymptomatic carriers via fractal-fractional differential operator of Mit-tag-Leffler kernel. Appl. Numer. Math. 2023, 185, 446–465. [Google Scholar] [CrossRef]
  19. Malik, A.; Alkholief, M.; Aldakheel, F.M.; Ali Khan, A.; Ahmad, Z.; Kamal, W.; Khalil Gatasheh, M.; Alshamsan, A. Sensitivity analysis of COVID-19 with quarantine and vaccination: A fractal-fractional model. Alex. Eng. J. 2022, 61, 8859–8874. [Google Scholar] [CrossRef]
  20. Addai, E.; Zhang, L.; Asamoah, J.K.K.; Preko, A.K.; Arthur, Y.D. Fractal–fractional age-structure study of omicron SARS-CoV-2 variant transmission dynamics. Partial. Differ. Equ. Appl. Math. 2022, 6, 100455. [Google Scholar] [CrossRef]
  21. Addai, E.; Zhang, L.; Ackora-Prah, J.; Gordon, J.F.; Asamoah, J.K.K.; Essel, J.F. Fractal-fractional order dynamics and nu-merical simulations of a Zika epidemic model with insecticide-treated nets. Phys. A 2022, 603, 127809. [Google Scholar] [CrossRef]
  22. Addai, E.; Zhang, L.; Asamoah, J.K.; Essel, J.F. A fractional order age-specific smoke epidemic model. Appl. Math. Model. 2023, 119, 99–118. [Google Scholar] [CrossRef]
  23. Atangana, A.; Owolabi, K.M. New numerical approach for fractional differential equations. Math. Model Nat. Phenom 2018, 13, 3. [Google Scholar] [CrossRef]
  24. Ahmed, E.; El-Sayed, A.M.A.; El-Saka, H.A. Equilibrium points, stability and numerical solutions of fractional-order predator–prey and rabies models. J. Math. Anal. Appl. 2007, 325, 542–553. [Google Scholar] [CrossRef]
  25. Ritchie, H.; Roser, M. “Smoking”. 2013. Available online: https://ourworldindata.org/smoking (accessed on 11 December 2022).
  26. Atangana, A.; Araz, S.I. New Numerical Scheme with Newton Polynomial: Theory, Methods, and Applications; Academic Press: Cambridge, MA, USA, 2021. [Google Scholar]
Figure 1. Numerical smoking epidemic trajectory under the different fractal dimension β and different fractional order α for Group-1. (a) Dynamics of potential smokers. (b) Dynamics of occasional smokers. (c) Dynamics of the smokers class. (d) Dynamics of the smokers who temporarily quit. (e) Dynamics of the smokers who permanently quit.
Figure 1. Numerical smoking epidemic trajectory under the different fractal dimension β and different fractional order α for Group-1. (a) Dynamics of potential smokers. (b) Dynamics of occasional smokers. (c) Dynamics of the smokers class. (d) Dynamics of the smokers who temporarily quit. (e) Dynamics of the smokers who permanently quit.
Fractalfract 07 00370 g001aFractalfract 07 00370 g001b
Figure 2. Numerical smoking epidemic trajectory under the different fractal dimension β and different fractional order α for Group-2. (a) Dynamics of potential smokers. (b) Dynamics of occasional smokers. (c) Dynamics of the smokers class. (d) Dynamics of the smokers who temporarily quit. (e) Dynamics of the smokers who permanently quit.
Figure 2. Numerical smoking epidemic trajectory under the different fractal dimension β and different fractional order α for Group-2. (a) Dynamics of potential smokers. (b) Dynamics of occasional smokers. (c) Dynamics of the smokers class. (d) Dynamics of the smokers who temporarily quit. (e) Dynamics of the smokers who permanently quit.
Fractalfract 07 00370 g002aFractalfract 07 00370 g002b
Figure 3. Dynamics of the smoking epidemic under the CF fractal-fractional operator ( α = β = 0.95 ) when the rate of government intervention varied for Group-1. (a) Dynamics of potential smokers. (b) Dynamics of occasional smokers. (c) Dynamics of the smokers class. (d) Dynamics of the smokers who temporarily quit. (e) Dynamics of the smokers who permanently quit.
Figure 3. Dynamics of the smoking epidemic under the CF fractal-fractional operator ( α = β = 0.95 ) when the rate of government intervention varied for Group-1. (a) Dynamics of potential smokers. (b) Dynamics of occasional smokers. (c) Dynamics of the smokers class. (d) Dynamics of the smokers who temporarily quit. (e) Dynamics of the smokers who permanently quit.
Fractalfract 07 00370 g003aFractalfract 07 00370 g003b
Figure 4. Dynamics of the smoking epidemic under the CF fractal-fractional operator ( α = β = 0.95 ) when the rate of government intervention varied for Group-2. (a) Dynamics of potential smokers. (b) Dynamics of occasional smokers. (c) Dynamics of the smokers class. (d) Dynamics of the smokers who temporarily quit. (e) Dynamics of the smokers who permanently quit.
Figure 4. Dynamics of the smoking epidemic under the CF fractal-fractional operator ( α = β = 0.95 ) when the rate of government intervention varied for Group-2. (a) Dynamics of potential smokers. (b) Dynamics of occasional smokers. (c) Dynamics of the smokers class. (d) Dynamics of the smokers who temporarily quit. (e) Dynamics of the smokers who permanently quit.
Fractalfract 07 00370 g004aFractalfract 07 00370 g004b
Table 1. Description of the variables and parameters.
Table 1. Description of the variables and parameters.
  VariablesDescription
N 1 , N 2 Total population size for Group 1 and Group 2
P 1 , P 2 Potential smokers for Group 1 and Group 2
Q 1 , Q 2 Occasional smokers for Group 1 and Group 2
S 1 , S 2 Smokers for Group 1 and Group 2
R 1 , R 2 Smokers who temporarily quit smoking for Group 1 and Group 2
U 1 , U 2 Smokers who permanently quit smoking for Group 1 and Group 2
  ParametersDescription
Λ 1 Recruitment rate of Group 1
Λ 2 Recruitment rate of Group 2
β 1 Effective contact rate between S 1 and P 1
β 2 Effective contact rate between S 2 and P 2
μ 1 , μ 2 Natural death rate for Group 1 and Group 2
λ 1 , λ 2 The rate at which occasional smokers become regular smokers for Group 1 and Group 2,
η 1 , η 2 The rate of quitting smoking for both Group 1 and Group 2
ν 1 , ν 2 The contact rate between smokers and temporary quitters who revert back to smoking
( 1 δ 1 ) The fraction of smokers who temporary quit smoking (at the rate η 1 )
( 1 δ 2 ) The fraction of smokers who temporary quit smoking (at the rate η 2 )
δ 1 , δ 2 The remaining fraction of smokers who permanently quit smoking
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Addai, E.; Adeniji, A.; Peter, O.J.; Agbaje, J.O.; Oshinubi, K. Dynamics of Age-Structure Smoking Models with Government Intervention Coverage under Fractal-Fractional Order Derivatives. Fractal Fract. 2023, 7, 370. https://doi.org/10.3390/fractalfract7050370

AMA Style

Addai E, Adeniji A, Peter OJ, Agbaje JO, Oshinubi K. Dynamics of Age-Structure Smoking Models with Government Intervention Coverage under Fractal-Fractional Order Derivatives. Fractal and Fractional. 2023; 7(5):370. https://doi.org/10.3390/fractalfract7050370

Chicago/Turabian Style

Addai, Emmanuel, Adejimi Adeniji, Olumuyiwa J. Peter, Janet O. Agbaje, and Kayode Oshinubi. 2023. "Dynamics of Age-Structure Smoking Models with Government Intervention Coverage under Fractal-Fractional Order Derivatives" Fractal and Fractional 7, no. 5: 370. https://doi.org/10.3390/fractalfract7050370

APA Style

Addai, E., Adeniji, A., Peter, O. J., Agbaje, J. O., & Oshinubi, K. (2023). Dynamics of Age-Structure Smoking Models with Government Intervention Coverage under Fractal-Fractional Order Derivatives. Fractal and Fractional, 7(5), 370. https://doi.org/10.3390/fractalfract7050370

Article Metrics

Back to TopTop