Asymptotic stability regions of non-linear dynamical systems are estimated by using homogeneous polynomial Lyapunov functions, which include quadratic Lyapunov functions as a special case. The non-linear system is represented as a convex hull of some linear systems in a specified region of the state space. It is shown that the estimation can be formulated as an LMI optimization problem in an augmented state space.