I have tried for a long time to find a Lyapunov function for the specific problem
$$\begin{align} x' &= -x - 2y + xy^2 \\ y' &= 3x - 3y + y^3 \end{align}$$
Do you know any Lyapunov function what is going to work to determine the stability of the origin? According to a phase portrait produced by pplane the stability of the origin is asymptotically stable.
Thanks in advance!