Home Dynamics of two-species delayed competitive stage-structured model described by differential-difference equations
Article Open Access

Dynamics of two-species delayed competitive stage-structured model described by differential-difference equations

  • Sufang Han , Yaqin Li , Guoxin Liu , Lianglin Xiong and Tianwei Zhang
Published/Copyright: May 16, 2019

Abstract

Overf the last few years, by utilizing Mawhin’s continuation theorem of coincidence degree theory and Lyapunov functional, many scholars have been concerned with the global asymptotical stability of positive periodic solutions for the non-linear ecosystems. In the real world, almost periodicity is usually more realistic and more general than periodicity, but there are scarcely any papers concerning the issue of the global asymptotical stability of positive almost periodic solutions of non-linear ecosystems. In this paper we consider a kind of delayed two-species competitive model with stage structure. By means of Mawhin’s continuation theorem of coincidence degree theory, some sufficient conditions are obtained for the existence of at least one positive almost periodic solutions for the above model with nonnegative coefficients. Furthermore, the global asymptotical stability of positive almost periodic solution of the model is also studied. The work of this paper extends and improves some results in recent years. An example and simulations are employed to illustrate the main results of this paper.

1 Introduction

In [1], Zeng et al. proposed the following two-species competitive model with stage structure:

x˙1(t)=a1(t)x1(t)+b1(t)x2(t),x˙2(t)=a2(t)x2(t)b2(t)x2(t)c(t)x22(t)β1(t)x2(t)x3(t),x˙3(t)=x3(t)a3(t)b3(t)x3(t)β2(t)x2(t), (1.1)

where x1 and x2 are immature and mature population densities of one species, respectively, and x3 represents the population density of another species. The competition is between x2 and x3. By means of the fixed point theory and Lyapunv functional, Zeng et al. studied the existence and uniqueness of globally attractive positive T-periodic solution of system (1.1).

In real world, almost periodicity is more realistic and more general than periodicity. Therefore, more and more scholars have focused on the study of almost periodic dynamics of non-linear ecosystem [2, 3, 4, 5, 6]. Moreover, population models with delays have attracted much attention in recent years. Time delays represent an additional level of complexity that can be incorporated in a more detailed analysis of a particular system. So, this article is to consider the following delayed two-species almost periodic competitive model with stage structure:

x˙1(t)=a1(t)x1(t)+b1(t)x2(tτ(t)),x˙2(t)=a2(t)x2(t)b2(t)x2(t)c(t)x22(t)β1(t)x2(t)x3(tδ(t)),x˙3(t)=x3(t)a3(t)b3(t)x3(t)β2(t)x2(tσ(t)), (1.2)

where ai, bi, βj, c, τ, δ and σ are all nonnegative almost periodic functions, i = 1,2,3, j = 1,2.

Let ℝ, ℤ and ℕ+ denote the sets of real numbers, integers and positive integers, respectively, C(𝕏,𝕐) and C1(𝕏,𝕐) be the space of continuous functions and continuously differential functions which map 𝕏 into 𝕐, respectively. Especially, C(𝕏): = C(𝕏,𝕏), C1(𝕏): = C1(𝕏,𝕏). In relation to a continuous bounded function f, we use the following notations:

f=infsRf(s),f+=supsRf(s),|f|=supsR|f(s)|,f¯=limT1T0Tf(s)ds.

Throughout this paper, we always make the following assumption for system (1.2):

  1. ā1 > 0, ā2 > 2, ā3 > 0, 1 > 0 and 3 > 0.

fThe main purpose of this article is to study the existence and global asymptotic stability of positive almost periodic solution of system (1.2) by using the coincidence degree theory and Lyapunov functional. Finally, a example and some simulations are also given to illustrate the main results.

2 Preliminaries

Definition 2.1

([7, 8]) xC(ℝ, ℝn) is called almost periodic, if for any ϵ > 0, it is possible to find a real number l = l(ϵ) > 0, for any interval with length l(ϵ), there exists a number τ = τ(ϵ) in this interval such that ∥x(t+τ) − x(t)∥ < ϵ, ∀ t ∈ ℝ, where ∥⋅∥ is arbitrary norm of ℝn. τ is called to the ϵ-almost period of x, T(x, ϵ) denotes the set of ϵ-almost periods for x and l(ϵ) is called to the length of the inclusion interval for T(x, ϵ). The collection of those functions is denoted by AP(ℝ, ℝn). Let AP(ℝ): = AP(ℝ, ℝ).

Lemma 2.1

([5]) If xAP(ℝ) is differentiable, then forϵ > 0, it follows:

  1. there exists ξϵ ∈ [0, +∞) such that x(ξϵ) ∈ [x+ϵ, x+] and x′(ξϵ) = 0;

  2. there exists ηϵ ∈ [0, +∞) such that fx(ηϵ) ∈ [x, x+ϵ] and x′(ηϵ) = 0.

Lemma 2.2

([6]) If xAP(ℝ) is differentiable, for any interval [a, b] with ba > 0, let ξ, η ∈ [a, b] and

I1={s[ξ,b]:x˙(s)0},I2={s[η,b]:x˙(s)0},

then

x(t)x(ξ)+I1x˙(s)ds,t[ξ,b],x(t)x(η)+I2x˙(s)ds,t[η,b].

Lemma 2.3

([6]) Assume that xAP(ℝ), there exist ξ ∈ [a, b], ξ ∈ (−∞, a] and ξ ∈ [b, +∞) such that

x(ξ_)=x(ξ¯)andx(ξ)x(s),s[ξ_,ξ¯].

Lemma 2.4

([6]) Assume that xAP(ℝ), there exist η ∈ [a, b], η ∈ (−∞, a] and η ∈ [b, +∞) such that

x(η_)=x(η¯)andx(η)x(s),s[η_,η¯].

Lemma 2.5

([7]) If fAP(ℝ) and f¯=m(f)=1T0Tf(s)ds>0, then we have

1Taa+Tf(s)dsf¯2,3f¯2,TT0,

where a is an arbitrary real constant and T0 > 0 is a constant independent of a.

The method to be used in this paper involves the applications of the continuation theorem of coincidence degree.

Let 𝕏 and 𝕐 be real Banach spaces, L : DomL ⊆ 𝕏 → 𝕐 be a linear mapping and N : 𝕏 → 𝕐 be a continuous mapping. The mapping L is called a Fredholm mapping of index zero if ImL is closed in 𝕐 and dimKerL = codimImL < +∞. If L is a Fredholm mapping of index zero and there exist continuous projectors P : 𝕏 → 𝕏 and Q: 𝕐 → 𝕐 such that ImP = KerL, Ker Q = ImL = Im (IQ). It follows that LP = L|DomL∩KerP : (IP)𝕏 → ImL is invertible and its inverse is denoted by KP. If Ω is an open bounded subset of 𝕏, the mapping N will be called L-compact on Ω if QN(Ω) is bounded and KP(IQ)N : Ω → 𝕏 is compact. Since Im Q is isomorphic to KerL, there exists an isomorphism J : ImQ → KerL.

Lemma 2.6

([9]) Let Ω ⊆ 𝕏 be an open bounded set, L be a Fredholm mapping of index zero and N be L-compact on Ω. If all the following conditions hold:

  1. Lxλ Nx, ∀ x∂Ω ∩ DomL, λ ∈ (0, 1);

  2. QNx ≠ 0, ∀ x∂Ω ∩ KerL;

  3. deg{JQN, Ω ∩ KerL, 0} ≠ 0, where deg(⋅, ⋅, ⋅) is the Brouwer degree.

Then Lx = Nx has a solution on Ω ∩ DomL.

For fAP(ℝ), we denote by

Λ(f)={ϖR:limT1T0Tf(s)eiϖsds0},mod(f)={j=1mnjϖj:njZ,mN,ϖjΛ(f),j=1,2,m}

the set of Fourier exponents and the module of f, respectively.

3 Almost periodic solution

Let

f2+:=ln3(a¯2b¯2)c¯+3(a¯2b¯2)ω2,f3+:=ln3a¯3b¯3+3a¯3ω2,μ(s)=a2(s)b2(s)β1(s)ef3+,ν(s)=a3(s)β2(s)ef2+,sR,

where ω is defined as that in (3.3).

Theorem 3.1

Assume that (H1) holds. Suppose further that

  1. μ > 0 and ν > 0,

then system (1.2) admits at least one positive almost periodic solution.

Proof

Under the invariant transformation (x1, x2, x3)T = (ey1, ey2, ey3)T, system (1.2) reduces to

y˙1(t)=a1(t)+b1(t)ey2(tτ(t))ey1(t):=F1(t),y˙2(t)=a2(t)b2(t)c(t)ey2(t)β1(t)ey3(tδ(t)):=F2(t),y˙3(t)=a3(t)b3(t)ey3(t)β2(t)ey2(tσ(t)):=F3(t). (3.0)

It is easy to see that if system (3.0) has one almost periodic solution (y1, y2, y3)T, then (x1, x2, x3)T = (ey1, ey2, ey3)T is a positive almost periodic solution of system (1.2). Therefore, to complete the proof it suffices to show that system (3.0) has one almost periodic solution.

Take 𝕏 = 𝕐 = 𝕍1 ⨁ 𝕍2, where

V1={z=(y1,y2,y3)TAP(R,R3):mod(yi)mod(Li),ϖΛ(yi),|ϖ|θ0,i=1,2,3},V2={z=(y1,y2,y3)T(k1,k2,k3)T,k1,k2,k3R},

where

L1=L1(t,φ)=a1(t)+b1(t)eφ2(τ(0))eφ1(0),L2=L2(t,φ)=a2(t)b2(t)c(t)eφ2(0)β1(t)eφ3(δ(0)),L3=L3(t,φ)=a3(t)b3(t)eφ3(0)β2(t)eφ2(σ(0)),

φ = (φ1, φ2, φ3)TC([−τ0, 0], ℝ3), τ0 = maxi = 1, 2{τ+, δ+, σ+}, θ0 is a given positive constant. Define the norm

z=max{supsR|y1(s)|,supsR|y2(s)|,supsR|y3(s)|},z=(y1,y2,y3)TX=Y,

then 𝕏 and 𝕐 are Banach spaces with the norm ∥⋅∥. Set

L:DomLXY,Lz=L(y1,y2,y3)T=(y1,y2,y3)T,

where DomL = {z = (y1, y2, y3)T ∈ 𝕏: y1, y2, y3C1(ℝ)} and

N:XY,Nz=Ny1(t)y2(t)y3(t)=F1(t)F2(t)F3(t).

With these notations, system (3.0) can be written in the form

Lz=Nz,zX.

It is not difficult to verify that KerL = 𝕍2, ImL = 𝕍1 is closed in 𝕐 and dim KerL = 3 = codim ImL. Therefore, L is a Fredholm mapping of index zero (see Lemma 2.12 in [6]). Now define two projectors P: 𝕏 → 𝕏 and Q: 𝕐 → 𝕐 as

Pz=Py1y2y3=m(y1)m(y2)m(y3)=Qz,z=y1y2y3X=Y.

Then P and Q are continuous projectors such that ImP = KerL and ImL = KerQ = Im(IQ).

Furthermore, through an easy computation we find that the inverse KP : ImL → KerP ∩ DomL of LP has the form

KPz=KPy1y2y3=0ty1(s)dsm0ty1(s)ds0ty2(s)dsm0ty2(s)ds0ty3(s)dsm0ty3(s)ds,z=y1y2y3ImL.

Then QN : 𝕏 → 𝕐 and KP(IQ)N : 𝕏 → 𝕏 read

QNz=QNy1y2y3=m(F1)m(F2)m(F3),

KP(IQ)Nz=KP(IQ)Ny1y2y3=f[y1(t)]Qf[y1(t)]f[y2(t)]Qf[y2(t)]f[y3(t)]Qf[y3(t)],zImL,

where f(x) is defined by f[x(t)] = 0t [Nx(s) − QNx(s)] ds.ąd’ Then N is L-compact on Ω (see Lemma 2.13 in [6]).

In order to apply Lemma 2.6, we need to search for an appropriate open-bounded subset Ω.

Corresponding to the operator equation Lz = λ z, λ ∈ (0, 1), we have

y˙1(t)=λa1(t)+b1(t)ey2(tτ(t))ey1(t),y˙2(t)=λa2(t)b2(t)c(t)ey2(t)β1(t)ey3(tδ(t)),y˙3(t)=λa3(t)b3(t)ey3(t)β2(t)ey2(tσ(t)). (3.1)

Suppose that z = (y1, y2, y3)T ∈ DomL ⊆ 𝕏 is a solution of system (3.1) for some λ ∈ (0, 1), where DomL = {z = (y1, y2, y3)T ∈ 𝕏 : yiC1(ℝ), i = 1, 2, 3}.

There must exist {αn(i):nN+} such that

yi(αn(i))yi+1n,yi+,yi+=supsRyi(s),i=1,2,3. (3.2)

From (H1) and Lemma 2.5, ∀ t0 ∈ ℝ, there exists a constant ω ∈ (0, +∞) independent of t0 such that

1Tt0t0+Tai(s)dsa¯i2,3a¯i2,1Tt0t0+Tbi(s)dsb¯i2,3b¯i2,1Tt0t0+T[a2(s)b2(s)]dsa¯2b¯22,3(a¯2b¯2)2,Tω,i=1,3. (3.3)

For ∀ n0 ∈ ℕ+, we consider αn0(i)ω,αn0(i) (i = 1, 2, 3), where ω is defined as that in (3.3). By Lemma 2.3, there exist ξiαn0(i)ω,αn0(i),ξ_i,αn0(i)ω and ξ¯iαn0(i),+ such that

yi(ξ_i)=yi(ξ¯i),yi(ξi)yi(s),sξ_i,ξ¯i,i=1,2,3. (3.4)

By (3.4), we obtain from system (3.1) that

0=ξ_1ξ¯1a1(s)+b1(s)ey2(sτ(s))ey1(s)ds,0=ξ_2ξ¯2a2(s)b2(s)c(s)ey2(s)β1(s)ey3(sδ(s))ds,0=ξ_3ξ¯3a3(s)b3(s)ey3(s)β2(s)ey2(sσ(s))ds. (3.5)

From the second equation of system (3.5), it follows from (3.3)-(3.4) that

c¯2(ξ¯2ξ_2)ey2(ξ2)ξ_2ξ¯2c(s)ey2(s)ds=ξ_2ξ¯2a2(s)b2(s)β1(s)ey3(sδ(s))dsξ_2ξ¯2[a2(s)b2(s)]ds3(a¯2b¯2)2(ξ¯2ξ_2),

which implies that

y2(ξ2)ln3(a¯2b¯2)c¯. (3.6)

Let Ii=sξi,αn0(i):y˙i(s)0(i=1,2,3). It follows from the first equation of system (3.1) and (3.3) that

I2y˙2(s)ds=I2λa2(s)b2(s)c(s)ey2(s)β1(s)ey3(sδ(s))dsI2a2(s)b2(s)dsαn0(2)ωαn0(2)a2(s)b2(s)ds3(a¯2b¯2)ω2. (3.7)

By Lemma 2.2, it follows from (3.6)-(3.7) that

y2(t)y2(ξ2)+I2y˙2(s)dsln3(a¯2b¯2)c¯+3(a¯2b¯2)ω2:=f2+,t[ξ2,αn0(2)],

which implies that

y2(αn0(2))f2+.

In view of (3.2), letting n0 → +∞ in the above inequality leads to

y2+=limn0+y2(αn0(2))f2+. (3.8)

Similar to the arguments as that in (3.6)-(3.7), we can obtain from the third equation of system (3.5) that

y3(ξ3)ln3a¯3b¯3andI3y˙3(s)ds3a¯3ω2. (3.9)

By a similar discussion as that in (3.8), it follows from (3.9) that

y3+ln3a¯3b¯3+3a¯3ω2:=f3+. (3.10)

On the other hand, from (H2) and Lemma 2.5, ∀ t0 ∈ ℝ, there exists a constant l ∈ [ω, +∞) independent of t0 such that

1Tt0t0+Tμ(s)dsμ¯2,3μ¯2,1Tt0t0+Tν(s)dsν¯2,3ν¯2,Tl. (3.11)

For ∀ n0 ∈ ℤ, by Lemma 2.4, there exist ηi ∈ [n0l, n0l+l], ηi ∈ (−∞, n0l] and ηi ∈ [n0l+l, +∞) such that

yi(η_i)=yi(η¯i),yi(ηi)yi(s),sη_i,η¯i,i=2,3. (3.12)

By (3.12), we obtain from system (3.1) that

0=η_2η¯2a2(s)b2(s)c(s)ey2(s)β1(s)ey3(sδ(s))ds,0=η_3η¯3a3(s)b3(s)ey3(s)β2(s)ey2(sσ(s))ds. (3.13)

By the first equation of system (3.13), we have from (3.3), (3.11) and (3.12) that

3c¯2(η¯2η_2)ey2(η2)η_2η¯2c(s)ey2(s)ds=η_2η¯2a2(s)b2(s)β1(s)ey3(sδ(s))dsη_2η¯2a2(s)b2(s)β1(s)ef3+ds=η_2η¯2μ(s)dsμ¯2(η¯2η_2),

which yields that

y2(η2)lnμ¯3c¯. (3.14)

Furthermore, by the second equation of system (3.1), it follows that

n0ln0l+l|y˙2(s)|ds=n0ln0l+lλa2(s)b2(s)c(s)ey2(s)β1(s)ey3(sδ(s))dsa2++b2++c+ef2++β1+ef3+l. (3.15)

Then

y2(t)y2(η2)n0ln0l+l|y˙2(s)|dslnμ¯3c¯a2++b2++c+ef2++β1+ef3+l:=f2,t[n0l,n0l+l]. (3.16)

Obviously, f2 is a constant independent of n0. So it follows from (3.16) that

y2=infsRy2(s)=infn0Zmins[n0l,n0l+l]y2(s)infn0Zf2=f2. (3.17)

Similar to the arguments as that in (3.14)-(3.15), we can obtain from the third equation of systems (3.1) and (3.13) that

y3(η3)lnν¯3b¯3 (3.18)

and

n0ln0l+l|y˙3(s)|ds=n0ln0l+lλa3(s)b3(s)ey3(s)β2(s)ey2(sσ(s))dsa3++b3+ef3++β2+ef2+l. (3.19)

It follows from (3.18)-(3.19) that

y3(t)y3(η3)n0ln0l+l|y˙3(s)|dslnν¯3b¯3a3++b3+ef3++β2+ef2+l:=f3,t[n0l,n0l+l]. (3.20)

Obviously, f3 is a constant independent of n0. So it follows from (3.20) that

y3=infsRy3(s)=infn0Zmins[n0l,n0l+l]y3(s)infn0Zf3=f3. (3.21)

Finally, there must exist {θn : n ∈ ℕ+} such that

y1(θn)y1,y1+1n, (3.22)

where y1 = infs ∈ ℝy1(s). For ∀ n0 ∈ ℕ+, we consider [θn0ω, θn0], where ω is defined as that in (3.3). By Lemma 2.4, there exist η1 ∈ [θn0ω, θn0], η1 ∈ (−∞, θn0ω] and η1 ∈ [θn0, +∞) such that

y1(η_1)=y1(η¯1),y1(η1)y1(s),sη_1,η¯1. (3.23)

In view of the first equation of system (3.1), we get

0=η_1η¯1a1(s)+b1(s)ey2(sτ(s))ey1(s)ds,

which implies from (3.2) and (3.23) that

3a¯12(η¯1η_1)η_1η¯1a1(s)ds=η_1η¯1b1(s)ey2(sτ(s))ey1(s)dsb¯1ef22ey1(η1)(η¯1η_1),

which yields that

y1(η1)lnb¯1ef23a¯1. (3.24)

Let J = {s ∈ [η1, θn0]: 1(s) ≤ 0}. Then we have from the first equation of system (3.1) and (3.2) that

Jy˙1(s)ds=Jλa1(s)+b1(s)ey2(sτ(s))ey1(s)dsJa1(s)dsθn0ωθn0a1(s)ds3a¯1ω2. (3.25)

By Lemma 2.2, we get from (3.24)-(3.25) that

y1(t)y1(η1)+Jy˙1(s)dslnb¯1ef23a¯13a¯1ω2:=f1,tη1,θn0.

So

y1(θn0)f1y1f11n0.

Letting n0 → +∞ in the last inequality leads to

y1f1. (3.26)

Furthermore, by the first equation of system (3.5), we have from (3.3)-(3.4) that

a¯12(ξ¯1ξ_1)ξ_1ξ¯1a1(s)ds=ξ_1ξ¯1b1(s)ey2(sτ(s))ey1(s)ds3b¯1ef2+2ey1(ξ1)(ξ¯1ξ_1),

which yields that

y1(ξ1)ln3b¯1ef2+a¯1. (3.27)

And

I1y˙1(s)ds=I1a1(s)+b1(s)ey2(sτ(s))ey1(s)dsαn0(1)ωαn0(1)b1(s)ey2(sτ(s))ey1(s)ds3b¯1ef2+ω2ef1. (3.28)

Similar to the argument as that in (3.8), we obtain from (3.27)-(3.28) that

y1+ln3b¯1ef2+a¯1+3b¯1ef2+ω2ef1:=f1+. (3.29)

Set C=i=13(|fi|+|fi+|)+1. Clearly, C is independent of λ ∈ (0, 1). Let Ω = {z ∈ 𝕏 : ∥z𝕏 < C}. Therefore, Ω satisfies condition (a) of Lemma 2.6.

Now we show that condition (b) of Lemma 2.6 holds, i.e., we prove that QNz ≠ 0 for all z = (y1, y2, y3)T∂Ω ∩ KerL = ∂Ω ∩ ℝ3. If it is not true, then there exists at least one constant vector z0=(y10,y20,y30)T ∂Ω such that

0=ma1(t)+b1(t)ey20ey10,0=ma2(t)b2(t)c(t)ey20β1(t)ey30,0=ma3(t)b3(t)ey30β2(t)ey20.

Similar to the argument as that in (3.8), (3.10), (3.17), (3.21), (3.26) and (3.29), it follows that

f1y10f1+,f2y20f2+,f3y30f3+.

Then z0Ω ∩ ℝ3. This contradicts the fact that z0∂Ω. This proves that condition (b) of Lemma 2.6 holds.

Finally, we will show that condition (c) of Lemma 2.6 is satisfied. Let us consider the homotopy

H(ι,z)=ιQNz+(1ι)Φz,(ι,z)[0,1]×R3,

where

Φz=Φy1y2y3=a¯1+b¯1ef2ey1a¯2b¯2c¯ey2a¯3b¯3ey3.

From the above discussion it is easy to verify that H(ι, z) ≠ 0 on ∂Ω ∩ KerL, ∀ ι ∈ [0, 1]. Further, Φ z = 0 has a solution:

(y1+,y2+,y3+)T=lnb¯1ef2a¯1,lna¯2b¯2c¯,lna¯3b¯3TΩ.

A direct computation yields

deg(Φ,ΩKerL,0)=signb¯1ef2ey1+000c¯ey2+000b¯3ey3+=1.

By the invariance property of homotopy, we have

deg(JQN,ΩKerL,0)=deg(QN,ΩKerL,0)=deg(Φ,ΩKerL,0)0,

where J is the identity mapping since ImQ = KerL. Obviously, all the conditions of Lemma 2.6 are satisfied. Therefore, system (3.0) has at least one almost periodic solution, that is, system (1.2) has at least one positive almost periodic solution. This completes the proof. □

Remark 3.1

In [2], Chen also obtained the existence of positive almost periodic solutions of system (1.1). However, the result in [2] required the following condition:

  1. ai>0,bi>0,βj>0 and c > 0, i = 1, 2, 3, j = 1, 2.

  2. a1>a2+m,c>b1++β2+ and b3>β1+, where m is some positive constant.

Obviously, (H1) in Theorem 3.1 is weaker than (F1). Further, by Theorem 3.1, it is easy to obtain the existence of positive almost periodic of system (1.1) without (F2). So the work of this paper extends and improves the result in [2].

Corollary 3.1

Assume that (H1)-(H2) hold. Suppose further that ai, bi, βj, c, τ, σ and δ of system (1.2) are continuous nonnegative periodic functions with periods αi, ζi, ηj, ξ, ρ, ϱ and ς, respectively, i = 1, 2, 3, j = 1, 2, then system (1.2) has at least one positive almost periodic solution.

By Corollary 3.1, we obtain

Corollary 3.2

Assume that (H1)-(H2) hold. Suppose further that ai, bi, βj, c, τ, σ and δ of system (1.2) are continuous nonnegative ω-periodic functions, i = 1, 2, 3, j = 1, 2, then system (1.2) has at least one positive ω-periodic solution.

Let

M1:=(a2b2)+c,M2:=a3+b3,μ0(s)=a2(s)b2(s)β1(s)M2,ν0(s)=a3(s)β2(s)M1,sR.

In order to complement or improve (H2) in Theorem 3.1, we give the following result:

Theorem 3.2

Assume that (H1) holds. Suppose further that

  1. a1>0,b3>0 and c > 0,

  2. μ0 > 0 and ν0 > 0,

then system (1.2) admits at least one positive almost periodic solution.

Proof

Proceeding as in the proof of Theorem 3.1, in order to use Lemma 2.6, it remains to search for an appropriate open and bounded subset Ω ⊆ 𝕏.

Consider the operator equations (3.1). By Lemma 2.1, for ∀ϵ ∈ (0, 1), there exist ξi = ξi(ϵ) such that

y˙i(ξi)=0,yi(ξi)[yi+ϵ,yi+],yi+=supsRyi(s),i=1,2,3. (3.30)

From system (3.1), it follows from (3.30) that

0=a1(ξ1)+b1(ξ1)ey2(ξ1τ(ξ1))ey1(ξ1),0=a2(ξ2)b2(ξ2)c(ξ2)ey2(ξ2)β1(ξ2)ey3(ξ2δ(ξ2)),0=a3(ξ3)b3(ξ3)ey3(ξ3)β2(ξ3)ey2(ξ3σ(ξ3)). (3.31)

By the second equation of system (3.31), we have from (3.30) that

cey2+ϵc(ξ2)ey2(ξ2)=a2(ξ2)b2(ξ2)β1(ξ2)ey3(ξ2σ(ξ2))(a2b2)+=supsR[a2(s)b2(s)],

which yields that

y2+ln(a2b2)+eϵc.

Letting ϵ → 0 in the above inequality, we get

y2+ln(a2b2)+c:=g2+. (3.32)

So we can obtain from the first equation of system (3.31) that

a1ey1+ϵa1(ξ1)ey1(ξ1)=b1(ξ1)ey2(ξ1τ(ξ1))b1+ey2+b1+(a2b2)+eϵc.

Letting ϵ → 0 in the above inequality, we get

y1+lnb1+(a2b2)+a1c:=g1+. (3.33)

Similar to the arguments as that in (3.32)-(3.33), it follows from the third equation of system (3.31) that

y3+lna3+b3:=g3+. (3.34)

By the similar discussions as that in (3.17), (3.21) and (3.26), we could find gi such that

yigi,i=1,2,3.

Set C0=i=13(|gi|+|gi+|)+1. Clearly, C0 is independent of λ ∈ (0, 1). Let Ω = {z ∈ 𝕏 : ∥z∥𝕏 < C0}. From the proof in Theorem 3.1, it is easy to verify that Ω satisfies conditions (a)-(c) of Lemma 2.6. Obviously, all the conditions of Lemma 2.6 are satisfied. Therefore, system (3.0) has one almost periodic solution, that is, system (1.2) has at least one positive almost periodic solution. This completes the proof. □

Corollary 3.3

Assume that (H1), (H3) and (H4) hold. Suppose further that ai, bi, βj, c, τ, σ and δ of system (1.2) are continuous nonnegative periodic functions with periods αi, ζi, ηj, ξ, ρ, ϱ and ς, respectively, i = 1, 2, 3, j = 1, 2, then system (1.2) has at least one positive almost periodic solution.

From Corollary 3.3, we have

Corollary 3.4

Assume that (H1), (H3) and (H4) hold. Suppose further that ai, bi, βj, c, τ, σ and δ of system (1.2) are continuous nonnegative ω-periodic functions, i = 1, 2, 3, j = 1, 2, then system (1.2) has at least one positive ω-periodic solution.

4 Global asymptotical stability

Theorem 4.1

Assume that (H1)-(H2) hold. Suppose further that

  1. τ, σ, δC1(ℝ), sups ∈ ℝ{τ̇(s), σ̇(s), δ̇(s)} < 1, a1 > 0 and

    infsRc(s)b1(ξ1(s))1τ˙(ξ1(s))β2(φ1(s))1σ˙(φ1(s))>0,infsRb3(s)β1(ψ1(s))1δ˙(ψ1(s))>0,

    where ξ−1, φ−1 and ψ−1 are the inverse functions of ξ(t) = tτ(t), φ(t) = tσ(t) and ψ(t) = tδ(t), respectively, ∀ t ∈ ℝ.

Then system (1.2) has a unique positive almost periodic solution, which is globally asymptotically stable.

Proof

By Theorem 3.1, we know that system (1.2) has at least one positive almost periodic solution (x1,x2,x3)T. Suppose that (x1, x2, x3)T is another positive solution of system (1.2).

From (H5), there must exist 0 < θ < a1 such that

infsRc(s)b1(ξ1(s))1τ˙(ξ1(s))β2(φ1(s))1σ˙(φ1(s))>θ,infsRb3(s)β1(ψ1(s))1δ˙(ψ1(s))>θ.

Define

V(t)=i=16Vi(t),

where

V1(t)=|x1(t)x1(t)|,Vi(t)=|lnxi(t)lnxi(t)|,i=2,3,

V4(t)=tτ(t)tb1(ξ1(s))1τ˙(ξ1(s))|x2(s)x2(s)|ds,V5(t)=tσ(t)tβ2(φ1(s))1σ˙(φ1(s))|x2(s)x2(s)|ds,V6(t)=tδ(t)tβ1(ψ1(s))1δ˙(ψ1(s))|x3(s)x3(s)|ds.

By calculating the upper right derivative of Vi(i = 1, 2, 3, 4, 5, 6) along the positive solution of system (1.2), it follows that

D+V1(t)=sign[x1(t)x1(t)]{a1(t)x1(t)+b1(t)x2(tτ(t))a1(t)x1(t)+b1(t)x2(tτ(t))}a1(t)|x1(t)x1(t)|+b1(t)|x2(tτ(t))x2(tτ(t))|, (4.1)

D+V2(t)=sign[x2(t)x2(t)]{a2(t)b2(t)c(t)x2(t)β1(t)x3(tδ(t))a2(t)b2(t)c(t)x2(t)β1(t)x3(tδ(t))}c(t)|x2(t)x2(t)|+β1(t)|x3(tδ(t))x3(tδ(t))|, (4.2)

D+V3(t)=sign[x3(t)x3(t)]{a3(t)b3(t)x3(t)β2(t)x2(tσ(t))a3(t)b3(t)x3(t)β2(t)x2(tσ(t))}b3(t)|x3(t)x3(t)|+β2(t)|x2(tσ(t))x2(tσ(t))|, (4.3)

D+V4(t)=b1(ξ1(t))1τ˙(ξ1(t))|x2(t)x2(t)|b1(t)|x2(tτ(t))x2(tτ(t))|, (4.4)

D+V5(t)=β2(φ1(t))1σ˙(φ1(t))|x2(t)x2(t)|β2(t)|x2(tσ(t))x2(tσ(t))| (4.5)

and

D+V6(t)=β1(ψ1(t))1δ˙(ψ1(t))|x3(t)x3(t)|β1(t)|x3(tδ(t))x3(tδ(t))|. (4.6)

Together with (4.1)-(4.6), it follows that

D+V(t)a1(t)|x1(t)x1(t)|{c(t)b1(ξ1(t))1τ˙(ξ1(t))β2(φ1(t))1σ˙(φ1(t))}|x2(t)x2(t)|{b3(t)β1(ψ1(t))1δ˙(ψ1(t))}|x3(t)x3(t)|θi=13|xi(t)xi(t)|,t0.

Therefore, V is non-increasing. Integrating the last inequality from 0 to t leads to

V(t)+θi=130t|xi(s)xi(s)|dsV(0)<+,t0,

that is,

i=130+|xi(s)xi(s)|ds<+,

which implies that

i=13lims+|xi(s)xi(s)|=0.

Thus, the almost periodic solution of system (1.2) is globally exponentially stable.

Next, we show that there is only one positive almost periodic solution of system (1.2). For any two positive almost periodic solutions (x1, x2, x3)T and (x1, x2, x3)T of system (1.2), we claim that xi(t) ≡ xi(t), ∀ t ∈ ℝ, i = 1, 2, 3. If not, without loss of generality, there must be at least one t0 ∈ ℝ such that x1(t0) ≠ x1(t0), i.e., |x1(t0)−x1(t0)|: = l > 0. The global asymptotical stability implies that there exists t1 > t0 such that

|x1(t)x¯1(t)|<l4,tt1. (4.7)

By the almost periodicity of x1 and x1, there must exist l1 > 0 and τ0 ∈ [t1t0, t1t0+l1] such that

|x1(t+τ0)x1(t)|<l4,|x¯1(t+τ0)x¯1(t)|<l4,tR. (4.8)

So we can easily know from (4.7)-(4.8) that

l=|x1(t0)x¯1(t0)||x1(t0)x1(t0+τ0)|+|x1(t0+τ0)x¯1(t0+τ0)|+|x¯1(t0+τ0)x¯1(t0)|<l4+l4+l4=3l4,

which is a contradiction. Then x1(t) ≡ x1(t), ∀ t ∈ ℝ. Similarly, we can prove xi(t) ≡ xi(t), ∀ t ∈ ℝ, i = 2, 3. Therefore, the almost periodic solution of system (1.2) is unique. This completes the proof. □

Together with Theorem 3.2, we can easily show that

Theorem 4.2

Assume that (H1), (H3), (H4) and (H5) hold, then system (1.2) has a unique positive almost periodic solution, which is globally asymptotically stable.

5 Example and simulations

Example 5.1

Considering the following delayed two-species competitive model with stage structure and different periods:

x˙1(t)=[1+|sin(2t)|]x1(t)+0.1|cos(3t)|x2(t1),x˙2(t)=x2(t)|cos(3t)|x2(t)x22(t)0.1sin2(3t)x2(t)x3(t1),x˙3(t)=x3(t)|cos(3t)|[1+0.99sin(2t)]x3(t)0.1|cos(3t)|x2(t2). (5.1)

Then system (5.1) has a unique positive almost periodic solution, which is globally asymptotically stable.

Proof

Corresponding to system (1.2), we have

a1(s)a2(s)a3(s)=1+|sin(2s)|1|cos(3s)|,b1(s)b2(s)b3(s)=0.1|cos(3s)||cos(3s)|1+sin(2s),β1(s)β2(s)=0.1sin2(3s)0.1|cos(3s)|,c(s)=1,sR.

Obviously, (H1) in Theorem 3.1 holds. Further, ā2 = c = 1, a¯3=b¯2=2π,b¯3=1+1.98π and β1+=β2+=0.1.

Let e1(t)=|sin(2t)| and e2(t)=|cos(3t)|,t>0. Note that e¯i=2π0.64, i = 1, 2. For ∀ t0 ∈ ℝ, we have

Ei(t)=t0t0+tei(s)dst1π,3π0.32,0.955,t1,i=1,2.

Therefore, we can choose ω = 1 so that (3.3) holds. By an easy calculation, we obtain that

f2+0.62,f3+1.12,μ¯>0.06>0,ν¯>0.472>0,

which implies that (H2) in Theorem 3.1 holds. Further, it is easy to verify that (H5) in Theorem 4.1 is satisfied. Therefore, all the conditions of Theorems 3.1-4.1 are satisfied. By Theorems 3.1-4.1, system (5.1) has a unique positive almost periodic solution, which is globally asymptotically stable (see Figures 2-5). This completes the proof. □

Figure 1 
Almost periodic oscillations of system (5.1)
Figure 1

Almost periodic oscillations of system (5.1)

Figure 2 
Global asymptotical stability of state variable x1 of system (5.1)
Figure 2

Global asymptotical stability of state variable x1 of system (5.1)

Figure 3 
Global asymptotical stability of state variable x2 of system (5.1)
Figure 3

Global asymptotical stability of state variable x2 of system (5.1)

Figure 4 
Global asymptotical stability of state variable x3 of system (5.1)
Figure 4

Global asymptotical stability of state variable x3 of system (5.1)

Remark 5.1

Clearly, system (5.1) is with some nonnegative coefficients. Thus, condition (F1) in Remark 3.1 is not satisfied. Further, corresponding to system (1.1), b3 = 0.01 < 0.1 = β1+, which implies that condition (F2) in Remark 3.1 is invalid. Therefore, by the main result obtained by paper [2], it is impossible to obtain the existence and global asymptotical stability of positive almost periodic solutions of system (5.1). So the work of this paper extends and improves the result in [2].

6 Conclusions

The stage-structured models have been studied extensively, and many important phenomena have been observed in recent years. In this paper we study an almost periodic nonautonomous delayed two-species competitive model with stage structure, and this motivation comes from a nonautonomous delayed two-species competitive model. We obtain easily verifiable sufficient criteria for the existence and globally asymptotic stability of positive almost periodic solutions of the above model. In order to obtain a more accurate description of the ecological system perturbed by human exploitation activities such as planting and harvesting and so on, we need to consider the impulsive differential equations. In this paper, we only studied system (1.2) without impulses. Whether system (1.2) with impulses can be discussed in the same methods or not is still an open problem. We will continue to study this problem in the future.

Acknowledgement

This work are supported by National Nature Science Foundation under Grant No. 11461082, 11601474 and 61472093, Key laboratory of numerical simulation of Sichuan Province under Grant No. 2017KF002.

References

[1] Aiello W.G., Freedman H.I., A time-delay model of single-specise growth with stage structure, Math. Biosci., 1990, 101, 139-153.10.1016/0025-5564(90)90019-USearch in Google Scholar

[2] Chen F.D., Almost periodic solution of the non-autonomous two-species competitive model with stage structure, Appl. Math. Comput., 2006, 181, 685-693.10.1016/j.amc.2006.01.055Search in Google Scholar

[3] Gopalsamy K., Global asymptotic stability in an almost periodic Lotka-Volterra system, J. Aust. Math. Soc. Ser. B, 1986, 27, 346-360.10.1017/S0334270000004975Search in Google Scholar

[4] Zhou H., Agarwal P., Existence of almost periodic solution for neutral Nicholson blowflies model, Advances in Difference Equations, 2017, 2017, No. 329.10.1186/s13662-017-1382-0Search in Google Scholar

[5] Zhang T.W., Multiplicity of positive almost periodic solutions in a delayed Hassell-Varley-type predator-prey model with harvesting on prey, Math. Meth. Appl. Sci., 2014, 37, 686-697.10.1002/mma.2826Search in Google Scholar

[6] Zhang T.W., Almost periodic oscillations in a generalized Mackey-Glass model of respiratory dynamics with several delays, Int. J. Biomath., 2014, 7, 1-22.10.1142/S1793524514500296Search in Google Scholar

[7] Fink A.M., Almost Periodic Differential Equation, Spring-Verlag, Berlin, Heidleberg, New York, 1974.10.1007/BFb0070324Search in Google Scholar

[8] He C.Y., Almost Periodic Differential Equations, Higher Education Publishing House, Beijing, 1992 (in Chinese).Search in Google Scholar

[9] Gaines R.E., Mawhin J.L., Coincidence Degree and Nonlinear Differential Equations, Springer, Berlin, 1977.10.1007/BFb0089537Search in Google Scholar

Received: 2018-07-04
Accepted: 2019-02-20
Published Online: 2019-05-16

© 2019 Liu et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Articles in the same Issue

  1. Regular Articles
  2. On the Gevrey ultradifferentiability of weak solutions of an abstract evolution equation with a scalar type spectral operator of orders less than one
  3. Centralizers of automorphisms permuting free generators
  4. Extreme points and support points of conformal mappings
  5. Arithmetical properties of double Möbius-Bernoulli numbers
  6. The product of quasi-ideal refined generalised quasi-adequate transversals
  7. Characterizations of the Solution Sets of Generalized Convex Fuzzy Optimization Problem
  8. Augmented, free and tensor generalized digroups
  9. Time-dependent attractor of wave equations with nonlinear damping and linear memory
  10. A new smoothing method for solving nonlinear complementarity problems
  11. Almost periodic solution of a discrete competitive system with delays and feedback controls
  12. On a problem of Hasse and Ramachandra
  13. Hopf bifurcation and stability in a Beddington-DeAngelis predator-prey model with stage structure for predator and time delay incorporating prey refuge
  14. A note on the formulas for the Drazin inverse of the sum of two matrices
  15. Completeness theorem for probability models with finitely many valued measure
  16. Periodic solution for ϕ-Laplacian neutral differential equation
  17. Asymptotic orbital shadowing property for diffeomorphisms
  18. Modular equations of a continued fraction of order six
  19. Solutions with concentration and cavitation to the Riemann problem for the isentropic relativistic Euler system for the extended Chaplygin gas
  20. Stability Problems and Analytical Integration for the Clebsch’s System
  21. Topological Indices of Para-line Graphs of V-Phenylenic Nanostructures
  22. On split Lie color triple systems
  23. Triangular Surface Patch Based on Bivariate Meyer-König-Zeller Operator
  24. Generators for maximal subgroups of Conway group Co1
  25. Positivity preserving operator splitting nonstandard finite difference methods for SEIR reaction diffusion model
  26. Characterizations of Convex spaces and Anti-matroids via Derived Operators
  27. On Partitions and Arf Semigroups
  28. Arithmetic properties for Andrews’ (48,6)- and (48,18)-singular overpartitions
  29. A concise proof to the spectral and nuclear norm bounds through tensor partitions
  30. A categorical approach to abstract convex spaces and interval spaces
  31. Dynamics of two-species delayed competitive stage-structured model described by differential-difference equations
  32. Parity results for broken 11-diamond partitions
  33. A new fourth power mean of two-term exponential sums
  34. The new operations on complete ideals
  35. Soft covering based rough graphs and corresponding decision making
  36. Complete convergence for arrays of ratios of order statistics
  37. Sufficient and necessary conditions of convergence for ρ͠ mixing random variables
  38. Attractors of dynamical systems in locally compact spaces
  39. Random attractors for stochastic retarded strongly damped wave equations with additive noise on bounded domains
  40. Statistical approximation properties of λ-Bernstein operators based on q-integers
  41. An investigation of fractional Bagley-Torvik equation
  42. Pentavalent arc-transitive Cayley graphs on Frobenius groups with soluble vertex stabilizer
  43. On the hybrid power mean of two kind different trigonometric sums
  44. Embedding of Supplementary Results in Strong EMT Valuations and Strength
  45. On Diophantine approximation by unlike powers of primes
  46. A General Version of the Nullstellensatz for Arbitrary Fields
  47. A new representation of α-openness, α-continuity, α-irresoluteness, and α-compactness in L-fuzzy pretopological spaces
  48. Random Polygons and Estimations of π
  49. The optimal pebbling of spindle graphs
  50. MBJ-neutrosophic ideals of BCK/BCI-algebras
  51. A note on the structure of a finite group G having a subgroup H maximal in 〈H, Hg
  52. A fuzzy multi-objective linear programming with interval-typed triangular fuzzy numbers
  53. Variational-like inequalities for n-dimensional fuzzy-vector-valued functions and fuzzy optimization
  54. Stability property of the prey free equilibrium point
  55. Rayleigh-Ritz Majorization Error Bounds for the Linear Response Eigenvalue Problem
  56. Hyper-Wiener indices of polyphenyl chains and polyphenyl spiders
  57. Razumikhin-type theorem on time-changed stochastic functional differential equations with Markovian switching
  58. Fixed Points of Meromorphic Functions and Their Higher Order Differences and Shifts
  59. Properties and Inference for a New Class of Generalized Rayleigh Distributions with an Application
  60. Nonfragile observer-based guaranteed cost finite-time control of discrete-time positive impulsive switched systems
  61. Empirical likelihood confidence regions of the parameters in a partially single-index varying-coefficient model
  62. Algebraic loop structures on algebra comultiplications
  63. Two weight estimates for a class of (p, q) type sublinear operators and their commutators
  64. Dynamic of a nonautonomous two-species impulsive competitive system with infinite delays
  65. 2-closures of primitive permutation groups of holomorph type
  66. Monotonicity properties and inequalities related to generalized Grötzsch ring functions
  67. Variation inequalities related to Schrödinger operators on weighted Morrey spaces
  68. Research on cooperation strategy between government and green supply chain based on differential game
  69. Extinction of a two species competitive stage-structured system with the effect of toxic substance and harvesting
  70. *-Ricci soliton on (κ, μ)′-almost Kenmotsu manifolds
  71. Some improved bounds on two energy-like invariants of some derived graphs
  72. Pricing under dynamic risk measures
  73. Finite groups with star-free noncyclic graphs
  74. A degree approach to relationship among fuzzy convex structures, fuzzy closure systems and fuzzy Alexandrov topologies
  75. S-shaped connected component of radial positive solutions for a prescribed mean curvature problem in an annular domain
  76. On Diophantine equations involving Lucas sequences
  77. A new way to represent functions as series
  78. Stability and Hopf bifurcation periodic orbits in delay coupled Lotka-Volterra ring system
  79. Some remarks on a pair of seemingly unrelated regression models
  80. Lyapunov stable homoclinic classes for smooth vector fields
  81. Stabilizers in EQ-algebras
  82. The properties of solutions for several types of Painlevé equations concerning fixed-points, zeros and poles
  83. Spectrum perturbations of compact operators in a Banach space
  84. The non-commuting graph of a non-central hypergroup
  85. Lie symmetry analysis and conservation law for the equation arising from higher order Broer-Kaup equation
  86. Positive solutions of the discrete Dirichlet problem involving the mean curvature operator
  87. Dislocated quasi cone b-metric space over Banach algebra and contraction principles with application to functional equations
  88. On the Gevrey ultradifferentiability of weak solutions of an abstract evolution equation with a scalar type spectral operator on the open semi-axis
  89. Differential polynomials of L-functions with truncated shared values
  90. Exclusion sets in the S-type eigenvalue localization sets for tensors
  91. Continuous linear operators on Orlicz-Bochner spaces
  92. Non-trivial solutions for Schrödinger-Poisson systems involving critical nonlocal term and potential vanishing at infinity
  93. Characterizations of Benson proper efficiency of set-valued optimization in real linear spaces
  94. A quantitative obstruction to collapsing surfaces
  95. Dynamic behaviors of a Lotka-Volterra type predator-prey system with Allee effect on the predator species and density dependent birth rate on the prey species
  96. Coexistence for a kind of stochastic three-species competitive models
  97. Algebraic and qualitative remarks about the family yy′ = (αxm+k–1 + βxmk–1)y + γx2m–2k–1
  98. On the two-term exponential sums and character sums of polynomials
  99. F-biharmonic maps into general Riemannian manifolds
  100. Embeddings of harmonic mixed norm spaces on smoothly bounded domains in ℝn
  101. Asymptotic behavior for non-autonomous stochastic plate equation on unbounded domains
  102. Power graphs and exchange property for resolving sets
  103. On nearly Hurewicz spaces
  104. Least eigenvalue of the connected graphs whose complements are cacti
  105. Determinants of two kinds of matrices whose elements involve sine functions
  106. A characterization of translational hulls of a strongly right type B semigroup
  107. Common fixed point results for two families of multivalued A–dominated contractive mappings on closed ball with applications
  108. Lp estimates for maximal functions along surfaces of revolution on product spaces
  109. Path-induced closure operators on graphs for defining digital Jordan surfaces
  110. Irreducible modules with highest weight vectors over modular Witt and special Lie superalgebras
  111. Existence of periodic solutions with prescribed minimal period of a 2nth-order discrete system
  112. Injective hulls of many-sorted ordered algebras
  113. Random uniform exponential attractor for stochastic non-autonomous reaction-diffusion equation with multiplicative noise in ℝ3
  114. Global properties of virus dynamics with B-cell impairment
  115. The monotonicity of ratios involving arc tangent function with applications
  116. A family of Cantorvals
  117. An asymptotic property of branching-type overloaded polling networks
  118. Almost periodic solutions of a commensalism system with Michaelis-Menten type harvesting on time scales
  119. Explicit order 3/2 Runge-Kutta method for numerical solutions of stochastic differential equations by using Itô-Taylor expansion
  120. L-fuzzy ideals and L-fuzzy subalgebras of Novikov algebras
  121. L-topological-convex spaces generated by L-convex bases
  122. An optimal fourth-order family of modified Cauchy methods for finding solutions of nonlinear equations and their dynamical behavior
  123. New error bounds for linear complementarity problems of Σ-SDD matrices and SB-matrices
  124. Hankel determinant of order three for familiar subsets of analytic functions related with sine function
  125. On some automorphic properties of Galois traces of class invariants from generalized Weber functions of level 5
  126. Results on existence for generalized nD Navier-Stokes equations
  127. Regular Banach space net and abstract-valued Orlicz space of range-varying type
  128. Some properties of pre-quasi operator ideal of type generalized Cesáro sequence space defined by weighted means
  129. On a new convergence in topological spaces
  130. On a fixed point theorem with application to functional equations
  131. Coupled system of a fractional order differential equations with weighted initial conditions
  132. Rough quotient in topological rough sets
  133. Split Hausdorff internal topologies on posets
  134. A preconditioned AOR iterative scheme for systems of linear equations with L-matrics
  135. New handy and accurate approximation for the Gaussian integrals with applications to science and engineering
  136. Special Issue on Graph Theory (GWGT 2019)
  137. The general position problem and strong resolving graphs
  138. Connected domination game played on Cartesian products
  139. On minimum algebraic connectivity of graphs whose complements are bicyclic
  140. A novel method to construct NSSD molecular graphs
Downloaded on 9.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/math-2019-0030/html
Scroll to top button