DDNSDiscrete Dynamics in Nature and Society1607-887X1026-0226Hindawi Publishing Corporation93179810.1155/2010/931798931798Research ArticlePermanence of a Discrete Model of Mutualism with Infinite Deviating ArgumentsLiXuepengYangWenshengZhangBinggenSchool of Mathematics and Computer ScienceFujian Normal UniversityFuzhouFujian 350007Chinafjnu.edu.cn201021012010201015072009130120102010Copyright © 2010This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We propose a discrete model of mutualism with infinite deviating arguments, that is x1(n+1)=x1(n)exp{r1(n)[(K1(n)+α1(n)s=0J2(s)x2(ns))/(1+s=0J2(s)x2(ns))x1(nσ1(n))]},x2(n+1)=x2(n)exp{r2(n)[(K2(n)+α2(n)s=0J1(s)x1(ns))/(1+s=0J1(s)x1(ns))x2(nσ2(n))]}. By some Lemmas, sufficient conditions are obtained for the permanence of the system.

1. Introduction

Chen and You  studied the following two species integro-differential model of mutualism:

dN1(t)dt=r1(t)N1(t)[K1(t)+α1(t)0J2(s)N2(t-s)ds1+0J2(s)N2(t-s)ds-N1(t-σ1(t))],dN2(t)dt=r2(t)N2(t)[K2(t)+α2(t)0J1(s)N1(t-s)ds1+0J1(s)N1(t-s)ds-N2(t-σ2(t))], where ri,Ki,αi, and σi,i=1,2 are continuous functions bounded above and below by positive constants: ai>Ki,  i=1,2;    JiC([0,+),[0,+)) and 0Ji(s)ds=1,  i=1,2. Using the differential inequality theory, they obtained a set of sufficient conditions to ensure the permanence of system (1.1). For more background and biological adjustments of system(1.1), one could refer to  and the references cited therein.

However, many authors  have argued that the discrete time models governed by difference equations are more appropriate than the continuous ones when the populations have nonoverlapping generations. Also, since discrete time models can also provide efficient computational models of continuous models for numerical simulations, it is reasonable to study discrete time models governed by difference equations. Another permanence is one of the most important topics on the study of population dynamics. One of the most interesting questions in mathematical biology concerns the survival of species in ecological models. It is reasonable to ask for conditions under which the system is permanent.

Motivated by the above question, we consider the permanence of the following discrete model of mutualism with infinite deviating arguments:

x1(n+1)=x1(n)exp{r1(n)[K1(n)+α1(n)s=0J2(s)x2(n-s)1+s=0J2(s)x2(n-s)-x1(n-σ1(n))]},x2(n+1)=x2(n)exp{r2(n)[K2(n)+α2(n)s=0J1(s)x1(n-s)1+s=0J1(s)x1(n-s)-x2(n-σ2(n))]}, where xi(n),  i=1,2 is the density of mutualism species i at the nth generation. For {ri(n)},{Ki(n)},{αi(n)},{Ji(n)}, and {σi(n)},  i=1,2 are bounded nonnegative sequences such that

0<rilriu,0<αilαiu,0<KilKiu,0<σilσiu,n=0Ji(n)=1. Here, for any bounded sequence {a(n)}, au=supnNa(n),  al=infnNa(n).

Let σ=supn{σi(n),  i=1,2}, we consider (1.2) together with the following initial condition:

xi(θ)=φi(θ)0,θN[-τ,0]={-τ,-τ+1,,0},φi(0)>0.

It is not difficult to see that solutions of (1.2) and (1.4) are well defined for all n0 and satisfy

xi(n)>0,for    nZ,i=1,2.

The aim of this paper is, by applying the comparison theorem of difference equation and some lemmas, to obtain a set of sufficient conditions which guarantee the permanence of system (1.2).

2. Permanence

In this section, we establish permanence results for system (1.2).

Following Comparison Theorem of difference equation is Theorem 2.6 of [13, page 241].

Lemma.

Let kNk0+={k0,k0+1,,k0+l,},r0. For any fixed k,g(k,r) is a non-decreasing function with respect to r, and for kk0, following inequalities hold: y(k+1)g(k,y(k)),u(k+1)g(k,u(k)). If y(k0)u(k0), then y(k)u(k) for all kk0.

Now let us consider the following single species discrete model:

N(k+1)=N(k)exp{a(k)-b(k)N(k)}, where {a(k)} and {b(k)} are strictly positive sequences of real numbers defined for kN={0,1,2,} and 0<alau,0<blbu. Similar to the proof of Propositions 1 and 3 in , we can obtain the following.

Lemma.

Any solution of system (2.1) with initial condition N(0)>0 satisfies mlimk+infN(k)limk+supN(k)M, where M=1blexp{au-1},m=albuexp{al-buM}.

Lemma 2.3 (see [<xref ref-type="bibr" rid="B13">14</xref>]).

Let x(n) and b(n) be nonnegative sequences defined on N, and c0 is a constant. If x(n)c+s=0n-1b(s)x(s),  for  nN, then x(n)cs=0n-1[1+b(s)],for  nN.

Lemma 2.4 (see [<xref ref-type="bibr" rid="B2">2</xref>]).

Let x:ZR be a nonnegative bounded sequences, and let H:NR be a nonnegative sequence such that n=0Ji(n)=1. Then limn+infx(n)limn+infs=-nH(n-s)x(s)limn+sups=-nH(n-s)x(s)limn+supx(n).

Proposition.

Let (x1(n),x2(n)) be any positive solution of system (1.2), then limn+supxi(n)Mi,i=1,2, where Mi=exp{2riu[Kiu+αiu]},i=1,2.

Proof.

Let (x1(n),x2(n)) be any positive solution of system (1.2), then from the first equation of system (1.2) we have x1(n+1)x1(n)exp{r1(n)[K1(n)+α1(n)s=0J2(s)x2(n-s)1+s=0J2(s)x2(n-s)]}=x1(n)exp{r1(n)[K1(n)1+s=0J2(s)x2(n-s)+α1(n)s=0J2(s)x2(n-s)1+s=0J2(s)x2(n-s)]}x1(n)exp{r1(n)[K1(n)1+α1(n)s=0J2(s)x2(n-s)s=0J2(s)x2(n-s)]}=x1(n)exp{r1(n)[K1(n)+α1(n)]}x1(n)exp{r1u[K1u+α1u]}. Let x1(n)=exp{u1(n)}, then u1(n+1)u1(n)+r1u[K1u+α1u]=r1u[K1u+α1u]+s=0nb(s)x(s), where b(s)={0,0sn-1,1,s=n. When u1(n) is nonnegative sequence, by applying Lemma 2.3, it immediately follows that u1(n+1)2r1u[K1u+α1u]. When u1(n) is negative sequence, (2.12) also holds. From (2.12), we have limn+supx1(n)exp{2r1u[K1u+α1u]}:=M1. By using the second equation of system (1.2), similar to the above analysis, we can obtain limn+supx2(n)exp{2r2u[K2u+α2u]}:=M2. This completes the proof of Proposition 2.5.

Now we are in the position of stating the permanence of system (1.2).

Theorem.

Under the assumption(1.3), system (1.2) is permanent, that is, there exist positive constants mi,Mi,i=1,2 which are independent of the solutions of system (1.2) such that, for any positive solution (x1(n),x2(n)) of system(1.2) with initial condition (1.4), one has milimn+infxi(n)limn+supxi(n)Mi,i=1,2.

Proof.

By applying Proposition 2.5, we see that to end the proof of Theorem 2.6 it is enough to show that under the conditions of Theorem 2.6limn+infxi(n)mi. From Proposition 2.5, For all ε>0, there exists a N1>0,  N1N,  For all n>N1, xi(n)Mi+ε. According to Lemma 2.4, from (2.13) and (2.14) we have limn+sups=0Ji(s)xi(n-s)=limn+supk=-nJi(n-k)xi(k)Mi,i=1,2. For above ɛ>0, according to (2.18), there exists a positive integer N2, such that, for all n>N2, s=0Ji(s)xi(n-s)Mi+ɛ,i=1,2. Thus, for all n>max{N1,N2}+σ, from the first equation of system(1.2), it follows that x1(n+1)x1(n)exp{r1(n)[K1l1+(M2+ε)-(M1+ε)]}x1(n)exp{r1lK1l1+(M2+ε)-r1u(M1+ε)}. It follows that, for nσ1(n), i=n-σ1(n)n-1x1(i+1)i=n-σ1(n)n-1x1(i)exp{r1lK1l1+(M2+ε)-r1u(M1+ε)}. Hence x1(n)x1(n-σ1(n))exp{r1lK1l1+(M2+ε)σ1l-r1u(M1+ε)σ1u}. In other words, x1(n-σ1(n))x1(n)exp{-r1lK1l1+(M2+ε)σ1l+r1u(M1+ε)σ1u}. From the first equation of system (1.2) and (2.23), for all n>max{N1,N2}+σ, it follows that x1(n+1)x1(n)exp{-r1lK1l1+(M2+ε)-r1uexp{-r1lK1l1+(M2+ε)σ1l+r1u(M1+ε)σ1u}x1(n)}. By applying Lemmas 2.1 and 2.2 to (2.24), it immediately follows thatlimn+infx1(n)r1lK1lr1u(1+(M2+ε))exp{r1lK1l1+(M2+ε)σ1l-r1u(M1+ε)σ1u}×exp{r1lK1l1+(M2+ε)-r1uexp{-r1lK1l1+(M2+ε)σ1l+r1u(M1+ε)σ1u}M1}. Setting ɛ0, it follows that limn+infx1(n)r1lK1lr1u(1+M2)exp{r1lK1l1+M2σ1l-r1uM1σ1u}×exp{r1lK1l1+M2-r1uexp{-r1lK1l1+M2σ1l+r1uM1σ1u}M1}. Similar to the above analysis, from the second equation of system (1.2), we have that limn+infx2(n)r2lK2lr2u(1+M1)exp{r2lK2l1+M1σ2l-r2uM2σ2u}×exp{r2lK2l1+M1-r2uexp{-r2lK2l1+M1σ2l+r2uM2σ2u}M2}. This completes the proof of Theorem 2.6.

ChenF. D.YouM. S.Permanence for an integrodifferential model of mutualismApplied Mathematics and Computation20071861303410.1016/j.amc.2006.07.085MR2316488ZBL1118.45006ChenF. D.Permanence in a discrete Lotka-Volterra competition model with deviating argumentsNonlinear Analysis: Real World Applications2008952150215510.1016/j.nonrwa.2007.07.001MR2441771ZBL1156.39300LiY. K.XuG. T.Positive periodic solutions for an integrodifferential model of mutualismApplied Mathematics Letters200114552553010.1016/S0893-9659(00)00188-9MR1832658ZBL0981.45002YangP.XuR.Global asymptotic stability of periodic solution in n-species cooperative system with time delaysJournal of Biomathematics1998136841846AgarwalR. P.Difference Equations and Inequalities: Theory, Method and Applications20002282ndNew York, NY, USAMarcel Dekkerxvi+971Monographs and Textbooks in Pure and Applied MathematicsMR1740241ChenF. D.Permanence and global attractivity of a discrete multispecies Lotka-Volterra competition predator-prey systemsApplied Mathematics and Computation2006182131210.1016/j.amc.2006.03.026MR2292013ChenX.ChenF. D.Stable periodic solution of a discrete periodic Lotka-Volterra competition system with a feedback controlApplied Mathematics and Computation200618121446145410.1016/j.amc.2006.02.039MR2270776ZBL1106.39003LiY. K.LuL. H.Positive periodic solutions of discrete n-species food-chain systemsApplied Mathematics and Computation2005167132434410.1016/j.amc.2004.06.082MR2170919MuroyaY.Persistence and global stability in Lotka-Volterra delay differential systemsApplied Mathematics Letters200417779580010.1016/j.aml.2004.06.009MR2072837ZBL1061.34057MuroyaY.Partial survival and extinction of species in discrete nonautonomous Lotka-Volterra systemsTokyo Journal of Mathematics200528118920010.3836/tjm/1244208288MR2149632ZBL1081.39014YangX. T.Uniform persistence and periodic solutions for a discrete predator-prey system with delaysJournal of Mathematical Analysis and Applications2006316116117710.1016/j.jmaa.2005.04.036MR2201755ZBL1107.39017ChenF. D.Permanence for the discrete mutualism model with time delaysMathematical and Computer Modelling2008473-443143510.1016/j.mcm.2007.02.023MR2378847ZBL1148.39017WangL.WangM. Q.Ordinary Difference Equation1991Xinjiang, ChinaXinjiang University PressTakeuchiY.Global Dynamical Properties of Lotka-Volterra Systems1996River Edge, NJ, USAWorld Scientificxii+302MR1440182