> # This is a sample worksheet for your homework assignment. Use commond `?dsolve[piecewise]' to find out how to modify this worksheet for your piecewise defined equation.

> with(DEtools):with(plots):

> mu:=10:

> sys:=D(x)(t)=mu*(y(t)-x(t)^3/3+x(t)),D(y)(t)=-x(t)/mu:# The van der Pol equation as a system of 1st order equations.

> solu:= dsolve({sys,x(0)=1,y(0)=.4},{x(t),y(t)},type=numeric):

> a:=odeplot(solu, [x(t),y(t)],0..20,numpoints=100):

> b:=plot([x, x^3/3-x, x=-3..3],-3..3,-1..1,color=black):

> display({a,b});

> xp:=odeplot(solu,[t,x(t)],0..50,color=red):yp:=odeplot(solu,[t,y(t)],0..50,color=green):

> display({xp,yp});

> phaseportrait([D(x)(t)=mu*(y(t)-x(t)^3/3+x(t)),D(y)(t)=-x(t)/mu], \
[x(t),y(t)],t=0..20,[[x(0)=1,y(0)=0]],stepsize=.05, \
scene=[x(t),y(t)],linecolour=sin(t*Pi/2)); # Phaseportrait.

> ode:=diff(y(t),t\$2)+mu*(y(t)^2-1)*diff(y(t),t)+y(t)=0: # You can solve it as a 2nd order equation as well.

> ic:=y(0)=1,D(y)(0)=0:

> solu2:=dsolve({ode,ic},y(t),type=numeric):

> solu2(2); # Check the solution as a given time t=2.

> odeplot(solu2,[t,y(t)],0..50,color=black);

>