
\documentclass[twoside]{article}
\usepackage{amssymb, amsmath} 
\pagestyle{myheadings}

\markboth{\hfil Positive periodic solutions \hfil EJDE--2002/71}
{EJDE--2002/71\hfil D. Jiang, J. Wei, \&  B. Zhang \hfil}
\begin{document}
\title{\vspace{-1in}\parbox{\linewidth}{\footnotesize\noindent
{\sc  Electronic Journal of Differential Equations},
Vol. {\bf 2002}(2002), No. 71, pp. 1--13. \newline
ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu
\newline ftp  ejde.math.swt.edu  (login: ftp)}
 \vspace{\bigskipamount} \\
 %
  Positive periodic solutions of functional
  differential equations and population models
 %
\thanks{ {\em Mathematics Subject Classifications:} 34K13, 92B05.
\hfil\break\indent
{\em Key words:} Functional differential equations, positive periodic
 solution, \hfil\break\indent
 population models.
\hfil\break\indent
\copyright 2002 Southwest Texas State University. \hfil\break\indent
Submitted July 5, 2002. Published July 30, 2002.
\hfill\break\indent
Jiang \& Wei were partially supported by the National
Natural Science Fundation of China.} }
\date{}
%
\author{Daqing Jiang, Junjie Wei, \& Bo Zhang}
\maketitle

\begin{abstract} 
  In this paper, we employ Krasnosel'skii's fixed point theorem
  for cones to study the existence of positive periodic solutions
  to a system of infinite delay equations,
  $$  x'(t) = A(t)x(t) + f(t,x_t). $$
  We prove two general theorems and establish new periodicity
  conditions for several  population growth models.
\end{abstract}

\newtheorem{theorem}{Theorem}[section]
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{lemma}[theorem]{Lemma}
\numberwithin{equation}{section}

\section{Introduction}

    In this paper we study the existence of positive periodic
solutions of the system of functional differential equations
\begin{equation}
\dot {x}(t) = A(t)x(t) + f(t,x_t)
\end{equation}
in which $A(t) = \mathop{\rm diag}[a_1(t),a_2(t),\dots, a_n(t)]$,
$a_j \in C(\mathbb{R},\mathbb{R})$ is $\omega$-periodic, $f(t, x_t)$ is a
function defined on $\mathbb{R}\times BC$, and $f(t,x_t)$ is $\omega$-periodic
whenever $x$ is $\omega$-periodic,
where $BC$ denotes the Banach space of bounded continuous functions
$\phi : \mathbb{R}\to \mathbb{R}^n$ with the norm
$\|\phi\| = \sup_{\theta\in \mathbb{R}}\sum^n_{j=1}|\phi_j(\theta)|$ where
$\phi = (\phi_1,\phi_2,\dots,\phi_n)^T,
\;\omega>0$ is a constant.
If $x\in BC$, then $x_t \in BC$ for any $t\in \mathbb{R}$ is defined by $x_t(\theta) =
x(t+\theta)$ for $\theta \in \mathbb{R}$.

One of the most used models, a prototype of (1.1), is the system of Volterra
integrodifferential equations
\begin{equation}
\dot {x}_i(t) = x_i(t) \Big[ a_i(t) - \sum^n_{j=1}b_{ij}(t)x_j(t)
-\sum^n_{j=1}\int^t_{-\infty}\!C_{ij}(t,s)g_{ij}(x_j(s))ds\Big]
\end{equation}
which governs the population growth of interacting species $x_j(t),\:
j=1,2,\dots, n$. The integral term here
specifies how much weight to
attach to the population at varies past times, in order to arrive at their
present effect on the resources availability. For an autonomous system (1.2),
there may be a stable equilibrium point of the population. When the potential
equilibrium point becomes unstable, there may exist a nontrivial periodic
solution.  Then the oscillation of solutions occurs. The existence of such
stable periodic solution is of quite fundamental importance biologically since
it concerns the long time survival of species. The study of such phenomena has
become an essential part of qualitative theory of differential equations. For
historical background, basic theory of periodicity, and discussions of
applications of (1.1) to
a variety of dynamical
models, we refer to the reader to, for example, the work of
Burton [1], Burton and Zhang [2], Cavani, Lizana, and Smith [3], Cushing [5],
Gopalsamy [7], Hadeler [8], Hatvani and Krisztin [9],
Hino, Murakami, Taito [10], Kuang [12], Makay [13], May [14],
May and Leonard [15], Stech [17], and the references therein.

It is our view that one of the most important problems in the study of
dynamical models and their applications is that of describing the nature of
the solutions for a large range of parameters involved.  From a numerical point
of view, the existence of a periodic solution of the approximation scheme must
also be studied. The usual approach to fulfill such requirements is to have a
set of test equations which are as general as possible and for which explicit
analytic conditions can be given. The next step is to characterize the
numerical methods which show the same existence result under the same
conditions when applied to the test equations.


In this part of our investigation, we provide a unified approach to the study
of existence of positive periodic solutions of system (1.1) under
general conditions and apply the results to some well-known models in
population dynamics. In Section 2, we prove the main results concerning equation (1.1).
Applications of these results to population models will be given in Section 3.


     Let $\mathbb{R}=(-\infty, +\infty), \mathbb{R}_+=[0,+\infty)$, and $\mathbb{R}_{-}=(-\infty,0]$
respectively. For each $x = (x_1,x_2,\dots,x_n)^T \in \mathbb{R}^n$,
the norm of $x$ is defined as $|x| = \sum^n_{j=1}|x_j|.\:$
$\mathbb{R}^n_+ = \{(x_1,x_2,\dots,x_n)^T \in \mathbb{R}^n : x_j \geq 0,\:
j=1,2,\dots,n\}$. We say that $x$ is ``positive" whenever $x \in
\mathbb{R}^n_+$.
$BC(X,Y)$ denotes the set of bounded continuous functions
$\phi: X\to Y$.



\section{Existence of Positive Periodic Solutions}


We establish the existence of positive periodic solutions of equation (1.1)
by applying Krasnosel'skii's fixed point theorem (see [6],[11]) on cones.
A compact operator will be constructed. It will be shown that
the operator has a fixed point, which corresponds to a periodic
solution of (1.1). We denote $f=(f_1,f_2,\dots, f_n)^T$ and assume
\begin{enumerate}
\item [(H1)]
$\int^\omega_0a_j(s)ds \not = 0$ for $j=1,2,\dots,n$.

\item [(H2)]  $f_j(t,\phi_t)\int^\omega_0a_j(s)ds\leq 0$ for all
$(t,\phi)\in \mathbb{R}\times BC(\mathbb{R},\mathbb{R}^n_+),\: j=1,2,\dots,n$.

\item [(H3)]
$f(t,x_t)$ is a continuous function of $t$ for each $x \in BC(\mathbb{R},\mathbb{R}^n_+)$.

\item [(H4)]
For any $L>0$ and $\varepsilon>0$, there exists $\delta >0$ such that \\
$[\phi,\psi \in BC,\: \|\phi\|\leq L,\: \|\psi\|\leq L,\:
\|\phi-\psi\| < \delta,\: 0\leq s\leq \omega] $ imply
\begin{equation}
|f(s,\phi_s) - f(s,\psi_s)| < \varepsilon.
\end{equation}
\end{enumerate}


\paragraph{Definition.} Let $X$ be a Banach space and $K$ be a closed, nonempty
subset of $X$. $K$ is a cone if
\begin{enumerate}
\item [(i)]
$\alpha u + \beta v \in K\;\;$ for all $u,v \in K$ and all $\alpha, \beta \geq
0$
\item [(ii)] $u, -u \in K$ imply $u=0$.
\end{enumerate}


\begin{theorem}[{Krasnosel'skii [11]}]  \label{thmA}
 Let $X$ be a Banach space, and let $K \subset X$ be a cone in $X$.
 Assume that $\Omega_1, \Omega_2$ are open subsets
of $X$ with $0 \in \Omega_1, \bar{\Omega}_1 \subset\Omega_2$, and let
\begin{equation*}
\Phi : K\cap (\bar {\Omega}_2\setminus\Omega_1) \to K
\end{equation*}
be a completely continuous operator such that either
\begin{enumerate}
\item [(i)]
$\|\Phi y\| \leq \|y\|\; \forall y \in K\cap \partial\Omega_1$ and
$\|\Phi y\| \geq \|y\|\; \forall y \in K\cap \partial\Omega_2$; or
\item [(ii)]
$\|\Phi y\| \geq \|y\|\; \forall y \in K\cap \partial\Omega_1$ and
$\|\Phi y\| \leq \|y\|\; \forall y \in K\cap\partial \Omega_2$.
\end{enumerate}
Then $\Phi$ has a fixed point in
$K\cap (\bar {\Omega}_2\setminus\Omega_1)$.
\end{theorem}

We now for $(t,s) \in \mathbb{R}^2,\:j=1,2,\dots,n$, we define
\begin{gather}
\sigma : = \min\Big\{ \exp(-2\int^\omega_0|a_j(s)|ds),
j=1,2,\dots,n\Big\} \\
G_j(t,s) = \frac{\exp(\int^t_s a_j(\nu)d\nu)}
{\exp(-\int^\omega_0 a_j(\nu)d\nu)-1}.
\end{gather}
We also define
\begin{equation*}
G(t,s) = \mathop{\rm diag}[G_1(t,s), G_2(t,s), \dots, G_n(t,s)].
\end{equation*}
It is clear that $G(t,s)=G(t+\omega,s+\omega)$ for
all $(t,s)\in \mathbb{R}^2$ and by (H2),
\begin{equation*}
G_j(t,s)f_j(u,\phi_u)\geq 0
\end{equation*}
for $(t,s)\in \mathbb{R}^2$ and $(u,\phi)\in \mathbb{R} \times BC(\mathbb{R},\mathbb{R}^n_+)$.

Next, we introduce two sets with $\omega$ and $\sigma$ as given above.
\begin{gather}
C_\omega = \{x\in C(\mathbb{R},\mathbb{R}^n):\; x(t+\omega) = x(t),\
; t \in \mathbb{R}\}, \nonumber \\
K = \{ x \in C_\omega: x_j(t)\geq \sigma \|x_j\|, t \in [0,\omega],
x = (x_1,x_2,\dots, x_n)^T\}.
\end{gather}
One may readily verify that $K$ is a cone.


Finally, we define an operator $\Phi : K\to K$ as
\begin{equation*}
(\Phi x)(t) = \int^{t+\omega}_tG(t,s)f(s,x_s)ds
\end{equation*}
for $x \in K,\: t\in \mathbb{R}$, where $G(t,s)$ is defined following (2.3). We denote
\begin{equation*}
(\Phi x) = \Big(\Phi_1x,\Phi_2x,\dots, \Phi_nx\Big)^T.
\end{equation*}


\begin{lemma} \label{lm2.1}
  $\Phi : K\to K$ is well-defined.
\end{lemma}

\paragraph{Proof.}
For each $x \in K$, since $f(t,x_t)$ is a continuous function of
$t$, we have $(\Phi x)(t)$ is continuous in $t$ and
\begin{eqnarray*}
(\Phi x)(t+\omega) & = & \int ^{t+2\omega}_{t+\omega}G(t+\omega,s)f(s,x_s)ds\\
& = & \int ^{t+\omega}_{t}G(t+\omega, v + \omega)f(v+\omega, x_{v+\omega})dv\\
& = & \int^{t+\omega}_tG(t,v)f(v,x_v)dv = (\Phi x)(t).
\end{eqnarray*}
Thus, $(\Phi x)\in C_\omega$.  Observe that
\begin{equation}
p_j := \frac{\exp(-\int^\omega_0 |a_j(\nu)|d\nu)}{\Big|
\exp(-\int^\omega_0 a_j(\nu)d\nu)-1\Big|}
\leq |G_j(t,s)| \leq
\frac{\exp(\int^\omega_0 |a_j(\nu)|d\nu)}
{|\exp(-\int^\omega_0 a_j(\nu)d\nu)-1|} =: q_j
\end{equation}
for all $s \in [t,t+\omega]$. Hence, for $x \in K$, we
have
\begin{equation}
\|\Phi_j x\| \leq q_j\int^\omega_0|f_j(s,x_s)|ds
\end{equation}
and
\begin{equation*}
(\Phi_j x)(t) \geq p_j\int^\omega_0|f_j(s,x_s)|ds\geq
\frac{p_j}{q_j}\|\Phi_j x\| \geq \sigma\|\Phi_j x\|.
\end{equation*}
Therefore, $(\Phi x) \in K$. This completes the proof. \hfill$\square$


\begin{lemma} \label{lm2.2}
$\Phi : K\to K$ is completely continuous, and
$x = x(t) $ is an $\omega$-periodic solution of (1.1) whenever
$x$ is a fixed point of $\Phi$.
\end{lemma}

\paragraph{Proof.}
 We first show that $\Phi$ is continuous. By (H4), for any $L>0$ and
$\varepsilon>0$, there exists a $\delta>0$ such that
$[\phi,\:\psi \in C_\omega,\:\|\phi\|\leq L,\: \|\psi\|\leq L$,
$\|\phi-\psi\|< \delta]$ imply
\begin{equation*}
\sup_{0\leq s\leq \omega} |f(s,\phi_s) - f(s,\psi_s)| <
\frac{\varepsilon}{q\omega}
\end{equation*}
where $q=\max_{1\leq j\leq n}q_j$.
If $x,\, y \in K$ with $\|x\|\leq L$, $\|y\|\leq L$, and $\|x - y\| < \delta$,
then
\begin{eqnarray*}
|(\Phi x)(t) - (\Phi y)(t)| &\leq & \int^{t+\omega}_t|G(t,s)||f(s,x_s) -
f(s,y_s)|ds\\
& \leq & q\int^\omega_0|f(s,x_s)-f(s,y_s)|ds< \varepsilon
\end{eqnarray*}
for all $t \in [0,\omega]$, where $|G(t,s)| = \max_{1\leq j\leq n}
|G_j(t,s)|$.
This yields
$\|(\Phi x) - (\Phi y)\| < \varepsilon$. Thus, $\Phi$ is continuous.


Next, we show that $f$ maps bounded sets into bounded sets. Indeed, let
$\varepsilon =1$. By (H4), for any $\mu>0$ there exists $\delta>0$ such
that
$[x,\:y \in BC$, $\|x\|\leq \mu$, $\|y\|\leq \mu$, $\|x-y\| <\delta]$ imply
$$ |f(s,x_s)-f(s,y_s)|< 1.
$$
Choose a positive integer
$N$ such that $\frac{\mu}{N}<\delta$. Let $x \in BC$ and define
$x^k(t) = x(t) \frac{k}{N}$ for $k=0,1,2,\dots, N$.
If $\|x\|\leq \mu$, then
$$\|x^k-x^{k-1}\| = \sup_{t\in \mathbb{R}}|x(t)\frac{k}{N} - x(t) \frac{k-1}{N}|\leq
\|x\|\frac{1}{N}\leq \frac{\mu}{N} < \delta.
$$
Thus,
\begin{equation*}
|f(s,x^k_s)-f(s,x^{k-1}_s)|<1
\end{equation*}
for all $s \in [0,\omega]$. This yields
\begin{eqnarray}
|f(s,x_s)| &\leq & \sum_{k=1}^N|f(s,x^k_s)-f(s,x^{k-1}_s)| + |f(s,0)|
\nonumber\\
& < & N + \sup_{s\in [0,\omega]}|f(s,0)| =: M_\mu.
\end{eqnarray}
It follows from (2.6) that for $t\in [0,\omega]$,
$$
\|\Phi x\| = \sup_{\theta\in \mathbb{R}}\sum^n_{j=1}|(\Phi_jx)(\theta)|
\leq\sum^n_{j=1} q_j\int^\omega_0|f_j(s,x_s)|ds
\leq q\omega M_\mu.
$$
Finally, for $t \in \mathbb{R}$ we have
\begin{eqnarray}
\frac{d}{dt}(\Phi x)(t) & = & G(t,t+\omega)f(t+\omega,x_{t+\omega})
- G(t,t)f(t,x_t) + A(t) (\Phi x)(t)\nonumber\\
& =& A(t)(\Phi x)(t) + [G(t,t+\omega) - G(t,t)]f(t,x_t)\nonumber \\
& = & A(t) (\Phi x)(t) + f(t,x_t).
\end{eqnarray}
Combine (2.6), (2.7), and (2.8) to obtain
\begin{eqnarray*}
|\frac{d}{dt}(\Phi x)(t)| & \leq & \|A\| q\int^\omega_0|f(s,x_s)|ds
+ |f(t,x_t)|\\
& \leq & \|A\|q\omega M_\mu + M_\mu
\end{eqnarray*}
where $\|A\| = \max_{1\leq j\leq n}\|a_j\|$.
Hence, $\{(\Phi x): x\in K, \|x\|\leq \mu\}$ is a family of uniformly bounded
and equicontinuous functions on $[0,\omega]$. By a theorem of Ascoli-Arzela
(Royden [16, p.169]), the function $\Phi$ is completely continuous.
It is clear from (2.8) that $x = x(t) $ is
an $\omega$-periodic solution of (1.1) whenever
$x$ is a fixed point of $\Phi$. This proves the lemma.

\begin{theorem} \label{thm2.1}
Assume (H1)-(H4) and that there are positive constants
$R_1$ and $R_2$ with $R_1<R_2$ such that
\begin{gather}
\sup_{\|\phi\|=R_1,\: \phi \in K} \int^\omega_0|f(s,\phi_s)|ds \leq
R_1/q\\
\inf_{\|\phi\|=R_2,\: \phi \in K} \int^\omega_0|f(s,\phi_s)|ds \geq
R_2/p
\end{gather}
for all $t\in [0,\omega]$ with $p=\min_{1\leq j\leq n}p_j$ and $q=\max_{1\leq
j\leq n}q_j$, where $p_j, q_j$ are defined in (2.5).
Then equation (1.1) has an $\omega$-periodic solution
$x$ with $R_1\leq \|x\|\leq R_2$.
\end{theorem}

\paragraph{Proof.}  Let $x \in K$ and $\|x\| = R_1$. By (2.9), we have
\begin{equation*}
|(\Phi x)(t)| \leq q\int^{t+\omega}_t|f(s,x_s)|ds
\leq q R_1/q = R_1.
\end{equation*}
This implies that $\|(\Phi x)\| \leq \|x\|$ for
$x \in K\cap \partial \Omega_1,\:
\Omega_1 = \{\psi \in C_\omega: \|\psi\|< R_1\}$. If $x \in K$ and $\|x\| =
R_2$, then
\begin{equation*}
(\Phi_j x)(t) \geq p_j\int^{t+\omega}_t|f_j(s,x_s)|ds.
\end{equation*}
Use (2.10) to obtain
$$|(\Phi x)(t)|\geq p\int^{t+\omega}_t|f(s,x_s)|ds \geq p(R_2/p)= R_2
$$
Thus, $\|(\Phi x)\| \geq \|x\|$ for $x \in K\cap \partial \Omega_2,\:
\Omega_2 = \{\psi \in C_\omega: \|\psi\|< R_2\}$.
By Theorem \ref{thmA}, $\Phi$ has a fixed point in $K\cap
(\bar{\Omega}_2\setminus \Omega_1)$. It follows from Lemma \ref{lm2.2} 
that (1.1) has an $\omega$-periodic solution $x$ with 
$R_1\leq \|x\|\leq R_2$. This completes the proof. \hfill$\square$

\begin{corollary} \label{coro2.1}
Assume (H1)-(H4) and
\begin{gather}
 \lim _{\phi \in K,\: \|\phi\|\to 0}
\frac{\int^\omega_0|f(s,\phi_s)|ds}{\|\phi\|} = 0,\\
 \lim_{\phi \in K,\: \|\phi\|\to \infty}
\frac{\int^\omega_0|f(s,\phi_s)|ds}{\|\phi\|} = \infty.
\end{gather}
Then (1.1) has an $\omega$-periodic solution.
\end{corollary}

\begin{theorem} \label{thm2.2}
Assume (H1)-(H4) and that there are positive constants
$R_1$ and $R_2$ with $R_1<R_2$ such that
\begin{gather}
\inf_{\|\phi\|=R_1,\: \phi \in K} \int^{\omega}_0|f(s,\phi_s)|ds \geq
R_1/p,\\
\sup_{\|\phi\|=R_2,\: \phi \in K}\int^{\omega}_0|f(s,\phi_s)|ds \leq
R_2/q.
\end{gather}
Then equation (1.1) has an $\omega$-periodic solution
$x$ with $R_1\leq \|x\|\leq R_2$.
\end{theorem}

The proof will be omitted here since it is similar to that of
Theorem \ref{thm2.1}.
We now state a corollary to the above theorem.

\begin{corollary} \label{coro2.2}
Assume (H1)-(H4) and that
\begin{gather}
\lim _{\phi \in K,\: \|\phi\|\to 0}
\frac{\int^\omega_0|f(s,\phi_s)|ds}{\|\phi\|} = \infty,\\
\lim_{\phi \in K,\: \|\phi\|\to \infty}
\frac{\int^\omega_0|f(s,\phi_s)|ds}{\|\phi\|} = 0.
\end{gather}
Then (1.1) has an $\omega$-periodic solution.
\end{corollary}

\section{Applications}

We now apply the main result obtained in the previous section to
some well-known models in population dynamics. First, we consider the
scalar integrodifferential equation
\begin{equation}
\dot {N}(t) = \alpha(t)N(t)\Big[1 -\frac{1}{N_0(t)}\int^0_{-
\infty}B(s)N(t+s)ds\Big]
\end{equation}
which governs the growth of the population $N(t)$ of a single species whose
members compete among themselves  for a limited amount of food and living
spaces. Equation (3.1) is a modification of the simple logistic equation
\begin{equation}
\dot {N}(t) = \alpha N(t)\Big[1 -\frac{N(t)}{N_0}\Big]
\end{equation}
where $\alpha$ is the intrinsic per capita growth rate and $N_0$ is the total
carrying capacity. We refer the reader to May [14] for a detailed model
construction from (3.2) to (3.1).
Suppose that
\begin{enumerate}
\item [(P1)] $\alpha(t),\: N_0(t)$ are real-valued and $\omega$-periodic
functions on $\mathbb{R}$ with $N_0(t)$ positive and
$\int^\omega_0\alpha(t)dt >0$.

\item [(P2)] $B(t)$ is nonnegative and piecewise continuous on $\mathbb{R}^-$ with
$\int^0_{-\infty}B(t)dt =1$.
\end{enumerate}


\begin{theorem} \label{thm3.1}
Under assumptions (P1) and (P2), equation (3.1) has at least one positive
$\omega$-periodic solution.
\end{theorem}

Let $a(t) = \alpha(t)$ and
$$f(t,x_t) =
-\frac{x(t)\alpha(t)}{N_0(t)}\int^0_{-\infty}B(s)x(t+s)ds.
$$
It is clear that $f(t,x_t)$ is $\omega$-periodic whenever $x$ is
$\omega$-periodic. We need to show that (H1)-(H4) hold. In fact,
$\int^\omega_0a(t)dt >0$ and $f(t,\phi_t)\leq 0$ for all $(t,\phi) \in
\mathbb{R}\times BC(\mathbb{R},\mathbb{R}_+)$. Thus, (H1) and (H2) are
satisfied. To verify (H3), let $x \in
BC(\mathbb{R},\mathbb{R}_+)$ with $\|x\|\leq B_1$. Notice that
$x(t)\alpha(t)/N_0(t)$ is continuous. By Lebesgue Dominant Convergence
Theorem, one can easily see that
$\int^0_{-\infty}B(s)x(t+s)ds$ is continuous for all $t \in \mathbb{R}$.
Thus, $f(t,x_t)$ is continuous in $t$.  Now let $x, y \in BC(\mathbb{R},\mathbb{R}_+)$ with
$\|x\|\leq L,\: \|y\|\leq L$ for some $L>0$. Then
\begin{eqnarray*}
\lefteqn{|f(t,x_t)-f(t,y_t)|} \\
&=& \Big|\frac{x(t)\alpha(t)}{N_0(t)}\int^0_{-
\infty}B(s)x(t+s)ds -
\frac{y(t)\alpha(t)}{N_0(t)}\int^0_{-\infty}B(s)y(t+s)ds\Big|\\
&\leq & \Big|\frac{x(t)\alpha(t)}{N_0(t)}\Big|\int^0_{-
\infty}B(s)|x(t+s)-y(t+s)| ds\\
& & +
\Big|\frac{(x(t)-y(t))\alpha(t)}{N_0(t)}\Big|\int^0_{-
\infty}B(s)|y(t+s)| ds \\
&\leq& \frac{L\|\alpha\|}{N_{0*}}\sup_{s \in \mathbb{R}_-}|x(t+s)-y(t+s)|
+\frac{|x(t)-y(t)|\|\alpha\|L}{N_{0*}},
\end{eqnarray*}
where $N_{0*}=\inf\{N(s): 0\leq s\leq \omega\}$. For any $\varepsilon>0$,
choose $\delta = \varepsilon N_{0*}/(2L\|a\|)$. If $\|x-y\|<\delta$, then
$$
|f(t,x_t)-f(t,y_t)|< L\|\alpha\|\delta/N_{0*} +
\delta\|\alpha\|L/N_{0*} = 2L\|\alpha\|\delta/N_{0*} = \varepsilon.
$$
This implies that (H4) holds.


We now verify (2.11) and (2.12) for equation (3.1).
For $\phi \in K$, we have $\phi(t)\geq \sigma\|\phi\|$ for all
$t \in [0,\omega]$. This yields
$$\frac{|f(t,\phi)|}{\|\phi\|} \leq \sup_{\tau\in
[0,\omega]}\frac{\alpha(\tau)}{N_0(\tau)}\int^0_{-\infty}B(s)ds
\|\phi\|\to 0$$
as $\|\phi\|\to 0$ and
$$\frac{|f(t,\phi)|}{\|\phi\|} \geq \inf_{\tau\in
[0,\omega]}\frac{\alpha(\tau)}{N_0(\tau)}\int^0_{-\infty}B(s)ds
(\sigma^2 \|\phi\|)\to +\infty$$
as $\|\phi\|\to \infty$. Thus, (2.11) and (2.12) are satisfied. By
Corollary \ref{coro2.1}, equation (3.1) has a positive $\omega$-periodic solution.

Next, consider a hematopoiesis model (Weng and Liang [19])
\begin{equation}
\dot{N}(t) = -\gamma(t)N(t) + \alpha(t) \int^{+\infty}_0B(s)
e^{-\beta(t)N(t-s)}ds
\end{equation}
where $N(t)$ is the number of red blood cells at time $t$, $\alpha, \beta,
\gamma \in C(\mathbb{R},\mathbb{R})$ are
$\omega$-periodic, and $B \in L^1(\mathbb{R}^+)$ is nonnegative and
piecewise continuous.  This is a generalized model
of the red cell system introduced by Wazewska-Czyzewska and Lasota [18]
\begin{equation}
\dot{n}(t) = -\gamma n(t) + \alpha e^{-\beta n(t-r)}
\end{equation}
where $\alpha, \beta, \gamma, r $ are constants with $r>0$.
Periodicity in equation (3.4) has been investigated extensively in the
literature (see Chow [4]).  Using Corollary \ref{coro2.2} with $a(t)=-\gamma(t)$ and
$$f(t,\phi_t)=\alpha(t)\int^{+\infty}_0B(s)e^{-\beta(t)\phi(t-s)}ds, $$
we obtain the following theorem whose proof is essentially the same as that of
Theorem \ref{thm3.1}, and therefore will be omitted.


\begin{theorem} \label{thm3.2}
 Suppose that $\alpha(t)\!>\!0,\: \beta(t)\geq 0$ for all
$t\in \mathbb{R}^+,\: \int^\omega_0\gamma(t)dt>0$.
$\int^{+\infty}_0B(s)ds >0$. Then
equation (3.3) has a positive $\omega$-periodic solution.
\end{theorem}

We now consider the Volterra equation mentioned in Section 1.1.
\begin{equation}
\dot {x}_i(t) = x_i(t) \Big[ a_i(t) - \sum^n_{j=1}b_{ij}(t)x_j(t)
-\sum^n_{j=1}\int^t_{-\infty}\!C_{ij}(t,s)g_{ij}(x_j(s))ds\Big]
\end{equation}
where $x_i(t)$ is the population of the $i$th species,
$a_i,\: b_{ij} \in C(\mathbb{R},\mathbb{R})$ are
$\omega$-periodic and $C_{ij}(t,s)$ is
piecewise continuous on $\mathbb{R}^2$.


\begin{theorem} \label{thm3.3}
Suppose that the following conditions hold for $i,j=1,2,\dots,n$.
\begin{enumerate}
\item [(i)]
 $\displaystyle \int^\omega_0a_i(s)ds \not = 0$,
\item [(ii)]
 $\displaystyle b_{ij}(t)\int^\omega_0a_i(s)ds \geq 0,\;\:
\displaystyle C_{ij}(t,s)\int^\omega_0a_i(\nu)d\nu \geq 0\:\:$ for all $(t,s) \in
\mathbb{R}^2$,
\item [(iii)]
 $g_{ij} \in C(\mathbb{R}^+,\mathbb{R}^+)$ is increasing with $g_{ij}(0)=0$,
\item [(iv)]
\hspace {.2cm} $\displaystyle \int^\omega_0b_{ii}(s)ds \not = 0$,
\item [(v)]
 $C_{ij}(t+\omega,s+\omega) = C_{ij}(t,s)$ for all $(t,s)\in
\mathbb{R}^2$ with  \\
$\sup_{t\in \mathbb{R}} \int^t_{-\infty}|C_{ij}(t,s)|ds < +\infty$.
\end{enumerate}
Then equation (3.5) has an $\omega$-periodic solution.
\end{theorem}

\paragraph{Proof.} For $x=(x_1,x_2, \dots, x_n)^T$, define
\begin{equation*}
f_i(t,x_t)= -x_i(t)\sum^n_{j=1}b_{ij}(t)x_j(t)
-x_i(t)\sum^n_{j=1}\int^t_{-\infty}\!C_{ij}(t,s)g_{ij}(x_j(s))ds
\end{equation*}
for $i=1,2,\dots,n$ and set $f = (f_1,f_2,\dots, f_n)^T$.
Then (3.5) can be written in the form of (1.1) with (H$_1)-$ (H4)
satisfied. Define
\begin{gather*}
b^* = \max\{ \|b_{ij}\|: i,j=1,2,\cdot,\cdot,\cdot,n\}\\
C^* = \max \Big\{\sup_{t\in \mathbb{R}} \sum^{n}_{j=1}
\int^t_{-\infty}|C_{ij}(t,s)|ds
: i=1,2,\cdot,\cdot,\cdot,n\Big\}\\
g^*(u) = \max \{g_{ij}(u): i,j=1,2,\cdot,\cdot,\cdot,n\}
\end{gather*}
Let $x \in K$, where $K$ is defined in (2.4).  Then
\begin{eqnarray*}
|f_i(t,x_t)| &\leq &
|x_i(t)|\Big[ b^*\|x\| + \sum^{n}_{j=1}
\int^t_{-\infty}|C_{ij}(t,s)|g_{ij}(\|x_j\|)ds\Big ]\\
&\leq & |x_i(t)| [b^*\|x\| + C^* g^*(\|x\|)]
\end{eqnarray*}
and
\begin{equation*}
\int^\omega_0|f(t,x_s)|ds \leq
\omega \|x\| [b^*\|x\| + C^* g^*(\|x\|)].
\end{equation*}
This implies
\begin{equation*}
\frac{\int^\omega_0|f(t,x_s)|ds}{\|x\|} \leq \omega[b^*\|x\| + C^*g^*(\|x\|)]
\to 0
\end{equation*}
as $\|x\|\to 0$. Since $x_i(t)\geq \sigma\|x_i\|$
for all $t \in \mathbb{R}$ whenever $x \in K$ and $b_{ij}(t),\:
C_{ij}(t,s)$ have the same sign, we obtain
\begin{eqnarray*}
\lefteqn{\int^\omega_0|f_i(t,x_s)|ds }\\
& =& \sum^n_{j=1}\int^\omega_0x_i(t)|b_{ij}(t)|x_j(t)dt
 + \sum^n_{j=1}\int^\omega_0x_i(t)\int^t_{-
\infty}|C_{ij}(t,s)|g_{ij}(x_j(s))dsdt\\
&\geq & \int^\omega_0|b_{ii}(t)|x_i(t)|^2dt
\geq \sigma^2\|x_i\|^2   \int^\omega_0|b_{ii}(t)|dt
\end{eqnarray*}
and
\[
\int^\omega_0|f(t,x_s)|ds
\geq \sigma^2\sum^n_{i=1}\|x_i\|^2  \min_{1\leq i \leq n}
\int^\omega_0|b_{ii}(t)|dt
\geq  \frac{\sigma^2}{n}\|x\|^2  \min_{1\leq i \leq n}
\int^\omega_0|b_{ii}(t)|dt.
\]
Here we have applied the inequality
$\Big(\sum^n_{i=1}\|x_i\|\Big)^2 \leq n \sum^n_{i=1}\|x_i\|^2$. Thus,
\begin{equation*}
\frac{\int^\omega_0|f(t,x_s)|ds}{\|x\|}
\to +\infty\;\; {\rm as}\;\;\|x\|\to +\infty.
\end{equation*}
By Corollary \ref{coro2.1}, equation (3.5) has an
$\omega$-periodic solution.


\paragraph{Remark.} % 3.1}.
 It is clear from the proof of Theorem \ref{thm3.3} that condition
(iv) can be replaced by
\begin{enumerate}
\item [(iv$^*$)] $\displaystyle \int^\omega_0\int^t_{-\infty}|C_{ii}(t,s)|dsdt
\not = 0$ and $g_{ii}(u) \to +\infty$ as $u\to +\infty$.
\end{enumerate}

We conclude this paper by investigating the following scalar equation
of advanced and delay type which is highly nonlinear and takes a quite
general form.
\begin{equation}
\dot{x}(t) = a(t)x(t) + \int^{+\infty}_{-\infty}C(s)g(t,x(t-\tau_0(t)),x(t-\tau_1(t)),
x(t+s))ds
\end{equation}
where $a, \tau_0, \tau_1 \in C(\mathbb{R},\mathbb{R})$ are $\omega$-periodic,
$g\in C(\mathbb{R}^4,\mathbb{R}^+)$, and $C\in L(\mathbb{R}_-,\mathbb{R})$
is piecewise continuous.


Let $u=(u_0,u_1,u_2)^T \in \mathbb{R}^3$ and define $g(t,u) = g(t,u_0,u_1,u_2)$. We
assume that $g(t+\omega,u) = g(t,u)$ for all $(t,y)\in \mathbb{R}^4$ and
\begin{enumerate}
\item [(P$^*_1$)] $\displaystyle \int^\omega_0a(s)ds \not =0,$

\item [(P$^*_2$)] $C(t)\displaystyle \int^\omega_0a(s)ds\leq 0\;$ for all $t \in \mathbb{R}$ and
$\;\displaystyle\int^{+\infty}_{-\infty}|C(u)|du = 1.$
\end{enumerate}

\begin{theorem} \label{thm3.4}
Under assumptions (P$^*_1$) and (P$^*_2$), equation (3.6) has at least
one positive $\omega$-periodic solution, provided one of the following
conditions hold
\begin{enumerate}
\item [(P$^*_3$)]
$$\lim_{|u|\uparrow 0}\; \max_{t\in [0,\omega]}\;\frac{g(t,u)}{|u|} = 0\;\;\;
{\rm and}\;\;\; \lim_{|u|\uparrow \infty,\; u_j\geq \sigma |u|,\: 0\leq j\leq
2}\;\min_{t\in [0,\omega]}\frac{g(t,u)}{|u|} =\infty,$$
\item [(P$^*_4$)]
$$\lim_{|u|\uparrow 0,\; u_j\geq \sigma |u|,\: 0\leq j\leq
2}\;\min_{t\in [0,\omega]}\frac{g(t,u)}{|u|} =\infty
\;\;\; {\rm and}\;\;
\lim_{|u|\uparrow +\infty}\; \max_{t\in [0,\omega]}\;\frac{g(t,u)}{|u|} = 0$$
\end{enumerate}
where $\sigma$ is given in (2.2) and $u=(u_0,u_1,u_2)^T$.
\end{theorem}

\paragraph{Proof}. We first show that under (P$^*_3$), conditions (2.11)
and (2.12) of Corollary \ref{coro2.1} hold. Let
$$
f(t,x_t) = \int^{+\infty}_{-\infty}C(s)g(t,x(t-\tau_0(t)),x(t-\tau_1(t)),
x(t+s))ds.
$$
By (P$^*_3$), for any $\varepsilon>0$, there exists $R_1>0$ such that
$|u|\leq R_1$ implies $|g(t,u)|\leq \varepsilon |u|$ for all $t\in
[0,\omega]$. Now let $x \in C_\omega$ with $\|x\|<R_1/3$. Then
$$|x(t-\tau_0(t))| + |x(t-\tau_1(t))| + |x(t+s)| \leq 3\|x\| < R_1.
$$
This yields
\begin{eqnarray*}
|f(t,x_t)| & \leq & \int^{+\infty}_{-\infty}|C(s)|g(t,x(t-\tau_0(t)),x(t-\tau_1(t)),
x(t+s))ds\\
& \leq & \int^{+\infty}_{-\infty}|C(s)|[\varepsilon (|x(t-\tau_0(t))|+
|x(t-\tau_1(t))| + |x(t+s)|)]ds\\
& \leq & \int^{+\infty}_{-\infty}|C(s)|ds (3\varepsilon \|x\|) = 3\varepsilon \|x\|
\end{eqnarray*}
and
$$ \frac{\int^\omega_0|f(t,x_t)|dt}{\|x\|} \leq 3\omega\varepsilon.
$$
Thus, (2.11) holds.  Next, it follows from the second part of (P$^*_3$)
that for any $M>0$ there
exists $R_2>0$ such that $[|u|\geq R_2,\: u_j\geq \sigma |u|,\: 0\leq j\leq
2]$
implies
\begin{equation}
g(t,u)\geq M|u|.
\end{equation}
Let $K$ be given in (2.4) and observe that if $x\in K$, then
$x(t-\tau_j(t))\geq \sigma \|x\|,\;|x(t+s)|\geq \sigma \|x\|$
for all $t\in \mathbb{R}$ and $j=0,1$, and
$$  |x(t-\tau_0(t))|+|x(t-\tau_1(t))| + |x(t+s)|\geq 3\sigma \|x\|
\geq 3\sigma\cdot R_2/(3\sigma) = R_2
$$
whenever $\|x\|\geq R_2/(3\sigma)$. By (3.7), we have
\begin{eqnarray*}
\lefteqn{g(t,x(t-\tau_0(t)),x(t-\tau_1(t)),x(t+s))}\\
& \geq & M(|x(t-\tau_0(t))|+|x(t-\tau_1(t))| + |x(t+s)|)
\geq M(3\sigma \|x\|).
\end{eqnarray*}
Using the above inequality, we get
\begin{eqnarray*}
\int^\omega_0|f(t,x_t)|dt &=&
\int^\omega_0\int^{+\infty}_{-\infty}|C(s)|g(t,x(t-\tau_0(t)),x(t-\tau_1(t)),
x(t+s))ds\,dt\\
&\geq & \int^\omega_0\int^{+\infty}_{-\infty}|C(s)|(3\sigma M\|x\|)ds\,dt
= 3\sigma\omega
M\|x\|.
\end{eqnarray*}
Thus,
$$\frac{\int^\omega_0|f(t,x_t)|dt}{\|x\|}\geq 3\sigma \omega M.$$
This proves (2.12) holds.  By Corollary \ref{coro2.1}, equation (3.6) has an
$\omega$-periodic solution. We omit the case when (P$^*_4$) holds. This
completes the proof.

\begin{thebibliography}{00} \frenchspacing

\bibitem{b1}
T. A. Burton, Volterra Integral and Differential equations,
Academic Press, New York, 1983.
\bibitem{b2}
T. A. Burton and B. Zhang, Uniform ultimate boundedness and periodicity
in functional differential equations, Tohoku Math J. 42(1990), 93-100.
\bibitem{c1}
M. Cavani, M. Lizana, and H. L. Smith, Stable periodic orbits for a predator-
prey model with delay, J. Math. Anal. Appl. 249(2000), 324-339.
\bibitem{c2}
S. -N. Chow, Existence of periodic solutions of autonomous functional
differential equations, 15(1974), 350-378.
\bibitem{c3}
J. M. Cushing, Integro-Differential Equations with Delay Models in Population
Dynamics, Lecture Notes in Biomath. 20, Springer-Verlag, Berlin, 1977.
\bibitem{d1}
K. Deimling, Nonlinear Functional Analysis, Springer-Verlag, New York, 1985.
\bibitem{g1}
K. Gopalsamy, Stability and Oscillations in Delay Differential Equations of
Population Dynamics, Kluwer Academic Publishers, London, 1992.
\bibitem{h1}
K. Hadeler, On the stability of the stationary state of a population growth
equation with time-lag, J. Math. Biology 3(1976), 197-201.
\bibitem{h2}
L. Hatvani and T. Krisztin, On the existence of periodic solutions for linear
inhomogeneous and quasilinear functional differential equations, J.
Differential equations 97(1992), 1-15.
\bibitem{h3}
Y. Hino, S. Murakami, and T. Naito, Functional Differential Equations with
Infinite Delay, Springer-Verlag, New York, 1991.
\bibitem{k1}
M. A. Krasnoselskii, Positive Solutions of Operator Equations, Gorninggen,
 Noordhoff, 1964.
\bibitem{k2}
Y. Kuang, Delay Differential Equations with Application in Population
Dynamics, Academic Press, New York, 1993.
\bibitem{m1}
G. Makay, Periodic solutions of dissipative functional differential equations,
Tohoku Math. J. 46(1994), 417-426.
\bibitem{m2}
R. M. May, Stability and Complexity in Model Ecosystems,
Princeton University Press, 1973.
\bibitem{m3}
R. M. May and W. J. Leonard, Nonlinear aspect of competition between three
species, SIAM J. Appl. Math. 29(1975), 243-253.
\bibitem{r1}
H. L. Royden, Real Analysis, Macmillan Publishing Company, New York, 1988.
\bibitem{s1}
W. Stech, The effect of time lags on the stability of the equilibrium state of
a population growth equation, J. Math. Biology 5(1978), 115-120.
\bibitem{w1}
M. Wazewska-Czyzewska and A. Lasota, Mathematical models of the red cell
system, Matematyta Stosowana 6(1976), 25-40.
\bibitem{w2}
P. Weng and M. Liang, The existence and behavior of periodic solutions of a
Hematopoiesis model, Mathematica Applicata 8(1995), 434-439.

\end{thebibliography}

\noindent\textsc{Daqing Jiang} (e-mail: daqingjiang@vip.163.com)\\
\textsc{Junjie Wei} (e-mail: weijj@nenu.edu.cn)\\[2pt]
Department of Mathematics, Northeast Normal University\\
Changchun, Jilin 130024,  China \smallskip

\noindent\textsc{Bo Zhang} \\
Department of Mathematics and Computer Science\\
Fayetteville State University \\
Fayetteville, NC 28301, USA \\
e-mail: bzhang@uncfsu.edu

\end{document}
