We consider a system dtdx=r1(t)G1(x)[∫h1(t)tf1(y(s))dsR1(t,s)−x(t)],dtdy=r2(t)G2(y)[∫h2(t)tf2(x(s))dsR2(t,s)−y(t)] with
increasing functions f1 and f2, which has at most one positive
equilibrium. Here the values of the functions ri,Gi,fi are positive for
positive arguments, the delays in the cooperative term can be distributed and
unbounded, both systems with concentrated delays and integro-differential
systems are a particular case of the considered system. Analyzing the relation
of the functions f1 and f2, we obtain several possible scenarios of the
global behaviour. They include the cases when all nontrivial positive solutions
tend to the same attractor which can be the positive equilibrium, the origin or
infinity. Another possibility is the dependency of asymptotics on the initial
conditions: either solutions with large enough initial values tend to the
equilibrium, while others tend to zero, or solutions with small enough initial
values tend to the equilibrium, while others infinitely grow. In some sense
solutions of the equation are intrinsically non-oscillatory: if both initial
functions are less/greater than the equilibrium value, so is the solution for
any positive time value. The paper continues the study of equations with
monotone production functions initiated in [Nonlinearity, 2013, 2833-2849].Comment: 18 pages, 3 figures, published in 2015 in Nonlinearit