By comparing with symplectic different methods, the quadratic element is an approximately symplectic method which can keep high accuracy approximate of symplectic structure for Hamiltonian chaos, and it is also energy conservative when there have chaos phenomenon. We use the quadratic finite element method to solve the H[Formula: see text]non–Heiles system, and this method was never used before. Combining with Poincar[Formula: see text] section, when we increase the energy of the systems, KAM tori are broken and the motion from regular to chaotic. Without chaos, three kinds of methods to calculate the Poincar[Formula: see text] section point numbers are the same, and the numbers are different with chaos. In long-term calculation, the finite element method can better keep dynamic characteristics of conservative system with chaotic motion.