Solving hyperchaotic systems using the spectral relaxation method

A new multistage numerical method based on blending a Gauss-Siedel relaxation method and Chebyshev pseudospectral method, for solving complex dynamical systems exhibiting hyperchaotic behavior, is presented. The proposed method, called the multistage spectral relaxation method (MSRM), is applied for the numerical solution of three hyperchaotic systems, namely, the Chua, Chen, and Rabinovich-Fabrikant systems. To demonstrate the performance of the method, results are presented in tables and diagrams and compared to results obtained using a Runge-Kutta-(4,5)-based MATLAB solver, ode45, and other previously published results.