# On a System Modelling a Population with Two Age Groups

Authors: Hongliang Gao; Ruyun Ma
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

A system of first order ordinary differential equations describing a population divided into juvenile and adult age groups is studied. The system is not cooperative but its linear part is, and this makes it possible to establish the existence and nonexistence results of positive solutions for the system in terms of the principal eigenvalue of the corresponding linearized system.


---

## Body

## 1. Introduction

In this paper, we will study the following problem:(1)
u
′
=
a
(
t
)
v
-
c
(
t
)
u
-
e
u
[
u
+
v
]
,
v
′
=
b
(
t
)
u
-
d
(
t
)
v
-
f
v
[
u
+
v
]
,
where a, b, c, and d
∈
C
(
R
,
(
0
,
∞
)
) are ω-periodic functions and e and f are positive constants.We are interested in dividing the individuals within a population into two age groups. The first group contains all newborns in addition to all young individuals who are unable to produce newborns; such group will be referred to as the juvenile group. The second group, which we will call the adult group, contains all individuals who can produce newborns in addition to old individuals who may not be able to produce newborns. The functionsu and v represent, respectively, the total number of individuals who belong to the juvenile and adult groups. As adults give birth to juveniles, the function a corresponds to the birth rate of the population. Juveniles are lost both through death and through becoming adults; the function c corresponds to this overall loss. The function b gives the rate at which juveniles become adults and the function d corresponds to the death rate of adult population. The terms -
e
u
[
u
+
v
] and -
f
v
[
u
+
v
] correspond to decrease in population size due to competition for limited resources.In natural environments the number of individuals of a population changes in time in different ways. Many observations show that the number of individuals of a population can have large oscillations in nature. In the earlier models the population is characterized by its size which is the total number of individuals within the population or the total biomass. One of such models is the Malthus model for the human population growth. P. F. Verhulst in 1838 introduced another model which is known as the realistic model; see [1, 2]. Models presenting qualitatively this type of behavior are density dependent unstructured population models; the most well-known model for interspecific competition has been proposed by Lotka and Volterra and has been studied extensively by Bonhoeffer, Borrelli, and Murray; see also [1–3]. In fact, up to the mid of the 20th century most models characterize the population by its size or total biomass. In such models the population is considered as homogeneous; that is, the models do not distinguish between the individuals within the population. Models that involve structured population are called structured population models. A structured model describes how individuals move in time among different groups and thus describes the dynamics of population groups and as a result it describes the dynamics of the whole population. For other related results, we refer the readers to [4–7] and the references therein.Recently, a model for the growth of a population of two age groups (adult and juvenile) in which there is competition for limited resources has been considered in [8], where the authors assumed that the population is homogenous with common birth rate, common death rate, and common inhibiting constants. They established a time-invariant structure under general conditions and discussed the stability of the equilibrium points. More concretely, at time t the net rate of change in the populations of the two groups is modeled by the system
(2)
u
′
=
a
v
-
c
u
-
e
u
[
u
+
v
]
,
v
′
=
b
u
-
d
v
-
f
v
[
u
+
v
]
;
the authors showed in [8] that if c
d
<
a
b, then (2) has a unique positive equilibrium point value in addition to the trivial equilibrium point u
=
v
=
0.Obviously, in [8], since a, b, c, and d are positive constants, the solutions of (2) can be explicitly given and some estimates can be carried out easily. However, when a, b, c, and d are positive functions, the method of [8] cannot be applied to deal with the system (1) any more. If a, b, c, and d are not constants, whether the system (1) has a positive solution or not is a natural question. Inspired by above considerations, in the present paper, we will first establish the lower and upper solutions method for more general system
(3)
u
′
+
p
(
t
)
u
=
f
(
u
,
v
)
,
v
′
+
q
(
t
)
v
=
g
(
u
,
v
)
,
where f and g
:
R
×
R
→
R are continuous functions and p, q are ω-periodic continuous functions, and then we will prove the existence of positive solutions for system (1) by applying above method.Our main results can be stated as below.Theorem 1.
Suppose that the functionsu
_
,
v
_
,
u
¯
,
v
¯
∈
X
∩
C
1
(
R
,
R
). (
u
_
,
v
_
) and (
u
¯
,
v
¯
) are ordered coupling lower and upper solutions of systems (3); the following condition is hold (H). There exists M
1,  N
1
>
0 such that, for any v
_
≤
v
≤
v
¯ and u
_
≤
u
≤
u
¯,
(4)
f
(
u
¯
,
v
)
-
f
(
u
,
v
)
≥
-
M
1
(
u
¯
-
u
)
,
f
(
u
,
v
)
-
f
(
u
_
,
v
)
≥
-
M
1
(
u
-
u
_
)
,
g
(
u
,
v
¯
)
-
g
(
u
,
v
)
≥
-
N
1
(
v
¯
-
v
)
,
g
(
u
,
v
)
-
g
(
u
,
v
_
)
≥
-
N
1
(
v
-
v
_
)
.
Then, the problem (3) has at least one solution (
u
*
,
v
*
) with u
_
≤
u
*
≤
u
¯,  v
_
≤
v
*
≤
v
¯.Theorem 2.
There exists a positive periodic solution of systems (1) if and only if λ
1
(
M
)
<
0.Remark 3.
To overcome the difficulties caused by the spatially heterogeneous, we discuss the system (3) by lower and upper solutions method established in Theorem 1 and obtain the necessary and sufficient conditions for the existence of positive periodic solutions of (1) in terms of the principal eigenvalue of the associated linear system. For other related results on the study of differential systems via lower and upper solutions method, we refer the readers here to [9–11] and the references listed therein.The rest of the paper is organized as follows. In Section2, we establish the lower and upper solutions methods for the system (3). In Section 3, we obtain the necessary and sufficient conditions for the existence of a positive periodic solution of (1).


## 2. Lower and Upper Solutions Method

In this section, we will develop lower and upper solutions method for system (3).LetX be a Banach space defined as
(5)
X
:
=
{
u
∈
C
(
R
,
R
)
∣
u
(
t
)
=
u
(
t
+
ω
)
,
t
∈
R
}
.Definition 4.
Assume that the functionsu
_, v
_, u
¯, v
¯
∈
X
∩
C
1
(
R
,
R
). Then, (
u
_
,
v
_
) and (
u
¯
,
v
¯
) are called ordered coupling lower and upper solutions of systems (3), respectively, if u
_
≤
u
¯ and v
_
≤
v
¯ satisfying
(6)
u
_
′
+
p
(
t
)
u
_
≤
f
(
u
_
,
v
)
,
v
_
≤
v
≤
v
¯
,
v
_
′
+
q
(
t
)
v
_
≤
g
(
u
,
v
_
)
,
u
_
≤
u
≤
u
¯
,
u
¯
′
+
p
(
t
)
u
¯
≥
f
(
u
¯
,
v
)
,
v
_
≤
v
≤
v
¯
,
v
¯
′
+
q
(
t
)
v
¯
≥
g
(
u
,
v
¯
)
,
u
_
≤
u
≤
u
¯
.Proof of Theorem1.
By the condition (H), there existM
≥
M
1 and  N
≥
N
1 such that ∫
0
ω
‍
[
p
(
t
)
+
M
]
d
t
>
0 and  ∫
0
ω
‍
[
q
(
t
)
+
N
]
d
t
>
0.
For anyφ, ψ
∈
X, we consider the following linear problem:
(7)
u
′
+
p
(
t
)
u
+
M
u
=
f
(
φ
,
ψ
)
+
M
φ
,
v
′
+
q
(
t
)
v
+
N
v
=
g
(
φ
,
ψ
)
+
N
ψ
.
It is well known that the system (7) is equivalent to the equation
(8)
(
u
(
t
)
,
v
(
t
)
)
=
(
∫
t
t
+
w
‍
G
1
(
t
,
s
)
[
f
(
φ
(
s
)
,
ψ
(
s
)
)
+
M
φ
(
s
)
]
d
s
,
kkkkk
∫
t
t
+
w
‍
G
2
(
t
,
s
)
[
g
(
φ
(
s
)
,
ψ
(
s
)
)
+
N
ψ
(
s
)
]
d
s
)
,
where
(9)
G
1
(
t
,
s
)
=
e
∫
t
s
‍
[
p
(
θ
)
+
M
]
d
θ
e
∫
0
ω
‍
[
p
(
θ
)
+
M
]
d
θ
-
1
,
G
2
(
t
,
s
)
=
e
∫
t
s
‍
[
q
(
θ
)
+
N
]
d
θ
e
∫
0
ω
‍
[
q
(
θ
)
+
N
]
d
θ
-
1
,
and, consequently,
(10)
(
u
,
v
)
=
T
(
φ
,
ψ
)
with
(11)
T
(
φ
,
ψ
)
=
(
∫
t
t
+
w
‍
G
1
(
t
,
s
)
[
f
(
φ
(
s
)
,
ψ
(
s
)
)
+
M
φ
(
s
)
]
d
s
,
kkkkki
∫
t
t
+
w
‍
G
2
(
t
,
s
)
[
g
(
φ
(
s
)
,
ψ
(
s
)
)
+
N
ψ
(
s
)
]
d
s
)
.
It is easy to see that T
:
X
2
→
X
2 is completely continuous.
LetP
=
{
(
u
,
v
)
∈
X
2
∣
u
_
≤
u
≤
u
¯
,
v
_
≤
v
≤
v
¯
}; then P is bounded closed convex subset of X
2. So, T
:
P
→
X
2 is also completely continuous. We will show that T
:
P
→
P.
Let(
φ
,
ψ
)
∈
P,  (
u
,
v
)
=
T
(
φ
,
ψ
), and w
=
u
¯
-
u
,
z
=
v
¯
-
v. Since u
_
≤
φ
≤
u
¯ and  v
_
≤
ψ
≤
v
¯, we have
(12)
u
¯
′
+
p
(
t
)
u
¯
+
M
u
¯
≥
f
(
u
¯
,
ψ
)
+
M
u
¯
,
v
¯
′
+
q
(
t
)
v
¯
+
N
v
¯
≥
g
(
φ
,
v
¯
)
+
N
v
¯
.
Hence, w, z satisfy
(13)
w
′
+
p
(
t
)
w
+
M
w
≥
f
(
u
¯
,
ψ
)
-
f
(
φ
,
ψ
)
+
M
(
u
¯
-
φ
)
≥
0
,
z
′
+
q
(
t
)
z
+
N
z
≥
g
(
φ
,
v
¯
)
-
g
(
φ
,
ψ
)
+
N
(
v
¯
-
ψ
)
≥
0
.
T is strongly positive and implies w
≥
0,  z
≥
0; that is, u
≤
u
¯,  v
≤
v
¯. By a similar method, we have u
≥
u
_ and v
≥
v
_. Consequently, T
:
P
→
P. By the Schauder fixed point theorem, T has a fixed point (
u
*
,
v
*
) in P. Therefore, the problem (3) has at least one solution (
u
*
,
v
*
) with u
_
≤
u
*
≤
u
¯ and  v
_
≤
v
*
≤
v
¯.


## 3.  Existence and Nonexistence of Positive Periodic Solutions

We consider the system(14)
u
′
=
A
(
t
)
u
+
B
(
t
)
v
,
v
′
=
C
(
t
)
u
+
D
(
t
)
v
,
where A, B, C, and D are ω periodic functions, B, C
∈
C
(
R
[
0
,
∞
)
).Lemma 5.
Suppose that there exist functionsu
0
,
v
0
∈
X
∩
C
1
(
R
,
R
) such that u
0
>
0
,
v
0
>
0, and
(15)
u
0
′
≥
A
(
t
)
u
0
+
B
(
t
)
v
0
,
v
0
′
≥
C
(
t
)
u
0
+
D
(
t
)
v
0
,
where equality does not hold in all of the equations in (15). Then, (14) satisfies the strong maximum principle; that is, if u, v
∈
X
∩
C
1
(
R
,
R
) such that
(16)
u
′
≥
A
(
t
)
u
+
B
(
t
)
v
,
v
′
≥
C
(
t
)
u
+
D
(
t
)
v
,
then either (i)  u, v
≡
0 or (ii)  u
>
0 and v
>
0.Proof.
Suppose that the result is false. Then, there existu
1, v
1 not both identically zero satisfying inequalities (16) but not satisfying (ii) in the conclusion of the theorem. For 0
≤
τ
≤
1 define u
τ
=
(
1
-
τ
)
u
0
+
τ
u
1 and v
τ
=
(
1
-
τ
)
v
0
+
τ
v
1. Then, there exists τ
0, 0
<
τ
0
≤
1, such that u
τ, v
τ
>
0 for 0
≤
τ
<
τ
0 and either u
τ
0 or v
τ
0 has a zero. We may assume without loss of generality that there exists t
1 such that u
τ
0
(
t
1
)
=
0. Then,
(17)
u
τ
0
′
≥
A
(
t
)
u
τ
0
+
B
(
t
)
v
τ
0
,
and so
(18)
u
τ
0
′
-
A
(
t
)
u
τ
0
≥
0
.
Moreover, K is chosen sufficiently large to ensure that ∫
0
ω
‍
[
K
-
A
(
t
)
]
d
t
>
0, so it follows from the strongly positive that u
τ
0
≡
0 or u
τ
0
>
0. Since u
τ
0
(
t
1
)
=
0, it follows that u
τ
0
≡
0. Hence, by (17), it also follows that v
τ
0
≡
0. Since both u
1
,
v
1 are not identically zero and τ
0
<
1, from u
τ
0
=
(
1
-
τ
0
)
u
0
+
τ
0
u
1
≡
0 and v
τ
0
=
(
1
-
τ
0
)
v
0
+
τ
0
v
1
≡
0, we have u
0
=
-
(
τ
0
/
(
1
-
τ
0
)
)
u
1 and v
0
=
-
(
τ
0
/
(
1
-
τ
0
)
)
v
1. This is impossible as u
0, v
0 satisfy (15) and u
1, v
1 satisfy (16) and so the proof is complete.Corollary 6.
SupposeA
(
t
)
+
B
(
t
)
<
0 and C
(
t
)
+
D
(
t
)
<
0. Then, the system (14) satisfies the maximum principle.Proof.
The result follows from Lemma5 by choosing u
0
=
v
0
=
K where K is any positive number.Lemma 7.
The system (14) has a principal eigenvalue; that is, there exists Λ
∈
R and functions u
,
v
∈
X
∩
C
1
(
R
,
R
) such that u
,
v
>
0 and
(19)
u
′
-
A
(
t
)
u
-
B
(
t
)
v
=
Λ
u
,
v
′
-
C
(
t
)
u
-
D
(
t
)
v
=
Λ
v
.Proof.
DefineL
:
[
X
∩
C
1
(
R
,
R
)
]
2
→
X
2 by
(20)
L
(
u
v
)
=
(
u
′
v
′
)
,
and define the matrix M
(
t
) by
(21)
M
(
t
)
=
(
A
(
t
)
B
(
t
)
C
(
t
)
D
(
t
)
)
.
By essentially the same argument as in [12, Lemma 12], if K
>
0 is sufficiently large, then L
-
M
+
K is an invertible operator such that (
L
-
M
+
K
)
-
1 is compact. If, moreover, K is chosen sufficiently large to ensure that A
(
t
)
+
B
(
t
)
-
K
<
0 and C
(
t
)
+
D
(
t
)
-
K
<
0, it follows from Corollary 6 that (
L
-
M
+
K
)
-
1 is strongly positive.
Since(
L
-
M
+
K
)
-
1 is compact and strongly positive, (
L
-
M
+
K
)
-
1 has a positive principal eigenvalue (μ). Thus, there exists U
0
=
(
u
0
v
0
) with u
0
,
v
0
>
0 such that (
L
-
M
+
K
)
-
1
U
0
=
μ
U
0. Hence, (
L
-
M
)
U
0
=
(
1
/
μ
-
K
)
U
0 and so L
-
M has a principal eigenvalue.We will denote the principal eigenvalue ofL
-
M by λ
1
(
M
).Corollary 8.
SupposeM
1
(
t
) and M
2
(
t
) are cooperative matrices (i.e., matrices with positive entries in the off-diagonal elements) such that M
1
(
t
)
≥
M
2
(
t
) (i.e., the (
i
,
j
)th element of M
1≥ the (
i
,
j
)th element of M
2 but M
1
≢
M
2). Then, λ
1
(
M
1
)
<
λ
1
(
M
2
).Proof.
There existsU
0
=
(
u
0
v
0
) such that u
0, v
0
>
0 and [
L
-
M
1
-
λ
1
(
M
1
)
]
U
0
=
0. Then,
(22)
[
L
-
M
2
-
λ
1
(
M
1
)
]
U
0
=
[
L
-
M
1
-
λ
1
(
M
1
)
]
U
0
+
(
M
1
-
M
2
)
U
0
=
(
M
1
-
M
2
)
U
0
≥
0
.
But (
M
1
-
M
2
)
U
0
≢
0 and so by Lemma 5 the system [
L
-
M
2
-
λ
1
(
M
1
)
]
U
=
0 satisfies the strong maximum principle. Hence, if γ denotes the principal eigenvalue for the system L
-
M
2
-
λ
1
(
M
1
)
I, it follows easily that γ
>
0. Clearly, L
-
M
2 has principal eigenvalue λ
1
(
M
1
)
+
γ
>
λ
1
(
M
1
) and so λ
1
(
M
1
)
<
λ
1
(
M
2
), and the proof is complete.System (1) can be rewritten as
(23)
L
(
u
v
)
=
M
(
t
)
(
u
v
)
-
N
(
u
v
)
,
where M
(
t
)
=
(
-
c
(
t
)
a
(
t
)
b
(
t
)
-
d
(
t
)
) and N
:
X
2
→
X
2 such that
(24)
N
(
u
v
)
=
(
e
u
[
u
+
v
]
f
v
[
u
+
v
]
)
.AlthoughM
(
t
) is a cooperative matrix, system (1) is not a cooperative system. We can give necessary and sufficient conditions for the existence of a positive solution.Proof of Theorem2.
Supposeλ
1
(
M
)
<
0. Then, there exists ϕ
1, ϕ
2
>
0 such that
(25)
(
L
-
M
)
(
ϕ
1
ϕ
2
)
=
λ
1
(
M
)
(
ϕ
1
ϕ
2
)
;
that is,(26)
ϕ
1
′
=
a
(
t
)
ϕ
2
-
c
(
t
)
ϕ
1
+
λ
1
(
M
)
ϕ
1
,
ϕ
2
′
=
b
(
t
)
ϕ
1
-
d
(
t
)
ϕ
2
+
λ
1
(
M
)
ϕ
2
.
Let(
u
_
v
_
)
=
ɛ
(
ϕ
1
ϕ
2
) and (
u
¯
v
¯
)
=
(
K
K
). We will show that (
u
_
v
_
) and (
u
¯
v
¯
) satisfy the hypotheses of Theorem 1 provided that ɛ
>
0 is chosen sufficiently small and K is chosen sufficiently large. Let K
=
max
⁡
{
max
⁡
a
(
t
)
/
e
,
max
⁡
b
(
t
)
/
f
}. Then, a
(
t
)
-
e
K
≤
0 and so
(27)
a
(
t
)
v
-
c
(
t
)
u
¯
-
e
u
¯
[
u
¯
+
v
]
=
[
a
(
t
)
-
e
K
]
v
-
c
(
t
)
K
-
e
K
2
≤
0
=
u
¯
′
whenever v
≥
0. Similarly,
(28)
v
¯
′
≥
b
(
t
)
u
-
d
(
t
)
v
¯
-
f
v
¯
[
u
+
v
¯
]
,
whenever u
≥
0.
Letɛ
0
=
min
⁡
{
min
⁡
a
(
t
)
/
e
max
⁡
ϕ
1
,
min
⁡
b
(
t
)
/
f
max
⁡
ϕ
2
}. Then, when ɛ
<
ɛ
0
,
a
(
t
)
-
ɛ
e
ϕ
1
≥
0 and b
(
t
)
-
ɛ
f
ϕ
2
≥
0.
Hence, whenɛ
<
ɛ
0,  u
_
=
ɛ
ϕ
1, and v
≥
ɛ
ϕ
2, we have
(29)
u
_
′
-
a
(
t
)
v
+
c
(
t
)
u
_
+
e
u
_
[
u
_
+
v
]
=
ɛ
[
ϕ
1
′
+
c
(
t
)
ϕ
1
-
a
(
t
)
ϕ
2
]
+
a
(
t
)
[
ɛ
ϕ
2
-
v
]
+
e
ɛ
ϕ
1
[
ɛ
ϕ
1
+
v
]
=
ɛ
λ
1
(
M
)
ϕ
1
+
ɛ
a
(
t
)
ϕ
2
-
[
a
(
t
)
-
ɛ
e
ϕ
1
]
v
+
ɛ
2
e
[
ϕ
1
]
2
≤
ɛ
λ
1
(
M
)
ϕ
1
+
ɛ
a
(
t
)
ϕ
2
-
[
a
(
t
)
-
ɛ
e
ϕ
1
]
ɛ
ϕ
2
+
ɛ
2
e
[
ϕ
1
]
2
=
ɛ
λ
1
(
M
)
ϕ
1
+
ɛ
2
e
ϕ
1
[
ϕ
1
+
ϕ
2
]
<
0
,
when ɛ is sufficiently small.
Similarly, whenv
_
=
ɛ
ϕ
2,  u
≥
ɛ
ϕ
1, and ɛ is sufficiently small, we have
(30)
v
_
′
-
b
(
t
)
u
+
d
(
t
)
v
_
+
f
v
_
[
u
+
v
_
]
<
0
.
Hence, by Theorem 1, there exists a positive solution of system (1).
Suppose now that system (1) has a solution (
u
0
v
0
) with u
0, v
0
>
0. Then, (
u
0
v
0
) is a solution of the system
(31)
u
′
+
[
c
(
t
)
+
e
q
(
t
)
]
u
-
a
(
t
)
v
=
0
,
v
′
-
b
(
t
)
u
+
[
d
(
t
)
+
f
q
(
t
)
]
v
=
0
,
where q
(
t
)
=
u
0
(
t
)
+
v
0
(
t
). Hence, (
u
0
v
0
) may be regarded as the principal eigenfunction corresponding to the principal eigenvalue λ
=
0 of the system (
L
-
M
q
(
t
)
)
U
=
0, where
(32)
M
q
(
t
)
=
(
-
c
(
t
)
-
e
q
(
t
)
a
(
t
)
b
(
t
)
-
d
(
t
)
-
f
q
(
t
)
)
,
U
=
(
u
v
)
.
Hence,λ
1
(
M
q
)
=
0. As M
q
≤
M but M
q
≢
M, it follows from Corollary 8 that λ
1
(
M
)
<
λ
1
(
M
q
)
=
0 and the proof is complete.Remark 9.
LetM
1
=
(
-
1
1
1
-
1
). Then λ
1
(
M
1
)
=
0. If M
>
M
1, then by Corollary 8 and Theorem 2 we know that system (1) has a positive solution. In fact, we can deduce from M
>
M
1 that a
b
>
c
d. Hence, our main results extend and complement the corresponding ones of [8] to some extent.Example 10.
Leta
(
t
)
=
2
+
sin
(
π
t
/
2
); b
(
t
)
=
e
t
/
2 and t
∈
[
0,2
] with b
(
t
+
4
)
=
b
(
t
) and  b
(
t
)
=
b
(
-
t
);  c
(
t
)
=
(
1
/
2
)
|
cos
⁡
(
π
t
/
4
)
|
+
1
/
2;  d
(
t
)
=
(
1
/
4
)
t
+
1
/
2 and  t
∈
[
0,2
] with d
(
t
+
4
)
=
d
(
t
) and  d
(
-
t
)
=
d
(
t
). We consider the following problem:
(33)
u
′
=
a
(
t
)
v
-
c
(
t
)
u
-
u
[
u
+
v
]
,
v
′
=
b
(
t
)
u
-
d
(
t
)
v
-
v
[
u
+
v
]
.
It is easy to verify that M
>
M
1; then λ
1
(
M
)
<
0 by Corollary 8. Therefore, we know that system (33) has a positive solution by Theorem 2.



---

*Source: 920348-2014-06-23.xml*