A fast sampling method for estimating the domain of attraction