\left(jE^{m}\right) _{M}^{N} | N=M=20 | N=M=40 | N=M=80 |
m=1 | 0.0068, j=6 | 0.0032, j=6 | 0.0016, j=6 |
m=2 | 0.0071, j=6 | 0.0033, j=6 | 0.0016, j=6 |
m=3 | 0.0073, j=6 | 0.0034, j=6 | 0.0017, j=6 |
In this study, we explore a stochastic age-dependent cooperative Lotka-Volterra (LV) system with an environmental noise. By applying the theory of M-matrix, we prove the existence and uniqueness of the global solution for the system. Since the stochastic age-dependent cooperative LV system cannot be solved explicitly, we then construct an Euler-Maruyama (EM) numerical solution to approach the exact solution of the system. The convergence rate and the pth-moment boundedness of the scheme have also been obtained. Additionally, numerical experiments have been conducted to verify our theoretical results.
Citation: Mengqing Zhang, Qimin Zhang, Jing Tian, Xining Li. The asymptotic stability of numerical analysis for stochastic age-dependent cooperative Lotka-Volterra system[J]. Mathematical Biosciences and Engineering, 2021, 18(2): 1425-1449. doi: 10.3934/mbe.2021074
[1] | Thomas Torku, Abdul Khaliq, Fathalla Rihan . SEINN: A deep learning algorithm for the stochastic epidemic model. Mathematical Biosciences and Engineering, 2023, 20(9): 16330-16361. doi: 10.3934/mbe.2023729 |
[2] | Guo Lin, Shuxia Pan, Xiang-Ping Yan . Spreading speeds of epidemic models with nonlocal delays. Mathematical Biosciences and Engineering, 2019, 16(6): 7562-7588. doi: 10.3934/mbe.2019380 |
[3] | Lin Zhang, Yongbin Ge, Zhi Wang . Positivity-preserving high-order compact difference method for the Keller-Segel chemotaxis model. Mathematical Biosciences and Engineering, 2022, 19(7): 6764-6794. doi: 10.3934/mbe.2022319 |
[4] | Xiao-Min Huang, Xiang-ShengWang . Traveling waves of di usive disease models with time delay and degeneracy. Mathematical Biosciences and Engineering, 2019, 16(4): 2391-2410. doi: 10.3934/mbe.2019120 |
[5] | Azmy S. Ackleh, Mark L. Delcambre, Karyn L. Sutton, Don G. Ennis . A structured model for the spread of Mycobacterium marinum: Foundations for a numerical approximation scheme. Mathematical Biosciences and Engineering, 2014, 11(4): 679-721. doi: 10.3934/mbe.2014.11.679 |
[6] | Zijuan Wen, Meng Fan, Asim M. Asiri, Ebraheem O. Alzahrani, Mohamed M. El-Dessoky, Yang Kuang . Global existence and uniqueness of classical solutions for a generalized quasilinear parabolic equation with application to a glioblastoma growth model. Mathematical Biosciences and Engineering, 2017, 14(2): 407-420. doi: 10.3934/mbe.2017025 |
[7] | Shouying Huang, Jifa Jiang . Global stability of a network-based SIS epidemic model with a general nonlinear incidence rate. Mathematical Biosciences and Engineering, 2016, 13(4): 723-739. doi: 10.3934/mbe.2016016 |
[8] | Abdennasser Chekroun, Mohammed Nor Frioui, Toshikazu Kuniya, Tarik Mohammed Touaoula . Global stability of an age-structured epidemic model with general Lyapunov functional. Mathematical Biosciences and Engineering, 2019, 16(3): 1525-1553. doi: 10.3934/mbe.2019073 |
[9] | Edward J. Allen . Derivation and computation of discrete-delayand continuous-delay SDEs in mathematical biology. Mathematical Biosciences and Engineering, 2014, 11(3): 403-425. doi: 10.3934/mbe.2014.11.403 |
[10] | Fang Wang, Juping Zhang, Maoxing Liu . Dynamical analysis of a network-based SIR model with saturated incidence rate and nonlinear recovery rate: an edge-compartmental approach. Mathematical Biosciences and Engineering, 2024, 21(4): 5430-5445. doi: 10.3934/mbe.2024239 |
In this study, we explore a stochastic age-dependent cooperative Lotka-Volterra (LV) system with an environmental noise. By applying the theory of M-matrix, we prove the existence and uniqueness of the global solution for the system. Since the stochastic age-dependent cooperative LV system cannot be solved explicitly, we then construct an Euler-Maruyama (EM) numerical solution to approach the exact solution of the system. The convergence rate and the pth-moment boundedness of the scheme have also been obtained. Additionally, numerical experiments have been conducted to verify our theoretical results.
Mathematical models have been used to understand and control the dynamics of the disease (influenza, covid 19, HIV etc.). After first SIR type mathematical model constructed by Karmic and MC-Kendrick, for different kind of models (SIR, SIS, SEIR, SI etc.) have been proposed and studied by many authors (see, e.g., [1,2,3,4]). One of the important virus that whole world fighting with is the Human Immunodeficiency Virus (HIV) that is a lent virus caused HIV infection [5]. HIV virus can be transmitted in many ways like sexual intercourse, direct contact with contaminated blood products, needle, or during birth or through breastfeeding (mother to child). Since there is no recovery after HIV, SI type models are more appropriate to comprehend the dynamics of the HIV. In [3], J.J.Wang et al. have been studied mother to child transmission of HIV. The system that obtained with the constructed model in [3] extended and the numerical solutions of the system of linear parabolic equations(PEs) is studied in [6]. In the paper [7], Ashyralyev, Hincal and Kaymakamzade investigated the boundedness of solution of the initial boundary value problem for the system of PEs of observing epidemic models with general nonlinear incidence rate. This model constructed in [8] and it is well-known that (see, [9]) such and many other initial boundary value problems for system of PEs can be reduced to the initial-value problem for system of ordinary differential equations
{dv1(t)dt+μv1(t)+Av1(t)=−F(t,v1(t),v2(t)),dv2(t)dt+(α+μ)v2(t)+Av2(t)=F(t,v1(t),v2(t))−G(t,v2(t)),dv3(t)dt+μv3(t)+Av3(t)=G(t,v2(t)),0<t<T,vk(0)=ψk,1≤k≤3 | (1.1) |
in a Hilbert.space.H.with a self-adjoint.positive.definite.operator A≥δI,δ>0.
Existence.and.uniqueness.theorems.of the.bounded.solution of linear and nonlinear systems was.established in the following theorem ([7], [9]).
Theorem 1. Assume.the following.hypotheses hold
1. ψn,1≤n≤3 belongs to D(A) and
‖ψn‖D(A)=M1. | (1.2) |
2. The.function.F:[0,T]×H×H⟶H be.continuous.function, that.is
‖F(t,w(t),z(t))‖H≤M2 | (1.3) |
in [0,T]×H×H and Lipschitz condition holds uniformly with respect to t
‖F(t,u,v)−f(t,z,w)‖H≤L1[‖u−z‖H+‖v−w‖H]. | (1.4) |
3. The.function.G:[0,T]×H⟶H be.continuous function, that is
‖G(t,v(t))‖H≤M3 | (1.5) |
in [0,T]×H and Lipschitz condition holds uniformly with respect to t
‖G(t,u)−G(t,z)‖H≤L2‖u−z‖H. | (1.6) |
Then, there.exists.a unique bounded solution v(t)=(v1(t),v2(t),v3(t))⊥ of problem (1.1).
In applications, theorems on the bounded solutions of several systems of nonlinear. PEs were given. Moreover the first order of accuracy DS
{v1k−v1k−1τ+μv1k+Av1k=−F(tk,v1k,v2k),v2k−v2k−1τ+(α+μ)v2k+Av2k=Ff(tk,v1k,v2k)−G(tk,v2k),v3k−v3k−1τ+μv3k+Av3k=G(tk,v2k),tk=kτ,1⩽k⩽N,Nτ=T,vn0=ψn,1≤n≤3 | (1.7) |
for the approximate solution of problem (1.1) was studied. The existence and uniqueness of a bounded solution of DS (1.7) uniformly with.respect.to time.step.τ was.established in the following theorem.
Theorem 2. If the assumptions (1.2)-(1.6) and μ+δ>2(L1+L2) hold, then there exists a unique solution vτ={vk}Nk=0 of DS (1.7) which is bounded uniformly with respect to τ.
Bounded solutions of several systems of nonlinear PEs and DSs for the approximate solution of these systems were constituted. Numerical results were provided.
In general, .it is not possible to get exact solution of nonlinear problems. Therefore, we are interested in constructing a uniformly bounded high order of accuracy DSs with respect to time step size for the approximate solutions initial value problem (1.1).
In this work, for the approximate solution.of problem (1.1), the second.order of accuracy.Crank-Nicholson DS is investigated. The existence and uniqueness theorems of bounded solution of Crank-Nicholson DS uniformly with respect to time step τ is proved. In practice, theoretical results are presented on four nonlinear systems of parabolic equations to explain how it works on one and multidimensional problems. Numerical results are provided.
In this section, it will be considered the second order of accuracy Crank-Nicholson DS
{v1k−v1k−1τ+μv1k+v1k−12+Av1k+v1k−12=−F(tk−τ2,v1k+v1k−12,v2k+v2k−12),v2k−v2k−1τ+(α+μ)v2k+v2k−12+Av2k+v2k−12=F(tk−τ2,v1k+v1k−12,v2k+v2k−12)−G(tk−τ2,v2k+v2k−12),v3k−v3k−1τ+μv2k+v2k−12+Av2k+v2k−12=G(tk−τ2,v2k+v2k−12),tk=kτ,1⩽k⩽N,Nτ=T,vn0=ψn,1≤n≤3 | (2.1) |
for the approximate solution.of the initial.value.problem (1.1).
It is interested to studied the existence and uniqueness of a bounded solution of Crank- Nicholson DS (2.1) uniformly with respect to time step τ under the assumptions of Theorem 2.
Unfortunately, we are not able to establish the Theorem 3 for the solution of Crank-Nicholson DS (2.1) under the assumption μ+δ>2(L1+L2) without restriction to T. Nevertheless, it could be established the such result under assumption 1+1/2τ(μ+δ)>2(L1+L2)T. It is more strong than μ+δ>2(L1+L2) that means it is under assumption with restriction for T. Thus such result we can proved under assumption with restriction for T. It is based on reducing this DS to an equivalent system of nonlinear equations and is used as an operator method to prove the main theorem on the existence and uniqueness of a bounded solution of DS (2.1) uniformly with respect to τ.
An equivalent system of nonlinear equations for the DS (2.1) is
{v1k=Bkψ1−∑km=1Bk−mRF(tm−τ2,v1m+v1m−12,v2m+v2m−12)τ,v2k=Bk1ψ2+∑km=1Bk−m1R1[F(tm−τ2,v1m+v1m−12,v2m+v2m−12)−G(tm−τ2v2m+v2m−12)]τ,v3k=Bkψ3+∑km=1Bk−mRG(tm−τ2v2m+v2m−12)τ,1≤k≤N | (2.2) |
in Cτ(H)×Cτ(H)×Cτ(H) and the using of successive approximations. Here.and.in.future B=(I−τ(μI+A)2)R,R=(I+τ(μI+A)2)−1,B1=(I−τ((μ+α)I+A)2)R1,R1=(I+τ((μ+α)I+A)2)−1 and Cτ(H)=C([0,T]τ,H) stands for the Banach space of the mesh functions vτ={vl}Nl=0 defined on [0,T]τ with values in H,equipped with the norm
∥vτ∥Cτ(H)=max0≤l≤N‖vl‖H. |
For the solution of DS (2.1), the recursive formula is
{jv1k−jv1k−1τ+μjv1k+jv1k−12+Ajv1k+jv1k−12=−F(tk−τ2,(j−1)v1k+(j−1)v1k−12,(j−1)v2k+(j−1)v2k−12),jv2k−jv2k−1τ+(α+μ)jv2k+jv2k−12+Ajv2k+jv2k−12=F(tk−τ2,(j−1)v1k+(j−1)v1k−12,(j−1)v2k+(j−1)v2k−12)−G(tk−τ2,(j−1)v2k+(j−1)v2k−12),jv3k−jv3k−1τ+μjv2k+jv2k−12+Ajv2k+jv2k−12=G(tk−τ2,(j−1)v2k+(j−1)v2k−12),tk=kτ,1⩽k⩽N,Nτ=T,jvn0=ψn,1≤n≤3,j=1,2,...,0vnk=Bkψn,n=1,3,0v2k=Bkψ2,0≤k≤N. | (2.3) |
From (2.2) and (2.3) it follows
{jv1k=Bkψ1−∑km=1Bk−mRF(tk−τ2,(j−1)v1k+(j−1)v1k−12,(j−1)v2k+(j−1)v2k−12)τ,jv2k=Bk1ψ2+∑km=1Bk−m1R1F(tk−τ2,(j−1)v1k+(j−1)v1k−12,(j−1)v2k+(j−1)v2k−12)τ−∑km=1Bk−m1R1G(tk−τ2,(j−1)v2k+(j−1)v2k−12)τ,jv3k=Bkψ3+∑km=1Bk−mRG(tk−τ2,(j−1)v2k+(j−1)v2k−12),(j−1)v2k)τ,1≤k≤N,j=1,2,...,0vmk=Bkψm,m=1,3,0v2k=Bkψ2,0≤k≤N. | (2.4) |
Theorem 3. Let the assumptions (1.2)-(1.6) be satisfied and 2(L1+L2)T<1+τ(μ+δ)2. Then, there exists a unique solution vτ={vk}Nk=0 of DS (2.1) that is bounded in Cτ(H)×Cτ(H)×Cτ(H) of uniformly wrt. τ.
Proof. Since v3k does not appear in equations for vnk−vnk−1τ,n=1,2, it is sufficient to analyze the behaviors of solutions v1k and v2k of (2.1). According to the method of recursive approximation (2.4), we get
vnk=0vnk+∞∑i=0[(i+1)vnk−ivnk],n=1,2, | (2.5) |
where
0vnk={Bkψn,n=1,3,Bk1ψ2,n=2. | (2.6) |
Applying formula (2.6), estimates
‖B‖H→H≤1,‖B1‖H→H≤1, | (2.7) |
we get
‖0vnk‖H≤‖ψn‖H≤M1. | (2.8) |
Applying formula (2.4), estimates (2.7) and
‖R‖H→H≤11+τ(μ+δ)2,‖R1‖H→H≤11+τ(μ+δ+α)2, | (2.9) |
we get
‖1v1k−0v1k‖H≤k∑m=1‖Bk−mR‖H→H‖f(tm−τ2,0v1m+0v1m−12,0v2m+0v2m−12)‖Hτ |
≤M2k∑m=1τ1+τ(μ+δ)2≤M2T1+τ(μ+δ)2, |
‖1v2k−0v2k‖H≤k∑m=1‖Bk−m1R1‖H→H[‖F(tm−τ2,0v1m+0v1m−12,0v2m+0v2m−12)‖H |
+‖G(tm−τ2,0v2m+0v2m−12)‖H]τ |
≤(M2+M3)k∑m=1τ1+τ(μ+δ+α)2≤(M2+M3)T1+τ(μ+δ+α)2 |
for any k=1,⋅⋅⋅,N. With using triangle inequality, is is obtained that
‖1v1k‖H≤M1+(M2+M3)T1+τ(μ+δ)2, |
‖1v2k‖H≤M1+(M2+M3)T1+τ(μ+δ)2 |
for any.k=1,⋅⋅⋅,N. Applying.formula (2.4), and estimates (2.7), (2.9), (1.4), (1.2) and (1.3), we get
‖2v1k−1v1k‖H≤τk∑m=1‖Bk−mR‖H→H |
×‖F(tm−τ2,1v1m+1v1m−12,1v2m+1v2m−12)−F(tm−τ2,0v1m+0v1m−12,0v2m+0v2m−12)‖H |
≤k∑m=1L1τ1+τ(μ+δ)2[‖1v1m+1v1m−12−0v1m+0v1m−12‖H+‖1v2m+1v2m−12−0v2m+0v2m−12‖H] |
≤2L1(M2+M3)T1+τ(μ+δ)2k∑m=1τ1+τ(μ+δ)2≤2(L1+L2)(M2+M3)T2(1+τ(μ+δ)2)2, |
‖2v2k−1v2k‖H≤τk∑m=1‖Bk−m1R1‖H→H |
×‖F(tm−τ2,1v1m+1v1m−12,1v2m+1v2m−12)−F(tm−τ2,0v1m+0v1m−12,0v2m+0v2m−12)‖H |
+τk∑m=1‖Bk−m1R1‖H→H‖G(tm−τ2,1v2m+1v2m−12)−G(tm−τ2,0v2m+0v2m−12)‖H |
≤L1k∑m=1τ1+τ(μ+δ+α)2 |
×[‖1v1m+1v1m−12−0v1m+0v1m−12‖H+‖1v2m+1v2m−12−0v2m+0v2m−12‖H] |
+L2k∑m=1τ1+τ(μ+δ+α)2‖1v2m+1v2m−12−0v2m+0v2m−12‖H |
≤(2L1+L2)(M2+M3)T1+τ(μ+δ)2k∑m=1τ1+τ(μ+δ+α)2≤2(L1+L2)(M2+M3)T2(1+τ(μ+δ)2)2 |
for any k=1,⋅⋅⋅,N. Then
‖2vnk‖H≤M1+(M2+M3)T1+τ(μ+δ)2+2(L1+L2)(M2+M3)T2(1+τ(μ+δ)2)2,n=1,2 |
for any k=1,⋅⋅⋅,N. Let
‖jvnk−(j−1)vnk‖H≤2j−1(L1+L2)j−1(M2+M3)Tj(1+τ(μ+δ)2)j,n=1,2. |
Applying formula (2.4), estimates (2.7), (1.4), (1.2) and (1.3), we get
‖(j+1)v1k−jv1k‖H≤τk∑m=1‖Bk−mR‖H→H |
×‖f(tm−τ2,jv1m+jv1m−12,jv2m+jv2m−12) |
−f(tm−τ2,(j−1)v1m+(j−1)v1m−12,(j−1)v2m+(j−1)v2m−12)‖H |
≤k∑m=1L1τ1+τ(μ+δ)2 |
×[‖jv1m+jv1m−12−(j−1)v1m+(j−1)v1m−12‖H |
+‖jv2m+jv2m−12−(j−1)v2m+(j−1)v2m−12‖H] |
≤2L1⋅2j−1(L1+L2)j−1(M2+M3)Tj(1+τ(μ+δ)2)jk∑m=1τ1+τ(μ+δ)2≤(2(L1+L2))j(M2+M3)Tj+1(1+τ(μ+δ)2)j+1, |
‖(j+1)v2k−jv2k‖H≤τk∑m=1‖Bk−m1R1‖H→H |
×‖F(tm−τ2,jv1m+jv1m−12,jv2m+jv2m−12) |
−F(tm−τ2,(j−1)v1m+(j−1)v1m−12,(j−1)v2m+(j−1)v2m−12)‖H |
+τk∑m=1‖Bk−m1R1‖H→H |
×‖G(tm−τ2,jv2m+jv2m−12)−F(tm−τ2,(j−1)v2m+(j−1)v2m−12)‖H |
≤k∑m=1L1τ1+τ(μ+δ+a)2 |
×[‖jv1m+jv1m−12−(j−1)v1m+(j−1)v1m−12‖H |
+‖jv2m+jv2m−12−(j−1)v2m+(j−1)v2m−12‖H] |
+k∑m=1L2τ1+τ(μ+δ+a)2‖jv2m+jv2m−12−(j−1)v2m+(j−1)v2m−12‖H |
≤(2L1+L2)2j−1(L1+L2)j−1(M2+M3)Tj(1+τ(μ+δ)2)jk∑m=1τ1+τ(μ+δ)2≤(2(L1+L2))j(M2+M3)Tj+1(1+τ(μ+δ)2)j+1 |
for any k=1,⋅⋅⋅,N. Then
‖(j+1)vnk‖H≤M1+(M2+M3)T1+τ(μ+δ)2 |
+2(L1+L2)(M2+M3)T2(1+τ(μ+δ)2)2+⋅⋅⋅+(2(L1+L2))j(M2+M3)Tj+1(1+τ(μ+δ)2)j+1,n=1,2 |
for any k=1,⋅⋅⋅,N. Therefore, for any j,j≥1, we have that
‖(j+1)vnk−jvnk‖H≤(2(L1+L2))j(M2+M3)Tj+1(1+τ(μ+δ)2)j+1,n=1,2, |
and
‖(j+1)vnk‖H≤M1+(M2+M3)T1+τ(μ+δ)2 |
+2(L1+L2)(M2+M3)T2(1+τ(μ+δ)2)2+⋅⋅⋅+(2(L1+L2))j(M2+M3)Tj+1(1+τ(μ+δ)2)j+1,n=1,2 |
by mathematical.induction. From that and formula (2.5) it is obtained
‖vnk‖H≤‖0vnk‖H+∞∑i=0‖(i+1)vnk−ivnk‖H |
≤M1+(M2+M3)T1+τ(μ+δ)2∞∑i=02i(L1+L2)iTi(1+τ(μ+δ)2)i,n=1,2 |
that proves the existence of a bounded solution of DS (2.1) which is bounded.in Cτ(H)×Cτ(H)×Cτ(H) of uniformly wrt. τ. Theorem 3 is proved.
A study of discretization, over time only, of the initial value problem also permits one to include general DSs in applications, if the differential operator A is replaced by the difference operator Ah that act in the Hilbert spaces and are uniformly self-adjoint positive definite in h for 0<h≤h0.
In this section it will be given considered some nonlinear partial differential equations(PDEs).
First, we consider the initial-boundary.value problem for one.dimensional system.of nonlinear PDEs.
{∂v1(t,x)∂t−(a(x)v1x(t,x))x+(δ+μ)v1(t,x)=−F(t,x;v1(t,x),v2(t,x)),∂v2(t,x)∂t−(a(x)v2x(t,x))x+(δ+μ+α)v2(t,x)=F(t,x;v1(t,x),v2(t,x))−G(t,x;v2(t,x)),∂v3(t,x)∂t−(a(x)v3x(t,x))x+(δ+μ)v3(t,x)=G(t,x;v2(t,x)),0<t<T,0<x<l,vn(0,x)=ψn(x),ψn(0)=ψn(l),φmx(0)=ψnx(l),x∈[0,l],n=1,2,3,vn(t,0)=vn(t,l),vnx(t,0)=vnx(t,l),0≤t≤T,n=1,2,3, | (3.1) |
where a(x),ψ(x) are given sufficiently smooth functions and δ>0 is the sufficiently large number. We will assume that a(x)≥a>0 and a(l)=a(0).
Assume the following.hypotheses hold
1. ψn,n=1,2,3 belongs to W22[0,l] and
‖ψn‖W22[0,l]≤M1. | (3.2) |
2. The function f:[0,T]×[0,l]×L2[0,l]×L2[0,l]→L2[0,l] be continuous function in t, that is
‖F(t,⋅,u(t,⋅),v(t,⋅))‖L2[0,l]≤M2 | (3.3) |
in [0,T]×[0,l]×L2[0,l]×L2[0,l] and Lipschitz condition holds uniformly with respect to t
‖G(t,⋅,u,v)−G(t,⋅,z,w)‖L2[0,l]≤L1[‖u−z‖L2[0,l]+‖v−w‖L2[0,l]]. | (3.4) |
3. The function g:[0,T]×[0,l]×L2[0,l]→L2[0,l] be continuous function in t, that is
‖G(t,⋅,u(t,⋅))‖L2[0,l]≤M3 | (3.5) |
in [0,T]×[0,l]×L2[0,l] and Lipschitz condition holds uniformly with respect to t
‖G(t,⋅,u)−G(t,⋅,z)‖L2[0,l]≤L2‖u−z‖L2[0,l]. | (3.6) |
Here and in future, Lm,m=1,2,Mm,m=1,2,3 are positive constants.
The discretization of problem (3.1) is carried out in two steps. In the first step, let us define the grid space
[0,l]h={x:xr=rh,0≤r≤K,Kh=l}. |
We introduce the Hilbert spaces L2h=L2([0,l]h) and W22h=W22([0,l]h)\ of the grid functions ψh(x)={ψr}K0 defined on [0,l]h, equipped with the norms
‖ψh‖L2h=(∑x∈[0,l]h|ψh(x)|2h)1/2 |
and
‖ψh‖W22h=‖ψh‖L2h+(∑x∈[0,l]h|(ψh)x¯x,j|2h)1/2 |
respectively. To the differential operator A generated by problem (3.1), we assign the difference operator Axh by the formula
Axhψh(x)={−(a(x)ψ¯x)x,r+δψr}K−11, | (3.7) |
acting in the space of grid functions ψh(x)={ψr}K0 satisfying the conditions ψ0=ψK, ψ1−ψ0=ψK−ψK−1. With the help of Axh, we arrive at the initial value problem
{dv1h(t,x)dt+μv1h(t,x)+Axhv1h(t,x)=−Fh(t,x;v1h(t,x),v2h(t,x)),dv2h(t,x)dt+(μ+α)v2h(t,x)+Axhv2h(t,x)=Fh(t,x;v1h(t,x),v2h(t,x))−Gh(t,x;v2h(t,x)),dv3h(t,x)dt+μv3h(t,x)+Axhv3h(t,x)=Gh(t,x;v2h(t,x)),0<t<T,x∈[0,l]h,vnh(0,x)=ψn(x),n=1,2,3,x∈[0,l]h | (3.8) |
for an infinite system of nonlinear ordinary differential equations (ODEs). In the second step, we replace problem (3.8) by DS (2.1)
{v1k−v1k−1τ+μv1k+v1k−12+Axhv1k+v1k−12=−Fh(tk−τ2,x,v1k+v1k−12,v2k+v2k−12),v2k−v2k−1τ+(α+μ)v2k+v2k−12+Axhv2k+v2k−12=Fh(tk−τ2,x,v1k+v1k−12,v2k+v2k−12)−Gh(tk−τ2,x,v2k+v2k−12),v3k−v3k−1τ+μv3k+v3k−12+Axhv3k+v3k−12=Gh(tk−τ2,x,v2k+v2k−12),tk=kτ,1⩽k⩽N,Nτ=T,x∈[0,l]h,vn0=ψn,n=1,2,3. | (3.9) |
Theorem 4. Let the assumptions (3.2)-(3.6) be satisfied and 2(L1+L2)T<1+τ(μ+δ)2. Then, there exists a unique solution vτ={vk}Nk=0 of DS (3.9) which is bounded in Cτ(L2h)×Cτ(L2h)×Cτ(L2h) of uniformly with respect to τ and h.
The proof of Theorem 4 is based on the abstract Theorem 3 and symmetry properties of the difference operator Axh defined by formula (3.7)[11].
Second, we consider the initial-boundary value problem for one dimensional system of nonlinear PDs with involution
{∂v1(t,x)∂t−(a(x)v1x(t,x))x−β(a(−x)vx(t,−x))x+(δ+μ)v1(t,x)=−F(t,x;v1(t,x),v2(t,x)),∂v2(t,x)∂t−(a(x)v2x(t,x))x−β(a(−x)vx(t,−x))x+(δ+μ+α)v2(t,x)=F(t,x;v1(t,x),v2(t,x))−G(t,x;v2(t,x)),∂v3(t,x)∂t−(a(x)v3x(t,x))x−β(a(−x)vx(t,−x))x+(δ+μ)v3(t,x)=G(t,x;v2(t,x)),0<t<T,−l<x<l,vn(0,x)=ψn(x),ψn(−l)=ψn(l)=0,x∈[−l,l],n=1,2,3,vn(t,−l)=vn(t,l)=0,0≤t≤T,n=1,2,3, | (3.10) |
where a(x),ψ(x) are given sufficiently smooth functions and δ>0 is the sufficiently large number. We will assume that a≥a(x)=a(−x)≥δ>0,δ−a|β|≥0.
Assume the following hypotheses:
1. ψn,n=1,2,3 belongs to W22[−l,l] and
‖ψn‖W22[−l,l]≤M1. | (3.11) |
2. The function F:[0,T]×[−l,l]×L2[−l,l]×L2[−l,l]→L2[−l,l] be continuous function in t, that is
‖F(t,⋅,u(t,⋅),v(t,⋅))‖L2[−l,l]≤M2 | (3.12) |
in [0,T]×[−l,l]×L2[−l,l]×L2[−l,l] and Lipschitz condition holds uniformly with respect to t
‖F(t,⋅,u,v)−F(t,⋅,z,w)‖L2[−l,l]≤L1[‖u−z‖L2[−l,l]+‖v−w‖L2[−l,l]]. | (3.13) |
3. The function G:[0,T]×[−l,l]×L2[−l,l]→L2[−l,l] be continuous function in t, that is
‖G(t,⋅,u(t,⋅))‖L2[0,l]≤M3 | (3.14) |
in [0,T]×[−l,l]×L2[−l,l] and Lipschitz condition holds uniformly with respect to t
‖G(t,⋅,u)−G(t,⋅,z)‖L2[−l,l]≤L2‖u−z‖L2[−l,l]. | (3.15) |
The discretization of problem (3.10) is carried out in two steps. In the first step, let us define the grid space
[−l,l]h={x:xr=rh,−K≤r≤K,Kh=l}. |
We introduce the Hilbert spaces L2h=L2([−l,l]h) and W22h=W22([−l,l]h)\ of the grid functions ψh(x)={ψr}K−K defined on [−l,l]h, equipped with the norms
‖ψh‖L2h=(∑x∈[−l,l]h|ψh(x)|2h)1/2 |
and
‖ψh‖W22h=‖ψh‖L2h+(∑x∈[−l,l]h|(ψh)x¯x,j|2h)1/2 |
respectively. To the differential operator A generated by problem (3.10), we assign the difference operator Axh by the formula
Axhψh(x)={−(a(x)ψ¯x(x))x,r−β(a(−x)ψ¯x(−x))x,r+δψr}K−1−K+1, | (3.16) |
acting in the space of grid functions ψh(x)={ψr}K−K satisfying the conditions ψ−K=ψK=0. With the help of Axh, we arrive at the initial value problem
{dv1h(t,x)dt+μv1h(t,x)+Axhv1h(t,x)=−Fh(t,x;v1h(t,x),v2h(t,x)),dv2h(t,x)dt+(μ+α)v2h(t,x)+Axhv2h(t,x)=Fh(t,x;v1h(t,x),v2h(t,x))−Gh(t,x;v2h(t,x)),dv3h(t,x)dt+μv3h(t,x)+Axhv3h(t,x)=Gh(t,x;v2h(t,x)),0<t<T,x∈[−l,l]h,vnh(0,x)=ψn(x),n=1,2,3,x∈[−l,l]h | (3.17) |
for an infinite system of nonlinear ODEs. In the second step, we replace problem (3.17) by DS (2.1)
{v1k−v1k−1τ+μv1k+v1k−12+Axhv1k+v1k−12=−Fh(tk−τ2,x,v1k+v1k−12,v2k+v2k−12),v2k−v2k−1τ+(α+μ)v2k+v2k−12+Axhv2k+v2k−12=Fh(tk−τ2,x,v1k+v1k−12,v2k+v2k−12)−Gh(tk−τ2,x,v2k+v2k−12),v3k−v3k−1τ+μv3k+v3k−12+Axhv3k+v3k−12=Gh(tk−τ2,x,v2k+v2k−12),tk=kτ,1⩽k⩽N,Nτ=T,x∈[−l,l]h,vn0=ψn,n=1,2,3. | (3.18) |
Theorem 5. Let the assumptions (3.11)-(3.15) be satisfied and 2(L1+L2)T<1+τ(μ+δ)2. Then, there exists a unique solution vτ={vk}Nk=0 of DS (3.18) which is bounded in Cτ(L2h)×Cτ(L2h)×Cτ(L2h) of uniformly with respect to τ and h.
The proof of Theorem 5 is based on the abstract Theorem 3 and symmetry properties of the difference operator Axh defined by formula (3.16)[12].
Third, let Ω⊂Rn be a bounded open domain with smooth boundary S, ¯Ω=Ω∪S. In [0,T]×Ω we consider the initial-boundary value problem for multidimensional system of nonlinear PDEs
{∂v1(t,x)∂t−n∑r=1(ar(x)v1xr)xr+(δ+μ)v1(t,x)=−F(t,x;v1(t,x),v2(t,x)),∂v2(t,x)∂t−n∑r=1(ar(x)v2xr)xr+(δ+μ+α)v2(t,x)=F(t,x;v1(t,x),v2(t,x))−G(t,x;v2(t,x)),∂v3(t,x)∂t−n∑r=1(ar(x)v3xr)xr+(δ+μ)v3(t,x)=G(t,x;v2(t,x)),0<t<T,x=(x1,...,xn)∈Ω,vm(0,x)=ψm(x),x∈¯Ω,m=1,2,3,vm(t,x)=0,0≤t≤T,x∈S,m=1,2,3, | (3.19) |
where ar(x) and ψm(x) are given sufficiently smooth functions and δ>0 is the sufficiently large number and ar(x)>0.
Assume the following hypotheses:
1. ψm,m=1,2,3 belongs to L2(¯Ω) and
‖ψm‖W22(¯Ω)≤M1. | (3.20) |
2. The function f:[0,T]×[0,l]×L2(¯Ω)×L2(¯Ω)→L2(¯Ω) be continuous function in t, that is
‖F(t,⋅,u(t,⋅),v(t,⋅))‖L2(¯Ω)≤M2 | (3.21) |
in [0,T]×[0,l]×L2(¯Ω)×L2(¯Ω) and Lipschitz condition holds uniformly with respect to t
‖F(t,⋅,u,v)−F(t,⋅,z,w)‖L2(¯Ω)≤L1[‖u−z‖L2(¯Ω)+‖v−w‖L2(¯Ω)]. | (3.22) |
3. The function G:[0,T]×[0,l]×L2(¯Ω)→L2(¯Ω) be continuous function in t, that is
‖G(t,⋅,u(t,⋅))‖L2(¯Ω)≤M3 | (3.23) |
in [0,T]×[0,l]×L2(¯Ω) and Lipschitz condition holds uniformly with respect to t
‖G(t,⋅,u)−G(t,⋅,z)‖L2(¯Ω)≤L2‖u−z‖L2(¯Ω). | (3.24) |
The discretization of problem (3.19) is also carried out in two steps. In the first step, let us define the grid sets
¯Ωh={x=xr=(h1r1,...,hmrm),r=(r1,...,rm), |
0≤rj≤Nj,hjNj=1,j=1,...,m}, |
Ωh=¯Ωh∩Ω,Sh=¯Ωh∩S. |
We introduce the Banach spaces L2h=L2(¯Ωh) and W22h=W22(¯Ωh)\ of the grid functions φh(x)={ψ(h1r1,...,hmrm)} defined on ¯Ωh, equipped with the norms
‖ψh‖L2h=(∑x∈¯Ωh|ψh(x)|2h1⋅⋅⋅hm)1/2 |
and
‖ψh‖W2h=‖ψh‖L2h+(∑x∈¯Ωhm∑r=1|(ψh)xr¯xr,jr|2h1⋅⋅⋅hm)1/2 |
respectively. To the differential operator A generated by problem (3.19), we assign the difference operator Axh by the formula
Axhvhx=−m∑r=1(ar(x)vh¯xr)xr,jr | (3.25) |
acting in the space of grid functions uh(x), satisfying the conditions vh(x)=0 for all x∈Sh. It is known that Axh is a self-adjoint positive definite operator in L2(¯Ωh). With the help of Axh, we arrive at the initial value problem
{dv1h(t,x)dt+μv1h(t,x)+Axhv1h(t,x)=−Fh(t,x;v1h(t,x),v2h(t,x)),dv2h(t,x)dt+(μ+α)v2h(t,x)+Axhv2h(t,x)=Fh(t,x;v1h(t,x),v2h(t,x))−Gh(t,x;v2h(t,x)),dv3h(t,x)dt+μv3h(t,x)+Axhv3h(t,x)=Gh(t,x;v2h(t,x)),0<t<T,x∈¯Ωh,vmh(0,x)=ψm(x),m=1,2,3,x∈¯Ωh | (3.26) |
for an infinite system of nonlinear ODEs. In the second step, we replace problem (3.26) by DS (2.1)
{v1k−v1k−1τ+μv1k+v1k−12+Axhv1k+v1k−12=−Fh(tk−τ2,x,v1k+v1k−12,v2k+v2k−12),v2k−v2k−1τ+(α+μ)v2k+v2k−12+Axhv2k+v2k−12=Fh(tk−τ2,x,v1k+v1k−12,v2k+v2k−12)−Gh(tk−τ2,x,v2k+v2k−12),v3k−v3k−1τ+μV3k+v3k−12+Axhv3k+v3k−12=Gh(tk−τ2,x,v2k+v2k−12),tk=kτ,1⩽k⩽N,Nτ=T,x∈¯Ωh,vm0=ψm,m=1,2,3. | (3.27) |
Theorem 6. Let the assumptions (3.20)-(3.24) be satisfied and 2(L1+L2)T<1+τ(μ+δ)2. Then, there exists a unique solution vτ={vk}Nk=0 of DS (3.27) which is bounded in Cτ(L2h)×Cτ(L2h)×Cτ(L2h) of uniformly with respect to τ and h.
The proof of Theorem 6 is based on the abstract Theorem 4 and symmetry properties of the difference operator Axh defined by formula (3.25) and the following theorem on coercivity inequality for the solution of the elliptic problem in L2h.
Theorem 7. For the solutions of the elliptic difference problem
{Axhvh(x)=gh(x), x∈Ωh,vh(x)=0, x∈Sh |
the following coercivity inequality
m∑r=1‖vhxr¯xr,jr‖L2h≤M‖gh‖L2h. |
holds (see [13]).
Fourth, in [0,T]×Ω we consider the initial-boundary value problem for multidimensional system of nonlinear PDEs
{∂v1(t,x)∂t−n∑r=1(ar(x)v1xr)xr+(δ+μ)v1(t,x)=−F(t,x;v1(t,x),v2(t,x)),∂v2(t,x)∂t−n∑r=1(ar(x)v2xr)xr+(δ+μ+α)v2(t,x)=f(t,x;v1(t,x),v2(t,x))−G(t,x;v2(t,x)),∂v3(t,x)∂t−n∑r=1(ar(x)v3xr)xr+(δ+μ)v3(t,x)=G(t,x;v2(t,x)),0<t<T,x=(x1,...,xn)∈Ω,vm(0,x)=ψm(x),x∈¯Ω,m=1,2,3,∂v∂→p(t,x)=0,0≤t≤T,x∈S,m=1,2,3, | (3.28) |
where ar(x) and ψm(x) are given sufficiently smooth functions and δ>0 is the sufficiently large number and ar(x)>0. Here, →p is the normal vector to Ω.
The discretization of problem (3.28) is also carried out in two steps. In the first step, to the differential operator A generated by problem (3.28), we assign the difference operator Axh by the formula
Axhvhx=−m∑r=1(ar(x)vh¯xr)xr,jr+δvh(x) | (3.29) |
acting in the space of grid functions vh(x), satisfying the conditions Dhvh(x)=0 for all x∈Sh. Here Dh is the approximation of operator ∂∂→p. It is known that Axh is a self-adjoint positive definite operator in L2(¯Ωh). With the help of Axh, we arrive at the initial value problem (3.26) for an infinite system of nonlinear ODEs. In the second step, we replace problem (3.26) by DS (2.1), we get DS (3.27).
Theorem 8. Let the assumptions (3.20)-(3.24) be satisfied and 2(L1+L2)T<1+τ(μ+δ)2. Then, there exists a unique solution vτ={vk}Nk=0 of DS (3.27) which is bounded in Cτ(L2h)×Cτ(L2h)×Cτ(L2h) of uniformly with respect to τ and h.
The proof of Theorem 8 is based on the abstract Theorem 3 and symmetry properties of the difference operator Axh defined by formula (3.29) and the following theorem on coercivity inequality for the solution of the elliptic problem in L2h.
Theorem 9. For the solutions of the elliptic difference problem
{Axhvh(x)=gh(x), x∈Ωh,Dhvh(x)=0, x∈Sh |
the following coercivity inequality holds (see [13]):
m∑r=1‖vhxr¯xr,jr‖L2h≤M‖gh‖L2h. |
In present section, we consider the initial-boundary value problem
{∂v1(t,x)∂t−λ+μv1(t,x)−β∂2v1(t,x)∂x2=−λ+(−1+μ+β)e−tsinx−sin(v1(t,x)v2(t,x)−e−2tsin2x),∂v2(t,x)∂t+(μ+α)v2(t,x)−d∂2v2(t,x)∂x2=(−1+μ+α+d)e−tsinx+sin(v1(t,x)v2(t,x)−e−2tsin2x)−cos(v2(t,x)−e−tsinx),∂v3(t,x)∂t+μv1(t,x)−γ∂2v1(t,x)∂x2=(−1+μ+γ)e−tsinx+cos(v2(t,x)−e−tsinx),0<t<1,0<x<π,vm(0,x)=sin(x),0≤x≤π,m=1,2,3,vm(t,0)=um(t,π)=0,0≤t≤1,m=1,2,3 | (4.1) |
for the system of nonlinear PDEs. The spatial factor, x, can be spatially discrete or spatially continuous. In either case, the spatial factor is used to describe the mobility of the population. This mobility can be due to travel and migration, and it could be between cities, towns or even countries, depending on the studied case. The exact solution of problem (4.1) is vm(t,x)=e−tsinx,m=1,2,3.
Numerical solutions of system (4.1) will be given for first and second order of DS. Firstly, we consider the first order of accuracy iterative DS
{jv1,kn−jv1,k−1nτ+μjv1,kn−βjv1,kn+1−2(jv1,kn)+jv1,kn−1h2=(−1+μ+β)e−tksinxn−sin((j−1)v1,kn(j−1)v2,kn−e−2tksin2xn),jv2,kn−jv2,k−1nτ+(α+μ)jv2,kn−djv2,kn+1−2(jv2,kn)+jv2,kn−1h2=(−1+μ+α+d)e−tksinxn+sin((j−1)v1,kn(j−1)v2,kn−e−2tksin2xn)−cos((j−1)v2,kn−e−tksinxn),jv3,kn−jv3,k−1nτ+μjv3,kn−γjv3,kn+1−2(jv3,kn)+jv3,kn−1h2=(−1+μ+γ)e−tksinxn+cos((j−1)v2,kn−e−tksinxn),tk=kτ,1≤k≤N,Nτ=1,xn=nh,1≤n≤M−1,Mh=π,jvm,0n=ψm(xn),jvm,k0=jum,kM=0,0≤k≤N,m=1,2,3,j=1,2,⋅⋅⋅,0vm,kn,0≤k≤N,0≤n≤M,m=1,2,3 is given | (4.2) |
and secondly, the second order of accuracy iterative Crank-Nicholson DS
\begin{equation} \left\{ \begin{array}{l} \frac{jv_{n}^{1,k}-jv_{n}^{1,k-1}}{\tau }+\mu \frac{ jv_{n}^{^{_{1,k}}}+jv_{n}^{^{_{1,k-1}}}}{2}-\beta \frac{_{j}v_{n+1}^{1,k}-2 \left( jv_{n}^{1,k}\right) +jv_{n-1}^{1,k}}{2h^{2}}-\beta \frac{ _{j}v_{n+1}^{1,k-1}-2\left( jv_{n}^{1,k-1}\right) +jv_{n-1}^{1,k-1}}{2h^{2}} \\ \\ = \left( -1+\mu +\beta \right) e^{-\left( t_{k}-\frac{\tau }{2}\right) }\sin x_{n} \\ \\ -\sin \left( \frac{\left( j-1\right) v_{n}^{^{_{1,k}}}+\left( j-1\right) v_{n}^{^{_{1,k-1}}}}{2}\frac{\left( j-1\right) v_{n}^{^{_{2,k}}}+\left( j-1\right) v_{n}^{^{_{2,k-1}}}}{2}-e^{-2\left( t_{k}-\frac{\tau }{2}\right) }\sin ^{2}x_{n}\right) , \\ \\ \frac{jv_{n}^{2,k}-jv_{n}^{2,k-1}}{\tau }+\left( \alpha +\mu \right) \frac{ jv_{n}^{^{_{2,k}}}+jv_{n}^{^{_{2,k-1}}}}{2}-d\frac{_{j}v_{n+1}^{2,k}-2\left( jv_{n}^{2,k}\right) +jv_{n-1}^{2,k}}{2h^{2}}-d\frac{_{j}v_{n+1}^{2,k-1}-2 \left( jv_{n}^{2,k-1}\right) +jv_{n-1}^{2,k-1}}{2h^{2}} \\ \\ = \left( -1+\mu +\alpha +d\right) e^{-\left( t_{k}-\frac{\tau }{2}\right) }\sin x_{n} \\ \\ +\sin \left( \frac{\left( j-1\right) v_{n}^{^{_{1,k}}}+\left( j-1\right) v_{n}^{^{_{1,k-1}}}}{2}\frac{\left( j-1\right) v_{n}^{^{_{2,k}}}+\left( j-1\right) v_{n}^{^{_{2,k-1}}}}{2}-e^{-2\left( t_{k}-\frac{\tau }{2}\right) }\sin ^{2}x_{n}\right) \\ \\ -\cos \left( \frac{\left( j-1\right) v_{n}^{^{_{2,k}}}+\left( j-1\right) v_{n}^{^{_{2,k-1}}}}{2}-e^{-\left( t_{k}-\frac{\tau }{2}\right) }\sin x_{n}\right) , \\ \\ \frac{jv_{n}^{3,k}-jv_{n}^{3,k-1}}{\tau }+\mu \frac{ jv_{n}^{^{_{3,k}}}+jv_{n}^{^{_{3,k-1}}}}{2}-\gamma \frac{_{j}v_{n+1}^{3,k}-2 \left( jv_{n}^{3,k}\right) +jv_{n-1}^{3,k}}{2h^{2}}-\gamma \frac{ _{j}v_{n+1}^{3,k-1}-2\left( jv_{n}^{3,k-1}\right) +jv_{n-1}^{3,k-1}}{2h^{2}} \\ \\ = \left( -1+\mu +\gamma \right) e^{-\left( t_{k}-\frac{\tau }{2}\right) }\sin x_{n}+\cos \left( \frac{\left( j-1\right) v_{n}^{^{_{2,k}}}+\left( j-1\right) v_{n}^{^{_{2,k-1}}}}{2}-e^{-\left( t_{k}-\frac{\tau }{2}\right) }\sin x_{n}\right) , \\ \\ t_{k} = k\tau ,1\leq k\leq N,N\tau = 1,x_{n} = nh,1\leq n\leq M-1,Mh = \pi , \\ \\ jv_{n}^{m,0} = \psi ^{m}(x_{n}),jv_{0}^{m,k} = jv_{M}^{m,k} = 0,0\leq k\leq N,m = 1,2,3,j = 1,2,\cdot \cdot \cdot , \\ \\ 0v_{n}^{m,k},0\leq k\leq N,0\leq n\leq M,m = 1,2,3\text{ is given} \end{array} \right. \end{equation} | (4.3) |
for the.approximate solution of the.initial-boundary.value.problem (4.1) for the system.of nonlinear. PEs. Here and in future j denotes the iteration.index and an.initial guess. _{0}u_{n}^{k}, k\geq 1, 0\leq n\leq M is to be made. For solving DS (4.3), the numerical.steps are.given.below. For 0\leq k < N, 0\leq n\leq M the algorithm.is as follows.[10] :
1. j = 1 .
2. _{j-1}v_{n}^{k} is.known.
3. _{j}v_{n}^{k} is.calculated.
4. If the.max absolute.error between _{j-1}v_{n}^{k} and _{j}v_{n}^{k} is greater.than the given.tolerance value \varepsilon = 10^{-8} , take j = j+1 and go.to.step 2. Otherwise, terminate.the iteration.process and take _{j}v_{n}^{k} as the result of the given problem. The errors are computed by
\left( jE^{m}\right) _{M}^{N} = \max\limits_{1\leq k\leq N, 1\leq n\leq M-1}\left\vert v^{m}(t_{k},x_{n})-\left( jv^{m}\right) _{n}^{k}\right\vert, m = 1,2,3 | (4.4) |
of the numerical solutions, where v^{m}(t_{k}, x_{n}) \;, m = 1, 2, 3 represents the exact solutions and \left(jv^{m}\right) _{n}^{k} \;, m = 1, 2, 3 represents.the numerical.solutions at (t_{k}, x_{n}) and the results of the first and second order of DS are given.in Table 1 and Table 2 respectively.
\left(jE^{m}\right) _{M}^{N} | N=M=20 | N=M=40 | N=M=80 |
m=1 | 0.0068, j=6 | 0.0032, j=6 | 0.0016, j=6 |
m=2 | 0.0071, j=6 | 0.0033, j=6 | 0.0016, j=6 |
m=3 | 0.0073, j=6 | 0.0034, j=6 | 0.0017, j=6 |
\left(jE^{m}\right) _{M}^{N} | N=M=20 | N=M=40 | N=M=80 |
m=1 | 5.5516e-5, j=7 | 1.3882e-5, j=7 | 3.4708e-6, j=7 |
m=2 | 8.7420e-5, j=7 | 2.1857e-5, j=7 | 5.4645e-6, j=7 |
m=3 | 1.1120e-4, j=7 | 2.7803e-5, j=7 | , 6.9510e-6j=7 |
According to Table 1 and Table 2, if N and M are doubled, the value of errors in the first order of accuracy DS decrease by a factor of 1/2, the errors in the second order of accuracy DS (4.3) decrease approximately by a factor of 1/4 . The errors presented in the tables indicate the stability of the DS and the accuracy of the results. Thus, the second order of accuracy DS increases faster than the first order of accuracy DS.
In the present paper, the.initial boundary value problem for the nonlinear system of PEs observing epidemic models with general nonlinear incidence rate is investigated. The main theorem on the existence and uniqueness of a bounded solution of Crank-Nicholson DS uniformly with respect to time step \tau is established. Applications of the theoretical results are presented for the four systems of one and multidimensional problems with different boundary conditions. Numerical results are given.
The publication has been prepared with the support of the "RUDN University Program 5-100" and published under target program BR05236656 of the Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan.
The author declares that there are no conflicts of interest.
[1] |
J. H. Vandermeer, D. H. Boucher, Varieties of mutualistic interaction in population models, J. Theor. Biol., 74 (1978), 549–558. doi: 10.1016/0022-5193(78)90241-2
![]() |
[2] |
W. J. Zuo, D. Q. Jiang, X. G. Sun, T. Hayat, A. Alsaedi, Long-time behaviors of a stochastic cooperative Lotka-Volterra system with distributed delay, Physica A., 506 (2018), 542–559. doi: 10.1016/j.physa.2018.03.071
![]() |
[3] | Q. Liu, Analysis of a stochastic non-autonomous food-limited Lotka-Volterra cooperative model, Appl. Math. Comput., 254 (2015), 1–8. |
[4] | E. P. Odum, Fundamentals of Ecology 5th Edition, Philadelphia, Saunders, 1971. |
[5] | J. D. Murray, Mathematical Biology, Springer, Berlin, 1989. |
[6] | M. E. Hale, The Biology of Lichens, Arnold, London, 1974. |
[7] |
J. Roughgarden, Evolution of marine symbiosis-A simple cost-benefit model, Ecology, 56 (1975), 1201–1208. doi: 10.2307/1936160
![]() |
[8] |
J. H. Huang, X. F. Zou, Traveling wavefronts in diffusive and cooperative Lotka-Volterra system with delays, J. Math. Anal. Appl., 271 (2002), 455–466. doi: 10.1016/S0022-247X(02)00135-X
![]() |
[9] |
G. C. Lu, Z. Y. Lu, Permanence for two-species Lotka-Volterra cooperative systems with delays, Math. Biosci. Eng., 5 (2008), 477–484. doi: 10.3934/mbe.2008.5.477
![]() |
[10] | G. C. Lu, Z. Y. Lu, X. Z. Lian, Delay effect on the permanence for Lotka-Volterra cooperative systems, Nonlinear. Anal-Real., 11 (2010), 2810–2816. |
[11] | Y. Nakata, Y. Muroya, Permanence for nonautonomous Lotka-Volterra cooperative systems with delays, Nonlinear. Anal-Real., 11 (2010), 528–534. |
[12] | L. Xu, J. Y. Liu, G. Zhang, Pattern formation and parameter inversion for a discrete Lotka-Volterra cooperative system, Chaos. Soliton. Fract., 110 (2018), 226–231. |
[13] | W. T. Li, Z. C. Wang, Traveling fronts in diffusive and cooperative Lotka-Volterra system with nonlocal delays, Z. Angew. Math. Phys., 58 (2007), 571–591. |
[14] | A. J. Lotka, Relation between birth rates and death rates, Science, 26 (1907), 121–130. |
[15] | X. R. Mao, Stochastic Differential Equations and Applications, Horwood, UK, 2007. |
[16] | X. M. Zhang, Z. H. Liu, Periodic oscillations in age-structured ratio-dependent predator-prey model with Michaelis-Menten type functional response, Physica D., 389 (2019), 51–63. |
[17] | F. J. Solis, R. A. Ku-Carrillo, Generic predation in age structure predator-prey models, Appl. Math. Comput., 231 (2014), 205–213. |
[18] |
M. Delgado, A. Suárez, Age-dependent diffusive Lotka-Volterra type systems, Math. Comput. Model, 45 (2007), 668–680. doi: 10.1016/j.mcm.2006.07.013
![]() |
[19] | N. Hritonenko, Y. Yatsenko, The structure of optimal time and age-dependent harvesting in the Lotka-McKendrik population model, Math. Biosci., 208 (2007), 48–62. |
[20] | D. W. Tudor, An age-dependent epidemic model with application to measles, Math. Biosci., 73 (1985), 131–147. |
[21] | D. Greenhalgh, Threshold and stability results for an epidemic model with an age-structured meeting rate, Math. Med. Biol., 5 (1988), 81–100. |
[22] | Q. M. Zhang, W. A. Liu, Z. K. Nie, Existence, uniqueness and exponential stability for stochastic age-dependent population, Appl. Math. Comput., 154 (2004), 183–201. |
[23] | Y. Li, M. Ye, Q. M. Zhang, Strong convergence of the partially truncated Euler-Maruyama scheme for a stochastic age-structured SIR epidemic model, Appl. Math. Comput., 362 (2019), 1–22. |
[24] | W. K. Pang, R. H. Li, M. Liu, Exponential stability of numerical solutions to stochastic age-dependent population equations, Appl. Math. Comput., 183 (2006), 152–159. |
[25] | R. H. Li, H. B. Meng, Q. Chang, Convergence of numerical solutions to stochastic age-dependent population equations, J. Comput. Appl. Math., 193 (2006), 109–120. |
[26] | R. H. Li, P. K. Leung, W. K. Pang, Convergence of numerical solutions to stochastic age-dependent population equations with Markovian switching, J. Comput. Appl. Math., 233 (2009), 1046–1055. |
[27] | R. H. Li, W. K. Pang, Q. H. Wang, Numerical analysis for stochastic age-dependent population equations with Poisson jumps, J. Math. Anal. Appl., 327 (2007), 1214–1224. |
[28] | W. K. Pang, R. H. Li, M. Liu, Convergence of the semi-implicit Euler method for stochastic age-dependent population equations, Appl. Math. Comput., 195 (2008), 466–474. |
[29] | X. Y. Li, Variational iteration method for nonlinear age-structured population models, Comput. Math. Appl., 58 (2009), 2177–2181. |
[30] | Y. Yang, C. F. Wu, Z. X. Li, Forced waves and their asymptotics in a Lotka-Volterra cooperative model under climate change, Appl. Math. Comput., 353 (2019), 254–264. |
[31] | Y. Zhao, S. L. Yuan, Q. M. Zhang, The effect of noise on the survival of a stochastic competitive model in an impulsive polluted environment, Appl. Math. Model., 40 (2016), 7583–7600. |
[32] | C. Q. Xu, S. L. Yuan, T. H. Zhang, Global dynamics of a predator-prey model with defence mechanism for prey, Appl. Math. Lett., 62 (2016), 42–48. |
[33] | Q. M. Zhang, Exponential stability of numerical solutions to a stochastic age-structured population system with diffusion, J. Comput. Appl. Math., 220 (2008), 22–33. |
[34] |
A. Chekroun, T. Kuniya, Global threshold dynamics of an infection age-structured SIR epidemic model with diffusion under the Dirichlet boundary condition, J. Diff. Equations, 269 (2020), 117–148. doi: 10.1016/j.jde.2020.04.046
![]() |
[35] | R. X. Lu, F. Y. Wei, Persistence and extinction for an age-structured stochastic SVIR epidemic model with generalized nonlinear incidence rate, Physica A., 513 (2019), 572–587. |
[36] | P. Yang, Y. S. Wang, Existence and properties of Hopf bifurcation in an age-dependent predation system with prey harvesting, Commun. Nonlinear. Sci. Numer. Simulat., 91 (2020), 105395. |
[37] | Z. X. Luo, Z. R. He, Optimal control for age-dependent population hybrid system in a polluted environment, Appl. Math. Comput., 228 (2014), 68–76. |
[38] | C. Burgos, J. C. Cortés, L. Shaikhet, R. J. Villanueva, A nonlinear dynamic age-structured model of e-commerce in spain: Stability analysis of the equilibrium by delay and stochastic perturbations, Commun. Nonlinear. Sci. Numer. Simulat., 64 (2018), 149–158. |
1. | Nezihal Gokbulut, Evren Hincal, Hasan Besim, Bilgen Kaymakamzade, Reducing the Range of Cancer Risk on BI-RADS 4 Subcategories via Mathematical Modelling, 2022, 133, 1526-1506, 93, 10.32604/cmes.2022.019782 |
\left(jE^{m}\right) _{M}^{N} | N=M=20 | N=M=40 | N=M=80 |
m=1 | 0.0068, j=6 | 0.0032, j=6 | 0.0016, j=6 |
m=2 | 0.0071, j=6 | 0.0033, j=6 | 0.0016, j=6 |
m=3 | 0.0073, j=6 | 0.0034, j=6 | 0.0017, j=6 |
\left(jE^{m}\right) _{M}^{N} | N=M=20 | N=M=40 | N=M=80 |
m=1 | 5.5516e-5, j=7 | 1.3882e-5, j=7 | 3.4708e-6, j=7 |
m=2 | 8.7420e-5, j=7 | 2.1857e-5, j=7 | 5.4645e-6, j=7 |
m=3 | 1.1120e-4, j=7 | 2.7803e-5, j=7 | , 6.9510e-6j=7 |
\left(jE^{m}\right) _{M}^{N} | N=M=20 | N=M=40 | N=M=80 |
m=1 | 0.0068, j=6 | 0.0032, j=6 | 0.0016, j=6 |
m=2 | 0.0071, j=6 | 0.0033, j=6 | 0.0016, j=6 |
m=3 | 0.0073, j=6 | 0.0034, j=6 | 0.0017, j=6 |
\left(jE^{m}\right) _{M}^{N} | N=M=20 | N=M=40 | N=M=80 |
m=1 | 5.5516e-5, j=7 | 1.3882e-5, j=7 | 3.4708e-6, j=7 |
m=2 | 8.7420e-5, j=7 | 2.1857e-5, j=7 | 5.4645e-6, j=7 |
m=3 | 1.1120e-4, j=7 | 2.7803e-5, j=7 | , 6.9510e-6j=7 |