Open Access

ARTICLE

# A Detailed Mathematical Analysis of the Vaccination Model for COVID-19

Abeer S. Alnahdi1,*, Mdi B. Jeelani1, Hanan A. Wahash2, Mansour A. Abdulwasaa3,4

1 Department of Mathematics and Statistics, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh, Saudi Arabia
3 Department of Statistics, Taiz University, Taiz, Yemen

* Corresponding Author: Abeer S. Alnahdi. Email:

Computer Modeling in Engineering & Sciences 2023, 135(2), 1315-1343. https://doi.org/10.32604/cmes.2022.023694

## Abstract

This study aims to structure and evaluate a new COVID-19 model which predicts vaccination effect in the Kingdom of Saudi Arabia (KSA) under Atangana-Baleanu-Caputo (ABC) fractional derivatives. On the statistical aspect, we analyze the collected statistical data of fully vaccinated people from June 01, 2021, to February 15, 2022. Then we apply the Eviews program to find the best model for predicting the vaccination against this pandemic, based on daily series data from February 16, 2022, to April 15, 2022. The results of data analysis show that the appropriate model is autoregressive integrated moving average ARIMA (1, 1, 2), and hence, a forecast about the evolution of theCOVID-19 vaccination in 60 days is presented. The theoretical aspect provides equilibrium points, reproduction number , and biologically feasible region of the proposed model. Also, we obtain the existence and uniqueness results by using the Picard-Lindel method and the iterative scheme with the Laplace transform. On the numerical aspect, we apply the generalized scheme of the Adams-Bashforth technique in order to simulate the fractional model. Moreover, numerical simulations are performed dependent on real data of COVID-19 in KSA to show the plots of the effects of the fractional-order operator with the anticipation that the suggested model approximation will be better than that of the established traditional model. Finally, the concerned numerical simulations are compared with the exact real available date given in the statistical aspect.

## Keywords

1  Introduction

Coronaviruses have three branches known as alpha, beta, and gamma. Recently, various strains of SARS-CoV-2 have emerged, including the most destructive and most dangerous delta variant, SARS-CoV-2 [1]. Human coronaviruses have been first distinguished during the 1960s. The first case reported of COVID-19 in Wuhan City in China’s Hubei Province on December 31, 2019, has been found to be contaminated with a new COVID-19 that has never been seen. As reported [2], this disease is believed to be transmitted from animal creatures to people, like SARS-CoV and MERS-CoV. Now the disease is sent from one individual to another. On March 23, 2022, according to the global statistics of this epidemic, the number of infected cases was estimated at 480,165,010, while the number of people recovered reached 414,597,310. By this date, there were 6,144,249 deaths related to the disease worldwide [3]. In this regard, an investigation of the patients who died found that most of them were elderly or patients who had been diagnosed with chronic diseases such as heart disease, diabetes, lung and kidney disease, etc.

Despite the different non-pharmaceutical control systems against COVID-19, a portion of the vaccines that have acquired crisis use approval (EUA) by the United States Centers for Disease Control and Prevention (CDC) are the Pfizer-BioNTech with 95% efficacy, the Moderna immunization with 94.5% viability, and the Janssen vaccine with 67% adequacy, and many others, see [4]. The vaccines have been successful against SARS-CoV-2 infections, including asymptomatic infections and symptomatic cases, severe COVID-19 disease, and deaths.

KSA began vaccinating people against the Coronavirus on December 17, 2020. On November 03, 2021, the use of the Pfizer vaccine was approved for the age groups 5–11 years. In the early stages of the vaccine, citizens and residents over 65 were targeted, which ended on February 18, 2021, whereas the second phase was launched on February 18, 2021, it targeted citizens and residents over 50 years old. The third phase targeted all citizens and residents wishing to receive the vaccine.

On the other hand, time series models are employed in the estimate process of variable behavior and other phenomena, as well as their future trends in diseases, epidemiology, climate sciences, economics, management, and other sciences. Lately, mathematical and statistical modeling have been used to forecast the behavior and patterns of some epidemics and diseases. The prediction method for these models includes four fundamental procedures, such as determining the model, estimation of unknown parameters, diagnostic process, and prediction process. The models (MA, AR, ARCH, ARMA, GARCH, ARIMA) are among the most commonly used time series models for forecasting [58].

In this context, the processes of making a mathematical model of a problem, interpreting the solution, validating the model, then making the model ready for utilization are not processes that can be conquered directly. Particularly these days, large numbers of the problems are complicated, nonlinear, have memory impact, or possess a stochastic construction; therefore, modeling methods and solutions specific to these problems must be developed.

Fractional calculus has shown gigantic improvement in applications to various real-world problems in different fields such as continuum mechanics, electromagnetic theory, and biological mathematics. Fractional calculus has become a substantial mathematical tool for the investigation of nonlinear derivative problems, see [913]. In modeling biological systems, fractional derivatives play a significant role in considering the nonlocality and memory impact properties that perfectly fit the test data of memory phenomena in various disciplines, such as mechanics, epidemiology, and psychology. The memory effect clarifies that the future condition of the fractional operator of a certain function relies upon its recorded behavior and present status. By using fractional derivatives, especially Caputo fractional derivatives, several real-world systems have been studied successfully in biomathematics and engineering [1417]. The fractional derivative’s definition has different methodologies with various kernels, so the researchers are keen on picking the best one. Specialists are drawn in by the nonsingular kernels as the classical singular kernels experienced difficulties in modeling some physical phenomena. Caputo et al. [18] overcome this trouble by expanding the Caputo fractional derivative to a nonsingular kernel. There are some interesting properties of the Caputo-Fabrizio operator that have been discussed by Losada et al. in [19]. ABC fractional derivative was given by Atangana et al. in [20], they considered the generalized Mittag-Leffler function as a nonsingular and nonlocal kernel. Some of the generalized characteristics of the ABC operator are introduced in [21,22]. This new operator commanded the attention of investigators because of its brilliant memory description [23,24]. Recent mathematical models on COVID-19 under the ABC operator show different approaches to manifest the transmission of the disease, see [2532]. In this regard, other papers dealt extensively with modeling and analysis of the transmission of the COVID-19 pandemic in order to limit the spread of this epidemic, see [3338]. Diagne et al. [39] formulated a mathematical model of the COVID-19 transmission mechanism with Vaccination and Treatment. A mathematical model of COVID-19 with the effect of vaccines was constructed by Gokbulut et al. [40]. They made the analysis of the model, then the idea of vaccination was used with a change in the vaccination rate among the population. Yavuz et al. [41] developed a mathematical model to reveal the effects of vaccine treatment on COVID-19 based on a system of ordinary differential equations.

The literature shows that very few works available on COVID-19 models are in the form of ABC fractional derivatives, even though they describe the dynamic behavior accurately compared to the classical derivatives. To our knowledge, there is no available literature on the COVID-19 fractional model with an expectation of the vaccination effect on KSA. Therefore, this is the main motivation behind our work. Also, the comparison of the real data with different fractional order simulations is one of our unique aims in this work. In addition, we extend and generalize the model studied by Yavuz et al. [41], based on a system of fractional differential equations involving ABC fractional derivatives.

This work is organized as follows. Section 2 provides a statistical analysis of COVID-19 vaccination in KSA. Some background material and auxiliary results of ABC fractional calculus are given in Section 3. The fractional mathematical model and its fundamental properties are described in Section 4. In Section 5, we prove the theoretical results of the given model. Sections 6 and 7 give numerical results and simulation results by estimating the parameters. In the last section, we provide brief conclusions of the study.

2  Statistical Analysis

In this section, we discuss some empirical results. First, we analyze and present the descriptive data of the vaccinated persons in KSA. Then we predict the number of people who will be fully vaccinated during the next period.

2.1 Descriptive Data Analysis

The latest statistics about the vaccination of people against the spread of the Coronavirus epidemic are collected from the official data published by the Saudi Ministry of Health as the fully vaccinated people reached 144,057 during the period January 06, 2021, to June 27, 2021, with an average of 8,401 per day. Then, data is collected on the number of fully vaccinated people daily during the period June 28, 2021, to February 15, 2022. Table 1 indicates that the number of vaccinated people in July 2021 reached 6,493,739, with an average of 209,476 people per day. Then it increased in August 2021 to reach 6,833,828 with an average of 220,446 people per day. Hence, we observe that the number of vaccinated people decreased continuously until the number of vaccinated people during the fifteen days of February 2022 reached 264,633. In general, the number of fully vaccinated people during the studied period reached 23,964,566, which represents 67.81% of the total population in KSA, which is 35,340,680 people. For more details, see Figs. 1 and 2.

Figure 1: Number of people fully vaccinated in the KSA during the period June 01, 2021 to February 15, 2022

Figure 2: Daily numbers of people completed vaccinations in KSA from June 28, 2021 to February 15, 2022

Table 2 summarizes the number of daily doses taken by people against coronavirus disease in KSA during the period June 28, 2021, to February 15, 2022. We find that the doses increased continuously from January 2021 to April 2021, then decreased in May 2021 and June 2021. After that, it increased dramatically in July 2021 and August 2021 to reach 19,252,318 doses with 33% of the total doses given during the study period. Then it decreased during the period from September 2021 to December 2021. The total doses given amounted to 57,964,264 during the studied period, with an average of 248,773 doses per day.

2.2 Prediction Strategies

In this portion, we are collecting data on the fully vaccinated people in KSA for the period June 01, 2021 to February 15, 2022. Moreover, the stability tests of these data are examined to use the prediction process as unit root tests and estimation of coefficients (ACF & APCF). Phillips-Peron’s and Dickey-Fuller’s tests show that the data series is unstable at the level, which means that there is a general trend in the series. As shown in Fig. 2. In this regard, we process the data by converting it to the logarithm and taking the first differences to remove the effect of the general trend from the series to be used in estimation and prediction, see Fig. 3. Furthermore, we repeat the unit root tests as in Table 3. Hence, we find that the calculated values are more significant than the critical values and Prob. = 0.0000, which is less than 5%, this indicates that the data series does not have a unit root and has become stable.

Figure 3: Transformation of the data of completely vaccinated people in KSA to the first difference

We apply ARIMA models for forecasting through the Eviews program on the data series about the number of vaccinated people completely in KSA. As well, we estimate ARIMA models by using the Ordinary Least Squares method, as shown in Table 4. On the other hand, the models that have no statistical significance are excluded and the comparison only in the statistically significant models, and we are choosing a proper model that achieves the best models in which the coefficient of determination is greater, less variance, less volatility, and less value to AIC indicator.

Moreover, we check the ARIMA (1, 1, 2) model by testing residuals and the shape of the autocorrelation and partial autocorrelation coefficients. It follows from Fig. 4 that the actual values match the estimated values, and this indicates that the differences were small, as well as the efficiency of the model and its suitability in the forecasting process.

Figure 4: Evaluation of the ARIMA model for actual and fitted values and residual limits

We see through the chosen model a decrease in the number of people who will be fully dosed to (927164) during the predictive period from February 16, 2022, to April 14, 2022. For more details, see Fig. 5 and Table 6.

Figure 5: Prediction of the number of fully vaccinated people in KSA with confidence intervals (95%) for the period February 16, 2022 to April 15, 2022

Finally, we estimate a linear model using the least-squares method to test the predictive ability of the model. Indeed, we take the actual values as a dependent variable and the estimated values as an independent variable. We conclude that the closer the estimated parameter to one, the more closely the estimated values are to the actual values. Through the results of Table 5, it is clear that the predicted parameter is close to one (0.955688), and this indicates the quality of the model in prediction and that there is a convergence of the predicted values from the actual values and is statistically significant. Additionally, we find that the value of F = 1287.919, whereas Prob. = 0.0000 is less than 0.05, which indicates the model is significant and good.

3  Auxiliary Results

Definition 3.1. [20] Let ρ[0,1], and ωH1(a,b), b > a. Then, the ABC fractional derivative of order ρ is defined by

ABCDa+ρω(t)=(ρ)1ρatEρ(ρρ1(tθ)ρ)ω(θ)dθ,t>0,(1)

The normalization function (ρ) satisfies (0)=(1)=1, where Eρ is the Mittag-Leffler function defined by

Eρ(z)=k=0zkΓ(kρ+1), Re(ρ)>0,zC.

Definition 3.2. [20] The correspondent AB fractional integral is given by

ABa+ρω(t)=1ρ(ρ)ω(t)+ρ(ρ)Γ(ρ)at(tθ)ρ1ω(θ)dθ,ρ[0,1].

Lemma 3.1. [22] Let ω(t) be a function defined on [a,b] and 0<ρ1. Then

ABABCa+ρABCDa+ρω(t)=ω(t)ω(a).

Lemma 3.2. [20] The Laplace transform of the ABC fractional derivative is defined by

LABCD0+ρω(t)=(ρ)1ρsρL{ω(t)}(s)sρ1ω(0)sρ+ρ1ρ.

Lemma 3.3. [20] The fractional differential system ABCD0+ρω(t)=f(t), ω(0)=ω0, gives the unique solution

ω(t)=ω0+1ρ(ρ)f(t)+ρ(ρ)Γ(ρ)at(tθ)ρ1f(θ)dθ,ρ[0,1].

4  Model Derivation in ABC Operator

Yavuz et al. [41] developed a mathematical model to reveal the effects of vaccine treatment on COVID-19 described by a system of ODEs as follows:

{S(ϰ)=Λ(αE(ϰ)+m+μ)S(ϰ),E(ϰ)=αE(ϰ)S(ϰ)+pV(ϰ)E(ϰ)(κI(ϰ)+c+μ)E(ϰ),I(ϰ)=κI(ϰ)E(ϰ)(z+μ+σ)I(ϰ),V(ϰ)=mS(ϰ)(pE(ϰ)+μ)V(ϰ),R(ϰ)=zI(ϰ)+cE(ϰ)μR(ϰ).(2)

With initial conditions

S(0)=S0,E(0)=E0,I(0)=I0,V(0)=V0,R(0)=R0.

In the present work, we make extend the model (2) by replacing the time derivative with the ABC fractional derivative. With this alteration, the right-and left-hand sides will not have the same dimensions. To conquer this affair, we add an auxiliary parameter K with the dimension of s, to change the nonlocal fractional operator so that the sides possess the same dimension [43]. Herewith, we reformulate the following fractional system:

{1K1ΘABCDΘS(ϰ)=Λ(αE(ϰ)+m+μ)S(ϰ),1K1ΘABCDΘE(ϰ)=αE(ϰ)S(ϰ)+pV(ϰ)E(ϰ)(κI(ϰ)+c+μ)E(ϰ),1K1ΘABCDΘI(ϰ)=κI(ϰ)E(ϰ)(z+μ+σ)I(ϰ),1K1ΘABCDΘV(ϰ)=mS(ϰ)(pE(ϰ)+μ)V(ϰ),1K1ΘABCDΘR(ϰ)=zI(ϰ)+cE(ϰ)μR(ϰ).(3)

With initial conditions

S(0)=S00,E(0)=E00,I(0)=I00,V(0)=V00,R(0)=R00,(4)

where ABCDΘ is the ABC fractional derivative of order Θ(0,1], the functions S(ϰ), E(ϰ), I(ϰ), V(ϰ), R(ϰ) and their fractional derivatives are continuous at ϰ ≥ 0, and (S,E,I,V,R)Λ5, where Λ:=C([0,T],R) is a Banach space endowed with the norm X:=supt[0,T]{|X(t)|}. Furthermore, we set Ω:=(Λ5,X) with the norm ||X||=||S||+||E||+||I||+||V||+||R||.

Here S(ϰ) is the class of Susceptible Individuals, E(ϰ) is the class of Exposed Individuals, I(ϰ) is the class of Infected individuals, V(ϰ) is the class of Individuals Vaccinated, R(ϰ) is the class of Recovered Individuals. The model parameter values and source are given in Table 7.

Total population at time ϰ, denoted by N(ϰ) and given by N(ϰ) = S(ϰ) + E(ϰ) + I(ϰ) + V(ϰ) + R(ϰ).

4.1 Positivity and Boundedness

In this portion, we show that the non-negative domain R5 is positively invariant region and the solutions of model are bounded, where R+5={ΦR5:Φ(ϰ)0}.

Theorem 4.1. There is a unique solution Φ(ϰ)=(S(ϰ),E(ϰ),I(ϰ),V(ϰ),R(ϰ))T for the models (3) and (4) on ϰ ≥ 0 in [0, T] (0<T<) and the solution will remain in R5. Furthermore, the solution is stated in the region ΔR+5, described by

Δ={(S(ϰ),E(ϰ),I(ϰ),V(ϰ),R(ϰ))R5:N(ϰ)Λμ}.

Proof. From the model (3), we have

ABCDΘS(ϰ)=K1Θ(Λ(αE(ϰ)+m+μ)S(ϰ)),ABCDΘE(ϰ)=K1Θ(αE(ϰ)S(ϰ)+pV(ϰ)E(ϰ)(κI(ϰ)+c+μ)E(ϰ)),ABCDΘI(ϰ)=K1Θ(κI(ϰ)E(ϰ)(z+μ+σ)I(ϰ)).ABCDΘV(ϰ)=K1Θ(mS(ϰ)(pE(ϰ)+μ)V(ϰ)).ABCDΘR(ϰ)=K1Θ(zI(ϰ)+cE(ϰ)μR(ϰ)).

The norm and all assumptions of the classical results are valid. It follows that:

{ABCDΘS(ϰ)K1Θ(αE(ϰ)+m+μ)S(ϰ),ABCDΘE(ϰ)K1Θ(κI(ϰ)+c+μ)E(ϰ),ABCDΘI(ϰ)K1Θ(z+μ+σ)I(ϰ)ABCDΘV(ϰ)K1Θ(pE(ϰ)+μ)V(ϰ)ABCDΘR(ϰ)K1ΘμR(ϰ)

For K>0 and ϰ > 0, we have

{S(ϰ)S0EΘ(ΘK1Θ(α||E||+m+μ)ϰΘ(Θ)(1Θ)K1Θ(α||E||+m+μ)),E(ϰ)E0EΘ(ΘK1Θ(κ||I||+c+μ)ϰΘ(Θ)(1Θ)K1Θ(κ||I||+c+μ)),I(ϰ)I0EΘ(ΘK1Θ(z+μ+σ)ϰΘ(Θ)(1Θ)K1Θ(z+μ+σ)),V(ϰ)V0EΘ(ΘK1Θ(p||E||+μ)ϰΘ(Θ)(1Θ)K1Θ(p||E||+μ))R(ϰ)R0EΘ(ΘK1ΘμϰΘ(Θ)(1Θ)K1Θμ).

This shows that if S0,E0,I0,V0, and R0 are positive, then S(ϰ), E(ϰ), I(ϰ), V(ϰ), and R(ϰ) are also positive. Now, we give a biologically feasible region for the suggested models (3) and (4). Consider summing all the equations in the model as follows

1K1ΘABCDΘN(ϰ)=1K1ΘABCDΘ(S(ϰ)+E(ϰ)+I(ϰ)+V(ϰ)+R(ϰ)).

It follows for the whole population that

ABCDΘN(ϰ)=K1Θ(ΛμN(ϰ)σI(ϰ))K1Θ(ΛμN(ϰ)),

which implies that

ABCDΘN(ϰ)K1Θ(ΛμN(ϰ)).

After applying the Laplace transform, we have

N(s)(K1Θ(1Θ)Λϑ+(Θ)ϑN(0))sΘ1(sΘ+λ)+K1ΘΘΛϑs1(sΘ+λ),

where ϑ:=[(Θ)+K1Θ(1Θ)μ], and λ=K1ΘΘμ[(Θ)+K1Θ(1Θ)μ]. Using the inverse Laplace, we get

N(ϰ)=(K1Θ(1Θ)Λϑ+(Θ)N(0)ϑ)EΘ,1(λϰΘ)+K1ΘΘΛϑϰΘEΘ,Θ+1(λϰΘ),

where Eα,β(z)=k=0zkΓ(kα+β), and L[ϰβ1Eα,β(±aϰα)]=sαβ(sαa). Subsequently, the solution of the model with the nonnegative conditions in Δ stays in Δ. In this way, the region Δ is positively invariant and attracts all the solutions in R5. This shows that the models (3) and (4) is epidemiologically well posed and the solutions.

4.2 Equilibrium Points and Reproduction Number

In this part, we will find the equilibrium points of the COVID-19 model. By equating each equation of model (3) to the zero, we can write

0=K1Θ(Λ(αE+m+μ)S),0=K1Θ(αES+pVE(κI+c+μ)E),0=K1Θ(κIE(z+μ+σ)I).0=K1Θ(mS(pE+μ)V).0=K1Θ(zI+cEμR).

For K>0, we have

0=Λ(αE+m+μ)S,0=αES+pVE(κI+c+μ)E,0=κIE(z+μ+σ)I.0=mS(pE+μ)V.0=zI+cEμR.

Hence, the disease free equilibrium (DFE) and endemic equilibrium (EE) are given by the following theorem.

Theorem 4.2. We have the next affirmations:

i. The model (3) has always DFE, E0=(S0,E0,I0,V0,R0), where S0=Λm+μ,E0=0,I0=0,V0=Λm(m+μ)μ, and R0=0.

ii. The model (3) has EE, E=(S,E,I,V,R), where

S=κΛαμ+ασ+κμ+κm+αzE=σ+μ+zκI=κ(Λ(αμ+mp)+αp(σ+μ+z))(κμ+pσ+pμ+pz)(κm+κμ+ασ+αμ+αz)c+μκV=κ2Λm(κμ+pσ+pμ+pz)(κm+κμ+ασ+αμ+αz)R=κ2(m(μ(c(μ+σ)μz)+Λpz)+μ(cμ(μ+σ)+z(αΛμ2)))κμ(κμ+pσ+pμ+pz)(κm+κμ+ασ+αμ+αz)+κ(μ+σ+z)+p(c(μ+σ)(μ+m)+z(αΛμ(μ+m)))κμ(κμ+pσ+pμ+pz)(κm+κμ+ασ+αμ+αz)+αμ(c(μ+σ)μz)+αp(σ+μ+z)2(c(μ+σ)μz)κμ(κμ+pσ+pμ+pz)(κm+κμ+ασ+αμ+αz).

By using the model (3) and next generation matrix method, the basic reproduction number R0 calculated as

R0=Λ(μα+mp)μ(μ+c)(μ+m).

For details, see [41].

4.3 Stability Analysis

Stability analysis of the equilibrium points is clarified in [41], exhaustively. In order of stability, the authors dealt with the next two theorems.

Theorem 4.3. The DFE point E0 of the epidemic model is locally asymptotically stable if R0<1, otherwise unstable.

Theorem 4.4. The EE point E of the epidemic model is locally asymptotically stable if R0>1, otherwise unstable.

Here we comment that R0 gives us data about the transmission of infectious sickness. For example in the normal utilization of contamination models, when R0>1, the disease infection will actually be able to begin transmitting in society, yet if R0<1 it implies that the disease infection continues declining in the general public. Overall, the bigger the worth of R0 the troublesome it is to control the pandemic from spreading in the public arena.

5  Existence Criteria

In this section, we apply the Picard-Lindel method and the Laplace transform to investigate the existence and uniqueness of solution for preventive and curative to fractional COVID-19 disease model.

5.1 Iterative Scheme with Laplace Transform

Theorem 5.1. For 0<Θ1, the following ABC-type FDE

1K1αABCDΘu(ϰ)=f(ϰ)(5)

has a unique solution, which is

u(ϰ)=u(0)+(1Θ)κ1Θ(Θ)f(ϰ)+Θκ1Θ(Θ)Γ(Θ)0ϰ(ϰν)Θ1f(ν)dν.(6)

Proof. Applying the Laplace transform on both sides of Eq. (5) we obtain

L{1K1ΘABCDΘu(ϰ)}(ϖ)=L{f(ϰ)}(ϖ),ϖ>0.

It follows from Theorem 3 in [20] that

(Θ)1ΘϖΘL(u(ϰ))(ϖ)ϖΘ1u(0)ϖΘ+Θ1Θ=K1ΘL{f(ϰ)}(ϖ),

which is equivalent to

L(u(ϰ))(ϖ)=1ϖu(0)+(1Θ)K1Θ(Θ)L{f(ϰ)}(ϖ)+ΘK1ΘϖΘ(Θ)L{f(ϰ)}(ϖ).

Applying the inverse Laplace transform give us

u(ϰ)=u(0)+(1Θ)K1Θ(Θ)f(ϰ)+L1{ΘK1ΘϖΘ(Θ)L{f(ϰ)}(ϖ)}(ϰ).

Now, we have

ΘK1Θ(Θ)(1ϖΘ)=ΘK1Θ(Θ)L{ϰΘ1Γ(Θ)}(ϖ).

Let F(ϖ)=ΘK1Θ(Θ)L{ϰΘ1Γ(Θ)}(ϖ), and G(ϖ)=L{f(ϰ)}(ϖ). It follows form convolution theorem that

L1{ΘK1ΘϖΘ(Θ)L{f(ϰ)}(ϖ)}(ϰ)=L1{F(ϖ)×G(ϖ)}(ϰ)=ΘK1Θ(Θ){ϰΘ1Γ(Θ)f(ϰ)}(ϖ)=ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰϖ)Θ1f(ϖ)dϖ.

Hence

u(ϰ)=u(0)+(1Θ)K1Θ(Θ)f(ϰ)+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰν)Θ1f(ν)dν.

By using Theorem 5.1, our model is equivalent to

{S(ϰ)S(0)=(1Θ)K1Θ(Θ){Λ(αE(ϰ)+m+μ)S(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{Λ(αE(r)+m+μ)S(r)}drE(ϰ)E(0)=(1Θ)K1Θ(Θ){αE(ϰ)S(ϰ)+pV(ϰ)E(ϰ)(κI(ϰ)+c+μ)E(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{αE(r)S(r)+pV(r)E(r)(κI(r)+c+μ)E(r)}drI(ϰ)I(0)=(1Θ)K1Θ(Θ){κI(ϰ)E(ϰ)(z+μ+σ)I(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{κI(r)E(r)(z+μ+σ)I(r)}drV(ϰ)V(0)=(1Θ)K1Θ(Θ){mS(ϰ)(pE(ϰ)+μ)V(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{mS(r)(pE(r)+μ)V(r)}drR(ϰ)R(0)=(1Θ)K1Θ(Θ){zI(ϰ)+cE(ϰ)μR(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{zI(r)+cE(r)μR(r)}dr(7)

The iterative scheme of the model (7) is given by

S0(ϰ)=S(0),E0(ϰ)=E(0),I0(ϰ)=I(0),V0(ϰ)=V(0),R0(ϰ)=R(0).

{Sn+1(ϰ)=(1Θ)K1Θ(Θ){Λ(αEn(ϰ)+m+μ)Sn(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{Λ(αEn(r)+m+μ)Sn(r)}drEn+1(ϰ)=(1Θ)K1Θ(Θ){αEn(ϰ)Sn(ϰ)+pVn(ϰ)En(ϰ)(κIn(ϰ)+c+μ)En(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{αEn(r)Sn(r)+pVn(r)En(r)(κIn(r)+c+μ)En(r)}drIn+1(ϰ)=(1Θ)K1Θ(Θ){κIn(ϰ)En(ϰ)(z+μ+σ)In(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{κIn(r)En(r)(z+μ+σ)In(r)}drVn+1(ϰ)=(1Θ)K1Θ(Θ){mSn(ϰ)(pEn(ϰ)+μ)Vn(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{mSn(r)(pEn(r)+μ)Vn(r)}drRn+1(ϰ)=(1Θ)K1Θ(Θ){zIn(ϰ)+cEn(ϰ)μRn(ϰ)}+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1{zIn(r)+cEn(r)μRn(r)}dr(8)

Taking the limit as n, we anticipate getting the exact solution.

5.2 Existence of a Unique Solution

In this portion, we apply the Picard-Lindel method to investigate the existence of solution for the fractional COVID-19 disease model, which its mathematical represente are presented by:

{1K1ΘABCDΘS(ϰ)=Λ(αE(ϰ)+m+μ)S(ϰ),1K1ΘABCDΘE(ϰ)=αE(ϰ)S(ϰ)+pV(ϰ)E(ϰ)(κI(ϰ)+c+μ)E(ϰ),1K1ΘABCDΘI(ϰ)=κI(ϰ)E(ϰ)(z+μ+σ)I(ϰ).1K1ΘABCDΘV(ϰ)=mS(ϰ)(pE(ϰ)+μ)V(ϰ).1K1ΘABCDΘR(ϰ)=zI(ϰ)+cE(ϰ)μR(ϰ).(9)

For sake of simpilicity, we define the functions Ψi, i = 1, 2,…, 5, as follows:

{Ψ1(ϰ,S(ϰ))=Λ(αE(ϰ)+m+μ)S(ϰ),Ψ2(ϰ,E(ϰ))=αE(ϰ)S(ϰ)+pV(ϰ)E(ϰ)(κI(ϰ)+c+μ)E(ϰ),Ψ3(ϰ,I(ϰ))=κI(ϰ)E(ϰ)(z+μ+σ)I(ϰ),Ψ4(ϰ,V(ϰ))=mS(ϰ)(pE(ϰ)+μ)V(ϰ),Ψ5(ϰ,R(ϰ))=zI(ϰ)+cE(ϰ)μR(ϰ).(10)

The model described in Eq. (9) becomes in the following form:

{1K1ΘABCDΘS(ϰ)=Ψ1(ϰ,S(ϰ)),1K1ΘABCDΘE(ϰ)=Ψ2(ϰ,E(ϰ)),1K1ΘABCDΘI(ϰ)=Ψ3(ϰ,I(ϰ)),1K1ΘABCDΘV(ϰ)=Ψ4(ϰ,V(ϰ)),1K1ΘABCDΘR(ϰ)=Ψ5(ϰ,R(ϰ)),(11)

By applying the operator ABΘ on Eq. (11), we obtain

{S(ϰ)S(0)=(1Θ)K1Θ(Θ)Ψ1(ϰ,S(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ1(r,S(r))dr,E(ϰ)E(0)=(1Θ)K1Θ(Θ)Ψ2(ϰ,E(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ2(r,E(r))dr,I(ϰ)I(0)=(1Θ)K1Θ(Θ)Ψ3(ϰ,I(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ3(r,I(r))dr,V(ϰ)V(0)=(1Θ)K1Θ(Θ)Ψ4(ϰ,V(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ4(r,V(r))dr,R(ϰ)R(0)=(1Θ)K1Θ(Θ)Ψ5(ϰ,R(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ5(r,R(r))dr.(12)

The kernels in Eq. (10) satisfies the Lipschitz condition for 0Li<1, i = 1, 2, …5, if and only if the nonlinear functions S(ϰ), E(ϰ), I(ϰ), V(ϰ) and R(ϰ) have an upper bound. Let S(ϰ) and S(ϰ) be two functions. Then

||Ψ1(ϰ,S(ϰ))Ψ1(ϰ,S(ϰ))||=||Λ(αE(ϰ)+m+μ)S(ϰ)Λ+(αE(ϰ)+m+μ)S(ϰ)||=||(αE(ϰ)+m+μ)(S(ϰ)S(ϰ))||(αsupϰ[0,T]E(ϰ)+m+μ)||S(ϰ)S(ϰ)||=L1||S(ϰ)S(ϰ)||,(13)

where L1:=α||E||+m+μ. Thus,

||Ψ1(ϰ,S(ϰ))Ψ1(ϰ,S(ϰ))||L1||S(ϰ)S(ϰ)||.(14)

Repeating the same procedure as in Eq. (13), we get

||Ψ2(ϰ,E(ϰ))Ψ2(ϰ,E(ϰ))||L2||E(ϰ)E(ϰ)||,||Ψ3(ϰ,I(ϰ))Ψ3(ϰ,I(ϰ))||L3||I(ϰ)I(ϰ)||,||Ψ4(ϰ,V(ϰ))Ψ4(ϰ,V(ϰ))||L4||V(ϰ)V(ϰ)||,||Ψ5(ϰ,R(ϰ))Ψ5(ϰ,R(ϰ))||L5||R(ϰ)R(ϰ)||,(15)

where L2:=α||S||+p||V||+κ||I||+c+μ, L3:=κ||E||+z+μ+σ, L4:=p||E||+μ, and L5:=μ.

Eq. (12) can be written in the recursive form given by

{Sn(ϰ)=S(0)+(1Θ)K1Θ(Θ)Ψ1(ϰ,Sn1(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ1(r,Sn1(r))dr,En(ϰ)+E(0)+(1Θ)K1Θ(Θ)Ψ2(ϰ,En1(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ2(r,En1(r))dr,In(ϰ)=I(0)+(1Θ)K1Θ(Θ)Ψ3(ϰ,In1(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ3(r,In1(r))dr,Vn(ϰ)=V(0)+(1Θ)K1Θ(Θ)Ψ4(ϰ,Vn1(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ4(r,Vn1(r))dr,Rn(ϰ)=R(0)+(1Θ)K1Θ(Θ)Ψ5(ϰ,Rn1(ϰ))+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ5(r,Rn1(r))dr.(16)

Now, we denote the difference between successive components by Φni, i = 1, 2, …, 5. Therefore,

Φn1(ϰ)=Sn(ϰ)Sn1(ϰ)=(1Θ)K1Θ(Θ)[Ψ1(ϰ,Sn1(ϰ))Ψ1(ϰ,Sn2(ϰ))]+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1[Ψ1(r,Sn1(r))Ψ1(r,Sn2(r))]dr,Φn2(ϰ)=En(ϰ)En1(ϰ)=(1Θ)K1Θ(Θ)[Ψ2(ϰ,En1(ϰ))Ψ2(ϰ,En2(ϰ))]+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1[Ψ2(r,En1(r))Ψ2(r,En2(r))]dr,Φn3(ϰ)=In(ϰ)In1(ϰ)=(1Θ)K1Θ(Θ)[Ψ3(ϰ,In1(ϰ))Ψ3(ϰ,In2(ϰ))]+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1[Ψ3(r,In1(r))Ψ3(r,In2(r))]dr,Φn4(ϰ)=Vn(ϰ)Vn1(ϰ)=(1Θ)K1Θ(Θ)[Ψ4(ϰ,Vn1(ϰ))Ψ4(ϰ,Vn2(ϰ))]+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1[Ψ4(r,Vn1(r))Ψ4(r,Vn2(r))]dr,Φn5(ϰ)=Rn(ϰ)Rn1(ϰ)=(1Θ)K1Θ(Θ)[Ψ5(ϰ,Rn1(ϰ))Ψ5(ϰ,Rn2(ϰ))]+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1[Ψ5(r,Rn1(r))Ψ5(r,Rn2(r))]dr,(17)

Taking into account that

{Sn(ϰ)=i=0nΦi1(ϰ),En(ϰ)=i=0nΦi2(ϰ),In(ϰ)=i=0nΦi3(ϰ),Vn(ϰ)=i=0nΦi4(ϰ),Rn(ϰ)=i=0nΦi5(ϰ).(18)

Then, we take the norm on both sides of Eq. (17). It follows from Eqs. (14) and (15) that

||Φn1(ϰ)||=(1Θ)K1ΘL1(Θ)||Φn11(ϰ)||+ΘK1ΘL1(Θ)Γ(Θ)0ϰ(ϰr)Θ1||Φn11(r)||dr,||Φn2(ϰ)||=(1Θ)K1ΘL2(Θ)||Φn12(ϰ)||+ΘK1ΘL2(Θ)Γ(Θ)0ϰ(ϰr)Θ1||Φn12(r)||dr,||Φn3(ϰ)||=(1Θ)K1ΘL3(Θ)||Φn13(ϰ)||+ΘK1ΘL3(Θ)Γ(Θ)0ϰ(ϰr)Θ1||Φn13(r)||dr,||Φn4(ϰ)||=(1Θ)K1ΘL4(Θ)||Φn14(ϰ)||+ΘK1ΘL4(Θ)Γ(Θ)0ϰ(ϰr)Θ1||Φn14(r)||dr,||Φn5(ϰ)||=(1Θ)K1ΘL5(Θ)||Φn15(ϰ)||+ΘK1ΘL5(Θ)Γ(Θ)0ϰ(ϰr)Θ1||Φn15(r)||dr.(19)

Next, we shall prove the main theorem based on the above results.

Theorem 5.2. The fractional models (3)--(4) has a unique solution for ϰ ∈ [0, T] if

((1Θ)K1Θ(Θ)+K1ΘTΘ(Θ)Γ(Θ))L1<1.(20)

Proof. Since S(t), E(t), I(t), V(t), and R(t) are bounded functions and satisfying the Lipschitz condition. Therefore, by virtue of Eq. (19), we obtain

||Φn1(ϰ)||||Sn(0)||[(1Θ)K1ΘL1(Θ)+TΘK1ΘL1(Θ)Γ(Θ)]n,||Φn2(ϰ)||||En(0)||[(1Θ)K1ΘL2(Θ)+TΘK1ΘL2(Θ)Γ(Θ)]n,||Φn3(ϰ)||||In(0)||[(1Θ)K1ΘL3(Θ)+TΘK1ΘL3(Θ)Γ(Θ)]n,||Φn4(ϰ)||||Vn(0)||[(1Θ)K1ΘL4(Θ)+TΘK1ΘL4(Θ)Γ(Θ)]n,||Φn5(ϰ)||||Rn(0)||[(1Θ)K1ΘL5(Θ)+TΘK1ΘL5(Θ)Γ(Θ)]n.(21)

Consequently, the sequences in Eq. (21) are exist and smooth, i.e., ||Φni(ϰ)||0, (i=1,2,,5) as n.

Now, we show that the functions in Eq. (21) are the solutions of the proposed model. We suppose

S(ϰ)S(0)=Sn(ϰ)ϝn1(ϰ),E(ϰ)E(0)=En(ϰ)ϝn2(ϰ),I(ϰ)I(0)=In(ϰ)ϝn3(ϰ),V(ϰ)V(0)=Vn(ϰ)ϝn4(ϰ),R(ϰ)R(0)=Rn(ϰ)ϝn5(ϰ).(22)

Then we show that the terms in Eq. (22) verify that ||ϝni(ϰ)||0, as n, (i = 1, 2, …, 5). Thus, we have

||ϝn1(ϰ)||=||(1Θ)K1Θ(Θ)[Ψ1(ϰ,S(ϰ))Ψ1(ϰ,Sn1(ϰ)]+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1[Ψ1(r,S(r))Ψ1(r,Sn1(r))]dr||(1Θ)K1Θ(Θ)||Ψ1(ϰ,S(ϰ))Ψ1(ϰ,Sn1(ϰ)||+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1||Ψ1(r,S(r))Ψ1(r,Sn1(r))||dr((1Θ)K1Θ(Θ)+K1ΘTΘ(Θ)Γ(Θ))L1||S(ϰ)Sn1(ϰ)||.

Repeating the procedure above, we get

||ϝn1(ϰ)||((1Θ)K1Θ(Θ)+K1ΘTΘ(Θ)Γ(Θ))n+1L1n+1M.

As n, ||ϝn1(ϰ)||0.By the same arguments above, it can be shown that ||ϝni(ϰ)||0, as n, (i = 2, 3, …, 5).

To prove the uniqueness result, we assume that the model (3) has another solution (S(ϰ),E(ϰ),I(ϰ),V(ϰ),R(ϰ)). Then

||S(ϰ)S(ϰ)||(1Θ)K1Θ(Θ)||Ψ1(ϰ,S(ϰ))Ψ1(ϰ,S(ϰ)||+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1||Ψ1(r,S(r))Ψ1(r,S(r))||dr((1Θ)K1Θ(Θ)+K1ΘTΘ(Θ)Γ(Θ))L1||S(ϰ)S(ϰ)||,

which implies

[1((1Θ)K1Θ(Θ)+K1ΘTΘ(Θ)Γ(Θ))L1]||S(ϰ)S(ϰ)||0.

It is clear that S(ϰ)=S(ϰ), if the following inequality is satisfied

[1((1Θ)K1Θ(Θ)+K1ΘTΘ(Θ)Γ(Θ))L1]>0.

This is satisfied by the hypothesis (20). Hence, S(ϰ)=S(ϰ).

Similarly, we obtain E(ϰ)=E(ϰ), I(ϰ)=I(ϰ), V(ϰ)=V(ϰ), and R(ϰ)=R(ϰ).

6  Numerical Results

This section gives the numerical solution of the COVID-19 model (3) under ABC fractional derivatives. The suggested fractional model is addressed numerically using a generalized Adams–Bashforth–Moulton technique [44,45].

By applying the operator ABΘ on the first equation of (11), we obtain

S(ϰ)S(0)=(1Θ)K1Θ(Θ)Ψ1(ϰ,S)+ΘK1Θ(Θ)Γ(Θ)0ϰ(ϰr)Θ1Ψ1(r,S)dr.

For ϰ=ϰn+1, n = 0, 1, 2, …, we get

S(ϰn+1)=S0+(1Θ)K1Θ(Θ)Ψ1(ϰn,S)+ΘK1Θ(Θ)Γ(Θ)0ϰn+1(ϰn+1r)Θ1Ψ1(r,S)dr=S0+(1Θ)K1Θ(Θ)Ψ1(ϰn,S)+ΘK1Θ(Θ)Γ(Θ)s=0nϰsϰs+1(ϰn+1r)Θ1Ψ1(r,S)dr.

Over [ϰs,ϰs+1], the function Ψ1(r,S) can be approximated by the interpolation polynomial

P(s,r)=rϰs1ϰsϰs1G(ϰs,Z(ϰs))rϰsϰsϰs1G(ϰs1,Z(ϰs1))=G(ϰs,Z(ϰs))h(rϰs1)G(ϰs1,Z(ϰs1))h(rϰs)G(ϰs,Zs)h(rϰs1)G(ϰs1,Zs1)h(rϰs),

which implies

S(ϰn+1)=S0+(1Θ)K1Θ(Θ)Ψ1(ϰn,S)+ΘK1Θ(Θ)Γ(Θ)s=0n{Ψ1(ϰs,S)hϰsϰs+1(ϰn+1r)Θ1(rϰs1)drΨ1(ϰs1,S)hϰsϰs+1(ϰn+1r)Θ1(rϰs)dr.

Set IΘ,1:=ϰsϰs+1(ϰn+1r)Θ1(rϰs1)dr and IΘ,2:=ϰsϰs+1(ϰn+1r)Θ1(rϰs)dr. Evaluating these integrals, we obtain

IΘ,1=hΘ+1(n+1s)Θ(ns+Θ+2)(ns)Θ(ns+2Θ+2)Θ(Θ+1),

IΘ,2=hΘ+1(n+1s)Θ(ns)Θ(ns+Θ+1)Θ(Θ+1).

Finally,

S(ϰn+1)=S0+(1Θ)K1Θ(Θ)Ψ1(ϰn,S)+ΘK1Θ(Θ)Γ(Θ)s=0n{hΘΨ1(ϰs,S)Γ(Θ+2)[(n+1s)Θ(ns+Θ+2)(ns)Θ(ns+2Θ+2)]hΘΨ1(ϰs1,S)Γ(Θ+2)[(n+1s)Θ(ns)Θ(ns+Θ+1)].(23)

In an analogous manner for the remainder of the equations of the model (3), we obtain the recursive formulae as below:

E(ϰn+1)=E0+(1Θ)K1Θ(Θ)Ψ1(ϰn,E)+ΘK1Θ(Θ)Γ(Θ)s=0n{hΘΨ1(ϰs,E)Γ(Θ+2)[(n+1s)Θ(ns+Θ+2)(ns)Θ(ns+2Θ+2)]hΘΨ1(ϰs1,E)Γ(Θ+2)[(n+1s)Θ(ns)Θ(ns+Θ+1)].(24)

I(ϰn+1)=I0+(1Θ)K1Θ(Θ)Ψ1(ϰn,I)+ΘK1Θ(Θ)Γ(Θ)s=0n{hΘΨ1(ϰs,I)Γ(Θ+2)[(n+1s)Θ(ns+Θ+2)(ns)Θ(ns+2Θ+2)]hΘΨ1(ϰs1,I)Γ(Θ+2)[(n+1s)Θ(ns)Θ(ns+Θ+1)].(25)

V(ϰn+1)=V0+(1Θ)K1Θ(Θ)Ψ1(ϰn,V)+ΘK1Θ(Θ)Γ(Θ)s=0n{hΘΨ1(ϰs,V)Γ(Θ+2)[(n+1s)Θ(ns+Θ+2)(ns)Θ(ns+2Θ+2)]hΘΨ1(ϰs1,V)Γ(Θ+2)[(n+1s)Θ(ns)Θ(ns+Θ+1)].(26)

R(ϰn+1)=R0+(1Θ)K1Θ(Θ)Ψ1(ϰn,R)+ΘK1Θ(Θ)Γ(Θ)s=0n{hΘΨ1(ϰs,R)Γ(Θ+2)[(n+1s)Θ(ns+Θ+2)(ns)Θ(ns+2Θ+2)]hΘΨ1(ϰs1,R)Γ(Θ+2)[(n+1s)Θ(ns)Θ(ns+Θ+1)].(27)

7  Simulations Results

The biological parameters estimated from the actual data reported in KSA for the period June 01, 2021 to February 15, 2022, and classified in Table 7 is used to acquire the simulation results. Further for initial data we use S(0) = 9.317558 millions, E(0) = 0, I(0) = 0.562300 million, V(0) = 24.000000 millions, R(0) = 0.580000 million, where N = 34.810000 millions is a total population. For the purposes of numerical simulations, the high estimate of 24,000,000 vaccinated population is considered to account for any uncertainty in the estimate for the current total number of vaccinated individuals population in KSA. We simulate the results for the given 300 days as in Figs. 610.

Figure 6: Graphical presentation for different fractional order of the susceptible class for the considered model

Figure 7: Graphical presentation for different fractional order of the exposed class for the considered model

Figure 8: Graphical presentation for different fractional order of the infected class for the considered model

Figure 9: Graphical presentation for different fractional order of the vaccinated class for the considered model

Figure 10: Graphical presentation for different fractional order of the recovered class for the considered model

We see that corresponding to different fractional order, the susceptible class is decreasing with various scenario in Fig. 6. Consequently the density of exposed class is also behave like susceptible at various fractional order, see Fig. 7. The infected class first increase then starts to decline (see Fig. 8) due to vaccination process until become stable as in Fig. 9. Thanking to vaccination the decline in infection will cause the increase in recovery class whose dynamics shows variation due to different fractional order as in Fig. 10. From these graphical presentations we observe the transmission dynamics in the presence of vaccine by using fractional order derivative. The fractional calculus helps us in better understanding of the transmission dynamics in a locality. Also the impact of vaccine is importance which understand from the recovery class presentation. Next, we compared the simulated plots of the population have been fully vaccinated in Fig. 11 with the graphs at different fractional order. We see that graphs are closely agreed which demonstrates the efficiency of our proposed model.

Figure 11: Comparison of real data with different fractional order simulation

8  Conclusion

In this research work, we have updated a COVID-19 model to a fractional order derivative of generalized type. On the statistical aspect, we have used some statistical analysis to collect data on vaccination in KSA for 300 days, and then the concerned statistical analysis has been shown. Consequently, the forecast about the evolution of the COVID-19 vaccination in 60 days has been presented. We have found, through the ARIMA(1, 1, 2) model, a decrease in the number of people who have been given full doses with (927164), and they constitute 2.6% of the total population in KSA. Data analysis showed that 67.81% of the population had been fully vaccinated during the study period. On the analytical aspect, we have established some adequate results for the existence and uniqueness of the solution through fixed point techniques. The respective results are important because a mathematical formulation should be preferentially checked for their existence. Finally, in terms of numerical aspects, we have extended the Adam-Bashforth method for the considered model to derive a scheme for numerical analysis. Moreover, we have then used the real data for parameters and some initial data of KSA to see the transmission dynamics of COVID-19 with the vaccinated class. Finally, the concerned numerical simulations have been compared with the exact real available date given in Section 2. We see that the real data plot and simulated data plots coincide very well. This phenomenon demonstrates the applicability of the method.

Authors’ Contributions: All authors contributed equally and significantly in writing this paper. All authors read and approved the final manuscript.

Acknowledgement: The authors would like to thank Imam Mohammad Ibn Saud Islamic University (IMSIU), for funding this research work.

Data Availability Statement: Data are available upon request.

Funding Statement: This research was supported by the Deanship of Scientific Research, Imam Mohammad Ibn Saud Islamic University (IMSIU), Saudi Arabia, Grant No. (21-13-18-069).

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

## References

1. Hu, B., Guo, H., Zhou, P., & Shi, Z. L. (2021). Characteristics of SARS-CoV-2 and COVID-19. Nature Reviews Microbiology, 19, 141-154. [Google Scholar] [CrossRef]
2. WHO (2021). World Health Organization Weakly Report. https://www.who.int/emergencies/diseases/novel-coronavirus-2019.
3. Worldometers.info (2021). Coronavirus cases-worldometer. https://www.worldometers.info/coronavirus.
4. FDA, U.S. FOOD & DRUG. COVID-19 Vaccines (2022). https://www.fda.gov/emergency-preparedness-and-response/coronavirus-disease-2019-covid-19/covid-19-vaccines.
5. Mathevet, T., Lepiller, M., & Mangin, A. (2004). Application of time series analyses to the hydrological functioning of an alpine karstic system: The case of Bange-L’Eua-Morte. Hydrology and Earth System Sciences, 8(6), 1051-1064. [Google Scholar] [CrossRef]
6. Khan, F., & Pilz, J. (2018). Modelling and sensitivity analysis of river flow in the upper Indus basin: Pakistan. International Journal of Water, 12(1), 1-21. [Google Scholar] [CrossRef]
7. Box, G., Jenkins, G. (1970). Time series analysis, forecasting and control. San Francisco: Holden-Day.
8. Maleki, M., Mahmoudi, M. R., Wraith, D., & Pho, K. H. (2020). Time series modelling to forecast the confirmed and recovered cases of COVID-19. Travel Medicine and Infectious Disease, 37, 101742. [Google Scholar] [CrossRef]
9. Podlubny, I. (1999). Fractional differential equations. San Diego: Academic Press.
10. Samko, S. G., Kilbas, A. A., Marichev, O. I. (1993). Fractional integrals and derivatives. Yverdon: Gordon & Breach.
11. Kilbas, A. A., Srivastava, H. M., Trujillo, J. J. (2006). Theory and applications of fractional differential equations. Amsterdam: North-Holland Mathematics Studies, Elsevier.
12. Ross, B. (2006). A brief history and exposition of the fundamental theory of fractional calculus. Fractional Calculus and Its Applications, 457, 1-36. [Google Scholar] [CrossRef]
13. Yang, X. J. (2019). General fractional derivatives: Theory, methods and applications. Chapman and Hall/CRC. DOI 10.1201/9780429284083. [CrossRef]
14. Toufik, M., & Atangana, A. (2017). New numerical approximation of fractional derivative with non-local and non-singular kernel: Application to chaotic models. The European Physical Journal Plus, 132(10), 1-16. [Google Scholar]
15. Koca, I. (2018). Modelling the spread of Ebola virus with Atangana-Baleanu fractional operators. The European Physical Journal Plus, 133(3), 100. [Google Scholar] [CrossRef]
16. Atangana, A., & Gomez-Aguilar, J. F. (2018). Decolonisation of fractional calculus rules: Breaking commutativity and associativity to capture more natural phenomena. The European Physical Journal Plus, 133(4), 166. [Google Scholar] [CrossRef]
17. Tateishi, A. A., Ribeiro, H. V., & Lenzi, E. K. (2017). The role of fractional time-derivative operators on anomalous diffusion. Frontiers in Physics, 52(5), 19. [Google Scholar] [CrossRef]
18. Caputo, M., & Fabrizio, M. (2015). A new definition of fractional derivative without singular kernel. Progress in Fractional Differentiation and Applications, 1(2), 73-85. [Google Scholar]
19. Losada, J., & Nieto, J. J. (2015). Properties of a new fractional derivative without singular kernel. Progress in Fractional Differentiation and Applications, 1(2), 87-92. [Google Scholar]
20. Atangana, A., & Baleanu, D. (2016). New fractional derivative with non-local and non-singular kernel. Thermal Science, 20(2), 757-763. [Google Scholar] [CrossRef]
21. Jarad, F., Abdeljawad, T., & Hammouch, Z. (2018). On a class of ordinary differential equations in the frame of Atangana-Baleanu fractional derivative. Chaos Solitons Fractals, 117, 16-20. [Google Scholar] [CrossRef]
22. Abdeljawad, T. (2017). A Lyapunov type inequality for fractional operators with nonsingular Mittag-Leffler kernel. Journal of Inequalities and Applications, 2017(1), 130. [Google Scholar] [CrossRef]
23. Atangana, A. (2018). Non validity of index law in fractional calculus: A fractional differential operator with Markovian and non-Markovian properties. Physica A: Statistical Mechanics and Its Applications, 505, 688-706. [Google Scholar] [CrossRef]
24. Atangana, A., & Gomez-Aguilar, J. F. (2018). Fractional derivatives with no-index law property: Application to chaos and statistics. Chaos Solitons Fractals, 114, 516-535. [Google Scholar] [CrossRef]
25. Abdo, M. S., Shah, K., Wahash, H. A., & Panchal, S. K. (2020). On a comprehensive model of the novel coronavirus (COVID-19) under Mittag-Leffler derivative. Chaos Solitons Fractals, 135, 109867. [Google Scholar] [CrossRef]
26. Almalahi, M. A., Panchal, S. K., Shatanawi, W., Abdo, M. S., & Shah, K. (2021). Analytical study of transmission dynamics of 2019-nCoV pandemic via fractal fractional operator. Results in Physics, 24, 104045. [Google Scholar] [CrossRef]
27. Jeelani, M. B., Alnahdi, A. S., Abdo, M. S., Abdulwasaa, M. A., & Shah, K. (2021). Mathematical modeling and forecasting of COVID-19 in Saudi Arabia under fractal-fractional derivative in caputo sense with power-law. Axioms, 10(3), 228. [Google Scholar] [CrossRef]
28. Abdulwasaa, M. A., Abdo, M. S., Shah, K., Nofal, T. A., & Panchal, S. K. (2021). Fractal-fractional mathematical modeling and forecasting of new cases and deaths of COVID-19 epidemic outbreaks in India. Results in Physics, 20, 103702. [Google Scholar] [CrossRef]
29. Atangana, A. (2020). Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination. Chaos Solitons Fractals, 136, 109860. [Google Scholar] [CrossRef]
30. Zhang, Z. (2020). A novel COVID-19 mathematical model with fractional derivatives: Singular and nonsingular kernels. Chaos Solitons Fractals, 139, 110060. [Google Scholar] [CrossRef]
31. Khan, A., Alshehri, H. M., Abdeljawad, T., Al-Mdallal, Q. M., & Khan, H. (2021). Stability analysis of fractional nabla difference COVID-19 model. Results in Physics, 22, 103888. [Google Scholar] [CrossRef]
32. Khan, Z. A., Khan, A., Abdeljawad, T., Khan, H. (2022). Computational analysis of fractional order imperfect testing infection disease model. Fractals. DOI 10.1142/S0218348X22401697. [CrossRef]
33. Ullah, I., Ahmad, S., Al-Mdallal, Q., Khan, Z. A., & Khan, H. (2020). Stability analysis of a dynamical model of tuberculosis with incomplete treatment. Advances in Difference Equations, 2020, 499. [Google Scholar] [CrossRef]
34. Shah, K., Khan, Z. A., Ali, A., Amin, R., & Khan, H. (2020). Haar wavelet collocation approach for the solution of fractional order COVID-19 model using caputo derivative. Alexandria Engineering Journal, 59(5), 3221-3231. [Google Scholar] [CrossRef]
35. Habenom, H., Aychluh, M., Suthar, D. L., Al-Mdallal, Q., & Purohit, S. D. (2022). Modeling and analysis on the transmission of COVID-19 pandemic in Ethiopia. Alexandria Engineering Journal, 61(7), 5323-5342. [Google Scholar] [CrossRef]
36. Bozkurt, F., Yousef, A., Abdeljawad, T., Kalinli, A., & Al Mdallal, Q. (2021). A fractional-order model of COVID-19 considering the fear effect of the media and social networks on the community. Chaos Solitons Fractals, 152, 111403. [Google Scholar] [CrossRef]
37. Sindhu, T. N., Shafiq, A., & Al-Mdallal, Q. M. (2021). On the analysis of number of deaths due to COVID-19 outbreak data using a new class of distributions. Results in Physics, 21, 103747. [Google Scholar] [CrossRef]
38. Shafiq, A., Lone, S. A., Sindhu, T. N., El Khatib, Y., & Al-Mdallal, Q. M. (2021). A new modified kies FrÃ©chet distribution: Applications of mortality rate of COVID-19. Results in Physics, 28, 104638. [Google Scholar] [CrossRef]
39. Diagne, M. L., Rwezaura, H., Tchoumi, S. Y., & Tchuenche, J. M. (2021). A mathematical model of COVID-19 with vaccination and treatment. Computational and Mathematical Methods in Medicine, 2021, [Google Scholar] [CrossRef]
40. Gokbulut, N., Kaymakamzade, B., Sanlidag, T., & Hincal, E. (2021). Mathematical modelling of COVID-19 with the effect of vaccine. AIP Conference Proceedings, 2325(1), 020065. [Google Scholar] [CrossRef]
41. Yavuz, M., CoşarÖ, F., Günay, F., & Özdemir, F. N. (2021). A new mathematical modeling of the COVID-19 pandemic including the vaccination campaign. Open Journal of Modelling and Simulation, 9(3), 299-321. [Google Scholar] [CrossRef]
42. covidvax.live (2021). http://covidvax.live/ar/location/sau.
43. Gómez-Aguilar, J. F., Rosales-García, J. J., Bernal-Alvarado, J. J., Córdova-Fraga, T., & Guzmán-Cabrera, R. (2012). Fractional mechanical oscillators. Revista Mexicana de Física, 58(4), 348-352. [Google Scholar]
44. Toufik, M., & Atangana, A. (2017). New numerical approximation of fractional derivative with non-local and non-singular kernel: Application to chaotic models. Chaos Solitons Fractals, 132, 444. [Google Scholar] [CrossRef]
45. Diethelm, K., & Freed, A. D. (1998). The FracPECE subroutine for the numerical solution of differential equations of fractional order. Forschung und Wissenschaftliches Rechnen, 1999, 57-71. [Google Scholar]

Alnahdi, A. S., Jeelani, M. B., Wahash, H. A., Abdulwasaa, M. A. (2023). A Detailed Mathematical Analysis of the Vaccination Model for COVID-19. CMES-Computer Modeling in Engineering & Sciences, 135(2), 1315–1343.

View

Like

## Related articles

• Haijian Shao, Xing Deng
• Jieren Cheng, Ruomeng Xu, Xiangyan...
• Baojia Wang, Pingzeng Liu, Zhang...
• Feng Guo, Pingzeng Liu, Wanming...
• Xiangyan Tang, Liang Wang, Jieren...