Motivated by the article [Getz and Mills. Host–parasitoid coexistence and egg-limited encounter rates. Am Nat. 1996;148:301–315], in this paper, we explore a discrete model involving a host and a parasitoid with search and egg limitations and arbitrary host-escape function. We added proportional refuge for the hosts to the model. We focus on the system's behavior at the equilibrium points and the nearby regions. In addition to the topological classification of these points, we examined local behavior. In the case of the extinction equilibrium, we obtain the global result. We describe the dynamical behavior scenarios in the neighborhood of the non-isolated exclusion equilibrium point (1:1 resonant). For the unique coexisting equilibrium, we prove the emergence of the Neimark–Sacker bifurcation and calculate the first Lyapunov exponent. This bifurcation can be either super or sub-critical. We have also established the occurrence of the Chenciner bifurcation. Our findings indicate that proportional refuge may or may not stabilize the system, and the choice of host-escape function plays a crucial role in shaping the system's dynamics. We also provide numerical examples to support our theoretical results.

This research delves into the generalized Beddington host–parasitoid model, which includes an arbitrary parasitism escape function. Our analysis reveals three types of equilibria: extinction, boundary, and interior. Upon examining the parameters, we discover that the first two equilibria can be globally asymptotically stable. The boundary equilibrium undergoes period-doubling bifurcation with a stable two-cycle and a transcritical bifurcation, creating a threshold for parasitoids to invade. Furthermore, we determine the interior equilibrium’s local stability and analytically demonstrate the period-doubling and Neimark–Sacker bifurcations. We also prove the permanence of the system within a specific parameter space. The numerical simulations we conduct reveal a diverse range of dynamics for the system. Our research extends the results in [Kapçak et al., 2013] and applies to a broad class of the generalized Beddington host–parasitoid model.

This paper examines the relationship between herbivores and plants with a strong Allee effect. When the plant reaches a particular size, the herbivore attacks it. We use the logistic equation to model plant growth and analyze its behavior without herbivores before investigating their interactions. Our study investigates the equilibrium points and their stability, discovering that different fixed points can become unstable due to various bifurcations such as transcritical, saddle-node, period-doubling, and Neimark–Sacker bifurcations. We have identified the Allee threshold, which, if exceeded, can cause both populations to become extinct below that level. However, we have discovered a coexistence equilibrium that is locally asymptotically stable for a range of parameter values above that threshold. Our additional numerical simulations suggest that this area of stability can be expanded. Our results indicate that this system is highly responsive to its parameters. We compare our findings to those of a system without strong Allee effects and conduct numerical simulations to verify our results. By including the Allee effect in the plant population, we enrich the local and global dynamics of the system.

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.

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.

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,…$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.

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.