Age distribution and the stability of simple discrete time population models

Differences in age specific demographic characteristics can considerably alter the behaviour of the population dynamics of a species or community of species. In this analysis techniques are developed which enable the stability of the equilibria of a set of models involving age structure to be investigated. The underlying model in all cases is a simple matrix representation of first order difference equations. The analysis enables the results of computer investigations of various population models of this type to be explained.     

Stability and bifurcation in age-dependent population dynamics

A general form for an integral cohort model reflecting Easterlin's hypothesis is constructed. It is found that the basic dynamics of the model can be determined by a knowledge of the slope of the net reproduction rate at the equilibrium. Several forms for the net reproduction rate are examined numerically, and the resulting bifurcation processes are compared with those of the related first order difference equation modeling a two age group population. The results suggest that the situation usually described as “chaos” is likely to be less prevalent than earlier studies of first and second order difference equations have suggested.     

Global stability in a delayed partial differential equation describing cellular replication

Here we consider the dynamics of a population of cells that are capable of simultaneous proliferation and maturation. The equations describing the cellular population numbers are first order partial differential equations (transport equations) in which there is an explicit temporal retardation as well as a nonlocal dependence in the maturation variable due to cell replication. The behavior of this system may be considered along the characteristics, and a global stability condition is proved.     

Properties of small neural networks

Networks containing neuronal models of the type considered in the previous paper can be described by a set of first order differential equations. Steady-state solutions and the stability of these solutions to small perturbations can be obtained. Networks of physiological interest which give rise to second, third and fourth order linear equations are analysed in detail. Conditions are derived under which such networks can be condensed into a single neuron of similar order. Simple mechanisms for memory storage, for the generation of oscillatory activity and for decision making in neural systems are suggested.     

Existence and stability of equilibria in age-structured population dynamics

The existence of positive equilibrium solutions of the McKendrick equations for the dynamics of an age-structured population is studied as a bifurcation phenomenon using the inherent net reproductive rate n as a bifurcation parameter. The local existence and uniqueness of a branch of positive equilibria which bifurcates from the trivial (identically zero) solution at the critical value n=1 are proved by implicit function techniques under very mild smoothness conditions on the death and fertility rates as functional of age and population density. This first requires the development of a suitable linear theory. The lowest order terms in the Liapunov-Schmidt expansions are also calculated. This local analysis supplements earlier global bifurcation results of the author. The stability of both the trivial and the positive branch equilibria is studied by means of the principle of linearized stability. It is shown that in general the trivial solution losses stability as n increases through one while the stability of the branch solution is stable if and only if the bifurcation is supercritical. Thus the McKendrick equations exhibit, in the latter case, a standard exchange of stability with regard to equilibrium states as they depend on the inherent net reproductive rate. The derived lower order terms in the Liapunov-Schmidt expansions yield formulas which explicitly relate the direction of bifurcation to properties of the age-specific death and fertility rates as functionals of population density. Analytical and numerical results for some examples are given which illustrate these results.     

On Gompertz growth model and related difference equations

Within the context of the dynamics of populations described by first order difference equations a datailed study of the Gompertz growth model is performed. This is mainly achieved by proving several theorems for a class of difference equations generalizing the Gompertz equation. Some interesting features of the discrete Gompertz model, not exhibited by other well known growth models, are finally pointed out.     

Stability and synchronization analysis of inertial memristive neural networks with time delays

This paper is concerned with the problem of stability and pinning synchronization of a class of inertial memristive neural networks with time delay. In contrast to general inertial neural networks, inertial memristive neural networks is applied to exhibit the synchronization and stability behaviors due to the physical properties of memristors and the differential inclusion theory. By choosing an appropriate variable transmission, the original system can be transformed into first order differential equations. Then, several sufficient conditions for the stability of inertial memristive neural networks by using matrix measure and Halanay inequality are derived. These obtained criteria are capable of reducing computational burden in the theoretical part. In addition, the evaluation is done on pinning synchronization for an array of linearly coupled inertial memristive neural networks, to derive the condition using matrix measure strategy. Finally, the two numerical simulations are presented to show the effectiveness of acquired theoretical results.     

Global stability in discrete population models with delayed-density dependence

We address the global stability issue for some discrete population models with delayed-density dependence. Applying a new approach based on the concept of the generalized Yorke conditions, we establish several criteria for the convergence of all solutions to the unique positive steady state. Our results support the conjecture stated by Levin and May in 1976 affirming that the local asymptotic stability of the equilibrium of some delay difference equations (including Ricker's and Pielou's equations) implies its global stability. We also discuss the robustness of the obtained results with respect to perturbations of the model.     

Stability results for delayed-recruitment models in population dynamics

This paper relates the stability properties of a class of delay-difference equations to those of an associated scalar difference equation. Simple but powerful conditions for testing global stability are presented which are independent of the length of the time delay involved. For models which do not have globally stable equilibria, estimates of stability regions are obtained. Some well known baleen whale models are used to illustrate the results.     

A new criterion for the global stability of simultaneous cell replication and maturation processes

 We analyze a population model of cells that are capable of simultaneous and independent proliferation and maturation. This model is described by a first order partial differential equation with a time delay and a retardation of the maturation variable, both due to cell replication. We provide a general criterion for global stability in such equations. Received: 26 August 1996 / Revised version: 22 March 1997     

The numerical solution of stiff differential equations

This paper first discusses the conditions in which a set of differential equations should give stable solutions, starting with linear systems assuming that these do not differ greatly in this respect from non-linear systems. Methods of investigating the stability of particular systems are briefly discussed. Most real biochemical systems are known from observation to be stable, but little is known of the regions over which stability persists; moreover, models of biochemical systems may not be stable, because of inaccurate choice of parameter values.The separate problem of stability and accuracy in numerical methods of approximating the solution of systems of non-linear equations is then treated. Stress is laid on the consistently unsatisfactory results given by explicit methods for systems containing "stiff" equations, and implicit multistep methods are particularly recommended for this class of problem, which is likely to include many biochemical model systems. Finally, an iteration procedure likely to give convergence both in multistep methods and in the steady-state approach is recommended, and areas in which improvement in methods is likely to occur are outlined.     

Simulations in a Plankton Model

A technique for analyzing the stability of a zooplankton, phytoplankton, nutrient interaction model is described. This is an extension of a two-dimensional predator-prey model, incorporating a nutrient food source. In this study difference equations, rather than differential equations, are used to simulate the system since the systems are more easily studied in this formulation. Using only a few reasonable constraints on the system, a remarkable stability is obtained.     

Transition state of the glycolytic pathway under FDP saturating conditions: experimental studies and a theoretical model

1. The transition state of the glycolytic pathway, under FDP saturating conditions, from no ADP to ADP-saturating levels, is studied in a metabolic model in vitro obtained from rat skeletal muscle. 2. When ADP is absent from the reaction mixture a steady state for NADH concentration is observed. After ADP addition, a new steady state is reached. The transition state from the first steady state to the second one shows a pulse of NADH. Both the profile and the size of this pulse depend on the enzyme concentration. 3. A kinetic model of the lower part of glycolysis (after PFK reaction) is proposed, and this is described by a set of first order coupled nonlinear differential equations. The results obtained through stability analysis and numerical integration of these equations agree with the experimental ones. 4. The possible role of the above mentioned transition state on the transmitter mechanism of glycolytic oscillations from PFK to the lower part of the glycolysis is discussed.     

具有连续变量的多时滞二阶非线性中立型差分方程的振动性

用分析的方法研究了一类具有连续变量的多时滞二阶非线性中立型差分方程解的振动性,给出了该类方程解振动和差分算子振动的几个充分条件.     

Dynamic reduction with applications to mathematical biology and other areas

In a difference or differential equation one is usually interested in finding solutions having certain properties, either intrinsic properties (e.g. bounded, periodic, almost periodic) or extrinsic properties (e.g. stable, asymptotically stable, globally asymptotically stable). In certain instances it may happen that the dependence of these equations on the state variable is such that one may (1) alter that dependency by replacing part of the state variable by a function from a class having some of the above properties and (2) solve the 'reduced' equation for a solution having the remaining properties and lying in the same class. This then sets up a mapping Τ of the class into itself, thus reducing the original problem to one of finding a fixed point of the mapping. The procedure is applied to obtain a globally asymptotically stable periodic solution for a system of difference equations modeling the interaction of wild and genetically altered mosquitoes in an environment yielding periodic parameters. It is also shown that certain coupled periodic systems of difference equations may be completely decoupled so that the mapping Τ is established by solving a set of scalar equations. Periodic difference equations of extended Ricker type and also rational difference equations with a finite number of delays are also considered by reducing them to equations without delays but with a larger period. Conditions are given guaranteeing the existence and global asymptotic stability of periodic solutions.     

含阻尼项二阶泛函差分方程的振动性质

从解的渐近状态入手,应用分类讨论方法和Raccati技巧,讨论了一类广泛的二阶泛函差分方程解的振动性与非振动性,建立了四个新的振动性定理,推广并改进了已有文献中的相关结果.     

The geometrical analysis of a predator-prey model with two state impulses

Using successor functions and Poincaré-Bendixson theorem of impulsive differential equations, the existence of periodical solutions to a predator-prey model with two state impulses is investigated. By stability theorem of periodic solution to impulsive differential equations, the stability conditions of periodic solutions to the system are given. Some simulations are exerted to prove the results.     

Boundedness,persistence and stability for classes of forced difference equations arising in population ecology

Journal of Mathematical Biology - Boundedness, persistence and stability properties are considered for a class of nonlinear, possibly infinite-dimensional, forced difference equations which arise...     

Non-linear bioconvection in a deep suspension of gyrotactic swimming micro-organisms

 The non-linear structure of deep, stochastic, gyrotactic bioconvection is explored. A linear analysis is reviewed and a weakly non-linear analysis justifies its application by revealing the supercritical nature of the bifurcation. An asymptotic expansion is used to derive systems of partial differential equations for long plume structures which vary slowly with depth. Steady state and travelling wave solutions are found for the first order system of partial differential equations and the second order system is manipulated to calculate the speed of vertically travelling pulses. Implications of the results and possibilities of experimental validation are discussed. Received: 26 May 1997 / Revised version: 10 May 1998     

On the Equivalence of Mathematical Models For Cell Proliferation Kinetics

Various types of mathematical models, such as partial differential equations, ordinary differential equations and difference equations, are available in the literature to describe the kinetics of cell proliferation, and different studies of cell kinetic phenomena have been conducted using these models. This paper discusses the equivalence between the different models identifying the conditions and approximations under which one type of models may be derived from another. Such an equivalence study is highly useful for an integration of the diverse results that have been obtained using different models in order to gain a more complete understanding of cell kinetic phenomena.