In this study, the performance of the Taylor’s decomposition method in the nine-dimensional Lorenz system is investigated. The proposed method is applied to the nine-dimensional Lorenz system for both chaotic and hyperchaotic cases. Phase portraits and time series graphs of the problem are drawn. The results obtained are compared with both the fourth order Runge-Kutta method and multi-domain compact finite difference relaxation method. Comparisons show that the obtained results are consistent with the results of other methods. Thus, the accuracy and efficiency of the method for nonlinear systems were emphasized by using the nine-dimensional Lorenz system.