Figure 1: Contour of Lyapunov function(colour), the domain D (red), the ROA (black) and the gradients (arrows)