Equilibrium Solutions

> EQUILeqns := { rhs( ode1 ) = 0, rhs( ode2 ) = 0 } ;

> EQUILsoln := solve( EQUILeqns, convert(VARS,set) ) ;

Notice that there are four equilibria - since there are four intersections of an - and a -nullcline, this is consistent with the nullclines

>

> display( [DIRplot, NULLplot], title=`Direction Field and Nullclines` );

>