Finding a region of attraction using sum of squares optimization

This is a question for @blegat