<p>The dynamics of tax evasion are modeled as a system of nonlinear equations, which is first transformed into an arbitrary order system of nonlinear equations in the Caputo sense. The transformed system is then solved using the Genocchi wavelet method. The Genocchi wavelet operational matrix is derived by utilizing block pulse functions. Comprehensive analyses are conducted, including convergence, error bounds, stability, existence, uniqueness, and boundedness of the solutions. Graphical visualizations of the solutions are presented in both two-dimensional and three-dimensional forms. The order of the fractional nonlinear system is varied, and a comparative analysis between the Genocchi wavelet method and the Adams-Bashforth-Moulton method is performed. Absolute error analyses are presented through graphs as well as tables to show the efficacy of the Genocchi wavelet approach.</p>