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

Abeer S.Alnahdi,Mdi B.Jeelani,Hanan A.Wahash and Mansour A.Abdulwasaa

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

2Department of Mathematics,Dr.Babasaheb Ambedkar Marathwada University,Aurangabad,India

3Department of Statistics,Taiz University,Taiz,Yemen

4Department of Statistics,Dr.Babasaheb Ambedkar Marathwada University,Aurangabad,India

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 the COVID-19 vaccination in 60 days is presented.The theoretical aspect provides equilibrium points,reproduction number R0,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 COVID-19; Eviews program; forecasting; ABC fractional derivative; Picard-Lindel method; Adams-Bashforth technique

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,SARSCoV-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[5–8].

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[9–13].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 [14–17].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 [25–32].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[33–38].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.

Table 1:The number of people fully vaccinated from June 01,2021 to February 15,2022 in KSA[42]

Figure 1:Number of people fully vaccinated in the KSA during the period June 01,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.

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

Table 2: New doses for people from January 06,2021 to February 15,2022 in KSA[42]

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.

Table 4: Test result of ARIMA(p,d,q)models

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.

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.

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

Table 5: Test results of the predictive ability of a linear completed vaccinations in KSA

Table 6: Test results of the predictive ability of a linear completed vaccinations in KSA

Table 6 (continued)DateForecastUCLLCLDateForecast UCLLCL 2/19/20221329533778 38673/21/202215650655081485 2/20/20221355835452 37453/22/202215724665241446 2/21/20221376936936 36293/23/202215797675401408 2/22/20221381237925 34853/24/202215872685571372 2/23/20221385138903 33483/25/202215946695751337 2/24/20221391839957 32233/26/202216021705921303 2/25/20221398741014 31073/27/202216096716101270 2/26/20221405342053 29973/28/202216172726291238 2/27/20221411843086 28923/29/202216248736481207 2/28/20221418444119 27923/30/202216324746681178 3/1/20221425145149 26983/31/202216400756881149 3/2/20221431846176 26094/1/202216477767091121 3/3/20221438547201 25244/2/202216555777311094 3/4/20221445348224 24424/3/202216632787531067 3/5/20221452049246 23654/4/202216711797751042 3/6/20221458950267 22914/5/202216789807991017 3/7/20221465751286 22204/6/20221686881822993 3/8/20221472652304 21534/7/20221694782847970 3/9/20221479553321 20884/8/20221702683872947 3/10/20221486454338 20264/9/20221710684898926 3/11/20221493455354 19674/10/20221718785924904 3/12/20221500456370 19104/11/20221726786951884 3/13/20221507557386 18554/12/20221734886951863 3/14/20221514658401 18024/13/20221743089007844 3/15/20221521759416 17524/14/20221751290035825 3/16/20221528860431 17034/15/20221759491065807 3/17/20221536061446 1656

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

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

Definition 3.2.[20]The correspondentABfractional integral is given by

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

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

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

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:

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 parameterKwith the dimension ofs,to change the nonlocal fractional operator so that the sides possess the same dimension[43].Herewith,we reformulate the following fractional system:

With initial conditions

whereABCDΘ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().

Table 7: Description of parameters and numerical values of the models(3)and(4)

4.1 Positivity and Boundedness

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

Theorem 4.1.There is a unique solution)=(S(),E(),I(),V(),R())Tfor the models(3)and(4)on≥0 in[0,T](0

Proof.From the model(3),we have

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

ForK >0 and0,we have

This shows that if S0,E0,I0,V0,and R0are 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

It follows for the whole population that

which implies that

ABCDΘN()≤K1?Θ(Λ?μN(yùn)()).

After applying the Laplace transform,we have

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

whereEα,β(z)=andSubsequently,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).

ForK >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:

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

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 E0of the epidemic model is locally asymptotically stable ifR0<1,otherwise unstable.

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

Here we comment thatR0gives us data about the transmission of infectious sickness.For example in the normal utilization of contamination models,whenR0>1,the disease infection will actually be able to begin transmitting in society,yet ifR0<1 it implies that the disease infection continues declining in the general public.Overall,the bigger the worth ofR0the 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

has a unique solution,which is

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

It follows from Theorem 3 in[20]that

which is equivalent to

Applying the inverse Laplace transform give us

Now,we have

Let F(?)=and G(?)=L{f()}(?).It follows form convolution theorem that

Hence

By using Theorem 5.1,our model is equivalent to

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

Taking the limit asn→∞,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:

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

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

By applying the operatorAB?Θon Eq.(11),we obtain

The kernels in Eq.(10)satisfies the Lipschitz condition for 0≤Li <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

whereL1:=α||E||+m+μ.Thus,

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

whereL2:=α||S||+p||V||+κ||I||+c+μ,L3:=κ||E||+z+μ+σ,L4:=p||E||+μ,andL5:=μ.

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

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

Taking into account that

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

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

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

Consequently,the sequences in Eq.(21)are exist and smooth,i.e.,→0,(i=1,2,...,5)asn→∞.

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

Then we show that the terms in Eq.(22)verify that→0,asn→∞,(i=1,2,...,5).Thus,we have

Repeating the procedure above,we get

Asn→∞,→0.By the same arguments above,it can be shown thatasn→∞,(i=2,3,...,5).

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

which implies

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

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 operatorAB?Θon the first equation of(11),we obtain

For=+1,n=0,1,2,...,we get

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

which implies

Finally,

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

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.6–10.

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

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 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

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.