Discrete Dynamics in Nature and Society
Volume 2010 (2010), Article ID 931798, 7 pages
doi:10.1155/2010/931798
Research Article

Permanence of a Discrete Model of Mutualism with Infinite Deviating Arguments

School of Mathematics and Computer Science, Fujian Normal University, Fuzhou, Fujian 350007, China

Received 15 July 2009; Accepted 13 January 2010

Academic Editor: Binggen Zhang

Copyright © 2010 Xuepeng Li and Wensheng Yang. This 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.

Abstract

We propose a discrete model of mutualism with infinite deviating arguments, that is 𝑥 1 ( 𝑛 + 1 ) = 𝑥 1 ( 𝑛 ) e x p { 𝑟 1 ( 𝑛 ) [ ( 𝐾 1 ( 𝑛 ) + 𝛼 1 ( 𝑛 ) 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) ) / ( 1 + 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) ) 𝑥 1 ( 𝑛 𝜎 1 ( 𝑛 ) ) ] } , 𝑥 2 𝑥 ( 𝑛 + 1 ) = 2 ( 𝑛 ) e x p { 𝑟 2 ( 𝑛 ) [ ( 𝐾 2 ( 𝑛 ) + 𝛼 2 ( 𝑛 ) 𝑠 = 0 𝐽 1 ( 𝑠 ) 𝑥 1 ( 𝑛 𝑠 ) ) / ( 1 + 𝑠 = 0 𝐽 1 ( 𝑠 ) 𝑥 1 ( 𝑛 𝑠 ) ) 𝑥 2 ( 𝑛 𝜎 2 ( 𝑛 ) ) ] } . By some Lemmas, sufficient conditions are obtained for the permanence of the system.

1. Introduction

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

𝑑 𝑁 1 ( 𝑡 ) 𝑑 𝑡 = 𝑟 1 ( 𝑡 ) 𝑁 1 𝐾 ( 𝑡 ) 1 ( 𝑡 ) + 𝛼 1 ( 𝑡 ) 0 𝐽 2 ( 𝑠 ) 𝑁 2 ( 𝑡 𝑠 ) 𝑑 𝑠 1 + 0 𝐽 2 ( 𝑠 ) 𝑁 2 ( 𝑡 𝑠 ) 𝑑 𝑠 𝑁 1 𝑡 𝜎 1 , ( 𝑡 ) 𝑑 𝑁 2 ( 𝑡 ) 𝑑 𝑡 = 𝑟 2 ( 𝑡 ) 𝑁 2 ( 𝐾 𝑡 ) 2 ( 𝑡 ) + 𝛼 2 ( 𝑡 ) 0 𝐽 1 ( 𝑠 ) 𝑁 1 ( 𝑡 𝑠 ) 𝑑 𝑠 1 + 0 𝐽 1 ( 𝑠 ) 𝑁 1 ( 𝑡 𝑠 ) 𝑑 𝑠 𝑁 2 𝑡 𝜎 2 ( , 𝑡 ) ( 1 . 1 ) where 𝑟 𝑖 , 𝐾 𝑖 , 𝛼 𝑖 , and 𝜎 𝑖 , 𝑖 = 1 , 2 are continuous functions bounded above and below by positive constants: 𝑎 𝑖 > 𝐾 𝑖 , 𝑖 = 1 , 2 ; 𝐽 𝑖 𝐶 ( [ 0 , + ) , [ 0 , + ) ) and 0 𝐽 𝑖 ( 𝑠 ) 𝑑 𝑠 = 1 , 𝑖 = 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 [14] and the references cited therein.

However, many authors [512] 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:

𝑥 1 ( 𝑛 + 1 ) = 𝑥 1 𝑟 ( 𝑛 ) e x p 1 𝐾 ( 𝑛 ) 1 ( 𝑛 ) + 𝛼 1 ( 𝑛 ) 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) 1 + 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) 𝑥 1 𝑛 𝜎 1 , 𝑥 ( 𝑛 ) 2 ( 𝑛 + 1 ) = 𝑥 2 𝑟 ( 𝑛 ) e x p 2 𝐾 ( 𝑛 ) 2 ( 𝑛 ) + 𝛼 2 ( 𝑛 ) 𝑠 = 0 𝐽 1 ( 𝑠 ) 𝑥 1 ( 𝑛 𝑠 ) 1 + 𝑠 = 0 𝐽 1 ( 𝑠 ) 𝑥 1 ( 𝑛 𝑠 ) 𝑥 2 𝑛 𝜎 2 , ( 𝑛 ) ( 1 . 2 ) where 𝑥 𝑖 ( 𝑛 ) , 𝑖 = 1 , 2 is the density of mutualism species 𝑖 at the 𝑛 th generation. For { 𝑟 𝑖 ( 𝑛 ) } , { 𝐾 𝑖 ( 𝑛 ) } , { 𝛼 𝑖 ( 𝑛 ) } , { 𝐽 𝑖 ( 𝑛 ) } , and { 𝜎 𝑖 ( 𝑛 ) } , 𝑖 = 1 , 2 are bounded nonnegative sequences such that

0 < 𝑟 𝑙 𝑖 𝑟 𝑢 𝑖 , 0 < 𝛼 𝑙 𝑖 𝛼 𝑢 𝑖 , 0 < 𝐾 𝑙 𝑖 𝐾 𝑢 𝑖 , 0 < 𝜎 𝑙 𝑖 𝜎 𝑢 𝑖 , 𝑛 = 0 𝐽 𝑖 ( 𝑛 ) = 1 . ( 1 . 3 ) Here, for any bounded sequence { 𝑎 ( 𝑛 ) } , 𝑎 𝑢 = s u p 𝑛 𝑁 𝑎 ( 𝑛 ) , 𝑎 𝑙 = i n f 𝑛 𝑁 𝑎 ( 𝑛 ) .

Let 𝜎 = s u p 𝑛 { 𝜎 𝑖 ( 𝑛 ) , 𝑖 = 1 , 2 } , we consider (1.2) together with the following initial condition:

𝑥 𝑖 ( 𝜃 ) = 𝜑 𝑖 [ ] ( 𝜃 ) 0 , 𝜃 𝑁 𝜏 , 0 = { 𝜏 , 𝜏 + 1 , , 0 } , 𝜑 𝑖 ( 0 ) > 0 . ( 1 . 4 )

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

𝑥 𝑖 ( 𝑛 ) > 0 , f o r 𝑛 𝑍 , 𝑖 = 1 , 2 . ( 1 . 5 )

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 2 4 1 ].

Lemma. Let 𝑘 𝑁 + 𝑘 0 = { 𝑘 0 , 𝑘 0 + 1 , , 𝑘 0 + 𝑙 , } , 𝑟 0 . For any fixed 𝑘 , 𝑔 ( 𝑘 , 𝑟 ) is a non-decreasing function with respect to 𝑟 , and for 𝑘 𝑘 0 , following inequalities hold: 𝑦 ( 𝑘 + 1 ) 𝑔 ( 𝑘 , 𝑦 ( 𝑘 ) ) , 𝑢 ( 𝑘 + 1 ) 𝑔 ( 𝑘 , 𝑢 ( 𝑘 ) ) . If 𝑦 ( 𝑘 0 ) 𝑢 ( 𝑘 0 ) , then 𝑦 ( 𝑘 ) 𝑢 ( 𝑘 ) for all 𝑘 𝑘 0 .

Now let us consider the following single species discrete model:

𝑁 ( 𝑘 + 1 ) = 𝑁 ( 𝑘 ) e x p { 𝑎 ( 𝑘 ) 𝑏 ( 𝑘 ) 𝑁 ( 𝑘 ) } , ( 2 . 1 ) where { 𝑎 ( 𝑘 ) } and { 𝑏 ( 𝑘 ) } are strictly positive sequences of real numbers defined for 𝑘 𝑁 = { 0 , 1 , 2 , } and 0 < 𝑎 𝑙 𝑎 𝑢 , 0 < 𝑏 𝑙 𝑏 𝑢 . Similar to the proof of Propositions 1 and 3 in [6], we can obtain the following.

Lemma. Any solution of system (2.1) with initial condition 𝑁 ( 0 ) > 0 satisfies 𝑚 l i m 𝑘 + i n f 𝑁 ( 𝑘 ) l i m 𝑘 + s u p 𝑁 ( 𝑘 ) 𝑀 , ( 2 . 2 ) where 1 𝑀 = 𝑏 𝑙 e x p { 𝑎 𝑢 𝑎 1 } , 𝑚 = 𝑙 𝑏 𝑢 𝑎 e x p 𝑙 𝑏 𝑢 𝑀 . ( 2 . 3 )

Lemma 2.3 (see [14]). Let 𝑥 ( 𝑛 ) and 𝑏 ( 𝑛 ) be nonnegative sequences defined on 𝑁 , and 𝑐 0 is a constant. If 𝑥 ( 𝑛 ) 𝑐 + 𝑛 1 𝑠 = 0 𝑏 ( 𝑠 ) 𝑥 ( 𝑠 ) , f o r 𝑛 𝑁 , ( 2 . 4 ) then 𝑥 ( 𝑛 ) 𝑐 𝑛 1 𝑠 = 0 [ ] 1 + 𝑏 ( 𝑠 ) , f o r 𝑛 𝑁 . ( 2 . 5 )

Lemma 2.4 (see [2]). Let 𝑥 𝑍 𝑅 be a nonnegative bounded sequences, and let 𝐻 𝑁 𝑅 be a nonnegative sequence such that 𝑛 = 0 𝐽 𝑖 ( 𝑛 ) = 1 . Then l i m 𝑛 + i n f 𝑥 ( 𝑛 ) l i m 𝑛 + i n f 𝑛 𝑠 = 𝐻 ( 𝑛 𝑠 ) 𝑥 ( 𝑠 ) l i m 𝑛 + s u p 𝑛 𝑠 = 𝐻 ( 𝑛 𝑠 ) 𝑥 ( 𝑠 ) l i m 𝑛 + s u p 𝑥 ( 𝑛 ) . ( 2 . 6 )

Proposition. Let ( 𝑥 1 ( 𝑛 ) , 𝑥 2 ( 𝑛 ) ) be any positive solution of system (1.2), then l i m 𝑛 + s u p 𝑥 𝑖 ( 𝑛 ) 𝑀 𝑖 , 𝑖 = 1 , 2 , ( 2 . 7 ) where 𝑀 𝑖 = e x p 2 𝑟 𝑢 𝑖 𝐾 𝑢 𝑖 + 𝛼 𝑢 𝑖 , 𝑖 = 1 , 2 . ( 2 . 8 )

Proof. Let ( 𝑥 1 ( 𝑛 ) , 𝑥 2 ( 𝑛 ) ) be any positive solution of system (1.2), then from the first equation of system (1.2) we have 𝑥 1 ( 𝑛 + 1 ) 𝑥 1 𝑟 ( 𝑛 ) e x p 1 𝐾 ( 𝑛 ) 1 ( 𝑛 ) + 𝛼 1 ( 𝑛 ) 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) 1 + 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) = 𝑥 1 𝑟 ( 𝑛 ) e x p 1 𝐾 ( 𝑛 ) 1 ( 𝑛 ) 1 + 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 + 𝛼 ( 𝑛 𝑠 ) 1 ( 𝑛 ) 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) 1 + 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) 𝑥 1 𝑟 ( 𝑛 ) e x p 1 𝐾 ( 𝑛 ) 1 ( 𝑛 ) 1 + 𝛼 1 ( 𝑛 ) 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) 𝑠 = 0 𝐽 2 ( 𝑠 ) 𝑥 2 ( 𝑛 𝑠 ) = 𝑥 1 ( 𝑟 𝑛 ) e x p 1 ( 𝐾 𝑛 ) 1 ( 𝑛 ) + 𝛼 1 ( 𝑛 ) 𝑥 1 𝑟 ( 𝑛 ) e x p 𝑢 1 𝐾 𝑢 1 + 𝛼 𝑢 1 . ( 2 . 9 ) Let 𝑥 1 ( 𝑛 ) = e x p { 𝑢 1 ( 𝑛 ) } , then 𝑢 1 ( 𝑛 + 1 ) 𝑢 1 ( 𝑛 ) + 𝑟 𝑢 1 𝐾 𝑢 1 + 𝛼 𝑢 1 = 𝑟 𝑢 1 𝐾 𝑢 1 + 𝛼 𝑢 1 + 𝑛 𝑠 = 0 𝑏 ( 𝑠 ) 𝑥 ( 𝑠 ) , ( 2 . 1 0 ) where 𝑏 ( 𝑠 ) = 0 , 0 𝑠 𝑛 1 , 1 , 𝑠 = 𝑛 . ( 2 . 1 1 ) When 𝑢 1 ( 𝑛 ) is nonnegative sequence, by applying Lemma 2.3, it immediately follows that 𝑢 1 ( 𝑛 + 1 ) 2 𝑟 𝑢 1 𝐾 𝑢 1 + 𝛼 𝑢 1 . ( 2 . 1 2 ) When 𝑢 1 ( 𝑛 ) is negative sequence, (2.12) also holds. From (2.12), we have l i m 𝑛 + s u p 𝑥 1 ( 𝑛 ) e x p 2 𝑟 𝑢 1 𝐾 𝑢 1 + 𝛼 𝑢 1 = 𝑀 1 . ( 2 . 1 3 ) By using the second equation of system (1.2), similar to the above analysis, we can obtain l i m 𝑛 + s u p 𝑥 2 ( 𝑛 ) e x p 2 𝑟 𝑢 2 𝐾 𝑢 2 + 𝛼 𝑢 2 = 𝑀 2 . ( 2 . 1 4 ) 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 𝑚 𝑖 , 𝑀 𝑖 , 𝑖 = 1 , 2 which are independent of the solutions of system (1.2) such that, for any positive solution ( 𝑥 1 ( 𝑛 ) , 𝑥 2 ( 𝑛 ) ) of system(1.2) with initial condition (1.4), one has 𝑚 𝑖 l i m 𝑛 + i n f 𝑥 𝑖 ( 𝑛 ) l i m 𝑛 + s u p 𝑥 𝑖 ( 𝑛 ) 𝑀 𝑖 , 𝑖 = 1 , 2 . ( 2 . 1 5 )

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.6 l i m 𝑛 + i n f 𝑥 𝑖 ( 𝑛 ) 𝑚 𝑖 . ( 2 . 1 6 ) From Proposition 2.5, For all 𝜀 > 0 , there exists a 𝑁 1 > 0 , 𝑁 1 𝑁 , For all 𝑛 > 𝑁 1 , 𝑥 𝑖 ( 𝑛 ) 𝑀 𝑖 + 𝜀 . ( 2 . 1 7 ) According to Lemma 2.4, from (2.13) and (2.14) we have l i m 𝑛 + s u p 𝑠 = 0 𝐽 𝑖 ( 𝑠 ) 𝑥 𝑖 ( 𝑛 𝑠 ) = l i m 𝑛 + s u p 𝑛 𝑘 = 𝐽 𝑖 ( 𝑛 𝑘 ) 𝑥 𝑖 ( 𝑘 ) 𝑀 𝑖 , 𝑖 = 1 , 2 . ( 2 . 1 8 ) For above 𝜀 > 0 , according to (2.18), there exists a positive integer 𝑁 2 , such that, for all 𝑛 > 𝑁 2 , 𝑠 = 0 𝐽 𝑖 ( 𝑠 ) 𝑥 𝑖 ( 𝑛 𝑠 ) 𝑀 𝑖 + 𝜀 , 𝑖 = 1 , 2 . ( 2 . 1 9 ) Thus, for all 𝑛 > 𝑚 𝑎 𝑥 { 𝑁 1 , 𝑁 2 } + 𝜎 , from the first equation of system(1.2), it follows that 𝑥 1 ( 𝑛 + 1 ) 𝑥 1 𝑟 ( 𝑛 ) e x p 1 𝐾 ( 𝑛 ) 𝑙 1 𝑀 1 + 2 𝑀 + 𝜀 1 + 𝜀 𝑥 1 𝑟 ( 𝑛 ) e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 + 𝜀 𝑟 𝑢 1 𝑀 1 . + 𝜀 ( 2 . 2 0 ) It follows that, for 𝑛 𝜎 1 ( 𝑛 ) , 𝑛 1 𝑖 = 𝑛 𝜎 1 ( 𝑛 ) 𝑥 1 ( 𝑖 + 1 ) 𝑛 1 𝑖 = 𝑛 𝜎 1 ( 𝑛 ) 𝑥 1 𝑟 ( 𝑖 ) e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 + 𝜀 𝑟 𝑢 1 𝑀 1 . + 𝜀 ( 2 . 2 1 ) Hence 𝑥 1 ( 𝑛 ) 𝑥 1 𝑛 𝜎 1 𝑟 ( 𝑛 ) e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 𝜎 + 𝜀 𝑙 1 𝑟 𝑢 1 𝑀 1 𝜎 + 𝜀 𝑢 1 . ( 2 . 2 2 ) In other words, 𝑥 1 𝑛 𝜎 1 ( 𝑛 ) 𝑥 1 𝑟 ( 𝑛 ) e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 𝜎 + 𝜀 𝑙 1 + 𝑟 𝑢 1 𝑀 1 𝜎 + 𝜀 𝑢 1 . ( 2 . 2 3 ) From the first equation of system (1.2) and (2.23), for all 𝑛 > 𝑚 𝑎 𝑥 { 𝑁 1 , 𝑁 2 } + 𝜎 , it follows that 𝑥 1 ( 𝑛 + 1 ) 𝑥 1 𝑟 ( 𝑛 ) e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 + 𝜀 𝑟 𝑢 1 𝑟 e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 𝜎 + 𝜀 𝑙 1 + 𝑟 𝑢 1 𝑀 1 𝜎 + 𝜀 𝑢 1 𝑥 1 . ( 𝑛 ) ( 2 . 2 4 ) By applying Lemmas 2.1 and 2.2 to (2.24), it immediately follows that l i m 𝑛 + i n f 𝑥 1 𝑟 ( 𝑛 ) 𝑙 1 𝐾 𝑙 1 𝑟 𝑢 1 𝑀 1 + 2 𝑟 + 𝜀 e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 𝜎 + 𝜀 𝑙 1 𝑟 𝑢 1 𝑀 1 𝜎 + 𝜀 𝑢 1 𝑟 × e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 + 𝜀 𝑟 𝑢 1 𝑟 e x p 𝑙 1 𝐾 𝑙 1 𝑀 1 + 2 𝜎 + 𝜀 𝑙 1 + 𝑟 𝑢 1 𝑀 1 𝜎 + 𝜀 𝑢 1 𝑀 1 . ( 2 . 2 5 ) Setting 𝜀 0 , it follows that l i m 𝑛 + i n f 𝑥 1 𝑟 ( 𝑛 ) 𝑙 1 𝐾 𝑙 1 𝑟 𝑢 1 1 + 𝑀 2 𝑟 e x p 𝑙 1 𝐾 𝑙 1 1 + 𝑀 2 𝜎 𝑙 1 𝑟 𝑢 1 𝑀 1 𝜎 𝑢 1 𝑟 × e x p 𝑙 1 𝐾 𝑙 1 1 + 𝑀 2 𝑟 𝑢 1 𝑟 e x p 𝑙 1 𝐾 𝑙 1 1 + 𝑀 2 𝜎 𝑙 1 + 𝑟 𝑢 1 𝑀 1 𝜎 𝑢 1 𝑀 1 . ( 2 . 2 6 ) Similar to the above analysis, from the second equation of system (1.2), we have that l i m 𝑛 + i n f 𝑥 2 𝑟 ( 𝑛 ) 𝑙 2 𝐾 𝑙 2 𝑟 𝑢 2 1 + 𝑀 1 𝑟 e x p 𝑙 2 𝐾 𝑙 2 1 + 𝑀 1 𝜎 𝑙 2 𝑟 𝑢 2 𝑀 2 𝜎 𝑢 2 𝑟 × e x p 𝑙 2 𝐾 𝑙 2 1 + 𝑀 1 𝑟 𝑢 2 𝑟 e x p 𝑙 2 𝐾 𝑙 2 1 + 𝑀 1 𝜎 𝑙 2 + 𝑟 𝑢 2 𝑀 2 𝜎 𝑢 2 𝑀 2 . ( 2 . 2 7 ) This completes the proof of Theorem 2.6.

References

  1. F. D. Chen and M. S. You, “Permanence for an integrodifferential model of mutualism,” Applied Mathematics and Computation, vol. 186, no. 1, pp. 30–34, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  2. F. D. Chen, “Permanence in a discrete Lotka-Volterra competition model with deviating arguments,” Nonlinear Analysis: Real World Applications, vol. 9, no. 5, pp. 2150–2155, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  3. Y. K. Li and G. T. Xu, “Positive periodic solutions for an integrodifferential model of mutualism,” Applied Mathematics Letters, vol. 14, no. 5, pp. 525–530, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  4. P. Yang and R. Xu, “Global asymptotic stability of periodic solution in n-species cooperative system with time delays,” Journal of Biomathematics, vol. 13, no. 6, pp. 841–846, 1998.
  5. R. P. Agarwal, Difference Equations and Inequalities: Theory, Method and Applications, vol. 228 of Monographs and Textbooks in Pure and Applied Mathematics, Marcel Dekker, New York, NY, USA, 2nd edition, 2000. View at MathSciNet
  6. F. D. Chen, “Permanence and global attractivity of a discrete multispecies Lotka-Volterra competition predator-prey systems,” Applied Mathematics and Computation, vol. 182, no. 1, pp. 3–12, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  7. X. Chen and F. D. Chen, “Stable periodic solution of a discrete periodic Lotka-Volterra competition system with a feedback control,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1446–1454, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  8. Y. K. Li and L. H. Lu, “Positive periodic solutions of discrete n-species food-chain systems,” Applied Mathematics and Computation, vol. 167, no. 1, pp. 324–344, 2005. View at Publisher · View at Google Scholar · View at MathSciNet
  9. Y. Muroya, “Persistence and global stability in Lotka-Volterra delay differential systems,” Applied Mathematics Letters, vol. 17, no. 7, pp. 795–800, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. Y. Muroya, “Partial survival and extinction of species in discrete nonautonomous Lotka-Volterra systems,” Tokyo Journal of Mathematics, vol. 28, no. 1, pp. 189–200, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  11. X. T. Yang, “Uniform persistence and periodic solutions for a discrete predator-prey system with delays,” Journal of Mathematical Analysis and Applications, vol. 316, no. 1, pp. 161–177, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  12. F. D. Chen, “Permanence for the discrete mutualism model with time delays,” Mathematical and Computer Modelling, vol. 47, no. 3-4, pp. 431–435, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  13. L. Wang and M. Q. Wang, Ordinary Difference Equation, Xinjiang University Press, Xinjiang, China, 1991.
  14. Y. Takeuchi, Global Dynamical Properties of Lotka-Volterra Systems, World Scientific, River Edge, NJ, USA, 1996. View at MathSciNet