>
BDH2-2.mws
This worksheet was developed using template-SYS.mws. Be sure you can use the template to create similar plots and obtain appropriate information from other problems.
>
Getting started
Example on pages 164-165
> ode1 := diff( x(t), t ) = 3*x(t)+y(t) ;
> ode2 := diff( y(t), t ) = x(t)-y(t) ;
> MODEL := { ode1, ode2 } ; # define the differential equation model
> IC := [ x(0) = -0.7, y(0) = 3 ] ; # specify an initial condition in the form h(0)=200, u(0)=5
> VAR := { x(t), y(t) } ; # identify the variables in the model in the form h(t), u(t)
>
Equilibrium Solutions
Direction Field, Phase Portrait, and Solution Curves ( DEplot )
Competing Species Model (pages 166-168)
> ode1 := diff( x(t), t ) = 2*x(t)*(1-x(t)/2) - x(t)*y(t) ;
> ode2 := diff( y(t), t ) = 3*y(t)*(1-y(t)/3) - 2*x(t)*y(t) ;
> MODEL := { ode1, ode2 } ; # define the differential equation model
> #IC := [ x(0) = -0.7, y(0) = 3 ] ; # specify an initial condition in the form h(0)=200, u(0)=5
> VAR := { x(t), y(t) } ; # identify the variables in the model in the form h(t), u(t)
>
Nullclines
Equilibrium Solutions
Direction Field, Phase Portrait, and Solution Curves ( DEplot )
> restart;