>

>	# (Anything after `#' is for comments. Maple does not treat them as input.) This is a tutorial to some common commands on solving differential equations.

>	with(plots): # Load up the graphic package. End any command with `:' if you don't want to see the output. Otherwise, end it with `;'.

Warning, the name changecoords has been redefined

>	with(DEtools): # Load up another package for differential equations.

>

?dsolve;

>	# Asking question in the form ``?[subject]'' is a very useful tool to learn new commands and techniques.

>	ode:=diff(y(t),t$2)+2*diff(y(t),t)+5*y(t)=cos(t);

>	dsolve(ode,y(t));

>	ic:=y(0)=10,D(y)(0)=-1;

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

>	a:=odeplot(solu,[t,y(t)],0..10,color=black):

>	display({a});# To plot two or more graphs in one plot, use ``display({a,b})'' with ``b'' any another assigned plot variable similar to ``a'' above. You may use plot package, not necessarily the same as ``odeplot'', to generate ``b''.

>	# The following example shows how to solve a system numerically and plot them together against the t-axis.

>	sys := {diff(x(t),t)=y(t),diff(y(t),t)=-x(t),x(0)=1,y(0)=2}; # Put the equations and initial conditions together.

>	syssolu:=dsolve(sys,{x(t),y(t)},type=numeric):

>	aa:=odeplot(syssolu,[t,x(t)],0..10,color=black):

>	bb:=odeplot(syssolu,[t,y(t)],0..10,color=red):

>	display({aa,bb});

>	?dfieldplot; # Find out how to plot direction field.

>	DEplot(diff(y(t),t)=y(t)*(1-y(t)),y(t),t=0..10,{[y(0)=0.5],[y(0)=2]}); #Slope field with solutions.

>	phaseportrait([D(x)(t)=-x(t)-y(t),D(y)(t)=2*x(t)-y(t)], [x(t),y(t)],t=0..20,[[x(0)=1,y(0)=0],[x(0)=3,y(0)=3]],x=-4..4,y=-4..4,stepsize=.05, scene=[x(t),y(t)],linecolour=sin(t*Pi/2)); # Here is another way to do vector field with solutions.

>