Phase Portrait
> IC := [ x(0) = 0.4, y(0) = 0.2 ],
> [ x(0) = 0.4, y(0) = 0.3 ],
> [ x(0) = 0.4, y(0) = 0.4 ],
> [ x(0) = 0.4, y(0) = 0.6 ],
> [ x(0) = 2, y(0) = 1 ],
> [ x(0) = 2, y(0) = 2 ],
> [ x(0) = 2, y(0) = 3 ],
> [ x(0) = 2, y(0) = 4 ] ;
>
> PPplot := DEplot( MODEL, VAR, DOMAIN, [ IC ], RANGE,
> scene=[ x, y ], arrows=NONE,
> title = `Phase Portrait for Competing Species` ):
> PPplot;
>
> display( [ DFplot, PPplot ], title=`Direction Field and Phase Portrait` );
> display( [ NULLplot, DFplot, PPplot ], title=`Nullclines, Direction Field, and Phase Portrait` );
>