
Our aim of this work is to approximate the fixed points of generalized αm-nonexpansive mappings employing AA-iterative scheme in the structure of modular spaces. The results of fixed points for generalized αm-nonexpansive mappings is proven in this context. Moreover, the stability of the scheme and data dependence results are given for m-contraction mappings. In order to demonstrate that the AA-iterative scheme converges faster than some other schemes for generalized αm-nonexpansive mappings, numerical examples are shown at the end.
Citation: Muhammad Waseem Asghar, Mujahid Abbas, Cyril Dennis Enyi, McSylvester Ejighikeme Omaba. Iterative approximation of fixed points of generalized αm-nonexpansive mappings in modular spaces[J]. AIMS Mathematics, 2023, 8(11): 26922-26944. doi: 10.3934/math.20231378
[1] | Sunisa Theswan, Sotiris K. Ntouyas, Jessada Tariboon . Coupled systems of ψ-Hilfer generalized proportional fractional nonlocal mixed boundary value problems. AIMS Mathematics, 2023, 8(9): 22009-22036. doi: 10.3934/math.20231122 |
[2] | Saima Rashid, Fahd Jarad, Khadijah M. Abualnaja . On fuzzy Volterra-Fredholm integrodifferential equation associated with Hilfer-generalized proportional fractional derivative. AIMS Mathematics, 2021, 6(10): 10920-10946. doi: 10.3934/math.2021635 |
[3] | Saima Rashid, Abdulaziz Garba Ahmad, Fahd Jarad, Ateq Alsaadi . Nonlinear fractional differential equations and their existence via fixed point theory concerning to Hilfer generalized proportional fractional derivative. AIMS Mathematics, 2023, 8(1): 382-403. doi: 10.3934/math.2023018 |
[4] | Sotiris K. Ntouyas, Bashir Ahmad, Jessada Tariboon . Nonlocal integro-multistrip-multipoint boundary value problems for ¯ψ∗-Hilfer proportional fractional differential equations and inclusions. AIMS Mathematics, 2023, 8(6): 14086-14110. doi: 10.3934/math.2023720 |
[5] | Thabet Abdeljawad, Pshtiwan Othman Mohammed, Hari Mohan Srivastava, Eman Al-Sarairah, Artion Kashuri, Kamsing Nonlaopon . Some novel existence and uniqueness results for the Hilfer fractional integro-differential equations with non-instantaneous impulsive multi-point boundary conditions and their application. AIMS Mathematics, 2023, 8(2): 3469-3483. doi: 10.3934/math.2023177 |
[6] | Kanagaraj Muthuselvan, Baskar Sundaravadivoo, Kottakkaran Sooppy Nisar, Suliman Alsaeed . Discussion on iterative process of nonlocal controllability exploration for Hilfer neutral impulsive fractional integro-differential equation. AIMS Mathematics, 2023, 8(7): 16846-16863. doi: 10.3934/math.2023861 |
[7] | Weerawat Sudsutad, Jutarat Kongson, Chatthai Thaiprayoon, Nantapat Jarasthitikulchai, Marisa Kaewsuwan . A generalized Gronwall inequality via ψ-Hilfer proportional fractional operators and its applications to nonlocal Cauchy-type system. AIMS Mathematics, 2024, 9(9): 24443-24479. doi: 10.3934/math.20241191 |
[8] | Ayub Samadi, Sotiris K. Ntouyas, Jessada Tariboon . Coupled systems of nonlinear sequential proportional Hilfer-type fractional differential equations with multi-point boundary conditions. AIMS Mathematics, 2024, 9(5): 12982-13005. doi: 10.3934/math.2024633 |
[9] | Weerawat Sudsutad, Chatthai Thaiprayoon, Sotiris K. Ntouyas . Existence and stability results for ψ-Hilfer fractional integro-differential equation with mixed nonlocal boundary conditions. AIMS Mathematics, 2021, 6(4): 4119-4141. doi: 10.3934/math.2021244 |
[10] | Lakhlifa Sadek, Tania A Lazǎr . On Hilfer cotangent fractional derivative and a particular class of fractional problems. AIMS Mathematics, 2023, 8(12): 28334-28352. doi: 10.3934/math.20231450 |
Our aim of this work is to approximate the fixed points of generalized αm-nonexpansive mappings employing AA-iterative scheme in the structure of modular spaces. The results of fixed points for generalized αm-nonexpansive mappings is proven in this context. Moreover, the stability of the scheme and data dependence results are given for m-contraction mappings. In order to demonstrate that the AA-iterative scheme converges faster than some other schemes for generalized αm-nonexpansive mappings, numerical examples are shown at the end.
Everywhere in the world there have been disasters like; tsunamis, earthquakes, floods, droughts and pandemics since the beginning of the time. These disasters are likely to happen at different times and in different areas. For instance, at various points in history, humans have been afflicted with a variety of infectious diseases, including Ebola virus, HIV, dengue virus, bird flu, malaria, tuberculosis, influenza, diarrhea, hepatitis C and B, rubella disease and many others. These highly contagious diseases afflicted both humans and animals. In addition to this, these contagious diseases also spread among plants such as; Moko diseases, Pine wilt, orange rust, sugar cane, Dutchelm and Karnal bunt etc. diseases [1]. These widespread diseases not only left the bad impact on world economic situation, but they also results in the death of a huge number of people.
Mathematical models play a key role in the study of infectious diseases. Due to COVID-19 pandemic recently, a lot of work has been done and some good results have been achieved so for. The research on contagious diseases is vital for the better understanding of spread of these diseases and its behavior on biological entities. Human population is rapidly expanding, which has a negative impact on the nature. The spreading of these diseases are far more different than that of the population; as some of the diseases spread more rapidly and some of them spread with a slow rate. The first human case on COVID-19 was reported on 31st December, 2019 in Wuhan, China. First symptoms of COVID-19 were thought to be pneumonia. The infected patient did not recover after being vaccinated and treated for pneumonia, and the medication was ineffective [2]. Second, it was discovered that this virus spread quickly from person to person. Moreover, the infected patients were not restricted to Wuhan but this infectious disease further rapidly spread geographically within China [3,4]. The prevalence of COVID-19 virus was not restricted to Wuhan, it spread to other parts of the globe as well. Thousands of people have died as this deadly virus has spread rapidly throughout China and the rest of the countries, for instance North America, Italy, France and Germany etc. This infectious virus also hits Asia, particularly India and Pakistan. It was identified in Pakistan on February 27th, 2020 [5]. Studies on other contagious diseases like HIV, chickenpox, hepatitis E, dengue fever, rubella disease, measles and tuberculosis are studied in [6,7,8,9]. For comprehensive description of the vaccination game for mathematical epidemiology by analytical and multi-agent simulation approaches, we refer the reader to [10].
Recently re-emerging disease Monkeypox, the Ebola virus disease (EVD) and the influence of quadratic Lévy noise on the dynamic of an SIC contagious was studied in detail in [11,12,13]. For the understanding of the dynamics of the these epidemic diseases, they used some nonstandard analytical methods and find that the noise intensities are correlated with some parameters. They use real data to confirmed their theoretical justifications. Stochastic delay differential equations modeling these pandemic diseases with a time-periodic delay is analyzed using the spectral element approach in [14]. A comprehensive study consist of numerical solutions to stochastic differential equations can be found in [15,16,17,18,19,20,21,22,23,24,25]. Some recent work on the analysis of these epidemic models can also be found in [26,27,28,29]. As these epidemic models play an important role in helping to understand the spread of diseases and to plan control policies for diseases. Some researcher used statistical analysis of diseases data using epidemic models, which may not be very helpful for the above purpose as data may not be correct. In addition to this, the mathematical descriptions of these epidemic models are coupled and nonlinear in nature and it is very hard or in some cases it is not possible to find the exact solution. For this reason, one must use some numerical techniques for the approximate solution of these models. To the best of our knowledge, most of the authors used low order methods to find the numerical solution, which may not be helpful to study the long time behavior of these models. To this end, we use a higher order numerical scheme based on spectral methods to study the effects of different parameters on the COVID-19 epidemic SEIR model. Spectral methods are well known for their exponential rate of convergence and one can get a very accurate solution for relatively very small collocations points.
The remaining structure of the paper is: Mathematical model is presented in section 2, followed by spectral method in section 3. The stability analysis along with some main results are provided in section 4, followed by the brief numerical simulations in section 5. Conclusions are drawn in section 6.
Let N(t) denotes the overall population, which is further classified into four subclasses. Susceptible, exposed, infected, and recovered individuals S(t), E(t), I(t) and R(t) respectively. According to proposed COVID-19 model, the recruitment rate denoted by Λ of susceptible individual, η the death rate, ξ denote the interaction rate which is equal to ξSI/1+aI2(t), where κ denotes the rate of incubation period of the exposed class and enter into the infected class, δ represents the recovery rate. In mathematical modeling of epidemic diseases, incident rate play an important rule. In many classical models the bilinear rate is most frequently used, especially when the susceptible population is huge in number. For this reason we also choose the incident rate to be bilinear. Another important aspect of the stochastic models is the the persistence time in the mean of the stochastic model, which is the time during which an epidemic population survived. One can think about this in the context of COVID-19 pandemic. Researcher trying to figure out how long a population model persist, which is consider in detail in [30]. The deterministic COVID-19 model along with bilinear incidence rate is given by:
dS(t)dt=Λ−ξS(t)I(t)1+aI2(t)−ηS(t),dE(t)dt=ξS(t)I(t)1+aI2(t)−(η+κ)E(t),dI(t)dt=κE(t)−(η+δ)I(t),dR(t)dt=δI(t)−ηR(t), | (2.1) |
where the initial values S(0)=S∗(t)≥0,E(0)=E∗(t)≥0,I(0)=I∗(t)≥0 and R(0)=R∗(t)≥0.
The model given in Eq (2.1) in stochastic form is given below:
S(t)=(Λ−ξS(t)I(t)1+aI2(t)−ηS(t))dt+ϵS(t)dB(t),E(t)=(ξS(t)I(t)1+aI2(t)−(η+κ)E(t))dt+ϵE(t)dB(t),I(t)=(κE(t)−(η+δ)I(t))dt+ϵI(t)dB(t),R(t)=(δI(t)−ηR(t))dt+ϵR(t)dB(t), | (2.2) |
where ϵ is the intensity of B(t), and B(t) is a Wiener process. The aim of the present research work is to apply a reliable and stable numerical technique to the systems given in Eqs (2.1) and (2.2), along with theoretical justifications.
Let Qn(τa), denotes the Legendre polynomials of nth order. The function g(τa) given on the interval [−1,1] is approximated by
g(τa)=n∑i=0giQi(τa), | (3.1) |
Xi=g(τai) indicates Legendre coefficients which are unknowns, τai,i=0,1,...,n are interpolating points −1=τa0<τa1<...<τan=1. The Legendre polynomial is given by:
Qi(τa)=12i[i2]∑a=0(−1)a(ia)(2i−2ai)τi−2a,τa∈[−1,1](i=0,1,2,...,n) | (3.2) |
i2={i2,ifiEven,i−12,ifiOdd. |
The proposed numerical scheme for model Eq (2.2) is based on the usage of Lagendre-Gauss quadrature with a weight function and Legendre-Gauss Lobatto points {τk}nk=0. Now, to solve the proposed model Eq (2.2), we begin by integration from [0,t] of Eq (2.2), to get:
S(t)=S(0)+∫t0{Λ−ξS(s)I(s)1+aI2(t)−ηS(s)}ds+∫t0{ϵS(s)dB(s)},E(t)=E(0)+∫t0{ξS(s)I(s)1+aI2(t)−(η+κ)E(s)}ds+∫t0{ϵE(s)dB(s)},I(t)=I(0)+∫t0{κE(s)−(η+δ)I(s)}ds+∫t0{ϵI(s)dB(s)},R(t)=R(0)+∫t0{δI(s)−ηR(s)}ds+∫t0{ϵR(s)dB(s)}, | (3.3) |
where S(0),E(0),I(0),R(0) are initial conditions for each class respectively. Since the proposed scheme is orthogonal on a standard interval [−1,1]. Therefore, by using simple transformations s=t2(ϑ+1), Eq (3.3) becomes:
S(t)=S(0)+t2∫1−1{Λ−ξS(t2(ϑ+1))I(t2(ϑ+1))1+aI2(t)−ηS(t2(ϑ+1))}dϑ+t2∫1−1{ϵS(t2(ϑ+1))}dB(ϑ),E(t)=E(0)+t2∫1−1{ξS(t2(ϑ+1))I(t2(ϑ+1))1+aI2(t)−(η+κ)E(t2(ϑ+1))}dϑ+t2∫1−1{ϵE(t2(ϑ+1))}dB(ϑ),I(t)=I(0)+t2∫1−1{κE(t2(ϑ+1))−(η+δ)I(t2(ϑ+1))}dϑ+t2∫1−1{ϵI(t2(ϑ+1))}dB(ϑ),R(t)=R(0)+t2∫1−1{δI(t2(ϑ+1))−ηR(t2(ϑ+1))}dϑ+t2∫1−1{ϵR(t2(ϑ+1))}dB(ϑ). | (3.4) |
The semi-discretised spectral system with weights functions of Eq (3.4) is given by:
S(t)=S(0)+t2N∑k=0{Λ−ξS(t2(ϑk+1))I(t2(ϑk+1))1+aI2(t)−ηS(t2(ϑk+1))}ωk,+t2N∑k=0ϵS(t2(ϑk+1))ϖk,E(t)=E(0)+t2N∑k=0{ξS(t2(ϑk+1))I(t2(ϑk+1))1+aI2(t)−(η+κ)E(t2(ϑk+1))}ωk,+t2N∑k=0ϵE(t2(ϑk+1))ϖk,I(t)=I(0)+t2N∑k=0{κE(t2(ϑk+1))−(η+δ)I(t2(ϑk+1))}ωk+t2N∑k=0ϵI(t2(ϑk+1))ϖk,R(t)=R(0)+t2N∑k=0{δI(t2(ϑk+1))−ηR(t2(ϑk+1))}ωk+t2N∑k=0ϵR(t2(ϑk+1))ϖk, | (3.5) |
where deterministic weighted function is:
ωk={2(L′N+1(sk))2(1−s2k)},k=0,1,...,N. |
While for stochastic terms the weighted function is given by:
ϖk={√ωk×randn(1,N)},k=0,1,...,N. |
Here the numerical solutions for each class using Eq (3.1), is given as follow:
S(t)=N∑n=0SnPn(t),E(t)=N∑n=0InPn(t),I(t)=N∑n=0InPn(t),R(t)=N∑n=0RnPn(t), | (3.6) |
where Sn, En, In, Rn are unknown Legendre coefficients. Using the above approximation in a system Eq (3.5), one get:
N∑n=0SnPn(t)=t2N∑k=0{Λ−ξ∑Nn=0SnPn(ξk)∑Nn=0InPn(ξk)1+aI2(t)−ηN∑n=0SnPn(ξk)}ωk+tϵ2N∑k=0N∑n=0SnPn(ξk)ϖk+N∑n=0SnPn(0),N∑n=0EnPn(t)=t2N∑k=0{ξ∑Nn=0SnPn(ξk)∑Nn=0InPn(ξk)1+aI2(t)−(η+κ)N∑n=0EnPn(ξk)}ωk+tϵ2N∑k=0N∑n=0EnPn(ξk)ϖk+N∑n=0EnPn(0),N∑n=0InPn(t)=t2N∑k=0{κN∑n=0EnPn(ξk)−(η+δ)N∑n=0InPn(ξk)}ωk+tϵ2N∑k=0N∑n=0InPn(ξk)ϖk+N∑n=0InPn(0),N∑n=0RnPn(t)=t2N∑k=0{δN∑n=0InPn(ξk)−ηN∑n=0RnPn(ξk)}ωk+tϵ2N∑k=0N∑n=0RnPn(ξk)ϖk+N∑n=0RnPn(0), | (3.7) |
for simplification, we use ξk=[t2(ϑk+1)]. Hence, Eq (3.7) consists of (4N+4) unknowns with 4N nonlinear equations. Now, the initial conditions are given by:
N∑n=0SnPn(0)=ν1,N∑n=0EnPn(0)=ν2,N∑n=0InPn(0)=ν3,N∑n=0RnPn(0)=ν4, | (3.8) |
therefore, Eq (3.7) along with Eq (3.8) forms a system of (4N+4) non-linear equations. We first find these unknown values from Eqs (3.7) and (3.8), then using it in Eq (3.6) gives the solution to the system given in Eq (2.2).
In this section, we examine the stability of the proposed scheme for the COVID-19 systems given in Eqs (2.1) and (2.2). Keeping in mind that the system's essential reproduction number of the system in Eq (2.1) represented by R0 calculated using next generation method:
R0={Λκξη(η+δ)(η+κ)}. |
Theorem 4.1. If R0<1, then system Eq (2.1) satisfy the disease free equilibrium, (S∗(t),E∗(t),I∗(t),R∗(t))→(Λη,0,0,0). If R0>1, then proposed system has stable endemic equilibrium E∗=(S∗(t),E∗(t),I∗(t),R∗(t)), where
S∗(t)=(Λ+ξR∗(t))(1+aI2(t))(ξI∗(t)+η(1+aI2(t))),E∗(t)=(η+δ)I∗(t)κ,I∗(t)=κE∗(t)(η+δ),R∗(t)=δI∗(t)η. |
Proof. Consider E∗ be the endemic equilibrium of Eq (2.1), then S∗(t),E∗(t),I∗(t),R∗(t) the equations in the system are satisfied:
Λ−ξS(t)I(t)1+aI2(t)−ηS(t)=0,ξS(t)I(t)1+aI2(t)−(η+κ)E(t)=0,κE(t)−(η+δ)I(t)=0,δI(t)−ηR(t)=0. | (4.1) |
Here we consider two cases, that is I∗(t)=0 and I∗(t)>0. First, if I∗(t)=0, from the last equation of system Eq (4.1), we are simply get the recovery class becomes zero therefore, R∗(t)=0, using this in second equation of proposed system Eq (4.1), the exposed class become zero as, E∗(t)=0. Moreover, using the first equation of system Eq (4.1), we found the susceptible class as S∗(t)=Λ/η. Now we can get the complete disease free equilibrium, (S∗(t),E∗(t),I∗(t),R∗(t))→(Λη,0,0,0). If I∗(t)>0, for such calculation we simply use Maple-13 and found suitable endemic equilibrium.
Lemma 4.2. Assume the given region Π={(S(0),E(0),I(0),R(0))∈R4+:S(t)>0,E(t)>0,I(t)>0,>0,R(t)>0; S(t)+E(t)+I(t)+R(t)≤Λη} is a positively invariant set for system Eq (2.1).
Proof. The observed total population are denoted by N(t) by Eq (2.1), we obtain
dN(t)dt=Λ−ηN(t). |
This implies that
N(t)=N(0)e−ηt+Λη≤Λη. |
As S(t)+E(t)+I(t)+R(t)≤Λη, thus the region Π is positively invariant.
Lemma 4.3. Each given initial values {S(0),E(0),I(0),R(0)∈R4+}, then solution {S(t),E(t),I(t),R(t)∈R4+} of Eq (2.2) has the properties below:
limt→+∞∫t0γS(s)dB(s)t=0, |
for each class. For the proof one can see Meng et al. [31]
Definition 4.4. Given stochastic system Eq (2.2), the infected class I(t) is called extinctive if limt→∞I(t)=0.
The basic reproduction number for the stochastic system Eq (2.2) denoted by ˉR0 is given by:
ˉR0=R0−{ϵ2Λ22η2(η+δ)(η+κ)}. |
Theorem 4.5. If ϵ2>max{ξηΛ,ξ22(η+δ)(η+κ)} or ϵ2<ξηΛ with ˉR0<1, then the infected class of stochastic system Eq (2.2), exponentially approaches to zero.
Proof.} Let us assume {S(t),E(t),I(t),R(t)} be the solution of stochastic system Eq (2.2) and satisfied the initial values {S(0),E(0),I(0),R(0)}. Then according to the Itˆo formula, we have the following form:
dlnI(t)=(ξS(t)−(η+δ)(η+κ)−ϵ2S2(t)2)dt+ϵS(t)dB(t). | (4.2) |
Integral from 0 to t of system Eq (4.2) we have
lnI(t)=lnI(0)+∫t0(ξS(t)−(η+δ)(η+κ)−ϵ2S2(t)2)dt+∫t0ϵS(t)dB(t). | (4.3) |
Now we discuss the following two different cases, if ϵ2>ξηΛ, then
lnI(t)≤lnI(0)+(ξ2ϵ2−(η+δ)(η+κ))t+∫t0ϵS(t)dB(t). | (4.4) |
Divide Eq (4.4) by a positive t, then
lnI(t)t≤1tlnI(0)−((η+δ)(η+κ)−ξ2ϵ2)+1t∫t0ϵS(t)dB(t). | (4.5) |
After this taking the limt→∞ and by Lemma 4.3, then Eq (4.5) converted to:
limt→∞lnI(t)t≤−((η+δ)(η+κ)−ξ2ϵ2)<0. |
Therefore when ϵ2>ξ22(η+δ)(η+κ), which show that limt→∞I(t)=0.
The second case if ϵ2<ξηΛ, then from Eq (4.3), we have
lnI(t)≤lnI(0)+(ξΛη−ϵ2Λ22η2−(η+δ)(η+κ))t+∫t0ϵS(t)dB(t). | (4.6) |
Dividing Eq (4.6) both sides by t>0, then we obtain
1tlnI(t)≤1tlnI(0)+(η+δ)(η+κ)(ξΛη(η+δ)(η+κ)−ϵ2Λ22η2(η+δ)(η+κ)−1)+1t∫t0ϵS(t)dB(t). | (4.7) |
Again by limt→∞ and use the Lemma 4.3, Eq (4.7) becomes:
limt→∞lnIh(t)t≤((η+δ)(η+κ))(ˉR0−1). |
For ˉR0<1, then
limt→∞lnI(t)t<0. |
This shows that limt→∞I(t)=0.
Theorem 4.6. The infected population of humans and birds, therefore I(t) respectfully are presents in Eq (2.2). For proof, we refer the reader [30].
In this section, we discuss the stochastic model numerically using the Legendre spectral collocation method (LSCM) for both deterministic and stochastic systems. Figures 1–8 demonstrate the approximate results to the proposed models. Using LSCM, we can precisely describe the behavior of the fundamental reproduction number R0. We established several parameters values for the current systems in order to demonstrate the simple and effectiveness of the suggested strategy for stochastic COVID-19. Our initial conditions are all set to 1. For the sake of simplification, we fix the noise intensity ϵ from 0 to 0.6 in all our simulations. This is also due to fact that noise intensity dramatically change the dynamics of stochastic models, which are used for mathematical description of epidemic diseases. If the intensity is large, it can shift the bifurcation points and produce noise-induced transitions and making it very hard to find the bifurcation boundaries.
Figure 1 is obtained using the parameters values Λ=4;a=1;η=1;ξ=1;κ=0.6;δ=1;ϵ=0, in deterministic system given in Eq (2.1). For above parameter values the reproduction numbers R0<1, thus using Theorem 4.1, we can surly see in Figure 1, that the proposed systems consists disease free equilibrium E0, that is (S∗(t),E∗(t),I∗(t),R∗(t))→(Λη,0,0,0), where the susceptible class is Λ/η=4. In Figure 2, the parameter values are Λ=4;a=1;η=1;ξ=1.7;κ=0.6;δ=1;ϵ=0, while he reproduction number for these parameters is R0=1.275>1. Now again applying Theorem 4.1, the system Eq (2.1), has stable endemic equilibrium E∗. We clearly see that all the classes are non-zero.
For Parameter values Λ=4;a=1;η=1;ξ=1;κ=0.6;δ=1;ϵ=0.6, we obtain Figure 3. Using this the stochastic reproduction number becomes ˉR0<1. Now by applying the Theorem 4.6, we can clearly observe from Figure 3, that the systems Eq (2.2) has stable disease free equilibrium, where the susceptible class is Λ/η=4 and all other classes are zeros. Using same parameter values given in Figure 4 and only choose ξ=1.7, for the solution to the system Eq (2.2). Simulation present that ˉR0>1, again by using Theorem 4.6, we see that given systems has stable endemic equilibrium and all classes of the system Eq (2.2) are non zero. Taking the parameters values Λ=4;a=1;η=1;ξ=1;κ=0.6;δ=1;ϵ=0.6, we get Figure 5. We compare the numerical solution for each class of the systems Eqs (2.1) and (2.2). By applying the Theorem 4.6 again to the parameter values R0<1 and ˉR0<1, we can see from Figure 5 that both systems have a stable disease-free equilibrium E0.
Figure 5 shows a comparison of the numerical solutions of deterministic and stochastic COVID-19 systems using the parameter values listed above. Using the above parameter values the simple calculation appears that and both the reproduction numbers are less then 1 clearly observe in figure 6. A stable endemic equilibrium E0 may be found utilizing the preceding results. For parameter values, Λ=4;a=1;η=1;ξ=1;κ=0.6;δ=1;ϵ=0.6 in Figure 7, we plot both humans and bird classes and also we compare the numerical solution for class to the systems Eqs (2.1) and (2.2) of COVID-19 model. With these parameters, both reproduction numbers are greater than 1, the system is in a stable endemic equilibrium, as seen in figure 7. We draw the comparison of the numerical solutions to Eqs (2.1) and (2.2), using same parameter values given in Figure 7. Since both reproduction numbers in these models are greater than 1, they have a stable infectious equilibrium E∗ as shown in Figure 8.
A mathematical concept called differentiation and integration is used to model the dynamic behavior of a disease, like COVID-19. These models can be used to forecast how a disease will propagate in a particular population like how many people will get sick by splitting the overall populations into several subgroups for this assignment. Two equilibrium points DFE and EE are derived for the proposed steady state models. The fundamental reproduction number R0 is determined using the next generation approach. We find through stability assessment that the provided model is locally asymptotically stable when R0<1. A number of numerical simulations are performed to confirm the theoretical justifications.
This work was financially supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [Project No. GRANT784].
The authors declare that they have no competing interests.
[1] |
M. Abbas, M. W. Asghar, M. De la Sen, Approximation of the solution of delay fractional differential equation using AA-iterative scheme, Mathematics, 10 (2022), 273. https://doi.org/10.3390/math10020273 doi: 10.3390/math10020273
![]() |
[2] | M. Abbas, T. Nazir, A new faster iteration process applied to constrained minimization and feasibility problems, Mat. Vesn., 66 (2014), 223–234. |
[3] |
A. A. N. Abdou, M. A. Khamsi, Fixed point theorems in modular vector spaces, J. Nonlinear Sci. Appl., 10 (2017), 4046–4057. http://doi.org/10.22436/jnsa.010.08.01 doi: 10.22436/jnsa.010.08.01
![]() |
[4] | R. P. Agarwal, D. O'Regan, D. R. Sahu, Iterative construction of fixed points of nearly asymptotically nonexpansive mappings, J. Nonlinear Convex Anal., 8 (2007), 61–79. |
[5] |
S. Banach, Surles opérations dans les ensembles abstraites et leurs applications aux équations intégrales, Fund. Math., 3 (1922), 133–181. https://doi.org/10.4064/fm-3-1-133-181 doi: 10.4064/fm-3-1-133-181
![]() |
[6] |
V. Berinde, Picard iteration converges faster than mann iteration for a class of quasi-contractive operators, Fixed Point Theory Appl., 2004 (2004), 716359. https://doi.org/10.1155/S1687182004311058 doi: 10.1155/S1687182004311058
![]() |
[7] | R. E. Bruck, Asymptotic behavior of nonexpansive mappings, Contemp. Math., 18 (1983), 1–47. |
[8] |
C. Byrne, Iterative oblique projection onto convex sets and the split feasibility problem, Inverse Probl., 18 (2002), 441. http://doi.org/10.1088/0266-5611/18/2/310 doi: 10.1088/0266-5611/18/2/310
![]() |
[9] |
C. Byrne, A unified treatment of some iterative algorithms in signal processing and image reconstruction, Inverse Probl., 20 (2003), 103. http://doi.org/10.1088/0266-5611/20/1/006 doi: 10.1088/0266-5611/20/1/006
![]() |
[10] |
Y. Censor, T. Bortfeld, B. Martin, A. Trofimov, A unified approach for inversion problems in intensity-modulated radiation therapy, Phys. Med. Biol., 51 (2006), 2353. http://doi.org/10.1088/0031-9155/51/10/001 doi: 10.1088/0031-9155/51/10/001
![]() |
[11] |
Y. Censor, T. Elfving, N. Kopf, T. Bortfeld, The multiple-sets split feasibility problem and its applications for inverse problems, Inverse Probl., 21 (2005), 2071. http://doi.org/10.1088/0266-5611/21/6/017 doi: 10.1088/0266-5611/21/6/017
![]() |
[12] |
P. L. Combettes, The convex feasibility problem in image recovery, Adv. Imag. Elect. Phys., 95 (1996), 155–270. https://doi.org/10.1016/S1076-5670(08)70157-5 doi: 10.1016/S1076-5670(08)70157-5
![]() |
[13] | O. Fouad, A. Radouane, M. Driss, A fixed point theorem on some multi-valued maps in modular spaces, Nonlinear Funct. Anal. Appl., 2022. |
[14] |
J. García-Falset, E. Llorens-Fuster, T. Suzuki, Fixed point theory for a class of generalized nonexpansive mappings, J. Math. Anal. Appl., 375 (2011), 185–195. https://doi.org/10.1016/j.jmaa.2010.08.069 doi: 10.1016/j.jmaa.2010.08.069
![]() |
[15] |
E. Hacıoǧlu, F. Gürsoy, A. R. Khan, Monotone generalized α-nonexpansive mappings on CATp(0) spaces, Ukrains'kyi Matematychnyi Zhurnal, 75 (2023), 970–986. https://doi.org/10.37863/umzh.v75i7.7188 doi: 10.37863/umzh.v75i7.7188
![]() |
[16] |
S. Ishikawa, Fixed points by a new iteration method, Proc. Amer. Math. Soc., 44 (1974), 147–150. https://doi.org/10.2307/2039245 doi: 10.2307/2039245
![]() |
[17] |
Z. J. Kadum, N. Y. Abdul-Hassan, New numerical methods for solving the initial value problem based on a symmetrical quadrature integration formula using hybrid functions, Symmetry, 15 (2023), 631. https://doi.org/10.3390/sym15030631 doi: 10.3390/sym15030631
![]() |
[18] |
W. Kassab, T. ţurcanu, Numerical reckoning fixed points of (ρE)-type mappings in modular vector spaces, Mathematics, 7 (2019), 390. https://doi.org/10.3390/math7050390 doi: 10.3390/math7050390
![]() |
[19] | S. H. Khan, Approximating fixed points of (λ, ρ)-firmly nonexpansive mappings in modular function spaces, arXiv, 2018. https://doi.org/10.48550/arXiv.1802.00681 |
[20] |
M. A. Khamsi, W. M. Kozlowski, S. Reich, Fixed point theory in modular function spaces, Nonlinear Anal., 14 (1990), 935–953. https://doi.org/10.1016/0362-546X(90)90111-S doi: 10.1016/0362-546X(90)90111-S
![]() |
[21] | M. A. Khamsi, W. M. Kozlowski, Fixed point theory in modular function spaces, Springer, 2015. https://doi.org/10.1007/978-3-319-14051-3 |
[22] |
W. M. Kozlowski, Notes on modular projections, Appl. Set-Valued Anal. Optim., 4 (2022), 337–348. https://doi.org/10.23952/asvao.4.2022.3.07 doi: 10.23952/asvao.4.2022.3.07
![]() |
[23] |
A. Latif, R. F. Al Subaie, M. O. Alansari, Fixed points of generalized multi-valued contractive mappings in metric type spaces, J. Nonlinear Var. Anal., 6 (2022), 123–138. https://doi.org/10.23952/jnva.6.2022.1.07 doi: 10.23952/jnva.6.2022.1.07
![]() |
[24] |
G. López, V. Martín-Márquez, F. Wang, H. K. Xu, Solving the split feasibility problem without prior knowledge of matrix norms, Inverse Probl., 28 (2012), 085004. https://doi.org/10.1088/0266-5611/28/8/085004 doi: 10.1088/0266-5611/28/8/085004
![]() |
[25] |
W. R. Mann, Mean value methods in iteration, Proc. Amer. Math. Soc., 4 (1953), 506–510. https://doi.org/10.1090/S0002-9939-1953-0054846-3 doi: 10.1090/S0002-9939-1953-0054846-3
![]() |
[26] | J. Musielak, Orlicz spaces and modular spaces, 1983, Springer. https://doi.org/10.1007/BFb0072210 |
[27] |
M. A. Noor, New approximation schemes for general variational inequalities, J. Math. Anal. Appl., 251 (2000), 217–229. https://doi.org/10.1006/jmaa.2000.7042 doi: 10.1006/jmaa.2000.7042
![]() |
[28] |
G. A. Okeke, D. Franci, J. K. Kim, New proofs of some fixed point thorems for mappings satisfying Reich type contraction in modular metric spaces, Nonlinear Funct. Anal. Appl., 2023. https://doi.org/10.22771/nfaa.2023.28.01.01 doi: 10.22771/nfaa.2023.28.01.01
![]() |
[29] |
A. M. Ostrowski, The round-off stability of iterations, ZAMM J. Appl. Math. Mech., 47 (1967), 77–81. https://doi.org/10.1002/zamm.19670470202 doi: 10.1002/zamm.19670470202
![]() |
[30] |
R. Pant, R. Shukla, Approximating fixed points of generalized α-nonexpansive mappings in Banach spaces, Numer. Funct. Anal. Optim., 38 (2017), 248–266. https://doi.org/10.1080/01630563.2016.1276075 doi: 10.1080/01630563.2016.1276075
![]() |
[31] |
B. Qu, N. Xiu, A note on the CQ algorithm for the split feasibility problem, Inverse Probl., 21 (2005), 1655. https://doi.org/10.1088/0266-5611/21/5/009 doi: 10.1088/0266-5611/21/5/009
![]() |
[32] |
A. Rahimi, A. Rezaei, B. Daraby, M. Ghasemi, A new faster iteration process to fixed points of generalized α-nonexpansive mappings in Banach spaces, Int. J. Nonlinear Anal. Appl., 2023. https://doi.org/10.22075/IJNAA.2021.23025.2460 doi: 10.22075/IJNAA.2021.23025.2460
![]() |
[33] |
H. J. Saeed, A. H. Ali, R. Menzer, A. D. Poţclean, H. Arora, New family of multi-step iterative methods based on homotopy perturbation technique for solving nonlinear equations, Mathematics, 11 (2023), 2603. https://doi.org/10.3390/math11122603 doi: 10.3390/math11122603
![]() |
[34] |
Ş. M. Şoltuz, T. Grosan, Data dependence for ishikawa iteration when dealing with contractive-like operators, Fixed Point Theory Appl., 2008 (2008), 242916. https://doi.org/10.1155/2008/242916 doi: 10.1155/2008/242916
![]() |
[35] |
B. S. Thakur, D. Thakur, M. Postolache, A new iteration scheme for approximating fixed points of nonexpansive mappings, Filomat, 30 (2016), 2711–2720. http://doi.org/10.2298/FIL1610711T doi: 10.2298/FIL1610711T
![]() |
[36] |
K. Ullah, M. Arshad, Numerical reckoning fixed points for suzuki's generalized nonexpansive mappings via new iteration process, Filomat, 32 (2018), 187–196. http://doi.org/10.2298/FIL1801187U doi: 10.2298/FIL1801187U
![]() |
[37] |
H. K. Xu, Iterative methods for the split feasibility problem in infinite-dimensional Hilbert spaces, Inverse Probl., 26 (2010), 105018. https://doi.org/10.1088/0266-5611/26/10/105018 doi: 10.1088/0266-5611/26/10/105018
![]() |
[38] |
Q. Yang, The relaxed CQ algorithm solving the split feasibility problem, Inverse Probl., 20 (2004), 1261. https://doi.org/10.1088/0266-5611/20/4/014 doi: 10.1088/0266-5611/20/4/014
![]() |
1. | Ricardo Almeida, Minimization Problems for Functionals Depending on Generalized Proportional Fractional Derivatives, 2022, 6, 2504-3110, 356, 10.3390/fractalfract6070356 | |
2. | Kishor D. Kucche, Ashwini D. Mali, Arran Fernandez, Hafiz Muhammad Fahad, On tempered Hilfer fractional derivatives with respect to functions and the associated fractional differential equations, 2022, 163, 09600779, 112547, 10.1016/j.chaos.2022.112547 | |
3. | Jessada Tariboon, Ayub Samadi, Sotiris K. Ntouyas, Nonlocal Boundary Value Problems for Hilfer Generalized Proportional Fractional Differential Equations, 2022, 6, 2504-3110, 154, 10.3390/fractalfract6030154 | |
4. | Surang Sitho, Sotiris K. Ntouyas, Chayapat Sudprasert, Jessada Tariboon, Integro-Differential Boundary Conditions to the Sequential ψ1-Hilfer and ψ2-Caputo Fractional Differential Equations, 2023, 11, 2227-7390, 867, 10.3390/math11040867 | |
5. | Sotiris Ntouyas, Bashir Ahmad, Jessada Tariboon, Nonlocal ψ-Hilfer Generalized Proportional Boundary Value Problems for Fractional Differential Equations and Inclusions, 2022, 2, 2673-9321, 377, 10.3390/foundations2020026 | |
6. | Weerawat Sudsutad, Chatthai Thaiprayoon, Bounmy Khaminsou, Jehad Alzabut, Jutarat Kongson, A Gronwall inequality and its applications to the Cauchy-type problem under ψ-Hilfer proportional fractional operators, 2023, 2023, 1029-242X, 10.1186/s13660-023-02929-x | |
7. | Ravi P. Agarwal, Snezhana Hristova, Boundary Value Problem for Multi-Term Nonlinear Delay Generalized Proportional Caputo Fractional Differential Equations, 2022, 6, 2504-3110, 691, 10.3390/fractalfract6120691 | |
8. | Sotiris K. Ntouyas, Bashir Ahmad, Jessada Tariboon, Nonlocal integro-multistrip-multipoint boundary value problems for ¯ψ∗-Hilfer proportional fractional differential equations and inclusions, 2023, 8, 2473-6988, 14086, 10.3934/math.2023720 | |
9. | Sunisa Theswan, Sotiris K. Ntouyas, Jessada Tariboon, Coupled systems of ψ-Hilfer generalized proportional fractional nonlocal mixed boundary value problems, 2023, 8, 2473-6988, 22009, 10.3934/math.20231122 | |
10. | Sotiris K. Ntouyas, Bashir Ahmad, Jessada Tariboon, Coupled Systems of Nonlinear Proportional Fractional Differential Equations of the Hilfer-Type with Multi-Point and Integro-Multi-Strip Boundary Conditions, 2023, 3, 2673-9321, 241, 10.3390/foundations3020020 | |
11. | Bashir Ahmad, Sotiris K. Ntouyas, 2024, Chapter 1, 978-3-031-62512-1, 1, 10.1007/978-3-031-62513-8_1 | |
12. | Najeeb Alam Khan, Oyoon Abdul Razzaq, Areeba Azam Bajwa, Muhammad Ayaz, Global dynamics and Impact of Gaussian noise intensity on the stochastic epidemic model with local fractional derivative, 2023, 98, 0031-8949, 064002, 10.1088/1402-4896/accf46 | |
13. | Ayub Samadi, Sotiris K. Ntouyas, Jessada Tariboon, Fractional Sequential Coupled Systems of Hilfer and Caputo Integro-Differential Equations with Non-Separated Boundary Conditions, 2024, 13, 2075-1680, 484, 10.3390/axioms13070484 | |
14. | Aurelian Cernea, On a Hilfer generalized proportional fractional integro-differential inclusion, 2024, 58, 24058963, 401, 10.1016/j.ifacol.2024.08.224 | |
15. | Sotiris K. Ntouyas, Phollakrit Wongsantisuk, Ayub Samadi, Jessada Tariboon, Hybrid System of Proportional Hilfer-Type Fractional Differential Equations and Nonlocal Conditions with Respect to Another Function, 2024, 12, 2227-7390, 1071, 10.3390/math12071071 | |
16. | Ayub Samadi, Sotiris K. Ntouyas, Jessada Tariboon, Coupled systems of nonlinear sequential proportional Hilfer-type fractional differential equations with multi-point boundary conditions, 2024, 9, 2473-6988, 12982, 10.3934/math.2024633 | |
17. | Weerawat Sudsutad, Jutarat Kongson, Chatthai Thaiprayoon, On generalized (k,ψ)-Hilfer proportional fractional operator and its applications to the higher-order Cauchy problem, 2024, 2024, 1687-2770, 10.1186/s13661-024-01891-x | |
18. |
James Viji, Velu Muthulakshmi, Pushpendra Kumar,
Oscillatory behavior of ΨΨ ‐Hilfer generalized proportional fractional initial value problems,
2024,
0170-4214,
10.1002/mma.10557
|
|
19. | V. M. Bulavatsky, On Some Generalizations of the Bi-Ordinal Hilfer’s Fractional Derivative, 2024, 60, 1060-0396, 541, 10.1007/s10559-024-00694-3 | |
20. | Weerawat Sudsutad, Jutarat Kongson, Chatthai Thaiprayoon, Nantapat Jarasthitikulchai, Marisa Kaewsuwan, A generalized Gronwall inequality via ψ-Hilfer proportional fractional operators and its applications to nonlocal Cauchy-type system, 2024, 9, 2473-6988, 24443, 10.3934/math.20241191 | |
21. | Ayub Samadi, Sotiris K. Ntouyas, Jessada Tariboon, Mixed Hilfer and Caputo Fractional Riemann–Stieltjes Integro-Differential Equations with Non-Separated Boundary Conditions, 2024, 12, 2227-7390, 1361, 10.3390/math12091361 | |
22. | Feryal Aladsani, Ahmed Gamal Ibrahim, Existence and Stability of Solutions for p-Proportional ω-Weighted κ-Hilfer Fractional Differential Inclusions in the Presence of Non-Instantaneous Impulses in Banach Spaces, 2024, 8, 2504-3110, 475, 10.3390/fractalfract8080475 | |
23. | Ayub Samadi, Sotiris K. Ntouyas, Asawathep Cuntavepanit, Jessada Tariboon, Hilfer proportional nonlocal fractional integro-multipoint boundary value problems, 2023, 21, 2391-5455, 10.1515/math-2023-0137 | |
24. | V.M. Bulavatsky, ON SOME GENERALIZATIONS OF THE BI-ORDINAL HILFER’S FRACTIONAL DERIVATIVE, 2024, 10195262, 36, 10.34229/KCA2522-9664.24.4.3 | |
25. | Weerawat Sudsutad, Chatthai Thaiprayoon, Aphirak Aphithana, Jutarat Kongson, Weerapan Sae-dan, Qualitative results and numerical approximations of the (k,ψ)-Caputo proportional fractional differential equations and applications to blood alcohol levels model, 2024, 9, 2473-6988, 34013, 10.3934/math.20241622 | |
26. |
Alaeddine Sahraoui, Zaid Laadjal, Qinghua Ma, Lamine Guedda,
On a Coupled System of ψψ ‐Liouville–Caputo‐Type Fractional Nonlinear Volterra–Fredholm Integro‐Differential Equations With Nonlocal Nonseparated Boundary Conditions,
2025,
0170-4214,
10.1002/mma.11082
|
|
27. | J. O. González-Cervantes, J. Bory-Reyes, B. Schneider, On a Proportional Fractional Holomorphic Functions of Two Complex Variables with Respect to Another Function, 2025, 46, 1995-0802, 416, 10.1134/S1995080224605198 | |
28. | Safoura Rezaei Aderyani, Reza Saadati, Donal O’Regan, Representation Formulas and Stability Analysis for Hilfer–Hadamard Proportional Fractional Differential Equations, 2025, 9, 2504-3110, 359, 10.3390/fractalfract9060359 |