Search

Abstract: This paper continues the work contained in two previous papers of the authors, devoted to the qualitative study of the dynamical system generated by a delay differential equation that models leukemia. The problem depends on five parameters and has two equilibria. As already known, at the non-zero equilibrium solution, for certain values of the parameters, Hopf bifurcation occurs. Our aim here is to investigate the Bautin bifurcation for the considered model. By using a second order approximation of the center manifold, we compute the first Lyapunov coefficient for Hopf bifurcation points. We explore the space of parameters in a zone of biological interest and find, by direct computation, points where the first Lyapunov coefficient equals zero. These are candidates for Bautin type bifurcation. For these points we compute the second Lyapunov coefficient, that determines the type of Bautin bifurcation. The computation of the second Lyapunov coefficient requires a fourth order approximation of the center manifold, that we determine. The points with null first Lyapunov coefficient obtained are given in tables and are also plotted on surfaces of Hopf bifurcation obtained by fixing two parameters. Next we vary two parameters around a point in the space of parameters with and numerically explore the behavior of the solution. The results confirm the Bautin type-bifurcation. [Copyright &y& Elsevier]

Abstract: In this paper, we study the influence of time delays on the dynamics of the classical Gompertz model. We consider the models with one discrete delay introduced in two different ways and the model with two delays which generalise those with one delay. We study asymptotic behaviour and bifurcations with respect to the ratio of delays . Our results show that in such model with two delays there is only one stability switch and for a threshold value of bifurcation parameter, Hopf bifurcation (HB) occurs. However, the type of HB, and therefore its stability (i.e. stability of periodic orbits arising due to it), strongly depends on the magnitude of . The function describing stability of HB is periodic with respect to . Within one period of length 4 five changes of HB stability are observed. We also introduce the second model with two delays which has a better biological interpretation than the first one. In that model several stability switches can occur, depending on the model parameters. We illustrate analytical results on the example of tumour growth model with parameters estimated on the basis of experimental data. [Copyright &y& Elsevier]