# Exponential Extinction of Nicholson's Blowflies System with Nonlinear Density-Dependent Mortality Terms

Authors: Wentao Wang
Journal: Abstract and Applied Analysis
 (2012)

---

## Abstract

This paper presents a new generalized Nicholson’s blowflies system with patch structure and nonlinear density-dependent mortality terms. Under appropriate conditions, we establish some criteria to guarantee the exponential extinction of this system. Moreover, we give two examples and numerical simulations to demonstrate our main results.


---

## Body

## 1. Introduction

To describe the population of the Australian sheep blowfly and agree well with the experimental date of Nicholson [1], Gurney et al. [2] proposed the following Nicholson's blowflies equation:
(1.1)N′(t)=-δN(t)+pN(t-τ)e-aN(t-τ).
Here, N(t) is the size of the population at time t,p is the maximum per capita daily egg production, (1/a) is the size at which the population reproduces at its maximum rate, δ is the per capita daily adult death rate, and τ is the generation time. There have been a large number of results on this model and its modifications (see, e.g., [3–8]). Recently, Berezansky et al. [9] pointed out that a new study indicates that a linear model of density-dependent mortality will be most accurate for populations at low densities and marine ecologists are currently in the process of constructing new fishery models with nonlinear density-dependent mortality rates. Consequently Berezansky et al. [9] presented the following Nicholson’s blowflies model with a nonlinear density-dependent mortality term
(1.2)N′(t)=-D(N(t))+PN(t-τ)e-aN(t-τ),
where P is a positive constant and function D might have one of the following forms: D(N)=aN/(N+b) or D(N)=a-be-N with positive constants a,b>0.Wang [10] studied the existence of positive periodic solutions for the model (1.2) with D(N)=a-be-N.  Hou et al. [11] investigated the permanence and periodic solutions for the model (1.2) with D(N)=aN/(N+b). Furthermore, Liu and Gong [12] considered the permanence for a Nicholson-type delay systems with nonlinear density-dependent mortality terms as follows:
(1.3)N1′(t)=-D11(t,N1(t))+D12(t,N2(t))+c1(t)N1(t-τ1(t))e-γ1(t)N1(t-τ1(t))N2′(t)=-D22(t,N2(t))+D21(t,N1(t))+c2(t)N2(t-τ2(t))e-γ2(t)N2(t-τ2(t)),
where
(1.4)Dij(t,N)=aij(t)Nbij(t)+NorDij(t,N)=aij(t)-bij(t)e-N,aij,bij,ci,γi:R→(0,+∞) are all continuous functions bounded above and below by positive constants, and τj(t)≥0 are bounded continuous functions, ri=supt∈Rτi(t)>0, and i,j=1,2.On the other hand, since the biological species compete and cooperate with each other in real world, the growth models given by patch structure systems of delay differential equation have been provided by several authors to analyze the dynamics of multiple species (see, e.g., [13–16] and the reference therein). Moreover, the extinction phenomenon often appears in the biology, economy, and physics field and the main focus of Nicholson's blowflies model is on the scalar equation and results on patch structure of this model are gained rarely [14, 16], so it is worth studying the extinction of Nicholson's blowflies system with patch structure and nonlinear density-dependent mortality terms. Motivated by the above discussion, we shall derive the conditions to guarantee the extinction of the following Nicholson-type delay system with patch structure and nonlinear density-dependent mortality terms:
(1.5)Ni′(t)=-Dii(t,Ni(t))+∑j=1,j≠in‍Dij(t,Nj(t))+∑j=1l‍cij(t)Ni(t-τij(t))e-γij(t)Ni(t-τij(t)),
where
(1.6)Dij(t,N)=aij(t)Nbij(t)+NorDij(t,N)=aij(t)-bij(t)e-N,aij,bij,cik,γik:R→(0,+∞) are all continuous functions bounded above and below by positive constants, and τik(t)≥0 are bounded continuous functions, ri=max1≤j≤l{supt∈Rτij(t)}>0, and i,j=1,2,…,n,k=1,2,…,l. Furthermore, in the case Dij(t,N)=aij(t)-bij(t)e-N, to guarantee the meaning of mortality terms we assume that aij(t)>bij(t) for t∈R and i,j=1,2,…,n. The main purpose of this paper is to establish the conditions ensuring the exponential extinction of system (1.5).For convenience, we introduce some notations. Throughout this paper, given a bounded continuous functiong defined on R, let g+ and g- be defined as
(1.7)g-=inft∈Rg(t),g+=supt∈Rg(t).LetRn(R+n) be the set of all (nonnegative) real vectors, we will use x=(x1,…,xn)T∈Rn to denote a column vector, in which the symbol (T) denotes the transpose of a vector. We let |x| denote the absolute-value vector given by |x|=(|x1|,…,|xn|)T and define ||x||=max1≤i≤n|xi|. Denote C=∏i=1n‍C([-ri,0],R)andC+=∏i=1n‍C([-ri,0],R+) as Banach spaces equipped with the supremum norm defined by ||φ||=sup-ri≤t≤0max1≤i≤n|φi(t)| for all φ(t)=(φ1(t),…,φn(t))T∈C (or ∈C+). If  xi(t) is defined on [t0-ri,ν) with t0,ν∈R and i=1,…,n, then we define xt∈C as xt=(xt1,…xtn)T where xti(θ)=xi(t+θ) for all θ∈[-ri,0] and i=1,…,n.The initial conditions associated with system (1.5) are of the form:
(1.8)Nt0=φ,φ=(φ1,…,φn)T∈C+,φi(0)>0,i=1,…,n.
We write Nt(t0,φ)(N(t;t0,φ)) for a solution of the initial value problem (1.5) and (1.8). Also, let [t0,η(φ)) be the maximal right-interval of existence of  Nt(t0,φ).Definition 1.1.
The system (1.5) with initial conditions (1.8) is said to be exponentially extinct, if there are positive constants M and κ such that |Ni(t;t0,φ)|≤Me-κ(t-t0),i=1,2…,n.  Denote it as Ni(t;t0,φ)=O(e-κ(t-t0)),i=1,2,…,n.
The remaining part of this paper is organized as follows. In Sections2 and 3, we shall derive some sufficient conditions for checking the extinction of system (1.5). In Section 4, we shall give two examples and numerical simulations to illustrate our results obtained in the previous sections.


## 2. Extinction of Nicholson's Blowflies System withDij(t,N)=aij(t)N/(bij(t)+N)(i,j=1,2,…,n)

Theorem 2.1.
Suppose that there exists positive constantK1 such that
(2.1)aii-bii++K1>∑j=1,j≠in‍aij+bij-+∑j=1l‍cij+γij-eK1,i=1,2,…,n.
Let
(2.2)E1={φ∣φ∈C+,φ(0)>0,0≤φi(t)<K1,∀t∈[-ri,0],i=1,2,…,n}.
Moreover, assume N(t;t0,φ) is the solution of (1.5) with φ∈E1 and Dij(t,N)=(aij(t)N/(bij(t)+N))(i,j=1,2,…,n). Then,
(2.3)0≤Ni(t;t0,φ)<K1,∀t∈[t0,η(φ)),i=1,2,…,n,η(φ)=+∞.Proof.
SetN(t)=N(t;t0,φ)forallt∈[t0,η(φ)). In view of φ∈C+, using Theorem 5.2.1 in [17, p. 81], we have Nt(t0,φ)∈C+ for all t∈[t0,η(φ)). Assume, by way of contradiction, that (2.3) does not hold. Then, there  exist t1∈[t0,η(φ)) and i∈{1,2,…,n} such that
(2.4)Ni(t1)=K1,0≤Nj(t)<K1∀t∈[t0-rj,t1),j=1,2,…,n.
Calculating the derivative of Ni(t), together with (2.1) and the fact that supu≥0ue-u=1/e and a(t)N/(b(t)+N)≤a(t)N/b(t) for all t∈R,N≥0, (1.5) and (2.4) imply that
(2.5)0≤Ni′(t1)=-Dii(t1,Ni(t1))+∑j=1,j≠in‍Dij(t1,Nj(t1))+∑j=1l‍cij(t1)Ni(t1-τij(t1))e-γij(t1)Ni(t1-τij(t1))≤-aii(t1)Ni(t1)bii(t1)+Ni(t1)+∑j=1,j≠in‍aij(t1)Nj(t1)bij(t1)+∑j=1l‍cij(t1)γij(t1)1e≤(-aii-bii++K1+∑j=1,j≠in‍aij+bij-+∑j=1l‍cij+γij-eK1)K1<0,
which is a contradiction and implies that (2.3) holds. From Theorem 2.3.1 in [18], we easily obtain η(φ)=+∞. This ends the proof of Theorem 2.1.Theorem 2.2.
Suppose that there exists positive constantK1 satisfying (2.1) and
(2.6)aii-bii++K1>∑j=1,j≠in‍aij+bij-+∑j=1l‍cij+,i=1,2,…,n.
Then the solution N(t;t0,φ) of (1.5) with φ∈E1 and Dij(t,N)=(aij(t)N/(bij(t)+N))(i,j=1,2,…,n) is exponentially extinct as t→+∞.Proof.
Define continuous functionsΓi(ω) by setting
(2.7)Γi(ω)=ω-aii-bii++K1+∑j=1,j≠in‍aij+bij-+∑j=1l‍cij+eωri,i=1,2,…,n.
Then, from (2.6), we obtain
(2.8)Γi(0)=-aii-bii++K1+∑j=1,j≠in‍aij+bij-+∑j=1l‍cij+<0,i=1,2,…,n.
The continuity of Γi(ω) implies that there exists λ>0 such that
(2.9)Γi(λ)=λ-aii-bii++K1+∑j=1,j≠in‍aij+bij-+∑j=1l‍cij+eλri<0,i=1,2,…,n.
Let
(2.10)yi(t)=Ni(t)eλ(t-t0),i=1,2,…,n.
Calculating the derivative of y(t) along the solution N(t) of system (1.5) with φ∈E1, we have
(2.11)yi′(t)=λyi(t)+eλ(t-t0)Ni′(t)=λyi(t)-aii(t)yi(t)bii(t)+Ni(t)+∑j=1,j≠in‍aij(t)yj(t)bij(t)+Nj(t)+∑j=1l‍cij(t)eλτij(t)yi(t-τij(t))e-γij(t)Ni(t-τij(t)),i=1,2,…,n.
Let M1 denote an arbitrary positive number and set
(2.12)M1>yi(t),∀t∈[t0-ri,t0],i=1,2,…,n.
We claim that
(2.13)yi(t)<M1,∀t∈[t0,+∞),i=1,2,…,n.
If this is not valid, there must exist t2∈(t0,+∞) and i∈{1,2,…,n} such that
(2.14)yi(t2)=M1,yj(t)<M1,∀t<t2,j=1,2,…,n.
Then, from (2.3) and (2.11), we have
(2.15)0≤yi′(t2)=λyi(t2)-aii(t2)yi(t2)bii(t2)+Ni(t2)+∑j=1,j≠in‍aij(t2)yj(t2)bij(t2)+Nj(t2)+∑j=1l‍cij(t2)eλτij(t2)yi(t2-τij(t2))e-γij(t2)Ni(t2-τij(t2))≤λM1-aii(t2)M1bii(t2)+K1+∑j=1,j≠in‍aij(t2)M1bij(t2)+∑j=1l‍cij(t2)eλriM1≤(λ-aii-bii++K1+∑j=1,j≠in‍aij+bij-+∑j=1l‍cij+eλri)M1<0.
This contradiction implies that (2.13) holds. Thus,
(2.16)Ni(t)=yi(t)e-λ(t-t0)≤M1e-λ(t-t0)∀t∈[t0-ri,+∞),i=1,2,…,n.
This completes the proof.


## 3. Extinction of Nicholson's Blowflies System withDij(t,N)=aij(t)-bij(t)e-N(i,j=1,2,…,n)

Theorem 3.1.
Suppose that there exists positive constantK2 such that
(3.1)aii->∑j=1,j≠in‍aij++(bii+-∑j=1,j≠in‍bij-)e-K2+∑j=1l‍cij+γij-e,i=1,2,…,n,(3.2)-aii(t)+bii(t)+∑j=1,j≠in‍(aij(t)-bij(t))≥0,i=1,2,…,n.
Let
(3.3)E2={φ∣φ∈C+,φ(0)>0,0≤φi(t)<K2,∀t∈[-ri,0],i=1,2,…,n}.
Moreover, assume N(t;t0,φ) is the solution of (1.5) with φ∈E2 and Dij(t,N)=aij(t)-bij(t)e-N(i,j=1,2,…,n). Then,
(3.4)0≤Ni(t;t0,φ)<K2,∀t∈[t0,η(φ)),i=1,2,…,n,(3.5)η(φ)=+∞.Proof.
SetN(t)=N(t;t0,φ)forallt∈[t0,η(φ)). Rewrite the system (1.5) as(3.6)N′(t)=f(t,Nt),
where f(t,ϕ)=(f1(t,ϕ),f2(t,ϕ),…,fn(t,ϕ))T and
(3.7)fi(t,ϕ)=-aii(t)+bii(t)e-ϕi(0)+∑j=1,j≠in‍(aij(t)-bij(t)e-ϕj(0))+∑j=1l‍cij(t)ϕi(-τij(t))e-γij(t)ϕi(-τij(t)),i=1,2,…,n,ϕ∈C.
In view of (3.2), whenever ϕ∈C satisfies ϕ≥0,ϕi(0)=0 for some i and t∈R, then
(3.8)fi(t,ϕ)=-aii(t)+bii(t)+∑j=1,j≠in‍(aij(t)-bij(t)e-ϕj(0))+∑j=1l‍cij(t)ϕi(-τij(t))e-γij(t)ϕi(-τij(t))≥-aii(t)+bii(t)+∑j=1,j≠in‍(aij(t)-bij(t))≥0.
Thus, using Theorem 5.2.1 in [17, p. 81], we have Nt(t0,φ)∈C+ for all t∈[t0,η(φ)) and φ∈E2⊂C+. Assume, by way of contradiction, that (3.4) does not hold. Then, there exist t3∈[t0,η(φ)) and i∈{1,2,…,n} such that
(3.9)Ni(t3)=K2,0≤Nj(t)<K2∀t∈[t0-rj,t3),j=1,2,…,n.
Calculating the derivative of Ni(t), together with (3.1) and the fact that supu≥0ue-u=1/e, (1.5) and (3.9) imply that
(3.10)0≤Ni′(t3)=-Dii(t3,Ni(t3))+∑j=1,j≠in‍Dij(t3,Nj(t3))+∑j=1l‍cij(t3)Ni(t3-τij(t3))×e-γij(t3)Ni(t3-τij(t3))≤-aii(t3)+bii(t3)e-K2+∑j=1,j≠in‍(aij(t3)-bij(t3)e-K2)+∑j=1l‍cij(t3)γij(t3)1e≤-aii-+∑j=1,j≠in‍aij++(bii+-∑j=1,j≠in‍bij-)e-K2+∑j=1l‍cij+γij-e<0,
which is a contradiction and implies that (3.4) holds. From Theorem 2.3.1 in [18], we easily obtain η(φ)=+∞. This ends the proof of Theorem 3.1.Theorem 3.2.
Let (3.1) and (3.2) hold. Moreover, suppose that there exist two positive constants λ~ and  M~ such that
(3.11)-aii(t)+bii(t)+∑j=1,j≠in‍(aij(t)-bij(t))≤M~e-λ~(t-t0),t∈R,i=1,2,…,n,(3.12)bii->1+K22bii++∑j=1,j≠in‍bij++∑j=1l‍cij+,i=1,2,…,n.
Then the solution N(t;t0,φ) of (1.5) with φ∈E2 and  Dij(t,N)=aij(t)-bij(t)e-N(i,j=1,2,…,n), is exponentially extinct as t→+∞.Proof.
Define continuous functionsΓi(ω) by setting
(3.13)Γi(ω)=ω-bii-+1+K22bii++∑j=1,j≠in‍bij++∑j=1l‍cij+eωri,i=1,2,…,n.
Then, from (3.12), we obtain
(3.14)Γi(0)=-bii-+1+K22bii++∑j=1,j≠in‍bij++∑j=1l‍cij+<0,i=1,2,…,n.
The continuity of Γi(ω) implies that there exists 0<μ<λ~ such that
(3.15)Γi(μ)=μ-bii-+1+K22bii++∑j=1,j≠in‍bij++∑j=1l‍cij+eμri<0,i=1,2,…,n.
Let
(3.16)xi(t)=Ni(t)eμ(t-t0),i=1,2,…,n.
Calculating the derivative of x(t) along the solution N(t) of system (1.5) with φ∈E2, in view of (3.4) and (3.11), we have
(3.17)xi′(t)=μxi(t)+eμ(t-t0)Ni′(t)=μxi(t)+eμ(t-t0)[-aii(t)+bii(t)e-Ni(t)+∑j=1,j≠in‍(aij(t)-bij(t)e-Nj(t))]+∑j=1l‍cij(t)eμτij(t)xi(t-τij(t))e-γij(t)Ni(t-τij(t))≤μxi(t)+eμ(t-t0)[∑j=1,j≠in‍(aij(t)-bij(t)(1-Nj(t)))-aii(t)+bii(t)(1-Ni(t)+12Ni2(t))+∑j=1,j≠in‍(aij(t)-bij(t)(1-Nj(t)))]+∑j=1l‍cij+eμrixi(t-τij(t))=μxi(t)+eμ(t-t0)[-aii(t)+bii(t)+∑j=1,j≠in‍(aij(t)-bij(t))]-bii(t)xi(t)+12bii(t)Ni(t)xi(t)+∑j=1,j≠in‍bij(t)xj(t)+∑j=1l‍cij+eμrixi(t-τij(t))≤μxi(t)+M~e(μ-λ~)(t-t0)-bii-xi(t)+K22bii+xi(t)+∑j=1,j≠in‍bij+xj(t)+∑j=1l‍cij+eμrixi(t-τij(t)).
Let M2 denote an arbitrary positive number and set
(3.18)M2>max{xi(t),M~}∀t∈[t0-ri,t0],i=1,2,…,n.
We claim that
(3.19)xi(t)<M2,∀t∈[t0,+∞),i=1,2,…,n.
If this is not valid, there must exist t4∈(t0,+∞) and i∈{1,2,…,n} such that
(3.20)xi(t4)=M2,xj(t)<M2,∀t<t4,j=1,2,…,n.
Then, from (3.15) and (3.17), we have
(3.21)0≤xi′(t4)≤μxi(t4)+M~e(μ-λ~)(t4-t0)-bii-xi(t4)+K22bii+xi(t4)+∑j=1,j≠in‍bij+xj(t4)+∑j=1l‍cij+eμrixi(t4-τij(t4))≤[μ+1-bii-+K22bii++∑j=1,j≠in‍bij++∑j=1l‍cij+eμri]M2<0.
This contradiction implies that (3.19) holds. Thus,
(3.22)Ni(t)=xi(t)e-μ(t-t0)≤M2e-μ(t-t0)∀t∈[t0-ri,+∞),i=1,2,…,n.
This completes the proof.


## 4. Numerical Examples

In this section, we give two examples and numerical simulations to demonstrate the results obtained in previous sections.Example 4.1.
Consider the following Nicholson's blowflies system with patch structure and nonlinear density-dependent mortality terms:(4.1)N1′(t)=-(25+|cos3t|)N1(t)5+|sin2t|+N1(t)+(1+|sin2t|)N2(t)3+|cos3t|+N2(t)+(1+|cos2t|)N3(t)3+|sin3t|+N3(t)+14(1+cos2t)N1(t-2|sint|)e-4N1(t-2|sint|)+14(1+sin2t)N1(t-2|cost|)e-4N1(t-2|cost|)N2′(t)=-(25+|sin3t|)N2(t)5+|cos2t|+N2(t)+(1+|cos2t|)N1(t)3+|sin3t|+N1(t)+(1+|sin2t|)N3(t)3+|cos3t|+N3(t)+14(1+sin2t)N2(t-2|cost|)e-4N2(t-2|cost|)+14(1+cos2t)N2(t-2|sint|)e-4N2(t-2|sint|)N3′(t)=-(25+|sin5t|)N3(t)5+|cos6t|+N3(t)+(1+|cos3t|)N1(t)3+|sin2t|+N1(t)+(1+|sin3t|)N2(t)3+|cos2t|+N2(t)+14(1+sin22t)N3(t-2|cos2t|)e-4N3(t-2|cos2t|)+14(1+cos22t)N3(t-2|sin2t|)e-4N3(t-2|sin2t|).
Obviously, aii-=25,bii+=6,(i=1,2,3),aij+=2,bij-=3,(i,j=1,2,3,i≠j),cij+=1/2,γij-=4,(i=1,2,3,j=1,2). Let  K1=e, then we have
(4.2)256+e=aii-bii++K1>∑j=1,j≠i3‍aij+bij-+∑j=12‍cij+γij-eK1=43+14e2,256+e=aii-bii++K1>∑j=1,j≠i3‍aij+bij-+∑j=12‍cij+=73.
Then (4.2) imply that the system (4.1) satisfies (2.1) and (2.6). Hence, from Theorems 2.1 and 2.2, the solution N(t) of system (4.1) with Dij(t,N)=aij(t)N/(bij(t)+N)(i,j=1,2,3) and φ∈E1={φ∣φ∈C+,φ(0)>0and0≤φi(t)<e,forall,t∈[-2,0],i=1,2,3} is exponentially extinct as t→+∞ and N(t)=N(t,0,φ)=O(e-κt),κ≈0.0001. The fact is verified by the numerical simulation in Figure 1.Figure 1
Numerical solutionN(t)=(N1(t),N2(t),N3(t))T of system (4.1) for initial value φ(t)≡(0.5,1.7,2.6)T.Example 4.2.
Consider the following Nicholson's blowflies system with patch structure and nonlinear density-dependent mortality terms:(4.3)N1′(t)=-(12+|sint|)+(11+|cost|)e-N1(t)+(1+12|sint|)-(12+12|cost|)e-N2(t)+(1+12|sint|)-(12+12|cost|)e-N3(t)+14(1+cos2t)N1(t-2|sint|)e-N1(t-2|sint|)+14(1+sin2t)N1(t-2|cost|)e-N1(t-2|cost|)N2′(t)=-(12+|cost|)+(11+|sint|)e-N2(t)+(1+12|cost|)-(12+12|sint|)e-N1(t)+(1+12|cost|)-(12+12|sint|)e-N3(t)+14(1+sin2t)N2(t-2|cost|)e-N2(t-2|cost|)+14(1+cos2t)N2(t-2|sint|)e-N2(t-2|sint|)N3′(t)=-(12+|sin2t|)+(11+|cos2t|)e-N3(t)+(1+12|sin2t|)-(12+12|cos2t|)e-N1(t)+(1+12|sin2t|)-(12+12|cos2t|)e-N2(t)+14(1+cos22t)N3(t-2|sin2t|)e-N3(t-2|sin2t|)+14(1+sin22t)N3(t-2|cos2t|)e-N3(t-2|cos2t|).
Obviously, aii-=12,bii-=11,bii+=12,(i=1,2,3),aij+=3/2,bij-=1/2,bij+=1,(i,j=1,2,3,i≠j),cij+=1/2,γij-=1,(i=1,2,3,j=1,2). Let K2=1, then we have
(4.4)12=aii->∑j=1,j≠i3‍aij++(bii+-∑j=1,j≠in‍bij-)e-K2+∑j=12‍cij+γij-e=3+12e,i=1,2,3,-aii(t)+bii(t)+∑j=1,j≠i3‍(aij(t)-bij(t))=0,i=1,2,3,11=bii->1+K22bii++∑j=1,j≠i3‍bij++∑j=12‍cij+=10,i=1,2,3.
Then (4.4) imply that the system (4.3) satisfies (3.1), (3.2), (3.11), and (3.12). Hence, from Theorems 3.1 and 3.2, the solution N(t) of system (4.1) with Dij(t,N)=aij(t)-bij(t)e-N(i,j=1,2,3) and φ∈E2={φ∣φ∈C+,φ(0)>0and0≤φi(t)<1,forall,t∈[-2,0],i=1,2,3} is exponentially extinct as t→+∞ and N(t)=N(t,0,φ)=O(e-κt),κ≈0.0001. The fact is verified by the numerical simulation in Figure 2.Figure 2
Numerical solutionN(t)=(N1(t),N2(t),N3(t))T of system (4.3) for initial value  φ(t)≡(0.1,0.5,0.9)T.Remark 4.3.
To the best of our knowledge, few authors have considered the problems of the extinction of Nicholson's blowflies model with patch structure and nonlinear density-dependent mortality terms. Wang [10] and Hou et al.  [11] have researched the permanence and periodic solution for scalar Nicholson's blowflies equation with a nonlinear density-dependent mortality term. Liu and Gong  [12] have considered the permanence for Nicholson-type delay systems with nonlinear density-dependent mortality terms and Takeuchi et al. [13] have investigated the global stability of population model with patch structure. Faria [14], Liu [15], and Berzansky et al. [16] have, respectively, studied the local and global stability of positive equilibrium for constant coefficients of Nicholson's blowflies model with patch structure. It is clear that all the results in [10–16] and the references therein cannot be applicable to prove the extinction of (4.1) and (4.3). This implies that the results of this paper are new.



---

*Source: 302065-2012-12-31.xml*