I don't think you will have much luck in finding a Lyapunov function for the general nonlinear case. You could try only if the known theorem about the eigenvalues of the linearized problem does not apply.

Which is not the case here. If we call and , the eigenvalues of are negative, and so

the origin is asymptotically stable.