Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

New insights on novel coronavirus 2019-nCoV/SARS-CoV-2 modelling in the aspect of fractional derivatives and fixed points

  • Correction on: Mathematical Biosciences and Engineering 19: 1588-1590
  • Extended orthogonal spaces are introduced and proved pertinent fixed point results. Thereafter, we present an analysis of the existence and unique solutions of the novel coronavirus 2019-nCoV/SARS-CoV-2 model via fractional derivatives. To strengthen our paper, we apply an efficient numerical scheme to solve the coronavirus 2019-nCoV/SARS-CoV-2 model with different types of differential operators.

    Citation: Sumati Kumari Panda, Abdon Atangana, Juan J. Nieto. New insights on novel coronavirus 2019-nCoV/SARS-CoV-2 modelling in the aspect of fractional derivatives and fixed points[J]. Mathematical Biosciences and Engineering, 2021, 18(6): 8683-8726. doi: 10.3934/mbe.2021430

    Related Papers:

    [1] Sumati Kumari Panda, Abdon Atangana, Juan J. Nieto . Correction: New insights on novel coronavirus 2019-nCoV/SARS-CoV-2 modelling in the aspect of fractional derivatives and fixed points. Mathematical Biosciences and Engineering, 2022, 19(2): 1588-1590. doi: 10.3934/mbe.2022073
    [2] Tahir Khan, Roman Ullah, Gul Zaman, Jehad Alzabut . A mathematical model for the dynamics of SARS-CoV-2 virus using the Caputo-Fabrizio operator. Mathematical Biosciences and Engineering, 2021, 18(5): 6095-6116. doi: 10.3934/mbe.2021305
    [3] A. D. Al Agha, A. M. Elaiw . Global dynamics of SARS-CoV-2/malaria model with antibody immune response. Mathematical Biosciences and Engineering, 2022, 19(8): 8380-8410. doi: 10.3934/mbe.2022390
    [4] Somayeh Fouladi, Mohammad Kohandel, Brydon Eastman . A comparison and calibration of integer and fractional-order models of COVID-19 with stratified public response. Mathematical Biosciences and Engineering, 2022, 19(12): 12792-12813. doi: 10.3934/mbe.2022597
    [5] Rahat Zarin, Usa Wannasingha Humphries, Amir Khan, Aeshah A. Raezah . Computational modeling of fractional COVID-19 model by Haar wavelet collocation Methods with real data. Mathematical Biosciences and Engineering, 2023, 20(6): 11281-11312. doi: 10.3934/mbe.2023500
    [6] A. M. Elaiw, Raghad S. Alsulami, A. D. Hobiny . Global dynamics of IAV/SARS-CoV-2 coinfection model with eclipse phase and antibody immunity. Mathematical Biosciences and Engineering, 2023, 20(2): 3873-3917. doi: 10.3934/mbe.2023182
    [7] Biplab Dhar, Praveen Kumar Gupta, Mohammad Sajid . Solution of a dynamical memory effect COVID-19 infection system with leaky vaccination efficacy by non-singular kernel fractional derivatives. Mathematical Biosciences and Engineering, 2022, 19(5): 4341-4367. doi: 10.3934/mbe.2022201
    [8] Adnan Sami, Amir Ali, Ramsha Shafqat, Nuttapol Pakkaranang, Mati ur Rahmamn . Analysis of food chain mathematical model under fractal fractional Caputo derivative. Mathematical Biosciences and Engineering, 2023, 20(2): 2094-2109. doi: 10.3934/mbe.2023097
    [9] M. Botros, E. A. A. Ziada, I. L. EL-Kalla . Semi-analytic solutions of nonlinear multidimensional fractional differential equations. Mathematical Biosciences and Engineering, 2022, 19(12): 13306-13320. doi: 10.3934/mbe.2022623
    [10] Noura Laksaci, Ahmed Boudaoui, Seham Mahyoub Al-Mekhlafi, Abdon Atangana . Mathematical analysis and numerical simulation for fractal-fractional cancer model. Mathematical Biosciences and Engineering, 2023, 20(10): 18083-18103. doi: 10.3934/mbe.2023803
  • Extended orthogonal spaces are introduced and proved pertinent fixed point results. Thereafter, we present an analysis of the existence and unique solutions of the novel coronavirus 2019-nCoV/SARS-CoV-2 model via fractional derivatives. To strengthen our paper, we apply an efficient numerical scheme to solve the coronavirus 2019-nCoV/SARS-CoV-2 model with different types of differential operators.



    In very recent years, the world has encountered with a deadly outbreak due to attack of a strange virus named novel coronavirus[1]. Initially, the outbreak started in Wuhan city of China, and rapidly spread to all parts of the world. The deadly infectious disease has already taken lives of thousands of people around the globe. It could be observed that the virus had hugely effected aged people and people with chronic health issues. It is believed that the spread is more aggressive in low temperatures; however, the science behind this has not been found yet. Despite many studies carried out on the origin and the spread of the virus, the researcher could not determine the exact cause of the spread of the virus including the finding up the life span of virus on surface. Some people believe that such a virus is from bats, pangolins, seafood, and others strongly believed the virus is man-made. The spread and fatality of the Coronavirus has made this world hostile place to live in. The spread of the strange virus within humans is a clear indication that mankind must have limits to interact with living beings. But the precautions taken against the virus prevented the severe spread of virus. Therefore, humans should take necessary care while interacting with other humans in day to day life. So far, there are many studies carried out to understand structure, the life span, the mode of transmission, effect on the human body, in terms of the affluence of temperature, age, etc[2,3,4]. However, to understand properly the spread, many researchers attempt to develop a mathematical model using the concept of differentiation and further some facts were derived through the observation. This observed factual information helps the mathematicians to come up with a suitable model with a solution that could show the predictable behaviour of the real-world problem as a function of time. Mathematicians can use this solution to understand, control, and predict the behavior of the real-world problem for a set of given data and parameters. In a latest study by Altaf and Atangana[5] suggested a mathematical model could depict the spread of the novel coronavirus, many factors are considered for the following model. The readers can refer to [5] to understand more about the variables in below mathematical model.

    dSpdt=ΠpapSpbpSp(Ip+ΨAp)NpbwSpM;dEpdt=bpSp(Ip+ΨAp)Np+bwSpM(1Υp)ωpEpΥpξpEpapEp;dIpdt=(1Υp)ωpEp(τp+ap)Ip;dApdt=ΥpξpEp(τϑp+ap)Ap;dRpdt=τpIp+τϑpApapRp;dMdt=cpIp+epApπM.} (1.1)

    Above all these studies suggested with help of the data sources, a reproductive number around 2.4, which was in good accordance with the value suggested by the World Health Organization(WHO).

    Further, fractional calculus plays a pivotal role in the mathematical modeling of many bio-medical and engineering disciplines. The following is detailed review of literature with references connected to the study.

    ● In [6,7,8], we can see that, new fractional derivatives with nonlocal and non-singular kernel and their properties.

    ● In [9,10,11,12], there is a vigorous study via numerical and/or computational methods for various fractional derivatives with nonlocal and non-singular kernel and also for ordinary differential equations.

    ● Addressed the integro-differential equations utilizing the Atangana-Baleanu fractional derivatives, that will have the generalised Mittag-Leffler functions in their kernels and have been studied extensively. They investigated the existence and uniqueness of solutions to such equations in Banach spaces using fixed point approach in [13,14].

    ● We can see valid discussions and new results on nonlocal functional integro-differential equations with different kind of approaches such as non-compactness, with impulsive conditions and Hilfer fractional derivative in [15,16,17].

    ● In [18,19,20], using the concept of fractional differentiation based on the generalized Mittag-Leffler function, they presented some interesting simulations.

    ● Applied various methods such as fractal-fractional Atangana-Baleanunu order, Atangana-Baleanu Caputo (ABC) derivative and stationary distribution extinction and optimal control in mathematical modelling for hepatitis B and drinking alcohol in [21,22,23,24,25].

    However, the analysis presented in their paper is based on numerical approximation, the numerical methods have been promising to be effective mathematical tools that could be used where analytical techniques failed. It is essential that, to note that, the exact solution of given model with its proof is important to be observed. Therefore in this paper, extended orthogonal spaces are adopted and established valid fixed point results. Afterward, using fractional derivatives, we present an overview of the existence of solutions for the novel coronavirus 2019-nCoV/SARS-CoV-2 model. An investigation of the proposed model's existence of solutions is presented. Furthermore, we use a numerical approach to address the proposed model using several forms of differential operators.

    The accessible collected works show that there is a clear connection between Metric Fixed Point Theory and many academic fields including but not limited to engineering, sciences, and technology. One of the clear applications of this is found in the theory of ordinary and partial differential equations. These equations are used to model real-world problems, thus their exact solutions are requested to be compared with experimental data. However on many occasions, exact solutions are impossible to be obtained, thus researchers relied on numerical approximation to provide a numerical solution that can be used for comparison with data. However, it is at least useful to provide conditions under which the model has an exact solution, in this case, metric fixed point theory is very useful. Several extensions or generalizations have been made and extensions have been done to capture problems that are more complex.

    In recent past, Kamran et al. [26] introduced a novel approach to the study of b-metric space as below:

    Definition 2.1 (see [26]). Let S be a non-empty set and μ:S×S[1,). A function D:S×S[0,) is called an extended b-metric if, for all u,v,wS, it satisfies

    (i). D(u,v)=0 if and only if u=v;

    (ii). D(u,v)=D(v,u);

    (iii). D(u,v)μ(u,v)[D(u,w)+D(w,v)].

    The pair (S,D) is called an extended b-metric space.

    His findings and/or methods also laid the groundwork for the observation of the role of extended dislocated b-metric[27], extended complex b-metric[28], controlled metric type[29], extended F-metric[30] and many more (For more info, reader can refer [31,32,33] and [34,35].

    Definition 2.2. Let S and S×S be a binary relation. If satisfies the below assertion:

    i.e., there exist u0 such that for all vS,vu0 or for all vS,u0v, then it is called an orthogonal set (clearly, O-set) and u0 is called an orthogonal element. We denote this O-set by (S,).

    Example 2.1. Let S=W. Define the binary relation on S by mn if there exist kW such that m=kn. It is easy to see that On for all nW. Hence (S,) is an O-set.

    Example 2.2. Let A={1,2,3}. Then the power set of A is P(A)={,{1},{2},{3},{1,2},{1,3},{2,3},{1,2,3}}. Define the binary relation on P(A) by CD such that =CD. It is easy to see that D for all DP(A). Hence (P(A),) is an O-set.

    Example 2.3. Let S be the set of all peoples in the world. Define the binary relation on S by uv if u can get novel coronavirus 2019-nCoV/SARS-CoV-2 infection by u. If u0 is a person such that he/she infected by novel coronavirus 2019-nCoV/SARS-CoV-2, then we have u0v for all vS. Thus (S,) is an O-set. In this regard O-set, u0 is not unique.

    More clearly, If we consider a social network of a particular person:

    The huge hub in the outline is a super-spreader, an individual who associates with numerous others, and along these lines can possibly spread the novel coronavirus 2019-nCoV/SARS-CoV-2 infection. Two spots are associated by a line on the off chance that they are in contact with one another, and the more contacts an individual has, the bigger their hub is on the network.

    Figure 1.  Reparable stage red hubs indicate virus opposite to blue hubs.

    Example 2.4. Let S be the set of all peoples who effected by any one of the listed diseases, so called: EBOV-Ebola virus, MERS-CoV, SARS-CoV and FECV-feline enteric conronavirus. Define the binary relation on S by uv if u and v will hit/target the disease module by drug combination for all u,vS. According to below given table, Antineoplastic drugs may furnish a synergistic repercussions in vigorously treating above mentioned EBOV-Ebola virus, MERS-CoV, SARS-CoV and FECV-feline enteric conronavirus diseases as in Table.1 (For more about eight potential repurposable drugs, the reader can refer [36]). Thus, (S,) is an O-set.

    Now we give the concepts of an extended O-sequence, an extended -continuous mapping, an extended O-complete orthogonal space, an extended -preserving mapping.

    Definition 2.3. The triplet (S,,D) is called an extended orthogonal space if (S,) is an O-set and (S) is an extended orthogonal space.

    Figure 2.  Still reparable but needs immunity as the red hubs are increasing in number.
    Figure 3.  Advanced stage or contaminating stage.

    Definition 2.4. Let (S,,D) be an extended orthogonal space. A sequence {un} is called an extended orthogonal sequence (shortly, EO-sequence) if for all nN,unun+1 (or) nN,un+1un

    Definition 2.5. Let (S,,D) be an extended orthogonal space. Then a mapping χ:SS is said to be an extended orthogonally continuous (shortly, E-continuous) in uS, if for each EO-sequence {un} in S with unu as n, we have χ(un)χ(u) as n. Also χ is said to be E-continuous on S if χ is E-continuous in each uS.

    Remark 1.1. Every continuous mapping is E-continuous and the converse is not true.

    Definition 2.6. Let (S,,D) be an extended orthogonal space. Then S is said to be an extended orthogonally complete (shortly, EO-complete) if every Cauchy EO-sequence is convergent.

    Remark 1.2. Every complete extended space is EO-complete space and the converse is need not to true.

    Definition 2.7. Let (S,) be an O-set. A mapping χ:SS is said to be -preserving if χ(u)χ(v) whenever uv. Also, χ:SS is said to be weakly -preserving if χ(u)χ(v) or χ(v)χ(u) whenever uv.

    Figure 4.  Repair is impossible but we can see how red exceeds blue in doubling size and number.
    Table 1.  Potentially Repurposable Drugs vs Reported Antivirals[36].
    Drug Name Structure Known Indications Reported Antiviral Profiles
    Toremifene Antineoplastic Ebola virus (shortly, EBOV)
    Camphor Antipruritic Swine Flu-H1N1
    Toremifene Antineoplastic Middle East respiratory syndrome-related coronavirus (shortly, MERS-CoV/EMC/2012)
    Equilin Estrogen Zaire Ebola virus Glycoprotein (shortly, ZEBOV-GP)
    Toremifene Antineoplastic Severe acute respiratory syndrome-related coronavirus (shortly, SERS-CoV)
    Emodin Anthroquinones Hepatitis-B-Virus (shortly, HBV)
    Toremifene Antineoplastic Feline Enteric Coronavirus (shortly, FECV)
    Oxymetholone Anabolic Steroid Human Immunodeficiency Virus (Shortly, HIV)

     | Show Table
    DownLoad: CSV

    Theorem 2.1. Let (S,,D) be a complete extended orthogonal space. Suppose that H:SS is an E-continuous mapping such that the following assertions hold:

    (i) there exists 0<k<1 such that D(Hu,Hv)kD(u,v) for all u,vS with uv and for each n0S, limn,mD(un,um)<1k, here un=Hnu0, n=1,2,3...

    (ii) H is E-preserving.

    Then H has a unique fixed point in S. Also the picard sequence {Hnu0} converges to the fixed point of H.

    Proof. From the definition of the orthogonality, it follows that u0H(u0) or H(u0)u0.

    Let,

    u1:=Hu0,u2:=Hu1=H2u0,
    u3:=Hu2=H2u1=H3u0,
    un:=Hun1=Hnu0, for all nN{0}.

    If un=un+1 for some nN{0}, then un is a fixed point of H and so the proof is completed. So, we may assume that unun+1 for all nN{0}. Thus we have D(Hun,Hun+1)>0 for all nN{0}. Since H is -preserving, we have unun+1 or un+1un for all nN{0}. This implies that {un} is an EO-sequence.

    Consider,

    D(Hun,Hun+1)=D(Hun1,Hun)kD(Hun1,Hun)k2D(Hun2,Hun1)knD(Hu0,Hu1), for all nN (2.1)

    By triangular inequality and Eq (2.1), we have,

    D(un,um)=μ(un,um)[D(Hun,Hun+1)+D(Hun+1,Hum)]μ(un,um)D(Hun,Hun+1)+μ(un,um)D(Hun+1,Hum)μ(un,um)D(Hun,Hun+1)+μ(un,um)μ(un+1,um)[D(Hun+1,Hun+2)+D(Hun+2,Hum)]μ(un,um)knD(Hu0,Hu1)+μ(un,um)μ(un+1,um)kn+1μ(u0,u1)       +.....+μ(un,um)μ(un+1,um)μ(un+2,um).....μ(um2,um)μ(um1,um)km1D(Hu0,Hu1)

    This will leads to,

    D(un,um)D(u0,u1)[μ(un,um)kn+μ(un,um)μ(un+1,um)kn+1       +.....+μ(un,um)μ(un+1,um)μ(un+2,um).....μ(um1,um)km1]D(u0,u1)[μ(u1,um)μ(u2,um)....μ(un1,um)μ(un,um)kn       +μ(u1,um)μ(u2,um).....μ(un,um)μ(un+1,um)kn+1+....       +μ(u1,um)μ(u2,um).....μ(um2,um)μ(um1,um)km1] (2.2)

    Since limn,mμ(un+1,um)k<1, the series Σn=1knni=1μ(ui,um) converges by ratio test for each mN.

    Let

    S=Σn=1knni=1μ(ui,um)
    Sn=Σj=1kjji=1μ(ui,um)

    Thus, for m>n, above inequality (2.2) implies

    D(un,um)D(u0,u1)[Sm1Sn].

    Letting n, it follows that {un} is a Cauchy EO-sequence in S. Since S is EO-complete, there exists uS such that unu as n. Since H is E-continuous, we have,

    Hu=H(limnHun)=limnun+1=u

    and so u is a fixed point of H.

    Next we prove the uniqueness of fixed point of H.

    Assume that u is an another fixed point of H. If unu as n, we have u=u. If un does not converge to u as n, there is a subsequence {unα} such that

    Hunαu,  for all αN. (2.3)

    By the choice of u0, we have u0u or uu0. Since H is E-preserving and Hnu=u for all nN, we have Hnu0u or uHnu0 for all nN.

    Consider,

    D(Hnαu0,u)=D(Hnαu0,Hnαu)knD(u0,u)  for all αN

    This implies that D(Hnαu0,u)0 as n since u0u, which is a contradiction to Eq (2.3).

    This yields unu as n. Hence u=u. Thus H has a unique fixed point.

    Example 2.5. Let S=[0,). Define D:S×SR+ and μ:S×S[1,) as D(u,v)=(uv)2, μ(u,v)=2u+3y+5 for all u,vS.

    Consider the sequence {Δn}nN defined as Δn=n6(2n+1)(n+1) for all nN{0}.

    Define a relation on S by uvuv{u,v}{Δn}. Thus (S,,D) is EO-complete space. Now define a mapping H:S×S by

    Hu={Δ0                              if Δ0uΔ1Δn1(Δn+1x)+Δn(uΔn)Δn+1Δn   if  Δ0uΔn+1 for each n1

    It is easy to see that H is E-continuous and H is E-preserving.

    Let u,vS with uv and D(Hu,Hv)>0. Without loss of generality we assume that u<v. This implies that u{Δ0,Δ1} and v=Δi for some iN{1}.

    Then,

    D(Hu,Hv)D(u,v)=D(Δ0,HΔi)D(Δ0,Δi)=D(Δ0,Δi1(Δi+1Δi)+Δi(ΔiΔi)Δi+1Δi)D(Δ0,Δi)=D(Δ0,Δi1)D(Δ0,Δi)=(Δ0Δi1)2(Δ0Δi)2=(Δi1)2(Δi)2=[(i1)2(2i1)2i236][(i+1)2(2i+1)2i236]=(i1)2(2i1)2(i+1)2(2i+1)2<k ; where k(0,1)for all iN{1} (2.4)

    Thus D(Hu,Hv)kD(u,v). Therefore, all the assertions managed by above theorem satisfied and so H has a unique fixed point u=Δ0.

    By applying various setting on above Theorem 2.1 we may get below corollaries. For suppose, if we consider μ(u,v)=b(1), then we can reduce above theorem to below corollary.

    Corollary 2.1. Let (S,,D) be a O-complete b-metric space. Suppose that H:SS is an -continuous mapping such that the following assertions hold:

    (i) there exists 0<k<1 such that D(Hu,Hv)kD(u,v) for all u,vS with uv and for each n0S, limn,mD(un,um)<1k, here un=Hnu0, n=1,2,3...

    (ii) H is -preserving.

    Then H has a unique fixed point in S.

    If we consider μ(u,v)=1, then we can reduce above theorem to below corollary.

    Corollary 2.2. Let (S,,D) be a O-complete metric space. Suppose that H:SS is an -continuous mapping such that the following assertions hold:

    (i) there exists 0<k<1 such that D(Hu,Hv)kD(u,v) for all u,vS with uv and for each n0S, limn,mD(un,um)<1k, here un=Hnu0, n=1,2,3...

    (ii) H is -preserving.

    Then H has a unique fixed point in S.

    Although humans have advanced their technologies and their medicine to protect the mankind, it is important to note that there are still diseases that continue to spread across the world that up to now humans have not found a permanent solution to eradicating these diseases. A few examples can be noted including HIV, malaria, Ebola, and other well-known infectious diseases that continue to take much life daily. In addition to these existing diseases, humans continue to be surprised by new ones that lead them to double their efforts to make their environment more saver, however, diseases are known to be the main source of death across the world. For example, information released by the WHO in 2001 indicated that diseases have added 26% of global mortality. From 2003, the world has witnessed a rise of different types of viruses for example SARS-CoV in 2003, Mers-CoV in 2012, Ebola in 2014, and now SARS-CoV-2 and others that are not listed here. Researchers from different backgrounds have investigated the biological structures, mode of transmission, and life span of these viruses. On the other hand, to protect life and livelihood medical bodies have developed several vaccines that could help prevent the spread of these diseases. Nevertheless, to have an idea of what could happen in a near future, mathematicians have developed some mathematical models that could be used first to provide the reproductive number of these diseases, second to predict the numbers of infected, recovery, and death cases. This was also valid with the breakout of the new virus known as covid-19, several investigations have been done, vaccines have been developed and some mathematical models have been proposed in the last two years with the aim to predict future behaviors of the spread. It's important to note that the use of mathematical modeling in the analysis of epidemiological disorders has risen significantly. In this work, we propose solutions for fractional-order models of the novel coronavirus 2019-nCoV/SARS-CoV-2.

    Let S={AC(I,R):A(t)>0 for all tI=[0,H],H>0}. Define the orthogonality relation on S by uvu(t)v(t)u(t) or u(t)v(t)v(t) for all tI.

    Define a mapping D:S×S[0,) and μ:S×S[0,) as D(u,v)=|uv| for all u,vS. Then (S,,D) is a O-complete metric space. Define a mapping H:S×S by,

    HSp(t)=Sp(0)+1ΛAB(Λ)A(t,Sp)+ΛAB(Λ)Γ(Λ)t0A(ϱ,Sp)(tϱ)Λ1dϱHEp(t)=Ep(0)+1ΛAB(Λ)B(t,Ep)+ΛAB(Λ)Γ(Λ)t0B(ϱ,Ep)(tϱ)Λ1dϱHIp(t)=Ip(0)+1ΛAB(Λ)C(t,Ip)+ΛAB(Λ)Γ(Λ)t0C(ϱ,Ip)(tϱ)Λ1dϱHAp(t)=Ap(0)+1ΛAB(Λ)D(t,Ap)+ΛAB(Λ)Γ(Λ)t0D(ϱ,Ap)(tϱ)Λ1dϱHRp(t)=Rp(0)+1ΛAB(Λ)E(t,Rp)+ΛAB(Λ)Γ(Λ)t0E(ϱ,Rp)(tϱ)Λ1dϱHMp(t)=Mp(0)+1ΛAB(Λ)F(t,Mp)+ΛAB(Λ)Γ(Λ)t0F(ϱ,Mp)(tϱ)Λ1dϱ} (3.1)

    Then, we see that H is -continuous. Now we show that H is -preserving. For each u,vS with uv and tI, we have,

    HSp(t)=Sp(0)+1ΛAB(Λ)A(t,Sp)+ΛAB(Λ)Γ(Λ)t0A(ϱ,Sp)(tϱ)Λ1dϱ1

    It follows that, [HSp(t)][HSq(t)][HSq(t)] and so [HSp(t)][HSq(t)]. Then H is -preserving.

    In this section, by using Corollary 2.2, we prove the existence and uniqueness of solutions of Eq (3.1) provided the following conditions satisfied:

    1) |A(t1,Sp)A(t2,Sp)||Sp(t1)Sp(t2)|;

    2) |A(k,Sp)Λ|<1;

    3) |tΛ1tΛ2||Sp(t1)Sp(t2)|;

    4) |(1Λ)Γ(Λ)+ΛAB(Λ)Γ(Λ)|<γ;  where 0<γ<1.

    Now consider,

    |HSp(t1)HSp(t2)|=|Sp(0)+1ΛAB(Λ)A(t1,Sp)+ΛAB(Λ)Γ(Λ)t10A(ϱ,Sp)(t1ϱ)Λ1dϱ      Sp(0)1ΛAB(Λ)A(t2,Sp)ΛAB(Λ)Γ(Λ)t20A(ϱ,Sp)(t2ϱ)Λ1dϱ||1ΛAB(Λ)A(t1,Sp)+ΛAB(Λ)Γ(Λ)t10A(ϱ,Sp)(t1ϱ)Λ1dϱ      1ΛAB(Λ)A(t2,Sp)ΛAB(Λ)Γ(Λ)t20A(ϱ,Sp)(t2ϱ)Λ1dϱ||1ΛAB(Λ)A(t1,Sp)1ΛAB(Λ)A(t2,Sp)|  +|ΛAB(Λ)Γ(Λ)t10A(ϱ,Sp)(t1ϱ)Λ1dϱΛAB(Λ)Γ(Λ)t20A(ϱ,Sp)(t2ϱ)Λ1dϱ||1ΛAB(Λ)||A(t1,Sp)A(t2,Sp)|      +|ΛAB(Λ)Γ(Λ)||t10A(ϱ,Sp)(t1ϱ)Λ1dϱt20A(ϱ,Sp)(t2ϱ)Λ1dϱ||1ΛAB(Λ)||A(t1,Sp)A(t2,Sp)|      +|ΛAB(Λ)Γ(Λ)||A(ϱ,Sp)[(t1ϱ)ΛΛ]t10A(ϱ,Sp)[(t2ϱ)ΛΛ]t20|

    This can be reduced in to:

    |HSp(t1)HSp(t2)||1ΛAB(Λ)||A(t1,Sp)A(t2,Sp)|      +|ΛAB(Λ)Γ(Λ)||A(ϱ,Sp)tΛ1ΛA(ϱ,Sp)tΛ2Λ||1ΛAB(Λ)||A(t1,Sp)A(t2,Sp)|      +|ΛAB(Λ)Γ(Λ)||A(ϱ,Sp)Λ(tΛ1tΛ2)||(1Λ)Γ(Λ)+ΛAB(Λ)Γ(Λ)||Sp(t1)Sp(t2)|k|Sp(t1)Sp(t2)|

    Therefore,

    |HSp(t1)HSp(t2)|k|Sp(t1)Sp(t2)|. (3.2)

    By following the same pattern as above, we get,

    |HEp(t1)HEp(t2)|k|Ep(t1)Ep(t2)|;
    |HIp(t1)HIp(t2)|k|Ip(t1)Ip(t2)|;
    |HAp(t1)HAp(t2)|k|Ap(t1)Ap(t2)|;
    |HRp(t1)HRp(t2)|k|Rp(t1)Rp(t2)|;
    |HMp(t1)HMp(t2)|k|Mp(t1)Mp(t2)|.

    From Eq (3.2), we can write, D(HSp(t1),HSp(t2))kD(Sp(t1),Sp(t2)). Similarly, we can get,

    D(HEp(t1),HEp(t2))kD(Ep(t1),Ep(t2));
    D(HIp(t1),HIp(t2))kD(Ip(t1),Ip(t2));
    D(HAp(t1),HAp(t2))kD(Ap(t1),Ap(t2));
    D(HRp(t1),HRp(t2))kD(Rp(t1),Rp(t2));
    D(HMp(t1),HMp(t2))kD(Mp(t1),Mp(t2)).

    Hence the system of novel corona virus model gratified all the assertions of Corollary 2.2. Hence Eq (3.1) has unique fixed point. Thus Eq (3.1) has a unique solution.

    In this section, we replace Atangana-Baleanu derivative with Caputo-Fabrizio fractional derivative and generalized model (1.1) in the following way:

    CFDSp=ΠpapSpbpSp(Ip+ΨAp)NpbwSpM;CFDEp=bpSp(Ip+ΨAp)Np+bwSpM(1Υp)ωpEpΥpξpEpapEp;CFDIp=(1Υp)ωpEp(τp+ap)Ip;CFDAp=ΥpξpEp(τϑp+ap)Ap;CFDRp=τpIp+τϑpApapRp;CFDM=cpIp+epApπM.} (3.3)

    By applying fractional integral operator, we convert model (3.3) into below integral form:

    Sp(t)Sp(0)=CFI[ΠpapSpbpSp(Ip+ΨAp)NpbwSpM];Ep(t)Ep(0)=CFI[bpSp(Ip+ΨAp)Np+bwSpM(1Υp)ωpEpΥpξpEpapEp];Ip(t)Ip(0)=CFI[(1Υp)ωpEp(τp+ap)Ip];Ap(t)Ap(0)=CFI[ΥpξpEp(τϑp+ap)Ap];Rp(t)Rp(0)=CFI[τpIp+τϑpApapRp];Mp(t)Mp(0)=CFI[cpIp+epApπM].} (3.4)

    Let S={AC(I,R):A(t)>0 for all tI=[0,H],H>0}. Define the orthogonality relation on S by uvu(t)v(t)u(t) or u(t)v(t)v(t) for all tI.

    Define a mapping D:S×S[0,) and μ:S×S[0,) as D(u,v)=|uv|2 and μ(u,v)=u(t)+v(t)+2 for all u,vS. Then (S,,D) is a complete extended orthogonal space. Define a mapping H:S×S by,

    HSp(t)=1ΛB(Λ)A(t,Sp(t))+ΛB(Λ)t0A(ϑ,Sp(ϑ))dϑHEp(t)=1ΛB(Λ)B(t,Sp(t))+ΛB(Λ)t0B(ϑ,Sp(ϑ))dϑHIp(t)=1ΛB(Λ)C(t,Sp(t))+ΛB(Λ)t0C(ϑ,Sp(ϑ))dϑHAp(t)=1ΛB(Λ)D(t,Sp(t))+ΛB(Λ)t0D(ϑ,Sp(ϑ))dϑHRp(t)=1ΛB(Λ)D(t,Sp(t))+ΛB(Λ)t0E(ϑ,Sp(ϑ))dϑHMp(t)=1ΛB(Λ)F(t,Sp(t))+ΛB(Λ)t0F(ϑ,Sp(ϑ))dϑ} (3.5)

    Then, we see that H is E-continuous. Now we show that H is E-preserving. For each u,vS with uv and tI, we have,

    HSp(t)=1ΛB(Λ)A(t,Sp(t))+ΛB(Λ)t0A(ϑ,Sp(ϑ))dϑ1

    It follows that, [HSp1(t)][HSq1(t)][HSq1(t)] and so [HSp1(t)][HSq1(t)]. Then H is E-preserving.

    In this section, by using Theorem 2.1, we prove the existence and uniqueness of solutions of Eq (3.5) provided the following conditions satisfied:

    1) |A(t,Sp1(t))A(t,Sp2(t))|2|Sp1(t)Sp2(t)|2

    2) (1Λ+σtB(Λ))2<ξ,  where  0<ξ<1

    Now consider,

    |HSp1(t)HSp2(t)|2=|1ΛB(Λ)A(t,Sp1(t))+ΛB(Λ)t0A(ϑ,Sp1(ϑ))dϑ     1ΛB(Λ)A(t,Sp2(t))ΛB(Λ)t0A(ϑ,Sp2(ϑ))dϑ|2=|1ΛB(Λ)(A(t,Sp1(t))A(t,Sp2(t)))     +ΛB(Λ)[t0(A(ϑ,Sp1(ϑ))A(ϑ,Sp2(ϑ)))dϑ]|2=|1ΛB(Λ)(A(t,Sp1(t))A(t,Sp2(t)))|2    +|ΛB(Λ)[t0(A(ϑ,Sp1(ϑ))A(ϑ,Sp2(ϑ)))dϑ]|2    +2|1ΛB(Λ)(A(t,Sp1(t))A(t,Sp2(t)))|                ×|ΛB(Λ)[t0(A(ϑ,Sp1(ϑ))A(ϑ,Sp2(ϑ)))dϑ]|=|1ΛB(Λ)|2|A(t,Sp1(t))A(t,Sp2(t))|2    +|ΛB(Λ)|2|A(ϑ,Sp1(ϑ))A(ϑ,Sp2(ϑ))|2|t0dϑ|2    +2|(1Λ)Λ(B(Λ))2||A(t,Sp1(t))A(t,Sp2(t))|2(t)=|1ΛB(Λ)|2|A(t,Sp1(t))A(t,Sp2(t))|2    +|ΛB(Λ)|2|A(ϑ,Sp1(ϑ))A(ϑ,Sp2(ϑ))|2t2    +2|(1Λ)Λ(B(Λ))2||A(t,Sp1(t))A(t,Sp2(t))|2(t)|1ΛB(Λ)|2|Sp1(t)Sp2(t)|2    +|ΛB(Λ)|2|Sp1(t)Sp2(t)|2t2    +2|(1Λ)Λ(B(Λ))2||Sp1(t)Sp2(t)|2t=[|1ΛB(Λ)|2+|ΛB(Λ)|2t2+2(1Λ)Λ(B(Λ))2t]|Sp1(t)Sp2(t)|2=(1ΛB(Λ)+ΛB(Λ)t)2|Sp1(t)Sp2(t)|2<ξ|Sp1(t)Sp2(t)|2

    Therefore,

    |HSp1(t)HSp2(t)|2ξ|Sp1(t)Sp2(t)|2. (3.6)

    By following the same pattern as above, we get,

    |HEp1(t)HEp2(t)|2ξ|Ep1(t)Ep2(t)|2;|HIp1(t)HIp2(t)|2ξ|Ip1(t)Ip2(t)|2;|HAp1(t)HAp2(t)|2ξ|Ap1(t)Ap2(t)|2;|HRp1(t)HRp2(t)|2ξ|Rp1(t)Rp2(t)|2;|HMp1(t)HMp2(t)|2ξ|Mp1(t)Mp2(t)|2. (3.7)

    From Eq (3.6) we can write, D(HSp1(t),HSp2(t))ξD(Sp1(t),Sp2(t)).

    Similarly,

    D(HEp1(t),HEp2(t))ξD(Ep1(t),Ep2(t));
    D(HIp1(t),HIp2(t))ξD(Ip1(t),Ip2(t));
    D(HAp1(t),HAp2(t))ξD(Ap1(t),Ap2(t));
    D(HRp1(t),HRp2(t))ξD(Rp1(t),Rp2(t));
    D(HMp1(t),HMp2(t))ξD(Mp1(t),Mp2(t)).

    Hence the system of novel corona virus model gratified all the assertions of Theorem 2.1. Hence Eq (3.5) has unique fixed point. Thus Eq (3.5) has a unique solution.

    Let S={AC(I,R):A(t)>0 for all tI=[0,H],H>0}.

    Define the orthogonality relation on S by uvu(t)v(t)u(t) or u(t)v(t)v(t) for all tI.

    Define a mapping D:S×S[0,) and μ:S×S[0,) as D(u,v)=|uv| for all u,vS. Then (S,,D) is a O-complete metric space. Define a mapping H:S×S by,

    HSp(t)=Sp(0)+1Γ(Λ)t0(tτ)Λ1A(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτHEp(t)=Ep(0)+1Γ(Λ)t0(tτ)Λ1B(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτHIp(t)=Ip(0)+1Γ(Λ)t0(tτ)Λ1C(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτHAp(t)=Ap(0)+1Γ(Λ)t0(tτ)Λ1D(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτHRp(t)=Rp(0)+1Γ(Λ)t0(tτ)Λ1E(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτHMp(t)=Mp(0)+1Γ(Λ)t0(tτ)Λ1F(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτ} (3.8)

    Then, we see that H is -continuous. Now we show that H is -preserving. For each u,vS with uv and tI, we have,

    HSp(t)=Sp(0)+1Γ(Λ)t0(tτ)Λ1A(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτ1

    It follows that, [HSp1(t)][HSq1(t)][HSq1(t)] and so [HSp1(t)][HSq1(t)]. Then H is -preserving.

    In this section, by using Corollary 2.2, we prove the existence and uniqueness of solutions of Eq (3.8) provided the following conditions satisfied:

    1) |A(τ,Sp1,Ep1,Ip1,Ap1,Rp1,Mp1)A(τ,Sp2,Ep2,Ip2,Ap2,Rp2,Mp2)|<|Sp1(t)Sp2(t)|

    2) tΛΓ(Λ+1)<ρ,  where  0<ρ<1

    Now consider,

    |HSp1(t)HSp2(t)|=|1Γ(Λ)t0(tτ)Λ1A(τ,Sp1,Ep1,Ip1,Ap1,Rp1,Mp1)dτ           1Γ(Λ)t0(tτ)Λ1A(τ,Sp2,Ep2,Ip2,Ap2,Rp2,Mp2)dτ|=|1Γ(Λ)(A(τ,Sp1,Ep1,Ip1,Ap1,Rp1,Mp1)A(τ,Sp2,Ep2,Ip2,Ap2,Rp2,Mp2))           ×t0(tτ)Λ1dτ||(A(τ,Sp1,Ep1,Ip1,Ap1,Rp1,Mp1)A(τ,Sp2,Ep2,Ip2,Ap2,Rp2,Mp2))|           ×|1Γ(Λ)t0(tτ)Λ1dτ|=|(A(τ,Sp1,Ep1,Ip1,Ap1,Rp1,Mp1)A(τ,Sp2,Ep2,Ip2,Ap2,Rp2,Mp2))||tΛΛΓ(Λ)|=tΛΓ(Λ+1)|(A(τ,Sp1,Ep1,Ip1,Ap1,Rp1,Mp1)A(τ,Sp2,Ep2,Ip2,Ap2,Rp2,Mp2))|<ρ|Sp1(t)Sp2(t)|

    Therefore,

    |HSp1(t)HSp2(t)|ρ|Sp1(t)Sp2(t)|. (3.9)

    By following the same pattern as above, we get,

    |HEp1(t)HEp2(t)|ρ|Ep1(t)Ep2(t)|;|HIp1(t)HIp2(t)|ρ|Ip1(t)Ip2(t)|;|HAp1(t)HAp2(t)|ρ|Ap1(t)Ap2(t)|;|HRp1(t)HRp2(t)|ρ|Rp1(t)Rp2(t)|;|HMp1(t)HMp2(t)|ρ|Mp1(t)Mp2(t)|. (3.10)

    From Eq.(3.6) we can write, D(HSp1(t),HSp2(t))ρD(Sp1(t),Sp2(t)). Similarly,

    D(HEp1(t),HEp2(t))ρD(Ep1(t),Ep2(t));
    D(HIp1(t),HIp2(t))ρD(Ip1(t),Ip2(t));
    D(HAp1(t),HAp2(t))ρD(Ap1(t),Ap2(t));
    D(HRp1(t),HRp2(t))ρD(Rp1(t),Rp2(t));
    D(HMp1(t),HMp2(t))ρD(Mp1(t),Mp2(t)).

    Hence the system of novel corona virus model gratified all the assertions of Corollary 2.2. Hence Eq (3.8) has unique fixed point. Thus Eq (3.8) has a unique solution.

    It has been recognized that dynamical spread exhibits nonlocal behaviors resembling processes like that of a power-law, fatigue, and crossover. These patterns are well described by some mathematical functions, for example, the power-law function is known to replicate power-law behaviors as time and power change. The fatigue effect can be well explained using an exponential decay function. The generalized Mittag-Leffler has been recognized to depict a passage from stretched fatigue effect to power-law. From these functions and the convolution, concepts of fractional differentiation and integration have been introduced and applied in many problems in the last decades. In this section, we investigate the model of 2019-nCoV/SARS-CoV-2 using differential operators based on the generalized Mittag-Leffler kernel.

    AB0DΛtSp(t)=A(t,Sp,Ep,Ip,Ap,Rp,Mp)AB0DΛtEp(t)=B(t,Sp,Ep,Ip,Ap,Rp,Mp)AB0DΛtIp(t)=C(t,Sp,Ep,Ip,Ap,Rp,Mp)AB0DΛtAp(t)=D(t,Sp,Ep,Ip,Ap,Rp,Mp)AB0DΛtRp(t)=E(t,Sp,Ep,Ip,Ap,Rp,Mp)AB0DΛtMp(t)=F(t,Sp,Ep,Ip,Ap,Rp,Mp)}

    The novel coronavirus 2019-nCoV/SARS-CoV-2 model further can be written as:

    HSp(t)=Sp(0)+1ΛAB(Λ)A(t,Sp)+ΛAB(Λ)Γ(Λ)t0A(ϱ,Sp)(tϱ)Λ1dϱHEp(t)=Ep(0)+1ΛAB(Λ)B(t,Ep)+ΛAB(Λ)Γ(Λ)t0B(ϱ,Ep)(tϱ)Λ1dϱHIp(t)=Ip(0)+1ΛAB(Λ)C(t,Ip)+ΛAB(Λ)Γ(Λ)t0C(ϱ,Ip)(tϱ)Λ1dϱHAp(t)=Ap(0)+1ΛAB(Λ)D(t,Ap)+ΛAB(Λ)Γ(Λ)t0D(ϱ,Ap)(tϱ)Λ1dϱHRp(t)=Rp(0)+1ΛAB(Λ)E(t,Rp)+ΛAB(Λ)Γ(Λ)t0E(ϱ,Rp)(tϱ)Λ1dϱHMp(t)=Mp(0)+1ΛAB(Λ)F(t,Mp)+ΛAB(Λ)Γ(Λ)t0F(ϱ,Mp)(tϱ)Λ1dϱ} (4.1)

    Hence we only consider the t=Δt(n+1)=tn+1,

    HSp(tn+1)=Sp(t0)+1ΛAB(Λ)A(t,Sp)+ΛAB(Λ)Γ(Λ)nσ=2tσ+1tσA(ϱ,Sp)(tn+1ϱ)Λ1dϱHEp(tn+1)=Ep(t0)+1ΛAB(Λ)B(t,Sp)+ΛAB(Λ)Γ(Λ)nσ=2tσ+1tσB(ϱ,Sp)(tn+1ϱ)Λ1dϱHIp(tn+1)=Ip(t0)+1ΛAB(Λ)C(t,Sp)+ΛAB(Λ)Γ(Λ)nσ=2tσ+1tσC(ϱ,Sp)(tn+1ϱ)Λ1dϱHAp(tn+1)=Ap(t0)+1ΛAB(Λ)D(t,Sp)+ΛAB(Λ)Γ(Λ)nσ=2tσ+1tσD(ϱ,Sp)(tn+1ϱ)Λ1dϱHRp(tn+1)=Rp(t0)+1ΛAB(Λ)E(t,Sp)+ΛAB(Λ)Γ(Λ)nσ=2tσ+1tσE(ϱ,Sp)(tn+1ϱ)Λ1dϱHMp(tn+1)=Mp(t0)+1ΛAB(Λ)F(t,Sp)+ΛAB(Λ)Γ(Λ)nσ=2tσ+1tσF(ϱ,Sp)(tn+1ϱ)Λ1dϱ} (4.2)

    Applying the method to our system, we get

    Sp(tn+1)=Sp(0)+1ΛAB(Λ)A(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)          +ΛAB(Λ)[(Δt)ΛΓ(Λ+1)nσ=2A(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                              ×[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[A(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                A(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[A(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2A(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +A(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]
    Ep(tn+1)=Ep(0)+1ΛAB(Λ)B(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)          +ΛAB(Λ)[(Δt)ΛΓ(Λ+1)nσ=2B(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                              ×[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[B(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                B(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[B(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2B(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +B(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]
    Ip(tn+1)=Ip(0)+1ΛAB(Λ)C(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)          +ΛAB(Λ)[(Δt)ΛΓ(Λ+1)nσ=2C(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                              ×[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[C(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                C(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[C(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2C(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +C(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]
    Ap(tn+1)=Ap(0)+1ΛAB(Λ)D(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)          +ΛAB(Λ)[(Δt)ΛΓ(Λ+1)nσ=2D(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                              ×[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[D(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                D(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[D(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2D(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +D(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]
    Rp(tn+1)=Rp(0)+1ΛAB(Λ)E(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)          +ΛAB(Λ)[(Δt)ΛΓ(Λ+1)nσ=2E(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                              ×[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[E(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                E(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[E(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2E(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +E(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]
    Mp(tn+1)=Mp(0)+1ΛAB(Λ)F(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)          +ΛAB(Λ)[(Δt)ΛΓ(Λ+1)nσ=2F(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                              ×[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[F(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                F(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[F(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2F(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +F(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]

    In this section, we presented the numerical solution of the suggested model by employing the newly established predictor-corrector numerical scheme. Since our equations are non-linear, we present the method for a general non-linear fractional differential equation with power-law and exponential kernal. We start with the following general equation.

    CF0DΛtv(t)=Λ(t,v(t))v(0)=v}

    Applying the associated integral operator, we obtain,

    v(t)v(0)=1ΛM(Λ)Λ(t,v(t))+ΛM(Λ)t0Λ(τ,v(τ))dτ

    at the point t=tn+1=Δt(n+1).

    Our equation becomes,

    v(tn+1)v(0)=1ΛM(Λ)Λ(tn+1,vpn+1)+ΛM(Λ)tn+10Λ(τ,v(τ))dτ.

    Then, between tn and tn+1, we approximate the function Λ(τ,v(τ)) using the Newton polynomial such that

    v(tn+1)v(tn)=1ΛM(Λ)[Λ(tn+1,vpn+1)Λ(tn,vn)]     +ΛM(Λ)tn+1tn[Λ(τn,vn1)+Λ(τn1,v(tn1))Λ(τn2,v(tn2))Δt(τtn2)     +Λ(τn,v(tn))2Λ(τn1,v(tn1))+Λ(τn2,v(tn2))2(Δt)2(τtn2)(τtn1)]dτ

    After integration, we obtain,

    v(tn+1)=v(tn)+1ΛM(Λ)[Λ(tn+1,vpn+1)Λ(tn,v(tn))]           +ΛM(Λ)[2312Λ(tn,v(tn)+512Λ(tn2,v(tn2)43Λ(tn1,v(tn1)]

    where, vpn+1=v0+(1Λ)ΛM(Λ)Λ(tn,v(tn))+ΛM(Λ)nk=0Λ(tk,v(tk))Δt.

    However, we need v1 to start the process, here,

    v(t1)=(1Λ)ΛM(Λ)Λ(t0,v(t0))+ΛM(Λ)Λ(t0,v(t0))Δt.

    To apply this method to our system, we transform the system to

    CF0DΛtSp(t)=A(t,Sp,Ep,Ip,Ap,Rp,Mp)CF0DΛtEp(t)=B(t,Sp,Ep,Ip,Ap,Rp,Mp)CF0DΛtIp(t)=C(t,Sp,Ep,Ip,Ap,Rp,Mp)CF0DΛtAp(t)=D(t,Sp,Ep,Ip,Ap,Rp,Mp)CF0DΛtRp(t)=E(t,Sp,Ep,Ip,Ap,Rp,Mp)CF0DΛtMp(t)=F(t,Sp,Ep,Ip,Ap,Rp,Mp)}

    Applying the method to our system, we get,

    Sp(tn+1)=Sp(tn)+(1Λ)ΛM(Λ)[A(tn+1,Spp,n+1,Epp,n+1,Ipp,n+1,App,n+1,Rpp,n+1,Mpp,n+1)            A(tn,Spp,n,Epp,n,Ipp,n,App,n,Rpp,n,Mpp,n)]            +(Λ)M(Λ)[2312A(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)            +512A(tn2,Sp,n2,Ep,n2,Ip,n2,Ap,n2,Rp,n2,Mp,n2)            43A(tn1,Sp,n1,Ep,n1,Ip,n1,Ap,n1,Rp,n1,Mp,n1)]
    Ep(tn+1)=Ep(tn)+(1Λ)ΛM(Λ)[B(tn+1,Spp,n+1,Epp,n+1,Ipp,n+1,App,n+1,Rpp,n+1,Mpp,n+1)            B(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)]            +(Λ)M(Λ)[2312B(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)            +512B(tn2,Sp,n2,Ep,n2,Ip,n2,Ap,n2,Rp,n2,Mp,n2)            43B(tn1,Sp,n1,Ep,n1,Ip,n1,Ap,n1,Rp,n1,Mp,n1)]
    Ip(tn+1)=Ip(tn)+(1Λ)ΛM(Λ)[C(tn+1,Spp,n+1,Epp,n+1,Ipp,n+1,App,n+1,Rpp,n+1,Mpp,n+1)            C(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)]            +(Λ)M(Λ)[2312C(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)            +512C(tn2,Sp,n2,Ep,n2,Ip,n2,Ap,n2,Rp,n2,Mp,n2)            43C(tn1,Sp,n1,Ep,n1,Ip,n1,Ap,n1,Rp,n1,Mp,n1)]
    Ap(tn+1)=Ap(tn)+(1Λ)ΛM(Λ)[D(tn+1,Spp,n+1,Epp,n+1,Ipp,n+1,App,n+1,Rpp,n+1,Mpp,n+1)            D(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)]            +(Λ)M(Λ)[2312D(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)            +512D(tn2,Sp,n2,Ep,n2,Ip,n2,Ap,n2,Rp,n2,Mp,n2)            43D(tn1,Sp,n1,Ep,n1,Ip,n1,Ap,n1,Rp,n1,Mp,n1)]
    Rp(tn+1)=Rp(tn)+(1Λ)ΛM(Λ)[E(tn+1,Spp,n+1,Epp,n+1,Ipp,n+1,App,n+1,Rpp,n+1,Mpp,n+1)            E(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)]            +(Λ)M(Λ)[2312E(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)            +512E(tn2,Sp,n2,Ep,n2,Ip,n2,Ap,n2,Rp,n2,Mp,n2)            43E(tn1,Sp,n1,Ep,n1,Ip,n1,Ap,n1,Rp,n1,Mp,n1)]
    Mp(tn+1)=Mp(tn)+(1Λ)ΛM(Λ)[F(tn+1,Spp,n+1,Epp,n+1,Ipp,n+1,App,n+1,Rpp,n+1,Mpp,n+1)            F(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)]            +(Λ)M(Λ)[2312F(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)            +512F(tn2,Sp,n2,Ep,n2,Ip,n2,Ap,n2,Rp,n2,Mp,n2)            43F(tn1,Sp,n1,Ep,n1,Ip,n1,Ap,n1,Rp,n1,Mp,n1)]

    where,

    Spp,n+1=Sp,0+(1Λ)ΛM(Λ)A(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)        +nk=0A(tk,Sp,k,Ep,k,Ip,k,Ap,k,Rp,k,Mp,k)Δt
    Epp,n+1=Ep,0+(1Λ)ΛM(Λ)B(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)        +nk=0A(tk,Sp,k,Ep,k,Ip,k,Ap,k,Rp,k,Mp,k)Δt
    Ipp,n+1=Ip,0+(1Λ)ΛM(Λ)C(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)        +nk=0C(tk,Sp,k,Ep,k,Ip,k,Ap,k,Rp,k,Mp,k)Δt
    App,n+1=Ap,0+(1Λ)ΛM(Λ)D(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)        +nk=0D(tk,Sp,k,Ep,k,Ip,k,Ap,k,Rp,k,Mp,k)Δt
    Rpp,n+1=Rp,0+(1Λ)ΛM(Λ)E(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)        +nk=0E(tk,Sp,k,Ep,k,Ip,k,Ap,k,Rp,k,Mp,k)Δt
    Mpp,n+1=Mp,0+(1Λ)ΛM(Λ)F(tn,Sp,n,Ep,n,Ip,n,Ap,n,Rp,n,Mp,n)        +nk=0F(tk,Sp,k,Ep,k,Ip,k,Ap,k,Rp,k,Mp,k)Δt

    And,

    Sp,1=Sp,0+(1Λ)ΛM(Λ)A(t0,Sp,0,Ep,0,Ip,0,Ap,0,Rp,0,Mp,0)        +ΛM(Λ)A(t0,Sp,0,Ep,0,Ip,0,Ap,0,Rp,0,Mp,0)
    Ep,1=Ep,0+B(t0,Sp,0,Ep,0,Ip,0,Ap,0,Rp,0,Mp,0)[(1Λ)ΛM(Λ)ΛΔtM(Λ)]
    Ip,1=Ip,0+C(t0,Sp,0,Ep,0,Ip,0,Ap,0,Rp,0,Mp,0)[(1Λ)ΛM(Λ)ΛΔtM(Λ)]
    Ap,1=Ap,0+D(t0,Sp,0,Ep,0,Ip,0,Ap,0,Rp,0,Mp,0)[(1Λ)ΛM(Λ)ΛΔtM(Λ)]
    Rp,1=Rp,0+E(t0,Sp,0,Ep,0,Ip,0,Ap,0,Rp,0,Mp,0)[(1Λ)ΛM(Λ)ΛΔtM(Λ)]
    Mp,1=Mp,0+F(t0,Sp,0,Ep,0,Ip,0,Ap,0,Rp,0,Mp,0)[(1Λ)ΛM(Λ)ΛΔtM(Λ)]

    Due to uncertainties associate with the etymology of the raise of infection, the initial time, within the Wuhan city, it is there for worth to include such uncertainty into a mathematical equation, additionally so far the world of research is very uncertain on when and how the curve of the spread will be flattened. To include such events into mathematical formulas, one should use a differential operator with a kernel that has similar properties, and the well-known differential operator with such is nothing else than the Caputo fractional derivative. In this section, we present the analysis with Caputo derivative

    C0DΛtSp(t)=A(t,Sp(t),Ep(t),Ip(t),Ap(t),Rp(t),Mp(t))C0DΛtEp(t)=B(t,Sp(t),Ep(t),Ip(t),Ap(t),Rp(t),Mp(t))C0DΛtIp(t)=C(t,Sp(t),Ep(t),Ip(t),Ap(t),Rp(t),Mp(t))C0DΛtAp(t)=D(t,Sp(t),Ep(t),Ip(t),Ap(t),Rp(t),Mp(t))C0DΛtRp(t)=E(t,Sp(t),Ep(t),Ip(t),Ap(t),Rp(t),Mp(t))C0DΛtMp(t)=F(t,Sp(t),Ep(t),Ip(t),Ap(t),Rp(t),Mp(t))}

    After integration, we have for the first equation

    Sp(t)=Sp(0)+1Γ(Λ)t0(tτ)Λ1A(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτ

    Hence we only consider the t=Δt(n+1)=tn+1,

    Sp(t+1)=Sp(0)+1Γ(Λ)nσ=0tσ+1tσ(tn+1τ)Λ1A(τ,Sp,Ep,Ip,Ap,Rp,Mp)dτ

    Approximating Ap(t)=Pp(t) within [tσ,tσ+1] following the work of Atangana and Seda[11], we obtain,

    Sp(tn+1)=Sp,0+(Δt)ΛΓ(Λ+1)nσ=2A(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[A(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                A(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[A(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2A(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +A(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]

    Also,

    Ep(tn+1)=Ep,0+(Δt)ΛΓ(Λ+1)nσ=2B(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[B(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                B(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[B(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2B(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +B(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]

    The same with Ip(tn+1),Ap(tn+1),Rp(tn+1) and Mp(tn+1).

    Ip(tn+1)=Ip,0+(Δt)ΛΓ(Λ+1)nσ=2C(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[C(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                C(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[C(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2C(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +C(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]
    Ap(tn+1)=Ap,0+(Δt)ΛΓ(Λ+1)nσ=2D(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[D(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                D(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[D(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2D(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +D(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]
    Rp(tn+1)=Rp,0+(Δt)ΛΓ(Λ+1)nσ=2E(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[E(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                E(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[E(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2E(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +E(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]
    Mp(tn+1)=Mp,0+(Δt)ΛΓ(Λ+1)nσ=2F(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)[(nσ+1)Λ(nσ)Λ]          +(Δt)ΛΓ(Λ+2)nσ=2[F(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                F(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ(nσ+3+2Λ)(nσ)Λ(nσ+3+3Λ)]       +(Δt)Λ2Γ(Λ+3)nσ=2[F(tσ,Sp,σ,Ep,σ,Ip,σ,Ap,σ,Rp,σ,Mp,σ)                                2F(tσ1,Sp,σ1,Ep,σ1,Ip,σ1,Ap,σ1,Rp,σ1,Mp,σ1)                                +F(tσ2,Sp,σ2,Ep,σ2,Ip,σ2,Ap,σ2,Rp,σ2,Mp,σ2)]                                  ×[(nσ+1)Λ{2(nσ)2+(3Λ+10)(nσ)+2Λ2+9Λ+12}                                  (nσ)Λ{2(nσ)2+(5Λ+10)(nσ)+6Λ2+18Λ+12}]

    In this section, we considered only Atangana-Baleanu sense and Caputo sense. We represent t numerical simulation of the suggested models using the software called Matlab. The numerical simulations are performed from the suggested numerical solutions obtained from the new suggested numerical schemes. The solutions are presented for different values of fractional orders including 1, 0.9, 0, 7, 0.5, 0.3 and 0.1 respectively. Starting with the model based on the Atangana-Baleanu fractional differential operator.

    Numerical simulation for Atangana-Baleanu based models are depicted in Figures 510 respectively.

    Figure 5.  Numerical solution with ABC derivative for Λ = 1.
    Figure 6.  Numerical solution with ABC derivative for Λ = 0.9.
    Figure 7.  Numerical solution with ABC derivative for Λ = 0.7.
    Figure 8.  Numerical solution with ABC derivative for Λ = 0.5.
    Figure 9.  Numerical solution with ABC derivative for Λ = 0.3.
    Figure 10.  Numerical solution with ABC derivative for Λ = 0.1.

    Numerical simulation for Caputo based models are depicted in Figures 1116 respectively.

    Figure 11.  Numerical solution with Caputo derivative for Λ = 1.
    Figure 12.  Numerical solution with Caputo derivative for Λ = 0.9.
    Figure 13.  Numerical solution with Caputo derivative for Λ = 0.7.
    Figure 14.  Numerical solution with Caputo derivative for Λ = 0.5.
    Figure 15.  Numerical solution with Caputo derivative for Λ = 0.3.
    Figure 16.  Numerical solution with Caputo derivative for Λ = 0.1.

    Extended orthogonal spaces are used, and valid fixed point results are obtained. Thereafter, we use fractional derivatives to provide a summary of the novel coronavirus 2019-nCoV/SARS-CoV-2 model's existence and unique solutions. Additionally, the model employed in this study is dubious because we did not include all information about the propagation into mathematical equations, and there are still various misunderstandings about the virus's behavior and spread. Our goal with this study was not to present a method to assist stop the spread, but rather to demonstrate that the model proposed by Altaf and Atangana[5] had unique solutions using a fixed-point technique using the initial conditions and theoretical values we utilized. In addition, to offer numerical solutions to the recommended models, a novel numerical scheme was applied. Most nonlinear systems require numerical simulations to analyze their behavior because their mathematical models are too complex to provide analytical answers. As a result of the numerical simulations produced, the virus is incredibly fast-spreading and can infect the entire world in a short period if precautions are not taken.

    The research of JJN has been partially supported by the Agencia Estatal de Investigacion (AEI) of Spain, projects MTM2016-75140-P and PID2020-113275GB-I00 and by the European Fund for Regional Development (FEDER). Also by Xunta de Galicia under grant ED431C 2019/02.

    The authors declare that they have no competing interests to influence the work reported in this paper.



    [1] Coronavirus disease (COVID-19) pandemic. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019.
    [2] S. Riddell, S. Goldie, A. Hill, D. Eagles, T. W. Drew, The effect of temperature on persistence of SARS-CoV-2 on common surfaces, Virology J., 17 (2020).
    [3] N. van Doremalen, D. H. Morris, M. G. Holbrook, A. Gamble, B. N. Williamson, A. Tamin, et al., Aerosol and surface stability of SARS-CoV-2 as compared with SARS-CoV-1, New Engl. J. Med., 382 (2020), 1564–1567. doi: 10.1056/NEJMc2004973
    [4] S. B. Kasloff, A. Leung, J. E. Strong, D. Funk, T. Cutts, Stability of SARS-CoV-2 on critical personal protective equipment, Sci. Rep., 11 (2021).
    [5] M. A. Khan, A. Atangana, Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative, Alexandria Eng. J., 59 (2020), 2379–2389. doi: 10.1016/j.aej.2020.02.033
    [6] A. Atangana, D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model, Therm. Sci., 20 (2016), 763-769. doi: 10.2298/TSCI160111018A
    [7] M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Prog. Fract. Differ. Appl., 1 (2015), 1–13.
    [8] J. Losada, J. J. Nieto, Properties of a new fractional derivative without singular kernel, Prog. Fract. Differ. Appl., 1 (2015), 87–92.
    [9] D. Kumar, F. Tchier, J. Singh, D. Baleanu, An efficient computational technique for fractal vehicular traffic flow, Entropy, 20 (2018), 259. doi: 10.3390/e20040259
    [10] E. F. D. Goufo, S. Kumar, S. B. Mugisha, Similarities in a fifth-order evolution equation with and with no singular kernel Chaos Solit. Fract., 130 (2020), 109467.
    [11] A. Atangana, S. I. Araz, RETRACTED: New numerical method for ordinary differential equations: Newton polynomial, J. Comput. Appl. Math., 372 (2020), 112622. doi: 10.1016/j.cam.2019.112622
    [12] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier Science B.V., Amsterdam, 2006.
    [13] C. Ravichandran, K. Logeswari, F. Jarad, New results on existence in the framework of Atangana-Baleanu derivative for fractional integro-differential equations, Chaos, Solitons Fractals, 125 (2019), 194–200. doi: 10.1016/j.chaos.2019.05.014
    [14] K. Logeswari, C. Ravichandran, A new exploration on existence of fractional neutral integro-differential equations in the concept of Atangana-Baleanu derivative, Phys. A: Stat. Mech. Appl., 544 (2020), 123454. doi: 10.1016/j.physa.2019.123454
    [15] M. A. Alqudah, C. Ravichandran, T. Abdeljawad, N. Valliammal, New results on Caputo fractional-order neutral differential inclusions without compactness, Adv. Differ. Equations, 2019 (2019).
    [16] R. Subashini, K. Jothimani, K. S. Nisar, C. Ravichandran, New results on nonlocal functional integro-differential equations via Hilfer fractional derivative, Alexandria Eng. J., 59 (2019), 2891–2899.
    [17] C. Ravichandran, K. Logeswari, S. K. Panda, K. S. Nisar, On new approach of fractional derivative by Mittag-Leffler kernel to neutral integro-differential systems with impulsive conditions, Chaos, Solitons Fractals, 139 (2020), 110012. doi: 10.1016/j.chaos.2020.110012
    [18] N. Valliammal, C. Ravichandran, K. S. Nisar, Solutions to fractional neutral delay differential nonlocal systems, Chaos, Solitons Fractals, 138 (2020), 109912. doi: 10.1016/j.chaos.2020.109912
    [19] D. Chalishajar, C. Ravichandran, S. Dhanalakshmi, R. Murugesu, Existence of fractional impulsive functional integro-differential equations in Banach spaces, Appl. Syst. Innovation, 2 (2019), 18. doi: 10.3390/asi2020018
    [20] S. K. Panda, Applying fixed point methods and fractional operators in the modelling of novel coronavirus 2019-nCoV/SARS-CoV-2, Results in Physics, 19 (2020), 103433. doi: 10.1016/j.rinp.2020.103433
    [21] A. Din, Y. Li, Stationary distribution extinction and optimal control for the stochastic hepatitis B epidemic model with partial immunity, Phys. Scr., 96 (2021), 074005. doi: 10.1088/1402-4896/abfacc
    [22] A. Din, Y. Li, Levy noise impact on a stochastic hepatitis B epidemic model under real statistical data and its fractal-fractional Atangana-Baleanunu order model, Phys. Scr., 96 (2021), 124008. doi: 10.1088/1402-4896/ac1c1a
    [23] A. Din, Y. Li, The extinction and persistence of a stochastic model of drinking alcohol, Results Phys., 28 (2021), 104649. doi: 10.1016/j.rinp.2021.104649
    [24] A. Din, Y. Li, F. M. Khan, Z. U. Khan, P. Liu, On analysis of fractional order mathematical model of Hepatitis B using Atangana-Baleanu Caputo (ABC) derivative, Fractals (2021), 2240017.
    [25] A. Din, Y. Li, A. Yusuf, A. I. Ali, Caputo type fractional operator applied to Hepatitis B system, Fractals, (2021), 2240023.
    [26] T. Kamran, M. Samreen, Q. U. Ain, A generalization of b-metric space and some fixed point theorems, Mathematics, 5 (2017).
    [27] S. K. Panda, E. Karapinar, A. Atangana, A numerical schemes and comparisons for fixed point results with applications to the solutions of Volterra integral equations in dislocated extended b-metric space, Alexandria Eng. J., 59 (2020), 815–827. doi: 10.1016/j.aej.2020.02.007
    [28] S. K. Panda, T. Abdeljawad, C. Ravichandran, A complex valued approach to the solutions of Riemann-Liouville integral, Atangana-Baleanu integral operator and non-linear Telegraph equation via fixed point method, Chaos, Solitons Fractals, 130 (2020), 109439. doi: 10.1016/j.chaos.2019.109439
    [29] N. Mlaiki, H. Aydi, N. Souayah, T. Abdeljawad, Controlled metric type spaces and the related contraction principle, Mathematics, 6 (2018), 194. doi: 10.3390/math6100194
    [30] S. K. Panda, T. Abdeljawad, C. Ravichandran, Novel fixed point approach to Atangana-Baleanu fractional and Lp-Fredholm integral equations, Alexandria Eng.J., 59 (2020), 1959–1970. doi: 10.1016/j.aej.2019.12.027
    [31] S. K. Panda, T. Abdeljawad, K. K. Swamy, New numerical scheme for solving integral equations via fixed point method using distinct (ωF)-contractions, Alexandria Eng. J., 59 (2020), 2015–2026. doi: 10.1016/j.aej.2019.12.034
    [32] T. Abdeljawad, R. P. Agarwal, E. Karapinar, P. S. Kumari, Solutions of the nonlinear integral equationand fractional differential equation using the technique of a fixed point with a numerical experiment in extended b-metric space, Symmetry, 11 (2019), 686. doi: 10.3390/sym11050686
    [33] E. Karapinar, P. S. Kumari, D. Lateef, A New Approach to the Solution of the Fredholm Integral Equation via a Fixed Point on Extended b-Metric Spaces, Symmetry, 10 (2018), 512. doi: 10.3390/sym10100512
    [34] P. S. Kumari, D. Panthi, Connecting various types of cyclic contractions and contractive self-mappings with Hardy-Rogers self-mappings, Fixed Point Theory Appl., 2016 (2016), 1–19. doi: 10.1186/s13663-015-0491-2
    [35] P. S. Kumari, I. R. Sarma, J. M. Rao, Metrization theorem for a weaker class of uniformities, Afrika Mat., 27 (2016), 667–672. doi: 10.1007/s13370-015-0369-9
    [36] Y. Zhou, Y. Hou, J. Shen, Y. Huang, W. Martin, F. Cheng, Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2, Cell Discovery, 6 (2020).
  • This article has been cited by:

    1. Hai-yang Xu, Heng-you Lan, Fan Zhang, General semi-implicit approximations with errors for common fixed points of nonexpansive-type operators and applications to Stampacchia variational inequality, 2022, 41, 2238-3603, 10.1007/s40314-022-01890-7
    2. Yan Wang, Rui Wu, Shanshan Gao, The Existence Theorems of Fractional Differential Equation and Fractional Differential Inclusion with Affine Periodic Boundary Value Conditions, 2023, 15, 2073-8994, 526, 10.3390/sym15020526
    3. Sumati Kumari Panda, Thabet Abdeljawad, Fahd Jarad, Chaotic attractors and fixed point methods in piecewise fractional derivatives and multi-term fractional delay differential equations, 2023, 46, 22113797, 106313, 10.1016/j.rinp.2023.106313
    4. Philip N. A. Akuka, Baba Seidu, C. S. Bornaa, Marko Gosak, Mathematical Analysis of COVID-19 Transmission Dynamics Model in Ghana with Double-Dose Vaccination and Quarantine, 2022, 2022, 1748-6718, 1, 10.1155/2022/7493087
    5. SUMATI KUMARI PANDA, ABDON ATANGANA, THABET ABDELJAWAD, EXISTENCE RESULTS AND NUMERICAL STUDY ON NOVEL CORONAVIRUS 2019-NCOV/ SARS-COV-2 MODEL USING DIFFERENTIAL OPERATORS BASED ON THE GENERALIZED MITTAG-LEFFLER KERNEL AND FIXED POINTS, 2022, 30, 0218-348X, 10.1142/S0218348X22402149
    6. Tahair Rasham, Separate families of fuzzy dominated nonlinear operators with applications, 2024, 70, 1598-5865, 4271, 10.1007/s12190-024-02133-0
    7. Sumati Kumari Panda, Thabet Abdeljawad, A. M. Nagy, On uniform stability and numerical simulations of complex valued neural networks involving generalized Caputo fractional order, 2024, 14, 2045-2322, 10.1038/s41598-024-53670-4
    8. Sumati Kumari Panda, Kumara Swamy Kalla, A.M. Nagy, Limmaka Priyanka, Numerical simulations and complex valued fractional order neural networks via (ɛ – μ)-uniformly contractive mappings, 2023, 173, 09600779, 113738, 10.1016/j.chaos.2023.113738
    9. Hammed Anuoluwapo Abass, Maggie Aphane, Morufu Oyedunsi Olayiwola, An inertial method for solving systems of generalized mixed equilibrium and fixed point problems in reflexive Banach spaces, 2023, 16, 1793-5571, 10.1142/S1793557123502078
    10. Yassine Adjabi, Fahd Jarad, Mokhtar Bouloudene, Sumati Kumari Panda, Revisiting generalized Caputo derivatives in the context of two-point boundary value problems with the p-Laplacian operator at resonance, 2023, 2023, 1687-2770, 10.1186/s13661-023-01751-0
    11. Nihar Kumar Mahato, Bodigiri Sai Gopinadh, , 2024, Chapter 15, 978-981-99-9545-5, 339, 10.1007/978-981-99-9546-2_15
    12. Sumati Kumari Panda, Velusamy Vijayakumar, Bodigiri Sai Gopinadh, Fahd Jarad, 2024, Chapter 6, 978-981-99-9545-5, 177, 10.1007/978-981-99-9546-2_6
    13. Sumati Kumari Panda, Ilyas Khan, Vijayakumar Velusamy, Shafiullah Niazai, Enhancing automic and optimal control systems through graphical structures, 2024, 14, 2045-2322, 10.1038/s41598-024-53244-4
    14. P K Santra, G S Mahapatra, Sanjoy Basu, Stability analysis of fractional epidemic model for two infected classes incorporating hospitalization impact, 2024, 99, 0031-8949, 065237, 10.1088/1402-4896/ad4692
    15. Anupam Das, Bipan Hazarika, Mohsen Rabbani, 2024, Chapter 10, 978-981-99-9206-5, 165, 10.1007/978-981-99-9207-2_10
    16. Tahair Rasham, Sumati Kumari Panda, Ghada Ali Basendwah, Aftab Hussain, Multivalued nonlinear dominated mappings on a closed ball and associated numerical illustrations with applications to nonlinear integral and fractional operators, 2024, 10, 24058440, e34078, 10.1016/j.heliyon.2024.e34078
    17. Sumati Kumari Panda, Velusamy Vijayakumar, Kottakkaran Sooppy Nisar, Applying periodic and anti-periodic boundary conditions in existence results of fractional differential equations via nonlinear contractive mappings, 2023, 2023, 1687-2770, 10.1186/s13661-023-01778-3
    18. Sumati Kumari Panda, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar, Bipan Hazarika, Solving existence results in multi-term fractional differential equations via fixed points, 2023, 51, 22113797, 106612, 10.1016/j.rinp.2023.106612
    19. Sumati Kumari Panda, Velusamy Vijayakumar, Results on finite time stability of various fractional order systems, 2023, 174, 09600779, 113906, 10.1016/j.chaos.2023.113906
    20. Sumati Kumari Panda, A.M. Nagy, Velusamy Vijayakumar, Bipan Hazarika, Stability analysis for complex-valued neural networks with fractional order, 2023, 175, 09600779, 114045, 10.1016/j.chaos.2023.114045
    21. Mokhtar Bouloudene, Fahd Jarad, Yassine Adjabi, Sumati Kumari Panda, Quasilinear Coupled System in the Frame of Nonsingular ABC-Derivatives with p-Laplacian Operator at Resonance, 2024, 23, 1575-5460, 10.1007/s12346-023-00902-z
    22. Sumati Kumari Panda, Vijayakumar Velusamy, Ilyas Khan, Shafiullah Niazai, Computation and convergence of fixed-point with an RLC-electric circuit model in an extended b-suprametric space, 2024, 14, 2045-2322, 10.1038/s41598-024-59859-x
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3109) PDF downloads(100) Cited by(22)

Figures and Tables

Figures(16)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog