{'fraction_non_alphanumeric': 0.05074608071523012, 'fraction_numerical': 0.027135931499087076, 'mean_word_length': 4.291139240506329, 'pattern_counts': {'":': 0, '<': 10, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 7, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A fractional‐order mathematical model for analyzing the pandemic trend of COVID‐19

Many countries worldwide have been affected by the outbreak of the novel coronavirus (COVID‐19) that was first reported in China. To understand and forecast the transmission dynamics of this disease, fractional‐order derivative‐based modeling can be beneficial. We propose in this paper a fractional‐order mathematical model to examine the COVID‐19 disease outbreak. This model outlines the multiple mechanisms of transmission within the dynamics of infection. The basic reproduction number and the equilibrium points are calculated from the model to assess the transmissibility of the COVID‐19. Sensitivity analysis is discussed to explain the significance of the epidemic parameters. The existence and uniqueness of the solution to the proposed model have been proven using the fixed‐point theorem and by helping the Arzela–Ascoli theorem. Using the predictor–corrector algorithm, we approximated the solution of the proposed model. The results obtained are represented by using figures that illustrate the behavior of the predicted model classes. Finally, the study of the stability of the numerical method is carried out using some results and primary lemmas.


INTRODUCTION
The WHO announced COVID-19 to be a pandemic epidemic on January 22, 2020. Nowadays, the pandemic is spreading across the world and affects almost every aspect of life. In addition to health problems, it undermines the global economic system and restricting people's contact. Researchers from diverse scientific fields have therefore committed Modeling and studied the pandemic trend of COVID-19 has been presented with fractional-order derivative. Using the generation matrix method, the basic reproduction number is calculated that located whether the disease would persist or disappears from the population. The equilibrium points for this system are calculated, where the graphical solutions show that the results of the model converge to their equilibrium points. Based on the derivatives of the basic reproduction number, the sensitivity of the parameters was analyzed. The existence and uniqueness of the solution to the proposed model have been proven using the fixed-point theorem and by helping the Arzela-Ascoli theorem.
themselves to the study of COVID-19. The objectives were to contribute to the enhancement of the comprehension, forecast, and interpretation from different points of view for this disease. Any new suggestions are a step forward in solving this health crisis. Mathematical models are a significant and efficient way to comprehend epidemic transmission dynamics. [1][2][3][4][5][6][7][8][9] These models can be helping predict disease transmission and thus help decision makers in planning and make necessary decisions.
The fractional-order models provide a more accurate fit than the integer-order models for the actual data for various diseases and other experiential studies around modeling and simulation. New mathematical models of fractional order are developed that can be used for simulation to forecast the disease outbreaks and flatten the infection and deaths curve. [10][11][12][13][14] The mathematicians have also put their efforts forward in the analysis of various nonlinear dynamics of problems related to the infection like epidemics. [15][16][17][18] Khan and Atangana 19 described the mathematical modeling and dynamics of a novel coronavirus (2019-nCoV), then formulated a new mathematical model for simulation of the dynamics COVID-19 with quarantine and isolation. 20 Atangana 21 proposed a model for the spread of COVID-19, with new fractal-fractional operators and taking into account the potential of transmission of COVID-19 from dead bodies to humans and the effect of lockdown. Yadav 27,28 analyzed the fractional-order COVID-19 model for simulating the potential transmission of epidemic and suggested some possible control strategies. Zhang et al 29 proposed and applied the time-dependent SEIR model to fit the time series of coronavirus evolution observed recently in China and then predict transmission. Matouk 30 analyzed complex dynamics in a model for COVID-19 and its fractional-order counterpart with multidrug resistance. Tang et al 31 adopted a deterministic model that was devised based on the clinical improvement to illuminate the dynamics for transmission of the novel coronavirus and determine the effect of the public health measures on the infection.
The paper is arranged as follows: In Section 2, we recalled some basic concepts and results used in the article. In Section 3, we developed the model pandemic trend Covid-19 with fractional order then calculated the basic reproduction number and equilibrium points along with sensitivity analysis for parameters. Using the fixed-point theory, the existence and uniqueness analysis of the model was performed in Section 4. In Section 5, we presented a numerical algorithm for solving the model based on the predictor-corrector method. In Section 6, we discussed the stability of the numerical method. In Section 7, we have displayed the numerical simulation results for the real data. A conclusion completes the paper.

BASIC CONCEPT
In this section, we present some of the fundamental definitions and rand results of fractional order derivatives which, can be used throughout this manuscript. Definition 1 (Samko et al. and Miller and Ross 32,33 ). For an integrable function f(t), the Caputo derivative of fractional order ∈ (0, 1) is given by The corresponding fractional integer of order with Re( ) > 0 is given by Lemma 1 (Li and Zeng 34 ). If 0 < < 1 and n≥0 is an integer, then there exists the positive constants ∁ ,1 and ∁ , 2 only dependent on , such that and (n + 2) +1 − 2(n + 1) +1 + n +1 ≤ ∁ ,2 (n + 1) −1 .
Lemma 3 (Diethelm 35 ). The function y ∈ C[0, T] is a solution of the following fractional differential equation: if and only if it is a solution of the nonlinear Volterra integral equation of the second kind

THEORETICAL APPROACH
In this section, we are considered to examine a classical model developed by Fatima et al 1 for the pandemic trend of 2019 coronavirus, which was first identified in the Chinese city of Wuhan in December 2019 and then spread quickly across the world. The population is divided into six subcategories: susceptible people denoted by S, exposed people E, infected people I, asymptomatic infectious people A, isolated or hospitalized people H, recovered or dead people R, and the reservoir for COVID-19 is denoted by W. The pandemic trend model of COVID-19 is given by the system of ordinary differential equations as follows: Subject to non-negative initial conditions.
where Υ = b × N, b is the birth rate, and N is the total number of people; the death rate of the population; 1 and 2 denote the rate of transmission infected of the susceptible people through sufficient contact with I and W, respectively. The parameters and denote the approximate transmissibility from A and H to I, respectively. ( ) −1 represented the transmission rate after the quarantine period and joined the class I and A. The parameters and (1 − ) are the moving rate from E to I and A, respectively. The rate of infected people who are hospitalized is 1 , and 2 is the recovery rate. 3 denotes the recovery rate of hospitalized patient. 1 and 2 are the infected and asymptomatically infected populations contributing to the reservoirs coronavirus, respectively. is the lifetime of the virus. Now, we modify Equations (9) and (10) by putting the Caputo fractional derivatives instead of the time derivatives; with this move, the dimensions on both sides can differ. To avoid this trouble, we use an auxiliary parameter , which has a sec. dimension to modify the fractional operator to have the same dimension for the sides. Equations (9) and (10) Subject to non-negative initial conditions

Equilibrium points and basic reproduction number
To find the disease-free equilibrium of the suggested fractional-order model (11), we solve the following system.
We get P 0 = (Υ∕ , 0, 0, 0, 0, 0, 0), which is the point where there is no disease. If the basic reproduction number  0 > 0, then the system (13) has a positive endemic equilibrium point, which is represented by , , , A vital indicator of the spread of infectious disease in the population is the basic reproduction number  0 , defined as the number of new infections caused by a typical infectious person in a disease-free equilibrium population. An epidemic will spread if  0 > 1. While that an outbreak will most likely not accrue if  0 < 1. To obtain  0 for the system (13), we use the generation method. 36 The infection components in the model system are E, I, A, H, and W, so we consider the following fractional system: where At point P 0 , the Jacobian matrix for F and V is given by The basic reproduction number is given by the spectral radius of J F (P 0 ) [J V (P 0 )] −1 as follows:
Thus for < 1,  0 is decreasing with , 2 , 3 , , and increasing with 1 , 2 , 1 , 2 , , , , Υ. Based on the result of derivatives, in general, cannot be commented on  0 sensitivity to other parameters 1 and , but for the parameter values in the proposed model, 1  0 < 0 and  0 < 0, thus  0 is decreasing with 1 and .

EXISTENCE AND UNIQUENESS ANALYSIS
In this section, we investigate the existence and uniqueness of the proposed model with the assistance of fixed-point theory. First, let us write Equation (11) as follows: By Lemma 3, the corresponding Volterra integral system of the second kind for Equation (16) is given by Theorem 1. The kernel of the proposed fractional model will be satisfying the Lipchitz condition if the following inequality hold: Proof. We will prove for the first kernel and similarly for the other. Consider function S(t) and S 1 (t), then By helping Equation (16), the proposed model Equation (11) can be written as follows: with initial condition (0) = 0 , where (t) = [

S(t) E(t) I(t) A(t) H(t) R(t) W(t)
Again, by Lemma 3, we have and assume that Ψ ∶ D → R 7 + be a continuous function; furthermore, define ∶= sup (t, )∈D ‖Ψ (t, (t))‖ then, for , which is a solution of Equation (18), where . It is clear that V is a non-empty set, since 0 ∈ V. We now define an operator on V. For each element ∈ V, Using that operator, Equation (19) can be written as = ; now we must prove that has a fixed-point. This will be done by using Schuder's second fixed-point theorem. First, we will show that ∈ V for ∈ V. For any ∈ V, we have Second, we show that is continuous. For every n , ∈ V, n = 1, 2, … , with lim n→∞ ‖ n − ‖ * = 0, we have lim n→∞ n (t) = (t), t ∈ [0, T], since all components of Ψ (t, (t)) are continuous on R, thus Ψ is continuous on On the other hand, Hence, ‖ n − ‖ * → 0 as n → ∞. Thus, is continuous. Finally, we have shown that (V) = { ∶ ∈ V} is relatively compact. By the Arzela-Ascoli theorem, it is enough to demonstrate that (V) is uniformly bounded and equicontinuous. Let ∈ (V), for all t ∈ [0, T], we have Thus, ‖( Since ≤ 1 and t 1 ≤ t 2 , then (t 1 − ) −1 ≥ (t 2 − ) −1 ; therefore, Hence, ‖( ) (t 1 ) − ( ) (t 2 )‖ * → 0 as t 2 → t 1 , and (V) is equicontinuous on [0, T]. Thus, (V) is relatively compact. Moreover, has a fixed-point, which is the required solution of Equation (18). This completes the proof.

NUMERICAL ALGORITHM
This section presents the numerical algorithm based on the predictor-corrector method. 37 Under the hypotheses of Theorem 2, there is a unique solution on [0, T]. Let t 0 = 0 < t 1 < … < t N = T be a uniformly divide of the interval [0, T] where t n = nh, n = 1, 2, … , N and h = T∕N. By Lemma 3, the solution of is equivalent to By Equation (30), we have Let S n+1 ≃ S (t n+1 ) , n = 0, 1, 2, … , N −1, by Lagrange interpolation, we approximate the kernel 1 (t, S(t)) over [t k , t k+1 ] as follows: Substituting Equation (32) into Equation (31), we get After the rearrangement of the summation on the right-hand side of Equation (33), we get where The quantity S n + 1 on the right-hand side of Equation (34) is predicted by S P n+1 applying the one-step Adams-Bashforth method to Equation (31) by replacing the function 1 ( , S( )) with the quantity 1 (t k , S k ) as follows: Similarly, we can obtain the numerical algorithm of the other equation of system Equation (16). Thus, the approximate solution is given by , , and where d k,n+1 = ( (n − k + 1) − (n − k) ) .  (36) and (37). Then, the fractional predictor-corrector method Equations (36),and (37) is conditional stable.

STABILITY ANALYSIS OF ITERATION METHOD
where, ∁ = max { 1, ( + 1) 2 1− } , and by using Lemma 2, we get where C is a positive constant. The initial values of infected people, death or recovery people, and hospitalized people as stated in the report of the WHO on February 4, 2020, are I 0 = 24545, R 0 = 907, and H 0 = 12627. 38 The initial values of A 0 , E 0 , and W 0 assumed as A 0 = 20000, E 0 = 80000, and W 0 = 50000. Since N = S 0 + E 0 + I 0 + A 0 + H 0 + R 0 , then S 0 = 7609967373.

Numerical simulation of the pandemic of COVID-19 in the world
Concerning the parameter values in Table 1, the basic reproduction number  0 > 1. That means that the pandemic will spread, and the equilibrium point of the model Equation (11) is positive. P * 0 = (2.02 × 10 7 , 4.957 × 10 9 , 7.43 × 10 5 , 2.66 × 10 6 , 2.7 × 10 6 , 1.37 × 10 10 , 3.4 × 10 2 ). All the approximate solutions are computed by using Wolfram Mathematica software with = 0.99. The graphical solutions of fractional system Equation (11) in the interval [0, 250] have been described in Figures 1-3, where the unit of time is days. Figures 1-3 show that the results of the model converge to their equilibrium point for different fractional-order derivatives and stable at that points. These figures indicate that the obtained plots have the same behavior pattern for different values of = 1, 0.9, 0.8, 0.7, 0.6.

Effect of parameters on the epidemic spread
The parameter values of the model play a significant role influences the spread of the epidemic. So, the single most effective way to restrict the spread of a disease is to create quarantine to decrease the mobility of individuals. If we choose the     Table 2, then the basic reproduction number  0 < 1 where Υ=0.00181. Thus, we get the results shown in Figures 4-6. Comparing the results with Fatima et al, 1 we note that the behavior of graphic solutions is the same. Figure 4-6 show that all variables in the pandemic model will be decreased and hit zero, indicating the system's stability.

CONCLUSION
In this paper, modeling and studied the pandemic trend of COVID-19 has been presented with fractional-order derivative.
Using the generation matrix method, the basic reproduction number is calculated that located whether the disease would persist or disappears from the population. The equilibrium points for this system are calculated, where the graphical solutions show that the results of the model converge to their equilibrium points. Based on the derivatives of the basic reproduction number, the sensitivity of the parameters was analyzed. The existence and uniqueness of the solution to the proposed model have been proven using the fixed-point theorem and by helping of the Arzela-Ascoli theorem. The fractional proposed model is solved using the predictor-corrector method in the sense of Caputo derivative. Using some essential lemmas, we proved that this method is conditionally stable. The results indicate that the disease will continue where that the basic reproduction number  0 > 1. We selected the same parameter values and same initial conditions in the classical model 1 for comparing the results of the fractional model with the classical model. We noted that both models have the same behavior. Moreover, parameters played a very significant role in limiting disease outbreaks.