In this work, the AutonomousSystems4D package is presented, which allows the qualitative analysis of non-linear differential equation systems in four dimensions, as well as drawing the phase surfaces by immersing R4 in R3. The package is programmed in the computational tool Octave. As a case study applied to the new Lorenz 4D System, sensitivity was found in the initial conditions, Lyapunov exponents, Kaplan Yorke dimension, a stable and unstable critical point, limit cycle, Hopf bifurcation, and hyperattractors. The package could be adapted to perform qualitative analysis and visualize phase surfaces to autonomous systems, e.g. Sprott 4D, Rossler 4D, etc. The package can be applied to problems such as: design, analysis, implementation of electronic circuits; to message encryption.