%% Created by Maple 7.00 (IBM INTEL LINUX)
%% Source Worksheet: PhasePortrait.mws
%% Generated: Mon Apr 28 16:13:21 2003
%\documentclass{article}
%\usepackage{maple2e}
% \def\emptyline{\vspace{12pt}}
%\DefineParaStyle{Maple Plot}
%\DefineParaStyle{Warning}
%\begin{document}
\pagestyle{empty}
\begin{maplegroup}
\begin{mapleinput}
\mapleinline{active}{1d}{restart:}{%
}
\end{mapleinput}

\end{maplegroup}
\begin{maplegroup}
\begin{mapleinput}
\mapleinline{active}{1d}{with(DEtools): with(plots): with(plottools):}{%
}
\end{mapleinput}

\mapleresult
\begin{maplettyout}
Warning, the name changecoords has been redefined
\end{maplettyout}

\begin{maplettyout}
Warning, the names arrow and translate have been redefined
\end{maplettyout}

\end{maplegroup}
\begin{maplegroup}
\begin{flushleft}
In this example, we are looking at x'=y,y'=-x(1+y), with two sets of
initial conditions: (x(0),y(0))=(1,0), (x(0),y(0))=(2,0).
\end{flushleft}

\end{maplegroup}
\begin{maplegroup}
\begin{mapleinput}
\mapleinline{active}{1d}{Digits := 10:}{%
}
\end{mapleinput}

\end{maplegroup}
\begin{maplegroup}
\begin{mapleinput}
\mapleinline{active}{1d}{plot1 :=
DEplot([diff(x(t),t)=x(t)-6*y(t),diff(y(t),t)=2*x(t)-7*y(t)],
[x(t),y(t)], t=-5..5, scene=[x(t),y(t)],
[[x(0)=1,y(0)=0],[x(0)=-4,y(0)=0],[x(0)=0,y(0)=1],[x(0)=4,y(0)=0]],
x=-4..4, y=-4..4, linecolor=[red,blue,green,black], arrows=SMALL,
method=rkf45, stepsize=0.05):}{%
}
\end{mapleinput}

\end{maplegroup}
\begin{maplegroup}
\begin{mapleinput}
\mapleinline{active}{1d}{display(\{plot1\}, title="Phase Portrait", titlefont = [HELVETICA,
BOLD, 12], labels = ["x","y"], labelfont = [HELVETICA, BOLD, 10]);}{%
}
\end{mapleinput}

\mapleresult
\begin{center}
\mapleplot{hw25_1301.eps}
\end{center}

\end{maplegroup}
\begin{maplegroup}
\begin{flushleft}
For a plot of x(t) only....
\end{flushleft}

\end{maplegroup}
\begin{maplegroup}
\begin{mapleinput}
\end{mapleinput}

\end{maplegroup}
%\end{document}
%% End of Maple 7.00 Output
