Abstract
Investigation of the Runge-Kutta s and Burlisch-Stoer s numerical methods for EDOs, in their applications on the chaotic behavior of the Van der Pol s and Henon-Heiles equations, which will also be studied with the Symplectic integrator scheme.