For a locally stable polynomial dynamical system its region of attraction can be estimated by a polynomial Lyapunov function level set. We formulate this problem in terms of global minimization of a polynomial function over single polynomial constraint. We describe a simple Matlab tool, which employs numerical algebraic geometry methods for computing all local solutions of the optimization problem and therefore its global solution.