Abstract

This manuscript investigates fractal-fractional order smoking models with relapse and harmonic mean type incidence rate under the Caputo derivative. We derive the existence and unique results about the solution for the considered model via fixed point theory. For the stability of the considered system, Ulam-Hyers (UH) approach is used. We compute the numerical solution by using fractional Adams-Bashforth method. For the simulation of the model, we consider different values of fractional order and fractal dimension by using some real values of the parameters. The proposed scheme is used to simulate the available data for some smoking community including potential, light, and quit smokers. Various graphical presentations are given to understand the dynamics of the model at various fractional orders.

1. Introduction

The first biological model that describes the dynamics of infectious disease was presented in 1927. Later on, scientists and researchers started to investigate different properties of the models such as the spreading behavior and trends of the diseases by studying the various aspects [14]. They have formulated several models for different diseases like pine wilt, HIV, viral disease including leishmania, TB, and COVID-19 [512].

Smoking is also similar to infectious diseases by spreading its behavior in the population. The ratio of diseases due to smoking is increasing day by day. Castillo-Garsow et al. [13] formulated for the first time a simple giving up smoking model with known spreading behavior of smoking in the community. The same authors modified and extended the work by adding another class of light smoking. The authors [14] focused on the control strategy of smoking epidemic by choosing optimal campaigns. Furthermore, some of the smokers may relapse because they may have frequent contacts with smokers, whereas some of them may cease smoking permanently. Rahman et al. [15] have been worked on a smoking model and included the relapse terms for the quit smokers.

The abovementioned models have been investigated under ordinary derivatives. During the last twenty years, fractional calculus (FC) has gained more interest from the researcher and been used in different fields of sciences. Mathematical models along with fractional differential equations (FDEs) have been proved for several smoking models. Compared with integer-order model, fractional-order models have better fitting degree with different experimental results in signal processing, mathematical biology and engineering [16, 1719]. In this regard, Mahdy et al. [20] found the approximate solution for a smoking model by utilizing the Sumudu transform with Caputo derivative. Sing et al. [21] has been introduced a giving up smoking dynamic fractional model with nonsingular kernel. Khan et al. [22] have been studied a biological model of smoking type with some iterative method. Mohamed et al. [23] used reduced differential transform method to solve the nonlinear smoking fractional-order model. Alrabaiah et al. [24] have been applied Adams-Bashforth-Moulton method to investigate the tobacco smoking fractional model order containing snuffing class. Therefore, for the past periods, to develop the real phenomena for a better degree of precision and accuracy, FDEs have been utilized very well. Many researchers have utilized several methods for studying the theoretical investigation of fractional-order mathematical models, (for instance see [2529]). For further detail, see [3035]. Adomian in 1980 introduced a useful decomposition method for the solution of nonlinear systems analytically. Later on, the abovementioned method has been slowly enforced as an actual tool for consideration semianalytical or estimated results to several systems of applied sciences. Mathematical models have been examined widely using the Homotopy method, decomposition method along with integral transforms, and difference methods, for details, see [30, 31]. Recently, many methods have been utilized to handle problems of fractional order (see details in [3638]).

Keeping in mind that derivative of noninteger can be defined in several ways. The first definition of fractional derivative was given by Riemann-Liouville. Later on in 1967, Caputo gave his own definition which has been increasingly used. The mentioned both definitions include singular kernels which often cause problem in numerical investigations. To overcome these difficulties, recently, Caputo and Fabrizio [39] have introduced a new definition. The said definition contains exponential function instead of singular kernel. In subsequent years, the said definition has been further generalized by Atangana and Baleanu [40] by replacing exponential function on Mittag-Leffler one. This fact has been proved that the concerned derivative also has interesting features (see [4146]).

Recently, the area involves fractal-fractional derivative has got much attention (see [47, 4951]). Motivated from the above work and from, we consider the model presented in [48] to fractal-fractional (FF) order in sense of Caputo operator which has various advantages. This model consists of four compartments, namely, people vulnerable to smoking , light smokers , smoker class , and quit smokers . This work also includes theoretical, practical analytical, and numerical results of smoking models with relapse and harmonic mean type incidence rate. Our considered model under Caputo operator for fractional-order and fractal dimension is as follows: with initial conditions where is the transmission rate that the potential smoker contact with the chain smoker, is the relapse rate, is the recruitment rate, is the natural death rate, and is the death rate induced by smoking. Also, is the conversion rate from light to chain smoker class. In same line, is the chain smokers rate when they quit smoking. We also discuss some stability results devoted to UH type. The mentioned stability has been recently investigated for various problems of FDEs (see [5557]).

The rest of the paper we organized is as follows: Section 2 is related to basic definitions and theorems. By using fixed point theorems, we show some suitable results for the uniqueness and existence in Section 3. With the help of famous AB technique, we find the numerical solution of the considered system in Section 4. Using the AB technique, we also perform the numerical simulation by using Matlab for getting the graphical representation for our analytic and briefly discuss the obtained results. Finally, we conclude our work in Section 5.

2. Basic Results

Definition 1 (see [47]). Let on be a continuous and differentiable function with order , then the FF order derivative can be defined as along-with , , where and .

Definition 2 (see [47]). Let be continuous on then the FF order integral of with order is defined as

Definition 3. The system (1) is UH stable if any real number such that for every and all the solutions , where , the inequality can be defined as is the unique solution for the considered model (1) such that Note: let us define a Banach space For the qualitative analysis , where with norm: .

3. Theoretical Results of Model (1)

Here in this section, we will investigate the model (1) for existence. Since the given integral is differentiable, so we can express the RHS of the model (1) as

In view of (7) and for , the proposed model may be written in the following form by changing with and using the integral of Riemann-Liouville, the solution of (8) will be where

Now, if we transform (1) to fixed point problem and let the operator can be defined by

To find the existence results of the considered model, we use the following theorem [54].

Theorem 4. If the operator be completely continuous and the set is bounded, then, the operator has at least one fixed point in .

Theorem 5. Suppose the operator is a continuous operator. Then, is compact.

Proof. First, we will show that which is defined in (12) is continuous. Consider is a bounded set in , then with , for all . Any , we have Hence, (14) implies that is uniformly bounded, where Beta function can be written as . Further, for equi-continuity of the operator , for any and , we obtain Hence, is equi-continuous and then the operator is bounded and continuous as well, therefore, by Arzelá-Ascoli theorem, the operator is relatively compact and so completely continuous. Furthermore, we use the following hypothesis:
(C) There exist constants such that, for each , , we have For existence uniqueness, we use fixed point approach as given in [54].

Theorem 6. Applying the hypothesis (C) and if , then, the model (1) has a unique solution if

Proof. Assume , such that We prove that , where and , we have Suppose the operator is defined in (12). Using the assumption and for every , , we obtain By this, is contraction by using (20). Therefore, equation (10) has one solution and so our model (1) has unique solution.
Now, we have to develop UH stability for the considered system (1), taking depending on the solution with . Then (i)(ii)

Lemma 7. The solution of perturbed equation satisfies the given relation

Theorem 8. With the assumption (C) and (22), the solution of the integral equation (1) is UH stable. Hence, the analytical results of the considered system are UH stable if , where is given in (17).

Proof. Suppose that be a unique solution and be any solution of (10), then using fractalfractional integration as in an equation (2), we have Which we have From (24), we can write as Thus, from the (25), we conclude that the solution of (10) is UH stable and therefore the proposed model (1) solution is UH stable.

4. Numerical Scheme

In this part of the paper, we are constructing the numerical algorithm for the considered model to perform numerical simulation. Here, for numerical method, the construction of equation (10) of the considered model goes to the following form

Now, we are presenting the numerical solution to the (26) and using the new approach . The first equation of the above system becomes

We obtained the approximate integral from the above equation as

Within the infinite interval in term of Lagrange interpolation polynomials the function along with , such that

putting (29) into (28), then, we can write (28) as

Simplifying the right side integrals of (30), we obtain the numerical iterative results for the class in (1) by using the FF derivatives in the Caputo form as:

Similarly, the remaining terms can be written as

4.1. Graphical Representations

In this section, we provide the numerical solution of our proposed model (1) using different values of parameters given in Table 1 for verification of the obtained scheme. We have taken two different sets of initial values of all the compartments in problem (1) for two different fractal dimension and fractional order .

Figures 1(a) and 1(b) show the dynamical behavior of potential smoker population at various fractal dimension and fractional order at two different initial values. On different six fractal-fractional values, the class increases and becomes stable which converges to the same point having two initial values. The increase occurs quickly at high order and slowly at low order and converges to the integer order as we increase the fractional order.

Figures 2(a) and 2(b) are the representation of the dynamical behavior of light smokers at different fractal dimension and fractional order of δ at two different initial values. On different six fractal-fractional values, the potential smoker class becomes stable and converges to the same point having two initial values.

Figures 3(a) and 3(b) show the dynamical behavior of smokers at various fractal dimension and fractional order at two different initial values. At six different fractional values, the class decreases and becomes stable which converges to the same converging point for two different initial approximations.

Figures 4(a) and 4(b) show the dynamical behavior of quit smoker at various fractal dimension and fractional order at two different starting values. On six different fractional values, the class declines quickly nearly at all fractional orders but then becomes stable which converges to the same point having two initial guesses.

5. Conclusion

In this manuscript, we have analyzed a giving up smoking model under the concept of fractal-fractional order derivative in Caputo sense. The considered model has been investigated for some theatrical analysis including existence theory and stability results. In this regard, sufficient results have been established for existence and uniqueness of solution by using Banach-contraction and Schauder’s theorems of nonlinear functional analysis. The Ulam-Hyers stability analysis has been developed by using the usual nonlinear analysis tools. Further, we have used fractional Adam Bash-forth method and developed an algorithm to compute numerical results. We have used various values of fractal dimensions and fractional orders to present the results graphically. From graphical presentation, one can observe that fractal and fractional calculus have the ability to present the dynamics of real-world problems more comprehensively.

Data Availability

Data sharing is not applicable to this article as no data sets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

All authors read and approved the final manuscript.

Acknowledgments

This research was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University through the Fast-track Research Funding Program to support publication in the top journal (Grant no. 42-FTTJ-70).