Hi. You have to analyze it's phase space right? Linearize it, calculate the eigenvalues, determine the type of fixed points whether they're sources, centers, saddles, sink, and then based on the initial conditions, conclude what the long term behavior is. For example, if I have initial conditions close to a sink, I'd expect the long term behavior to be drawn into the sink and reach equilibrium. You can just do a few runs: choose any old initial conditions, run it through a numerical integrator and compare the numerical results to it's phase space portrait.
I'll do the equilibrium point at the origin:
The Jacobian is
The eigenvalue equation then is:
The characteristic equation then is
The eigenvalues are thus 4 and 9. Both positive. This means the origin is a source: Initial points near it will move away from it.
Now need to to the others and see how they affect the dynamcis.