In this manuscript, we solve the ordinary model of nonlinear smoking mathematically by using the second kind of shifted Chebyshev polynomials. The stability of the equilibrium point is calculated. The schematic of the model illustrates our proposition. We discuss the optimal control of this model, and formularize the optimal control smoking work through the necessary optimality cases. A numerical technique for the simulation of the control problem is adopted. Moreover, a numerical method is presented, and its stability analysis discussed. Numerical simulation then demonstrates our idea. Optimal control for the model is further discussed by clarifying the optimal control through drawing before and after control. Fractional request differential equations (FDEs) are usually used to display frameworks that have memory and exist in a few thermoelasticity models and organic standards. FDEs show the realistic biphasic decline of infection of diseases but at a slower rate. FDEs are more suitable than integer order ones in modeling complex systems, such as biological systems.
Shifted second kind chebyshevsmoking modelstabilityhamiltonianlagrange multipliersoptimal controlIntroduction
Infectious diseases have a tremendous influence on human life. Every year billions of people suffer from or die due to various infectious diseases. Mathematical modeling is of considerable importance in epidemiology because it may provide understanding of the underlying mechanisms that influence the spread of a disease and thus may be used to offer control strategies. Many scientists explored the detection of illnesses and pests [1–72], wherein the idea that mathematical modeling helps in discovering diseases’ spread [17]. In 1766, Islamic scholars presented for the first time what was considered the beginning of modern epidemiology discovery to be since followed. Model smoking is one of the most famous and important model commonly used by researchers [8–14], [18–21], [22–36]. According to the World Health Organization [14], the global tobacco pandemic killed several million people. The cost of health care is increasing and its budgets are decreasing as developing economies are shrinking [17]. Smoking negatively affects overall health as well [1–4]. According to [1–3], more than five million individuals are killed every year from smoking tobacco, as such it is the main cause of preventable deaths. Serious illnesses caused by it include asthma, lung cancer, heart disease, and mouth ulcers. In [8–9], we discussed the smoking dynamics of the integer {of the integer what}. Moreover, [13] analyzed the above using the Homotopy analysis method. Four kinds of polynomials, Chebyshev [35], [57–59]: Chebyshev's model includes four kinds of polynomials. The outcomes of each and their applications can be found in several books. Moreover, many researchers use several of these polynomials in a single model [33,35]. However, the shifted Chebyshev polynomial of the second kind Zn∗(y) still needs further study [35]. The model in [56] has proven the equilibrium points.
In the present manuscript, we discuss the analytical solutions for alpha = 1, for the following order of smoking models [4,5,56]: DP=a(1−P)−bPS,DL=−aL+bPL−cLS,DS=−(a+d)S+cLS+fQ,DQ=−(a+f)Q+d(1−e)S,DR=−aR+edS.
The initial conditions are as follows:P0=0.55,L0=0.2,S0=0.17272,R0=0.01028,f=0.25,e=0.4,d=0.2,c=0.3,a=0.04,b=0.23,D=ddt.
In Fig. 1, a schematic proposition graph is a tool that is applied to display the interrelation between model states and enables the application of graph-theoretic tools to discover novel features of the model.
Schematic of the proposed model
Tab. 1 describes the parameters of the model.
Parameter definitions
T(t)
Total population with time
P(t)
Probable smokers
L(t)
Occasional smokers
S(t)
Heavy smokers
Q(t)
Temporary quitters
R(t)
Smokers who quit permanently
a
Rate of natural death
b
Contact rate between probable and occasional smokers
c
Contact rate between occasional and smokers heavy
d
Average number of [x] who quit smoking
e
The fraction of remaining smokers
f
Contact rate between smokers and temporary quitters
The organization of the remainder of this paper is as follows. The stability of equilibrium points is studied in Section 2. Optimal control for the smoking model is discussed in Section 3. In Section 4, the second kind of Chebyshev polynomials and their properties are given. In Section 5, the numerical implementation is given. Finally, a conclusion is presented in Section 6.
Stability of Equilibrium Points
For the stability behavior of this model at E0=(1, 0, 0, 0, 0), we use transformations [36]:
U = 1-P, P = 1-U, L=L, S=S, Q=Q, R=R,
U*= −P*, L*=L*, S*=S*, Q*=Q*, R*=R*.
It follows thatdUdt=−aU+bS−bSU,dLdt=−aL+bL−bLU−cLS,dSdt=−aS−dS+cLS+fQ,dQdt=−aQ−fQ+ds−deS,dRdt=−aR+edS.
where, U*=1-P*, U**=1-P**, U***=1-P*** and U****=1-P****.
Through Taylor approximation, the linearized form of the model is[X′Y′Z′V′W′]=B[XYZVW],
where X=U-0, Y=L-0, Z=S-0, V=Q-0, W=R-0, andB=[-a 0 b 0 00 -a+b 0 0 00 -a-d f 0 0 0 0 d-de -a-f 00 0 ed 0 -a].
The stability of (2.1) is the same as the linearized stability (2.3). The stability of (2.3) depends on eigenvalues B. The point P1 is asymptotically locally stable if eigenvalues are real negative parts; conversely, P1 is unstable if at least one eigenvalue of B is a nonnegative real part. For details on stability see [56,63].
Optimal Control for Smoking Model
Consider the state presented (1.2), in R5, with control functions admissible [28–32]:Ω={(uV(.),uR(.))∈(L∞(0,Tf)2)|0≤uV(.),uR(.)≤1,∀t∈[0,Tf]},
where Tf is the final time, and uV(.) and uR(.) are functions controls.
The objective function is defined asJ(uV(.),uR(.))=∫0Tf[AL(t)+BuV2(t)+CuR2(t)]dt,
where A, B, and C represent the number of occasional smokers, the rate of contact between smokers and quitters, and temporarily who regain support to smoking, and rate of smoking quitting.
We minimize the objective function as follows [28–32]:J(uV,uR)=∫0Tfηdt,
which is subjected to the constraintDP=ξ1,DL=ξ2,DS=ξ3,DQ=ξ4,DR=ξ5,ξk=ξ,k=1,2,3,4,5.
The following initial conditions are satisfied:P=P0,L=L0,S=S0,Q=Q0,R=R0.
OCP is defined, and we consider the following modified objective (cost) function:J¯=∫0Tf[H−∑k=15λkξk]dt,
where the Hamiltonian and control smoking objective functions are defined asH=η+∑k=15λkξk,H=AL+BuE2+CuM2+λ1∗(a(1−P)−bPS)+λ2∗(−aL+bPL−cLS)+λ3∗(−(a+d)S+cLS+fQ)+λ4∗(−(a+f)Q+d(1−e)S)+λ5∗(−aR+edS).
From (3.5) and (3.7), the conditions, necessary and sufficient for OPC areDλ1∗=∂H∂P,Dλ2∗=∂H∂L,Dλ3∗=∂H∂S,Dλ4∗=∂H∂Q,Dλ5∗=∂H∂R,
Therefore, the function controls uV∗,uR∗ are given byuV∗=Q∗[λ4∗−λ3∗]2B,
uR∗=S∗[λ3∗−λ4∗−eλ5∗]2C,
uV∗=min{1,max{0,Q∗[λ4∗−λ3∗]2B}},
uR∗=min{1,max{0,S∗[λ3∗−λ4∗−eλ5∗]2C}}.
For more on optimal controls for solving models, see [29–32,44,47].
Properties of the Second Kind of Chebyshev PolynomialsSecond Kind of Chebyshev Polynomials
Chebyshev polynomials Zn(y) of the second kind are rectangular polynomials of stage n in x presented on the interval [−1, 1] [34,35,57–59]:
Zn(y)=sin(n+1)θsinθ, where y = cosθ and θ ∈ [0, π].
Polynomials have rectangular with rating to the products indoor 〈Zn(y),Zm(y)〉=∫−111−y2Zn(y)Zm(y)dy={0,n≠m,π2,n=m,
where 1−y2 is a weight function.
Zn(y) can be generated by recurring relations
Zn(y) = 2y Zn−1(y) − Zn−2(y), n = 2, 3, …, n, with Z0(y) = 1. Z1(y) = 2y.
The analytical form Zn(y) of stage n is given by [34]:Zn(y)=∑i=0⌈n2⌉(−1)i(in−i)(2y)n−2i=∑i=0⌈π2⌉(−1)i2n−2iΓ(n−i+1)yn−2iΓ(i+1)Γ(n−2i+1),n>0.Here, ⌈n2⌉ is part integral n/2.
Second Kind of Shifted Chebyshev Polynomials
On the interval y ∈ [0, 1], Zn∗(y)is defined via substituting a change variable z = 2y − 1. Hence, Zn∗(y)is defined as [33–35] Zn∗(y)=Zn(2y−1), and satisfies the following relation: 2yZn−1∗(y2)=Z2n−1(y).
The following inner product is orthogonal on the interval [0, 1]:〈Zn∗(y),Zm∗(y)〉=∫01y−y2Zn∗(y)Zm∗(y)dy={0,n≠m,π8,n=m,
with weight function y−y2.Zn∗(y) can be produced via recurring relations
Zn∗(y)=2(2y−1)Zn−1∗(y)−Zn−2∗(y),n=2,3,…,n, with Z0∗(y)=1,Z0∗(y)=4y−2.
The following formula represents the analyticalZr∗(y)=∑j=0r(−1)j22r−2jΓ(2r−j+2)yr−jΓ(j+1)Γ(2r−2j+2),r>0,
The solution of this model can be written as Z∗(y).
Let g(y) be a square integral in [0, 1], and then the second kind of shifted Chebyshev polynomials can be represented as follows:g(y)=∑i=0∞aiZi∗(y),
The coefficients aj, j = 0, 1, …, are expressed asai=2π∫−11g(y+12)1−y2Zi(y)dy,
orai=8π∫01g(y)y−y2Zi∗(y)dy,
We use only the (r + 1) terms. Then,gr(y)=∑i=0raiZi∗(y).
Using the practice to construct an integral collocation style then z(y)as follows [33]:dkz(y)dxk≅∑n=0ranZn∗(y)=∑n=0ranwnk(y).
By integrating (4.8), we can obtain the following:dk−1z(y)dyn−1≅∑n=0ranwnk−1(y)+c1
From (4.8) and (4.12), we then havewn(k)(y)=∑i=0n(−1)i22n−2iΓ(2n−i+2)yn−iΓ(i+1)Γ(2n−2i+2),wn(k−1)(y)=∫wn(k)(y)dy=∑i=0n(−1)i22n−2iΓ(2n−i+2)yn−i+1Γ(i+1)Γ(2n−2i+2)(n−i+1),wn(k−2)(y)=∫wn(k−1)(y)dx=∑i=0n(−1)i22n−2iΓ(2n−i+2)yn−i+2Γ(i+1)Γ(2n−2i+2)(n−i+1)(n−i+2),wn(0)(y)=∫wn(1)(y)dx=∑i=0n(−1)i22n−2iΓ(2n−i+2)yi+kΓ(i+1)Γ(2n−2i+2)…(i+k−1)(i+k).
We presently register Eqs. (4.9)–(4.13) at (r + 1) points yp, p = 0, 1, …, r as follows:dkz(yp)dyk=Ω(k)E^,dk−1z(yp)dyk−1=Ω(k−1)E^,…,dz(yp)dy=Ω(1)E^,z(yp)=Ω(0)E^,where E^=[a0,a1,...,am,c1,c2,...cn]T, and Ω(r), Ω(r−1), …, Ω(0) are combined matrices.
Style Integral Collocation for Resolving the Nonlinear Smoking Model
In this section, we present the model's implementation through the following steps:
We first approximate the function using Eqs. (4.8)–(4.14) with m = 5 , as follows:dPdt≅∑n=05Anwn(1)(t),dLdt≅∑n=05Bnwn(1)(t),dSdt≅∑n=05Cnwn(1)(t),dQdt≅∑n=05Dnwn(1)(t),dRdt≅∑n=05Enwn(1)(t),P(t)≅∑n=05Anwn(0)(t)+c1,L(t)≅∑n=05Bnwn(0)(t)+c2,S(t)≅∑n=05Cnwn(0)(t)+c3,Q(t)≅∑n=05Qnwn(0)(t)+c4,R(t)≅∑n=05Enwn(0)(t)+c5,
where wn(0)(t) and wn(1)(t) are expressed aswn(0)(t)=∑i=0n(−1)i22n−2iΓ(2n−i+2)ti+kΓ(i+1)Γ(2n−2i+2)…(i+k−1)(i+k),wn(1)(t)=∑i=0n(−1)i22n−2iΓ(2n−i+2)tn−iΓ(i+1)Γ(2n−2i+2).
Then, the nonlinear smoking model (1.1) is transformed to∑n=05Anwn(1)(t)=a(1−∑n=05Anwn(0)(t)+c1)−b(∑n=05Anwn(0)(t)+c1)(∑n=05Bnwn(0)(t)+c2),∑n=05Bnwn(1)(t)=−a(∑n=05Bnwn(0)(t)+c2)+b(∑n=05Anwn(0)(t)+c1)(∑n=05Bnwn(0)(t)+c2)−c(∑n=05Bnwn(0)(t)+c2)(∑n=05Cnwn(0)(t)+c3),∑n=05Cnwn(1)(t)=−(a+d)(∑n=05Cnwn(0)(t)+c3)+c(∑n=05Bnwn(0)(t)+c2)(∑n=05Cnwn(0)(t)+c3)+f(∑n=05Dnwn(0)(t)+c4),∑n=05Dnwn(1)(t)=−(a+f)(∑n=05Dnwn(0)(t)+c4)+(d−de)(∑n=05Cnwn(0)(t)+c3),∑n=05Enwn(1)(t)=−a(∑n=05Enwn(0)(t)+c5)+ed(∑n=05Cnwn(0)(t)+c3).
Now, we collocate Eq. (5.2) at (r + 1 = 6) points ty, y = 0 − 5; as follows:∑n=05Anwn(1)(ty)=a(1−∑n=05Anwn(0)(ty)+c1)−b(∑n=05Anwn(0)(ty)+c1)(∑n=05Bnwn(0)(ty)+c2),∑n=05Bnwn(1)(ty)=−a(∑n=05Bnwn(0)(ty)+c2)+b(∑n=05Anwn(0)(ty)+c1)(∑n=05Bnwn(0)(ty)+c2)−c(∑n=05Bnwn(0)(ty)+c2)(∑n=05Cnwn(0)(ty)+c3),∑n=05Cnwn(1)(tp)=−(a+d)(∑n=05Cnwn(0)(tp)+c3)+c(∑n=05Bnwn(0)(tp)+c2)(∑n=05Cnwn(0)(tp)+c3)+f(∑n=05Dnwn(0)(tp)+c4),∑n=05Dnwn(1)(tp)=−(a+f)(∑n=05Dnwn(0)(tp)+c4)+(d−de)(∑n=05Cnwn(0)(tp)+c3),∑n=05Enwn(1)(tp)=−a(∑n=05Enwn(0)(tp)+c5)+ed(∑n=05Cnwn(0)(tp)+c3).
The roots of shifted Chebyshev polynomial Z6∗(t).
(ii) By putting the initial Eqs. (1.2) into (5.1), we can obtain five equations.
Eqs. (5.1) and (5.3) obtained in step (ii) represent nonlinear algebraic system equations.
(iii) We use Newton's iteration to resolve the system and solve for the unknowns.
Figs. 2–4 show the nonlinear smoking model's behavior before control.
The approximate solution of variables about (ICSM) at r = 5 before control
The relationship between L(t) and Q(t) before control. The relationship between L(t) and R(t) before control
Relationship between L(t) and S(t) before control. Relationship between P(t) and L(t) before control
Figs. 5–8, show the show the nonlinear smoking model's behavior after control.
Approximate solution of variables about (ICSM) at r = 5 after control
Relationship between L(t) and Q(t) after control. Relationship between L(t) and R(t) after control
Relationship on L(t) and S(t) after control. Relationship on P(t) and L(t) after control
Relationship between P(t) and Q(t) after control. Relationship between P(t) and R(t) after control
Conclusions
In this manuscript, a mathematical nonlinear smoking model is studied. The optimal control of this smoking model is discussed. The stability of the equilibrium point is calculated, and a schematic of the proposed model is presented. Moreover, the integral collocation style is used to obtain the approximate solutions of the model. ICM using the shifted Chebyshev polynomials of the second kind is a new technique for solving these problems. Under the application with the necessary optimality conditions, we have studied the problem control with numerical techniques for the simulation. Moreover, a numerical method and its stability is discussed.
The authors are grateful for Taif University. Taif University researchers support project number (TURSP-2020/160), Taif University, Taif, Saudi Arabia.
Funding Statement: This manuscript has been funded through “Taif University's Researchers Supporting Project Number (TURSP-2020/160), Taif University, Taif, Saudi Arabia”.
Conflicts of Interest: The authors declare that there are no conflicts of interest.
ReferencesA. H.Mokdad, J. S.Marks, D. F.Stroup and J. L.Gerberding, “Actual causes of death in the United States,” , vol. 291, pp. 1238–1245, 2004.O. K.Ham, “Stages and processes of smoking cessation among adolescents,” , vol. 29, pp. 301–315, 2007.V. S.Erturk. G.Zaman and S.Momani, “A numericl analytic method for approximating a giving up smoking model containing fractional derivatives,” , vol. 64, no. 2, pp. 3068–3074, 2012.J.Singh, D.Kumar, M.Al-Qurashi and D.Baleanu, “A new fractional model for giving up smoking dynamics,” , vol. 2017, pp. 1–16, 2017.F.Haq, K.Shah. G.Rahman and M.Shahzad, “Numerical solution of fractional order smoking model via laplace adomian decomposition method,” , vol. 57, no. 2, pp. 1061–1069, 2018.A. M. A.El-Sayed and S. M.Salman, “On a discretization process of fractional order riccati's differential equation,” , vol. 4, pp. 251–259, 2013.A. A.Elsadany and A. E.Matouk, “Dynamical behaviors of fractional-order lotka-voltera predator-prey model and its discretization,” , vol. 49, pp. 49, pp. 269–283, 2015.M.Khalid, F. S.Khan and A.Iqbal, “Perturbation iteration algorithm to solve fractional giving up smoking mathematical model,” , vol. 142, no. 9, pp. 1–6, 2016.G.Zaman, “Optimal campaign in the smoking dynamics,” , vol. 2011, pp. 1–9, 2011.G.Zaman, “Qualitative behavior of giving up smoking models,” , vol. 34, no. 2, pp. 403–415, 2011.J. L.Lubin and N. E.Cporaso, “Cigarette smoking and lung cancer: Modeling total exposure and intensity cancer epidemiology, “, vol. 15, no. 3, pp. 517–523, 2006.Y. A.Amer, A. M. S.Mahdy and E. S. M.Youssef, “Solving systems of fractional differential equations using sumudu transform method,” , vol. 7, no. 2, pp. 1–15, 2017.A.Zeb., I.Chohan and G.Zaman, “The homotopy analysis method for approximating of giving up smoking model in fractional order,” , vol. 3, pp. 914–919, 2012.Z.Alkhudhari, S.Al-Sheikh and S.Al-Tuwairqi, “Global dynamics of mathematical model on smoking,” , vol. 2014, pp. 914–919, 2014.A. A. M.Arafa, S. Z.Rida and H. M.Ali, “Generalized mittag-leffler function method for solving lorenz system,” , vol. 3, pp. 105–111, 2013.A. M. S.Mahdy, M. S.Mohamed, K.Lotfy, M.Kadry and A.El-Bary, “Plasma-elastic-thermal propagation of a rotator semiconductor magneto-electric medium during microtemperature and photothermal excitation processes subjected to mechanical ramp type,” , pp. 1–21, 2020.S. Z.Rida and A. A. M.Arafa, “New method for solving linear fractional differential equations,” , vol. 2011, pp. 1–8, 2011.I.Podlubny, “,” San Diego, CA: Academic Press, 1999.A. M. S.Mahdy, K. A.Gepreel, Kh.Lotfy and A. A.El-Bary, “A numerical method for solving the rubella ailment disease model,” , vol. 32, no. 7, pp. 1–15, 2021.F.Haq, K.Shah, G.Rahman and M.Shahzad, “Numerical solution of fractional order smoking model via laplace adomian decomposition method,” , vol. 57, no. 2, pp. 1061–1069, 2018.S.Momani and Z.Odibat, “Numerical approach to differential equations of fractional order,” , vol. 207, pp. 96–110, 2007.S.Deniz and N.Bildik, “Comparison of adomian decomposition method and taylor matrix method in solving different kinds of partial differential equations,” , vol. 4, no. 4, pp. 292–298, 2014.H.Bulut, H. M.Baskonus and F. B. M.Belgacem, “The analytical solutions of some fractional ordinary differential equations by sumudu transform method,” , vol. 2013, pp. 1–6, 2013.S.Rathore, D.Kumar, J.Singh and S.Gupta, “Homotopy analysis sumudu transform method for nonlinear equations,” , vol. 4, no. 4, pp. 1–13, 2012.J.Singh and D.Kumar, “Homotopy perturbation sumudu transform method for nonlinear equations,” , vol. 4, no. 4, pp. 165–175, 2011.A. M. S.Mahdy, Y. A.Amer, M. S.Mohamed and E.Sobhy, “General fractional financial models of awareness with caputo–Fabrizio derivative,” , vol. 12, no. 11, pp. 1–9, 2020.H.Schmid and A.Huber, “Analysis of switched-capacitor circuits using driving-point signal-flowgraphs,” , vol. 96, pp. 495–507, 2018.A. M. S.Mahdy and M.Higazy, “Numerical different methods for solving the nonlinear biochemical reaction model,” , vol. 5, no. 6, pp. 1–17, 2019.N. H.Sweilam and S. M.Al-Mekhlafi, “Optimal control for a nonlinear mathematical model of tumor under immune suppression a numerical approach,” , vol. 39, pp. 1581–1596, 2018.N. H.Sweilam, S. M.Al-Mekhlafi and D.Baleanu, “Optimal control for a fractional tuberculosis infection model including the impact of diabetes and resistant strains,” , vol. 17, pp. 125–137, 2019.N. H.Sweilam, O. M.Saad and D. G.Mohamed, “Fractional optimal control in transmission dynamics of west Nile model with state and control time delay a numerical approach,” , vol. 2019, no. 210, pp. 1–25, 2019.N. H.Sweilam, O. M.Saad and D. G.Mohamed, “Numerical treatments of the tranmission dynamics of west Nile virus and it's optimal control,” , vol. 7, no. 2, pp. 9–38, 2019.A. M. S.Mahdy and N. A. H.Mukhtar, “Second kind shifted chebyshev polynomials for solving the model nonlinear ODEs,” , vol. 7, no. 4, pp. 391–401, 2017.N. H.Sweilam, A. M.Nagy and A. A.Sayed, “Second kind shifted chebyshev polynomials for solving space fractional order diffusion equation,” chaos,” , vol. 73, pp. 141–147, 2015.J. C.Mason and D. C.Handscomb, “,” Chapman and Hall, New York, NY, CRC, Boca Raton, 2003.I. A.Moneim and G. A.Mosa, “Modelling the hepatitis C with different types of virus genome,” , vol. 7, no. 1, pp. 3–13, 2006.J. G.Liu, X. J.Yang, Y. Y.Feng and P.Cui, “On the (N+1)-dimensional local fractional reduced differential transform method and its applications,” , vol. 43, no. 15, pp. 8856–8866, 2020.J. G.Liu, X. J.Yang and Y. Y.Feng, “Analytical solutions of some integral fractional differential–difference equations,” , vol. 34, no. 1, pp. 205009, 2020.J. G.Liu, X. J.Yang, Y. Y.Feng and M.Iqbal, “Group analysis to the time fractional nonlinear wave equation,” , vol. 31, no. 4, pp. 2050029, 2020.A. M. S.Mahdy, “Numerical studies for solving fractional integro-differential equations,” , vol. 3, no. 2, pp. 127–132, 2018.H.Khan, R.Shah, P.Kumam, D.Baleanu and M.Arif, “Laplace decomposition for solving nonlinear system of fractional order partial differential equations,” , vol. 2020, no. 1, pp. 1–18, 2020.A.Ghaffar, A.Ali, S.Ahmed, S.Akram, D.Baleanuet al., “A novel analytical technique to obtain the solitary solutions for nonlinear evolution equation of fractional order,” , vol. 2020, no. 1, pp. 1–15, 2020.D.Baleanu, H.Mohammadi and S.Rezapour, “A fractional differential equation model for the COVID-19 transmission by using the caputo–Fabrizio derivative,” , vol. 2020, no. 1, pp. 1–27, 2020.N. H.Sweilam, S. M.Al-Mekhlafi, T.Assiri and A.Atangana, “Optimal control for cancer treatment mathematical model using atangana–Baleanu–Caputo fractional derivative,” , vol. 2020, no. 1, pp. 1–21, 2020.N. H.Sweilam, D. M.El-Sakout and M. M.Muttardi, “Compact finite difference method to numerically solving a stochastic fractional advection-diffusion equation,” , vol. 2020, no. 1, pp. 1–20, 2020.K. A.Gepreel, M.Higazy and A. M. S.Mahdy, “Optimal control, signal flow graph, and system electronic circuit realization for nonlinear anopheles mosquito model,” , vol. 31, no. 9, pp. 1–18, 2020.A. M. S.Mahdy, H.Higazy, K. A.Gepreel and A. A. A.El-dahdouh, “Optimal control and bifurcation diagram for a model nonlinear fractional SIRC,” , vol. 59, no. 5, pp. 3481–3501, 2020.Y. A.Amer, A. M. S.Mahdy and H. A. R.Namoos, “Reduced differential transform method for solving fractional-order biological systems,” , vol. 13, no. 20, pp. 8489–8493, 2018.Y. A.Amer, A. M. S.Mahdy, R. T.Shwayaa and E. S. M.Youssef, “Laplace transform method for solving nonlinear biochemical reaction model and nonlinear emden-fowler system,” , vol. 13, no. 17, pp. 7388–7394, 2018.M. M.Khader, N. H.Sweilam and A. M. S.Mahdy, “The chebyshev collection method for solving fractional order klein-gordon equation,” , vol. 13, pp. 31–38, 2014.A. M. S.Mahdy, K.Lotfy, M. H.Ahmed, A.El-Bary and E. A.Ismail, “Electromagnetic hall current effect and fractional heat order for micro temperature photo-excited semiconductor medium with laser pulses,” , vol. 17, pp. 1–9, 2020.A. K.Khamis, K.Lotfy, A. A.El-Bary, A. M. S.Mahdy and M. H.Ahmed, “Thermal-piezoelectric problem of a semiconductor medium during photo-thermal excitation,” , pp. 1–15, 2020.A. M. S.Mahdy, K.Lotfy, E. A.Ismail, A. A.El-Bary, M.Ahmedet al. “Analytical solutions of time-fractional heat order for a magneto-photothermal semiconductor medium with thomson effects and initial stress,” , vol. 18, pp. 1–11, 2020.A. M. S.Mahdy, K.Lotfy, W.Hassan and A. A.El-Bary, “Analytical solution of magneto-photothermal theory during variable thermal conductivity of a semiconductor material due to pulse heat flux and volumetric heat source,” , pp. 1–19, 2020.A. M. S.Mahdy and E. S. M.Youssef, “Numerical solution technique for solving isoperimetric variational problems,” , vol. 32, no. 1, pp. 1–14, 2021.A. M. S.Mahdy, M. S.Mohamed, K. A.Gepreel, A.AL-Amiri and M.Higazy, “Dynamical characteristics and signal flow graph of nonlinear fractional smoking mathematical model,” , vol. 141, pp. 1–16, 2020.E. H.Doha, A. H.Bhrawy and S. S.Ezz-Eldien, “Efficient chebyshev spectral methods for solving multi-term fractional orders differential equations,” , vol. 35, pp. 5662–5672, 2011.E. H.Doha, A. H.Bhrawy and S. S.Ezz-Eldien, “A chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order,” , vol. 62, pp. 2364–2373, 2011.A. H.Bhrawy and M. A.Alshomrani, “Shifted legendre spectral method for fractional-order multi-point boundary value problems,” , vol. 2012, pp. 1–19, 2012.A. M. S.Mahdy, “Numerical solutions for solving model time-fractional fokker–Planck equation,” , vol. 37, no. 2, pp. 1120–1135, 2021.A.Khan, H.Khan, J. F.Gómez-Aguilar and T.Abdeljawad, “Existence and hyers-ulam stability for a nonlinear singular fractional differential equations with mittag-leffler kernel,” , vol. 127, pp. 422–427, 2019.A.Khan, J. F.Gómez-Aguilar, T. S.Khan and H.Khan, “Stability analysis and numerical solutions of fractional order HIV/AIDS model,” , vol. 122, pp. 119–128, 2019.I.Ullah, S.Ahmad, Q.Al-Mdallal, Z. A.Khan, H.Khanet al., “Stability analysis of a dynamical model of tuberculosis with incomplete treatment,” , vol. 2020, no. 1, pp. 1–14, 2020.K. A.Gepreel, A. M. S.Mahdy, M. S.Mohamed and A.Al-Amiri, “Reduced differential transform method for solving nonlinear biomathematics models,” , vol. 61, no. 3, pp. 979–994, 2019.Y. A.Amer, A. M. S.Mahdy and E. S. M.Youssef, “Solving fractional integro-differential equations by using sumudu transform method and hermite spectral collocation method,” , vol. 54, no. 2, pp. 161–180, 2018.K. A.Gepreel, M. S.Mohamed, H.Alotaibi1 and A. M. S.Mahdy, “Dynamical behaviors of nonlinear coronavirus (COVID-19) model with numerical studies,” , vo. vol. 67, no. 1, pp. 675–686, 2021.M. I. A.Othman and A. M. S.Mahdy, “Numerical studies for solving a free convection boundary–layer flow over a vertical plate,” , vol. 22, no. 1, pp. 41–48, 2018.H. H.Abdel-Halim, M. I. A.Othman and A. M. S.Mahdy, “Variational iteration method for solving twelve order boundary value problems,” , vol. 3, no. 15, pp. 719–730, 2009.A. M. S.Mahdy, N. H.Sweilam and M.Higazy, “Approximate solutions for solving nonlinear fractional order smoking model,” , vol. 59, no. 2, pp. 739–752, 2020.I. H. A.Hassan, M. I. A.Othman and A. M. S.Mahdy, “Variational iteration method for solving: Twelve order boundary value problems,” , vol. 3, no. 13–16, pp. 719–730, 2009.A. M. S.Mahdy, M.Higazy and M. S.Mohamed, “Optimal and memristor-based control of a nonlinear fractional tumor-immune model,” , vol. 67, no. 3, pp. 3463–3486, 2021.A. M. S.Mahdy, M. S.Mohamed, K.Lotfy, M.Alhazmi, A. A.El-Baryet al. “Numerical solution and dynamical behaviors for solving fractional nonlinear rubella ailment disease model,” , vol. 24, pp. 1–10, 2021.