# 1

> ode1 := diff( x(t), t ) = -3*x(t) :

> ode2 := diff( y(t), t ) = x(t) - 3 * y(t) :

> MODEL := [ ode1, ode2 ] ;

> VAR := { x(t), y(t) } ;

>

c) Direction Field

d) Phase Plane

e) x- and y-solution curves for the IVP