Solving a system of nonlinear difference equations with bilinear dynamics

This paper presented a comprehensive study of a three-dimensional nonlinear system of difference equations, which can be reduced to a two-dimensional bilinear system. The system monitored the evolution of three sequences $\left(P_m\right),$ $\left(Q_m\right),$ $\left(R_m\right)$, governed by recursive relations. We investigated the solvability of this system and provided general closed-form solutions for various parameter conditions. Furthermore, the simulations provided valuable insights into the dynamic behavior of animals, modeled using recursive difference equations. The model encapsulated essential behavioral metrics, represented by the variables $P$, $Q$, and $R$, which corresponded to individual actions, social interactions, and environmental stressors, respectively. These variables adapted dynamically in response to internal and external influences, illustrating the system's sensitivity to various behavioral and environmental conditions.