Stability analysis is a critical part of the study of a dynamical system that includes: finding the attractors of the system, determining their stability or instability, and identifying basins of attraction (BoAs) of the attractors. There are some classical methods for finding parts of BoAs, like Lyapunov-based and Non-Lyapunov-based methods; however, due to the limitations of the current methods, the identification of a complex boundary of a BoA is not feasible in addition to the high computational load of these methods. The concept of Lyapunov exponents is a powerful tool for performing stability analysis of highly nonlinear systems. In this paper, a new and intelligent approach for identifying the boundary of the entire BoA using Lyapunov exponents, Monte Carlo techniques, and support vector machine learning algorithm is developed. The proposed approach can identify the boundary of BoAs due to the intelligent sampling brought by Monte Carlo technique, which also dramatically reduces the computational load. Three illustrative examples are presented to demonstrate the effectiveness of the proposed approach.