This paper considers the solution behavior and dynamical properties of the variable-order fractional Newton–Leipnik system defined via Liouville–Caputo derivatives of variable order. In contrast to integer-order models, the presence of variable-order fractional operators in the Newton–Leipnik structure enriches the model by providing memory-dependent effects that vary with time; hence, it is capable of a broader and more flexible range of nonlinear responses. Numerical simulations have been conducted to study how different order functions influence the trajectory and qualitative dynamics: clear transitions in oscillatory patterns have been identified by phase portraits, time-series profiles, and three-dimensional state evolution. The work goes further by considering the development of bifurcations and chaotic regimes and stability shifts and confirms the occurrence of several phenomena unattainable in fixed-order and/or integer-order formulations. Analysis of Lyapunov exponents confirms strong sensitivity to the initial conditions and further details how the memory effects either reinforce or prevent chaotic oscillations according to the type of order function. The results, in fact, show that the variable-order fractional Newton–Leipnik framework allows for more expressive and realistic modeling of complex nonlinear phenomena and points out the crucial role played by evolving memory in controlling how the system moves between periodic, quasi-periodic, and chaotic states.
Saadeh et al. (Fri,) studied this question.