Equilibrium Point #1:

> RANGE1 := x=-0.1..0.1, y=-0.1..0.1;

> DEplot( MODEL, VARS, DOMAIN, RANGE1, arrows=MEDIUM, title=`Direction Field Near (0,0)` );

>

The matrix for the linearized system at this equilibrium point is:

> A := subs( x=0, y=0, evalm(J) );

>

The eigenvalues of this matrix are:

> eigenvals( A );

Since both eigenvalues are positive, the origin is an unstable source for both the linearized and nonlinear system. This is consistent with the direction field that focuses on the portion of the phase portrait close to this equilibrium solution. (Can you see the two straight line solutions?)

>