Solutions and local stability of the Jacobsthal system of difference equations

We presented a comprehensive theory for deriving closed-form expressions and representations of the general solutions for a specific case of systems involving Riccati difference equations of order $m+1$, as discussed in the literature. However, our focus was on coefficients dependent on the Jacobsthal sequence. Importantly, this system of difference equations represents a natural extension of the corresponding one-dimensional difference equation, uniquely characterized by its theoretical solvability in a closed form. Our primary objective was to demonstrate a direct linkage between the solutions of this system and Jacobsthal and Lucas-Jacobsthal numbers. The system's capacity for theoretical solvability in a closed form enhances its distinctiveness and potential applications. To accomplish this, we detailed offer theoretical explanations and proofs, establishing the relationship between the solutions and the Jacobsthal sequence. Subsequently, our exploration addressed key aspects of the Jacobsthal system, placing particular emphasis on the local stability of positive solutions. Additionally, we employed mathematical software to validate the theoretical results of this novel system in our research.