This paper is motivated by the series of research papers that consider parasitoids’ external input upon the host–parasitoid interactions. We explore a class of host–parasitoid models with variable release and constant release of parasitoids. We assume that the host population has a constant rate of increase, but we do not assume any density dependence regulation other than parasitism acting on the host population. We compare the obtained results for constant stocking with the results for proportional stocking. We observe that under a specific condition, the release of a constant number of parasitoids can eventually drive the host population (pests) to extinction. There is always a boundary equilibrium where the host population extinct occurs, and the parasitoid population is stabilized at the constant stocking level. The constant and variable stocking can decrease the host population level in the unique interior equilibrium point; on the other hand, the parasitoid population level stays constant and does not depend on stocking. We prove the existence of Neimark–Sacker bifurcation and compute the approximation of the closed invariant curve. Then we consider a few host–parasitoid models with proportional and constant stocking, where we choose well-known probability functions of parasitism. By using the software package Mathematica we provide numerical simulations to support our study.

This paper is motivated by the series of research papers that consider parasitoids’ external input upon the host–parasitoid interactions. We explore a class of host–parasitoid models with variable release and constant release of parasitoids. We assume that the host population has a constant rate of increase, but we do not assume any density dependence regulation other than parasitism acting on the host population. We compare the obtained results for constant stocking with the results for proportional stocking. We observe that under a specific condition, the release of a constant number of parasitoids can eventually drive the host population (pests) to extinction. There is always a boundary equilibrium where the host population extinct occurs, and the parasitoid population is stabilized at the constant stocking level. The constant and variable stocking can decrease the host population level in the unique interior equilibrium point; on the other hand, the parasitoid population level stays constant and does not depend on stocking. We prove the existence of Neimark–Sacker bifurcation and compute the approximation of the closed invariant curve. Then we consider a few host–parasitoid models with proportional and constant stocking, where we choose well-known probability functions of parasitism. By using the software package Mathematica we provide numerical simulations to support our study.

In this paper, we explore the dynamics of a certain class of Beddington host-parasitoid models, where in each generation a constant portion of hosts is safe from attack by parasitoids, and the Ricker equation governs the host population. Using the intrinsic growth rate of the host population that is not safe from parasitoids as a bifurcation parameter, we prove that the system can either undergo a period-doubling or a Neimark–Sacker bifurcation when the unique interior steady state loses its stability. Then, we apply the new theory to the following well-known cases: May’s model, [Formula: see text]-model, Hassel and Varley (HV)-model, parasitoid-parasitoid (PP) model and [Formula: see text] model. We use numerical simulations to confirm our theoretical results.

AbstractWe consider the second-order rational difference equation xn+1=γ+δxnxn−12,$$ {x_{n+1}=\gamma +\delta \frac{x_{n}}{x^{2}_{n-1}}}, $$ where γ, δ are positive real numbers and the initial conditions x−1$x_{-1}$ and x0$x_{0}$ are positive real numbers. Boundedness along with global attractivity and Neimark–Sacker bifurcation results are established. Furthermore, we give an asymptotic approximation of the invariant curve near the equilibrium point.

We consider the second-order rational difference equation xn+1=γ+δxnxn−12,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {x_{n+1}=\gamma +\delta \frac{x_{n}}{x^{2}_{n-1}}}, $$\end{document} where γ, δ are positive real numbers and the initial conditions x−1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$x_{-1}$\end{document} and x0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$x_{0}$\end{document} are positive real numbers. Boundedness along with global attractivity and Neimark–Sacker bifurcation results are established. Furthermore, we give an asymptotic approximation of the invariant curve near the equilibrium point.

Abstract A certain class of a host–parasitoid models, where some host are completely free from parasitism within a spatial refuge is studied. In this paper, we assume that a constant portion of host population may find a refuge and be safe from attack by parasitoids. We investigate the effect of the presence of refuge on the local stability and bifurcation of models. We give the reduction to the normal form and computation of the coefficients of the Neimark–Sacker bifurcation and the asymptotic approximation of the invariant curve. Then we apply theory to the three well-known host–parasitoid models, but now with refuge effect. In one of these models Chenciner bifurcation occurs. By using package Mathematica, we plot bifurcation diagrams, trajectories and the regions of stability and instability for each of these models.

In this paper, we consider the dynamics of a certain class of host-parasitoid models, where some hosts are completely free from parasitism either with or without a spatial refuge and the host population is governed by the Beverton–Holt equation. We assume that, in each generation, a constant portion of the host population may find a refuge and be safe from the attack by parasitoids. We derive some criteria for the Neimark–Sacker bifurcation. Then, we apply the developed theory to the three well-known cases: [Formula: see text] model, Hassel and Varley model, and parasitoid–parasitoid model. Intensive numerical calculations suggest that the last two models undergo a supercritical Neimark–Sacker bifurcation.

By using KAM theory we investigate the stability of equilibrium points of the class of difference equations of the form xn+1=f(xn)xn−1,n=0,1,…\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$x_{n+1}=\frac{f(x _{n})}{x_{n-1}}, n=0,1,\ldots $\end{document} , f:(0,+∞)→(0,+∞)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$f:(0,+\infty )\to (0,+\infty )$\end{document}, f is sufficiently smooth and the initial conditions are x−1,x0∈(0,+∞)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$x_{-1}, x _{0}\in (0,+\infty )$\end{document}. We establish when an elliptic fixed point of the associated map is non-resonant and non-degenerate, and we compute the first twist coefficient α1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\alpha _{1}$\end{document}. Then we apply the results to several difference equations.

By using KAM theory we investigate the stability of equilibrium points of the class of difference equations of the form xn+1=f(xn)xn−1,n=0,1,…$x_{n+1}=\frac{f(x _{n})}{x_{n-1}}, n=0,1,\ldots $ , f:(0,+∞)→(0,+∞)$f:(0,+\infty )\to (0,+\infty )$, f is sufficiently smooth and the initial conditions are x−1,x0∈(0,+∞)$x_{-1}, x _{0}\in (0,+\infty )$. We establish when an elliptic fixed point of the associated map is non-resonant and non-degenerate, and we compute the first twist coefficient α1$\alpha _{1}$. Then we apply the results to several difference equations.

We investigate global dynamics of the following second order rational difference equation x n + 1 = x n x n − 1 + α x n + β x n − 1 a x n x n − 1 + b x n − 1 , where the parameters α , β , a , b are positive real numbers and initial conditions x − 1 and x 0 are arbitrary positive real numbers. The map associated to the right-hand side of this equation is always decreasing in the second variable and can be either increasing or decreasing in the first variable depending on the corresponding parametric space. In most cases, we prove that local asymptotic stability of the unique equilibrium point implies global asymptotic stability.

We investigate the local stability and the global asymptotic stability of the difference equation , with nonnegative parameters and initial conditions such that , for all . We obtain the local stability of the equilibrium for all values of parameters and give some global asymptotic stability results for some values of the parameters. We also obtain global dynamics in the special case, where , in which case we show that such equation exhibits a global period doubling bifurcation.