d) Phase Plane
> DOMAIN := t = 0 .. 1 ;
> WINDOW := x = -1 .. 1,
> y = -1 .. 1 ;
> IC := [ seq( seq( [x(0)=i, y(0)=j/3 ], j=-3..3 ), i=[-1,1] ),
> seq( seq( [x(0)=i/3,y(0)=j ], j=[-1,1] ), i=-3..3 ) ];
> PHASEplot := DEplot( MODEL, VAR, DOMAIN, WINDOW, IC, arrows=THIN, arrows = MEDIUM,
> title = `Phase Plane for #1` ):
> PHASEplot;
>