key: cord-0910860-3x15oevv
authors: Shah, Kamal; Sher, Muhammad; Rabai'ah, Hussam; Ahmadian, Ali; Salahshour, Soheil; Pansera, Bruno A.
title: Analytical and qualitative investigation of COVID‐19 mathematical model under fractional differential operator
date: 2021-08-22
journal: Math Methods Appl Sci
DOI: 10.1002/mma.7704
sha: 72df3fdb3e3db369f731f35c41bb758106864990
doc_id: 910860
cord_uid: 3x15oevv

In the current article, we aim to study in detail a novel coronavirus (2019‐nCoV or COVID‐19) mathematical model for different aspects under Caputo fractional derivative. First, from analysis point of view, existence is necessary to be investigated for any applied problem. Therefore, we used fixed point theorem's due to Banach's and Schaefer's to establish some sufficient results regarding existence and uniqueness of the solution to the proposed model. On the other hand, stability is important in respect of approximate solution, so we have developed condition sufficient for the stability of Ulam–Hyers and their different types for the considered system. In addition, the model has also been considered for semianalytical solution via Laplace Adomian decomposition method (LADM). On Matlab, by taking some real data about Pakistan, we graph the obtained results. In the last of the manuscript, a detail discussion and brief conclusion are provided.

In the last of 2019, a dangerous out break was confirmed in the big city of China called Wuhan which was due by coronavirus. There are many theories behind the virus outbreak in present-day literature. Some researchers have studied that it originated from bats to humans owing to the complete unlawful transmission of animals on a Wuhan seafood market. The virus involved has been identified both in dogs and in pangolin. Hence, it has been considered that several infected cases claimed that they had been working in a local fish and wild animal market in Wuhan from where they got the infection of COVID-19. Afterward, the researchers have confirmed that the disease's widespread nature is caused by direct contact between people. The mentioned disease was then named COVID-19. Nearly 100 000 people were infected in China up to March 2020. Then, the virus was transmitted from person to person very rapidly and spread up to the end of March 2020 in the whole globe. This pandemic is ongoing, and more than 7.38 million people almost in every country of the globe have been infected. Nearly 0.41 million people died from this pandemic. WHO announced that it is a global pandemic where scientists and researchers have not found success in controlling or curing this disease. 1, 2 Here, we remark that one way to open up this disease is due to immigration of the infected population from one area to another which affects the healthy population and spreads this disease very rapidly. In fact, it is like the flue which has killed millions people in past. Therefore, globally, all the countries of the entire world have reduced their activities to some restricted domain. Every country of the world restricts the movement of the people from one place to other. 3 Currently, researchers work on this outbreak. 4, 5 Various models in mathematical form have been developed for these diseases. For some study, we refer previous works [6] [7] [8] and work cited therein. All the diseases are dangerous, but infectious diseases greatly effect the health as well as economy of a state. In the past for many diseases, controlling procedure and its future predictions are made by using mathematical modeling. 9 Considering the threat of current pandemic, researchers have studied the COVID-19 disease from different aspects. [10] [11] [12] [13] [14] Mathematical models in this regard can play a vital role to restrict the disease spreading. For this need very recently, some models have been studied for COVID-19; see detail in other works. [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] It is often desirable to describe the behavior of some real-life system or phenomenon, whether physical, sociological, or even economic, in mathematical terms. The mathematical description of a system of phenomenon is called a mathematical model and is constructed with certain goals in mind. Put it simply, mathematical modeling should become part of the toolbox of public health research and decision making. One way to study problem dynamics in the real world is to use mathematical modeling. Mathematical modeling is the ability to convert real-world problems into mathematical formulations whose theoretical and numerical analysis can provide useful insights and guidance for native applications. In the field of mathematics, many mathematical models have been built according to the assumptions of researchers. The purpose of mathematical modeling is to represent different types of real world situation in the language of mathematics. To find out the different dynamics of a disease and therefore to overcome it at an early stage, mathematical modeling plays an important role there. Thus simply, we can say that mathematical modeling has important role to study a disease and control of the disease specially when vaccination is absent. A lot of mathematical models exist in the literature on stability, existence theory, and optimization of biological models; for example, see (other works [25] [26] [27] ). Currently, for COVID-19, some models have been constructed to study its different aspects (see Lai et al. and Fanelli and Piazza 28, 29 ) . In this way, Lin and his coauthors in their study 15 model COVID-19 under integer-order derivative as

is the susceptible populations, E(t) is exposed populations, I(t) is infectious populations, R(t) is removed population (recovered and dead), N(t) is total populations, D(t) is mimicking the public perception of risk regarding the number of severe and critical cases and deaths, and C(t) is the cumulative cases (reported and not reported). The authors establish this model on Wuhan city for the month of March 2020.

Detail of various parameters of model (1) Public reaction mean duration

As we know that most of biological models are based on the classical approach to get a system of nonlinear autonomous first-order differential equations. Therefore, there is still room to develop such a mathematical models by the tools of advanced fractional calculus. Several novel studies, including numerous biological models, have shown to be extremely helpful and more accurate than their counterparts. For example, Ullah et al 30 proposed a tuberculosis (TB) infection mathematical model for the Khyber Pakhtunkhwa province of Pakistan involving Caputo fractional-order derivative and tested it by the real data of the mentioned province and shows the advantages of fractional-order model. Khan et al 31 introduced pine wilt disease model involving fractional-order derivative of Caputo-Fabrizio type and confirmed its productiveness by establishing its unique solution. Qureshi and Atangana 32 analyzed ordinary and fractional-order models of dengue outbreak in Cape Verde islands during 2009 and showed that fractional-order operator in the sense of Caputo-Fabrizio has the smallest squared sum of errors. Muhammad and Atangana 33 analyzed the mathematical model of infectious disease called Ebola via Caputo, Caputo-Fabrizio, and Atanagana-Baleanu operators. The advanced studies show that the fractal-fractional-order operators are the best tools to analyze the mathematical models for real-world data; for detail, see previous studies. [34] [35] [36] The beauty of fractional calculus is that fractional derivatives (and integrals) are not a local (or point) property (or quantity). Thereby, this considers the history and nonlocal distributed effects. In other words, perhaps this subject translates the reality of nature better. Therefore, to make this subject available as popular subject to science and engineering community, it adds another dimension to understand or describe basic nature in a better way.

Existence theory is one of the interesting areas of differential equations. This area has been studied very well in the last couple of decades. In the present literature, various approaches were utilized to show existence and uniqueness of solution of differential and integral equations. Krassnoselskii's fixed point theorem, Lery-Schauder fixed point theorem, Schaefer's fixed point theorem, and degree theory are commonly used to investigate the solutions for fractional-order differential equations; see other studies. [37] [38] [39] [40] [41] [42] Because most of nonlinear problems cannot be solved for exact solution, we need powerful numerical or some analytical techniques. For good numerical results, one needs stable algorithms and methods. For such need, stability theory was founded. In the literature, there are different types of stability. But the most important type is Ulam-Hyers type stability which was introduced by Ulam 43 in 1940 and further studied by Hyers 44 in 1941. This stability answers the question when is it true that a function which approximately satisfies a functional equation must be close to an exact solution of this functional equation? If the problem accepts a solution, we say that the functional equation is stable. The mentioned stabilities have been studied for different physical problem; for example, see previous works. 42, 45, 46 Inspired of the work of Lin et al 15 and realistic nature of fractional calculus, we consider model (1) for existence, Ulam-Hyers stability, and for semianalytical solution under Caputo fractional-order derivative as

Initial conditions assume as follows:

First, we will develop qualitative theory for the above model, then we will discuss Ulam-Hyers stability and, finally, an approximate solution of the considered model (2) to be discussed. In our analysis, we used "Banach and Schaefer's fixed point theorem" to study the existence and stability of the mentioned model. At last, a general scheme is established to obtain series type solution. Also, we testify the approximate results for the model (2) on the real data of Pakistan for the last 60 days taken from a source. 47 At the end, graphical results are given with detail discussion.

Some basic results are recall as follows.

Kilbas et al. 27 The fractional integral of order ∈ R + , for a function W(t) ∈ L 1 ([a, b], R) is defined as

where 0 < ≤1 and integral on the right is assumed to converge.

Kilbas et al. 27 The Caputo fractional-order derivative of a function W(t) ∈ C n [a, b] is defined as

here, 0 < ≤1, n = [ ] + 1, and [ ] shows integer part of .

is given by

for c i ∈ R, i=0,1,2, … ,n−1.

Podlubny. 48 The Laplace transform of Caputo fractional-order derivative is defined as follows:

For a Banach's space

we have the following results.

is bounded, then T has a fixed point.

First of all, we will consider model (2) for qualitative theory. For this, we put the right hand side of model (2) as follows:

Using Equation (9), model (2) can be expressed as follows:

where 

On Lemma 2.3, Equation (10) can be converted to an equivalent integral form as follows:

The following assumptions are key to our analysis:

(H2) If there exist constants L Q > 0, and for each W,W ∈ X, such that

Define the map T : X → X as follows: Proof. With the help of Schaefer's fixed point theorem, the proof of the results will be present in four steps.

Step 1: First, we have to show T is continuous. Assume W i is continuous for i = 1, 2, 3, .., 7. Thus, Q(t, W(t)) is continuous. Let W j , W ∈ X such that W j → W, and we must have TW j → TW.

Consider

Since Q is continuous, thus TW j → TW; as a result ,T is continuous.

Step 2: The operator T is bounded. For any W ∈ X, the function T enjoys the growth condition as follows:

Let S be any bounded subset of X. We need to show T(S) is bounded. For any W ∈ S, since S is bounded, therefore, there exists K≥0, such that

Hence, for any W ∈ S using above growth conditions, one has

Hence, T(S) is bounded.

Step 3: For equi-continuity, let t 1 , t 2 ∈ [0, ], such that t 1 ≥t 2 , then

Hence, by Arzelá-Ascoli theorem, T(S) is relatively compact.

Step 4: Finally, to show that the set

is bounded, let W ∈ E, then for each t ∈ [0, ], we have

which shows E is bounded. Hence, by Schaefer's theorem, T possesses a fixed point; consequently, our considered problem (10) has at least one solution. Proof. Thanks to "Banach's contraction theorem," let W,W ∈ X, then

Hence, T has unique fixed point. Therefore, our problem (10) has unique solution.

Now in this part of the manuscript, we will discuss the Ulam-Hyers type stability which was introduced by Ulam, 43 in 1940 and studied further by Hyers. 44 The stated stability generalized further by Rassiass, 50 to Ulam-Hyer-Rassiass stability. For the last two decades this stability was studied very well. 45, 46, [50] [51] [52] First, we will present different definition of Ulam-Hyers type stability which can be found in Ali et al. 52 Let H : X → X be an operator satisfying HW = W, for W ∈ X.

Definition 4.1. Problem (16) is Ulam-Hyers type stable for > 0 and assume W ∈ X represents any solution of

There exists unique solutionW of problem (16) such that C q > 0 satisfying

For Υ ∈ C(R, R) with Υ(0) = 0, for any solution W of (17), and letW be at most one solution of (16) such that

then problem (16) is generalized Ulam-Hyers type stable.

The following relation is needed in future work. The perturb problem corresponding to Equation (10) is given by Proof. For any solution W ∈ X of problem (20) andW ∈ X represents unique solution to problem (10), then

Thus, problem (10) is Ulam-Hyers and generalized Ulam-Hyers type stable. 

there exists unique solutionW to problem (16) with C q > 0 satisfying

For Ω ∈ C[[0, ], R] if there exist C q, Ω and for > 0, let W represent any solution to problem (23) and W represent unique solution to problem (16) such that

then problem (16) is generalized Ulam-Hyers-Rassias type stable. 

Proof. Using Remark (4.8) and Lemma 2.3, one may get the required result.

Proof. Let W ∈ X be any solution andW ∈ X represent unique solution to problem (10), then

Thus, problem (10) is Ulam-Hyers-Rassias stable and generalized Ulam-Hyers-Rassias type stable. 

In this part of the article, a general algorithm for series solution to model (2) is developed. Using Laplace transform, we have Assuming solution in the form of series as

Decomposing the nonlinear terms S(t)I(t) by Adomian polynomials as By system (29) and (30), system (28) gives

Compare like terms on each sides of (31), we get

,

, For j≥0, we generalized the terms as

,

,

By application of inverse Laplace transform to systems (32)-(35), we have

] ,

] ,

] ,

] ,

] ,

] , Thus, for j≥0, the general term is given by

Similarly, we obtain series solution as

For justification of developed results, values have been taken about COVID-19 in Pakistan and also presented in Table  1 47 :

Case-I: G = 0, u = 0.02205∕da , a = 0.4239∕da , b 0 = 0.5944∕da In Figures 1-7 , we draw the graph of the solution for the given values of parameter defined as above. Using Matlab-16, the series solutions defined in Equation (40) have been plotted. From Figure 1 , we see that the susceptible population was decreasing for the initial population at given time in different orders. In same line as in Figure 2 , exposed class was increasing with the decrease of susceptive class, and hence, the number of infection is increasing as in Figure 3 because people are exposed to infection. This was due to the fact that people were not taking the matter seriously in Pakistan so the number of infections grew with proper speed. Initially, in first 15-20 days, the numbers of recovered class were increasing as people initially were getting ride easily due to the reason that cases were limited and treatment was proper as in Figure 4 . After that, the recovered rate is becoming slow now as numbers of infection increased, and also there is no proper cure yet introduced. Also now the infection has attacked mostly on aged people whose recovery will take time. From Figure 5 , we see that total population of healthy people is gradually decreasing with different rate due to fractional order derivative. Similarly, in Figure 6 , the people initially not taking the matter seriously but now the people are taking it seriously and so the Mimicking class is decreasing with different rates because of fractional order. The class of cumulative cases is also increasing with different rate as in Figure 7 . From Figures 1 to 7 , we see that fractional calculus provides a global dynamics of the novel coronavirus infection model. The smaller the fractional order, the faster the concerned decay or growth and vice versa. Case-II: G = 100, u = 0.000002205∕da , a = 0.0004239∕da , b 0 = 0.00003345∕da . In Figures 8-14 , we draw the graph of the solution for the given values of parameter defined above. Using Matlab-16, the series solutions defined in Equation (40) have been plotted.

From Figure 8 , we see that as susceptible papulation is decreasing in the presence of zoonotic cases and immigration for the given initial papulation at given time at different order, then in the same line as in Figure 9 , exposed class was also decreasing with the decrease of susceptive class, and hence, the number of infection was increasing as in Figure 10 , because more people were exposed to catch infection. This is due to the reason that people first were not taking the matter serious in Pakistan so the number of infection rate raised up with proper fast speed. Initially, in first 15-20 days, the numbers of recovered class were increasing as people initially were getting ride easily due to the reason that cases were limited and treatment was proper as in Figure 11 . The recovery rate is also gradually increased due to the increase of people getting ride from infection or going to die. From Figure 12 , we see that total papulation of healthy people is gradually decreasing with different rate due to fractional order derivative. Similarly, in Figure 13 , the people initially not taking the matter seriously but now the people are taking it seriously and so the Mimicking class is decreasing with different rates because of fractional order. The class of cumulative cases is also increasing with different rates as in Figure 14 . From  Figures 8-14 , we see that fractional calculus provides a global dynamics of the novel coronavirus-19 infection model. The smaller the fractional order, the faster the concerned decay or growth and vice versa.

We have examined a novel model of coronavirus-19 under fractional order derivative. We first have established the existence theory of the model along with stability results of Ulam type via the use of nonlinear analysis. After that by a coupled method of Laplace transform and Adomian decomposition, we have established a general algorithm for series type solution to the considered model. By using Matlab, we have plotted the graphs against some real data of Pakistan, we observed that permitting immigration and not taking the matter seriously by the people in first sixty days the infection rate has been raised roughly in the country. But now the people are taking the matter seriously, but in the previous sixty days proper, great numbers of people have been coughing throughout the country. By fractional calculus approach, we have examined the global dynamics of the current novel disease. We see that fractional calculus approach globally describes the dynamics of all the compartments of the considered model more comprehensively.

We are thankful to the reviewers for useful comments which have improved this work very well. No funding source is available.

There exists no competing interest regarding this research work.

Authors have equal contribution in this paper.

World Health Organization

Situation Report -72

Outbreak of pneumonia of unknown etiology in Wuhan, China: the mystery and the miracle

Data model to predict prevalence of COVID-19 in Pakistan

Analysis of the worldwide corona virus (COVID-19) pandemic trend; a modelling study to predict its spread

An efficient technique for a time fractional model of lassa hemorrhagic fever spreading in pregnant women

A new fractional SIRS-SI malaria disease model with application of vaccines, antimalarial drugs, and spraying

Semi-analytical study of Pine Wilt Disease model with convex rate under Caputo-Febrizio fractional order derivative

Mathematical Biology: An Introduction

Computational and theoretical modeling of the transmission dynamics of novel COVID-19 under Mittag-Leffler power law

Potent binding of 2019 novel coronavirus spike protein by a SARS coronavirus-specific human monoclonal antibody

Analysis of Caputo fractional-order model for COVID-19 with lockdown

Non pharmaceutical interventions for optimal control of COVID-19

A mathematical model of Coronavirus Disease (COVID-19) containing asymptomatic and symptomatic classes

A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action

An epidemic prediction from analysis of a combined HIV-COVID-19 co-infection model via ABC-fractional operator

COVID-19 Infected Lung Computed Tomography Segmentation and Supervised Classification Approach. CMC-COMPU Mate Conti

A stochastic numerical analysis based on hybrid NAR-RBFs networks nonlinear SITR model for novel COVID-19 dynamics

On nonlinear classical and fractional order dynamical system addressing COVID-19

Investigation of fractal-fractional order model of COVID-19 in Pakistan under Atangana-Baleanu Caputo (ABC) derivative

GraphCovidNet: A graph neural network based model for detecting COVID-19 from CT scans and X-rays of chest

An emotion care model using multimodal textual analysis on COVID-19

Prediction modelling of COVID using machine learning methods from B-cell dataset

Identification of dominant risk factor involved in spread of COVID-19 using hesitant fuzzy MCDM methodology

A fractional calculus approach to the dynamic optimization of biological reactive systems. Part I: Fractional models for biological reactions

An introduction to the fractional calculus and fractional differential equations

Theory and Applications of Fractional Differential Equations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges

Analysis and forecast of COVID-19 spreading in China, Italy and France

A fractional model for the dynamics of TB virus

A fractional order pine wilt disease model with Caputo-Fabrizio derivative

Mathematical analysis of dengue fever outbreak by novel fractional operators with field data

Dynamics of Ebola disease in the framework of different fractional derivatives

Modeling and analysis of the polluted lakes system with various fractional approaches

Analysis and numerical simulation of fractional model of bank data with fractal-fractional Atangana-Baleanu derivative

Fractional investigation of bank data with fractal-fractional Caputo derivative

A nonlinear alternative in Banach algebras with applications to functional differential equations

Fixed point theorems in ordered Banach algebras and applications

A fixed point theorem in Banach algebras involving three operators with applications

Qualitative analysis of multi-terms fractional order delay differential equations via the topological degree theory

On qualitative theory of fractional order delay evolution equation via the prior estimate method

Study of evolution problem under Mittag-Leffler type fractional order derivative

Problems in Modern Mathematics (Chapter VI)

On the stability of the linear functional equation

Applicability of topological degree theory to evolution equation with proportional delay

Study of time fractional order problems with proportional delay and controllability term via fixed point approach

Official Updates Coronavirus COVID-19 in Pakistan: Ministry of National Health Services Regulations & Coordination

Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications

Fixed Point Theory

On the stability of the linear mapping in Banach spaces

On Ulam's stability for a coupled systems of nonlinear implicit fractional differential equations

Ulam stability to a toppled systems of nonlinear implicit fractional order boundary value problem. Bound Value Prob

Analytical and qualitative investigation of COVID-19 mathematical model under fractional differential operator

ORCID Kamal Shah https://orcid.org/0000-0002-8851-4844 Hussam Rabai'ah https://orcid.org/0000-0003-4597-5787 Ali Ahmadian https://orcid.org/0000-0002-0106-7050