{VERSION 3 0 "IBM INTEL NT" "3.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 }{CSTYLE "" 0 21 "" 0 1 0 0 0 1 0 0 0 0 2 0 0 0 0 }{CSTYLE "
" -1 256 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 257 "" 0 1 0 
0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 0 1 0 0 0 0 
0 0 0 }{CSTYLE "" -1 259 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{PSTYLE "No
rmal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 
-1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 
-1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 }1 0 0 -1 -1 -1 0 0 0 0 
0 0 -1 0 }{PSTYLE "Warning" 2 7 1 {CSTYLE "" -1 -1 "" 0 1 0 0 255 1 0 
0 0 0 0 0 1 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Outpu
t" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 256 1 {CSTYLE "" -1 -1 "" 0 1 0 
0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 
0 257 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 260 1 {CSTYLE "" -1 -1 "" 0 1 0 
0 0 0 0 1 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 
0 261 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }0 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 262 1 {CSTYLE "" -1 -1 "" 0 1 0 
0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 
0 263 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 33 "Local Stability Worked Exa
mple 1." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 29 "Consider the nonlinear system" }}{PARA 
256 "" 0 "" {OLE 1 4096 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N:F:nyyy
yy]::yyyyyy:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
::::::::::::fyyyyya:nYf::G:jy;::::::::::::::::::::::::::::::::::::::::
:::::::::::::::::::::::::::::::::::::::JcvGYMt>^:fBWMtNHm=;:::::::n<Ej
yyiy=C:Qj;>:;JZC:bKB::>fmytvIu<oq;V:>j>Z:j:vCSmlJ::::::::::OJ;@jyyyyyy
;jysy;Z:::::::^<>:F::]KRnC=MtFGgml>:;::::::JGN:ry:>:<::::::=J:^J>:F:Al
qfG[maNFO=;::::::::_J;vyyuy:>:<::::::AZ:^E>:nYN:<j>J<JyK=j=<Z:>;;J?>:w
AQ:S:UJ:n;v;;JB<Z:F<nYvY::::::::::::::::::::::::::::::::::::::::::::::
:::::::::::::::::::::::::::=Z<B:<::::::::::;^:^Z;r:Eb<CR:=R:_sZke;?J;T
rXDk;]Q;TrODk;aQ;TrPDk;WZIj?@k:DJ>E=dj:Dm;aQ=D:NZ@h?N:J:b=D:@Jb<r[@PqH
pll@Ij=>:;:::::=B:wyyAb<vYxy?V:>:;P:<:TNC>:UTRcETcTX[US;SK]UW=EWMuUWm>
j\\@A:::cYH_WV;J:<::::::::yayY:^Z:j:>:yayQ:>:neCJuM:_h:F:MZ=Fg;J;f<>:E
J:F[<J:<J<B:?R:=j<j:Djyyyy=j<j>Z:^:n_;>=E:]c:=Z:f:V[<Z:f:^[<>j>>jBn;n^
@v;_jysy;B:K:;xyiwy;D<iC;Q:G;Sj`@@J:>XEZ:F;NZ:vCS=[LsfFaMR>@>Z::::::::
kJ;@:;J;>Z:vYxY:>Z::::::j<B:oi:B:N:Y<>D_=a=[;;B:::::::JF>:yay=J:B:::::
:nYyA<::::::::::::jysy:>:<:::::::::::::::::::vYxI:;Z:::::::::]:qi:;fyB
:>l;B:DZ<>ZJfcN?^??inOh[SIC:US:F[:>Z:N`DFl;>H=B:<JM<J:=j>J:^[:^IHN;sJB
:;::::::::N@]x:gm:JMjDB:qQ:[:JBAj:J:D:cTTUUSaEBWTSiEB_tUUURWmtAJ<jPF:C
:[Y:F;;JSfm;^C=Z:>IOo>JSJKEj:>:k=JSjdMj:>:mMO;_;oi;=J:>I:_khIJfF:FI:_;
ev:=B:;Jv:N@>s?F:FI><<jw;<:[V:;Z<b:<kc`a\\wgfc?C:UK:^:>X=j>>:_c<uNh@j:
<J:<jXVTM:_;qv:=:Qkc;_;UX;=:UC:a@N@uPn@j:jX:N@nm<F:;jP:N@nk>F:V?f<VY;s
yB:>l;Z:b:^DP@::jP@j:^:>X;j>>:_krwqAF:<Z:>HOB:<JMT:OjyyiyI:>:[Z:VY[j=B
:;JXE:<j:DJ;Dlc`qsLqlp`h_:f?;J<JrG:AB:<J:jS>HOJ:n>jysy=Z:>:[Z:VYv:<J:^
QB:=:::Uk:^:>X;j;<Z:>:?bBa<:f?:[AJ::::::::::::::::::::::::::::::::::::
:::::::::::5:" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 105 "
1) Find the equilibrium points of the system and classify the equilibr
ium points of their linearisations." }}{PARA 0 "" 0 "" {TEXT -1 70 "2)
 What does ths tell you about the stability of the nonlinear system." 
}}{PARA 0 "" 0 "" {TEXT -1 90 "3) Verify the answer to (2) by looking \+
at the phase plane diagram of the nonlinear system." }}{PARA 0 "" 0 "
" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 256 31 "To find the equilibriu
m points." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(linalg):" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "with(DEtools):" }}{PARA 7 "" 1 "" 
{TEXT -1 35 "Warning, new definition for adjoint" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 11 "f:=x+y-x^2;" }{TEXT -1 0 "" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 8 "g:=-x+y;" }{TEXT -1 0 "" }}{PARA 0 "" 0 "
" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "solve(\{f
=0,g=0\},\{x,y\});" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 40 "The equilib
rium points are (0,0), (2,2)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT 257 27 "To find the Jacobian matrix" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "v:=[f,g];" }{TEXT -1 0 "" }}
{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 21 "J:=jacobian(v,[x,y]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
0 "" }}}{EXCHG {PARA 260 "" 0 "" {TEXT -1 18 "At the point (0,0)" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "J0:=subs(\{x=0,y=0\},evalm(J
));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 
257 "" 0 "" {OLE 1 4608 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N:F:nyyy
yy]::yyyyyy:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
::::::::::::fyyyyya:nYf::G:I:wAyA:::::::::::::::::::::::::::::::::::::
::::::::::::::::::::::::::::::::::::::::::NDYmq^H;C:ELq^H_mvJ::::::::g
j<vyyuyJ<j?V:;J:>R<:T><::;uviHymhJwW:A:;f;<Z:j:vCSmlJ::::::::::OJ;@jyy
yyyy;jysy;Z:::::::^<>:F::]KRnC=MtFGgml>:;::::::JGN:ry:>:<::::::=J:^V>:
F:AlqfG[maNFO=;::::::::_J;vyyuy:>:<::::::AZ:^E>:nYN:<j>J<JyK=j=<Z:>;;J
?>:a:QJ:^;f;;JAjA>:[B:<jBJC>:c:wAeJ:n<;jEJFjFJyky;::::::::::::::::::::
:::::::::::::::::::::::::::::::::::::::::::::::=Z<B:<::::::::::;a:^Z;r
:Eb<=b>AFW_rZ_u:=J:@JCHR_DJ>=]<nlbNZPH??b[uc:KfF?bZSc:AFS?rZEe:AFWAR::
:::N\\=`hLbV:::HjAr<]R:?B\\_C:dJ;TrPDJ>MM;Hr:nyyYZDjysyQj;J:>R<:TNC>Z:
fc[_hb_ds?h_?^?GhoGfnGgioGM:ETV:::^xsNpkK:<::::::::yayY:^Z:j:<:yayQ:>:
FpYJlE:cE;F:MZ=Na<J;<jDJ:f:F[<J:<J<B:?R:=j<j:Djyyyy=j<j>Z:^:n_;>=E:]c:
=:E:Qb:B:E:Sb:;g:;HAW:Gc;YJCvYxY:<J>JZyyNyY:mk^Dj?B:G;Sj`@@J:>FEZ:F;;J
;<ja^G>D_mlVH[KRJ:<:::::::>=?R:>:?J:<jysy;J:<::::::EZ:^m>Z:J;vCJbNHVH>
@>Z::::::::kJ:vYxI:;Z::::::JywYB:::::::::::::yay=J:B::::::::::::::::::
:jysy:>:<::::::::f<<jwEJ:uI<:[V:=:;b:<Kd`ppPps<LaPpu<Lcxpp@Pq^iP_:f_;J
?<Z:N`D>v;^@?JlDj;Q`=`h=J:@JC:DJE;?:GVEB:pneHoY>:;::::ceVQmDG;;B:_;WY:
=Z:NCOO@>Z:N@fr@V:^f__H;POWGM:_;oT<=:O<O:_;]c=?:c=:_;Qy==:O<SZ:N@vpNV:
;Jl:Ooo:_;St?=:O<e:qi:Z:JBAJ:<Z:b:;B>:Zi;eJ<B:Uk:^:>x;j><J:N@>r<FZ:>Z:
>@;p>JSjwPj:Z:>RZDJPdj:JR:N@vvFF:B::d:uG>=:[;JSjrIk:JT>RJSJI\\k:JR:N@^
pUF:^@f<;jw?^y]:JBA:<Z<Jd::Jx:UK:^:>X=j>>:_kywyKF:;B:<Jr><<jw?Zw]:JBA:
<Z<B:GTSWuREuVTK<jPF:C:[Q:F;;JSdjvAj\\@j:>Z:NFwm>JSNy;Vi=F:vFw=N@mqR]j
:Jk:N@oqmej:jm:N@mQ[xj:Jk:N@oQv@k:J:::::::::::jm:N@mqKUk:>:?=JSNYmH@=:
I=[Z:VY;><:[V:b:DZJVdscRYEUTK<jP>:C:[q:F;;JSdJjgc>F:;B:ukcG;N@FqEF:fA:
_;sh>=J:f_:VTJSjOak:jX:N@FvTF:f?f<VYZ:JBA:<Z<JdRS:x@f<C:US:=J<Jr?:MJ:N
@Fb@F:<JP>FM:_;UW==:S;O:G;OjyyiyI:>:[Z:VY[j=B:;JXE:;B:=b:?bBaTXaEWEUUt
P<jP>:C:[q:VZ:B:;:::::::::::::::::::::::::::::::::::::::::::::::::::::
:::::::::::::::::5:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 99 "The eigenv
alues are imaginary with positive real parts and the linearisation is \+
an unstable spiral." }}{PARA 0 "" 0 "" {TEXT -1 119 "By the linearisat
ion theorem the nonlinear system mimics this behaviour and the point (
0,0) is a repeller and unstable." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }
{TEXT 258 0 "" }{TEXT -1 0 "" }}}{EXCHG {PARA 261 "" 0 "" {TEXT -1 18 
"At the point (2,2)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "J1:=
subs(\{x=2,y=2\},evalm(J));" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}
{EXCHG {PARA 262 "" 0 "" {OLE 1 4608 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j:
:>:B>N:F:nyyyyy]::yyyyyy::::::::::::::::::::::::::::::::::::::::::::::
:::::::::::::::::::::::::fyyyyya:nYf::G:I:wAyA::::::::::::::::::::::::
:::::::::::::::::::::::::::::::::::::::::::::::::::::::NDYmq^H;C:ELq^H
_mvJ::::::::gj<vyyuyJ<j?V:;J:>R<:T><::;uviHymhJwW:A:;f;<Z:j:vCSmlJ::::
::::::OJ;@jyyyyyy;jysy;Z:::::::^<>:F::]KRnC=MtFGgml>:;::::::JGN:ry:>:<
::::::=J:^V>:F:AlqfG[maNFO=;::::::::_J;vyyuy:>:<::::::AZ:^E>:nYN:<j>J<
JyK=j=<Z:>;;j?>:OJ:nY^;f;;JAjA>:[B:<jBJC>:a:c:e:gJ:v<>=F=nYvY:::::::::
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::j:b:<Z::::::
::::JjC:CR:Hj<dj:Dk;MQCHRsHj:>Z;N\\=Pc<>kkb:gFD?b_MK;TrXDJ=qL;\\RSDJ;m
LCDRlDK;QPCDRtHj::::::HRt:::Z=v;hjB@J;\\RS<ZDNZ@H:wyyAb<vYxy?V:>:;`:Z@
O<;B:UTRcETcTX[US;SK]UW=EWMuUWm>j\\@A:::cYH_WV>Z::::::::jysyA:CB:FZ:jy
sy?:;:?EB^v<J>M:=j>r:Cd:N:e:;j<>:Mb:>Z:^:<J;@j:f:FZ<vyyyy:f:F;B:C:WS:k
j<jRDj:j<j?D:<j<J@DJJEJjpKAJMTjAN<yayAZ:>;;JZyyNyYJpk^Dj?JMJ@fc[;>:;m<
B:M:?B:YLpJbNHEms>@[;;B:::::::JFNZ;J:N:;B:yayA:;B::::::f:<JRM:<:?ja:[L
sjsJRJ:<:::::::>=;jysy:>:<::::::wqy[:::::::::::::vYxI:;Z::::::::::::::
::::::yay=J:B::::::::jDB:qi:;fyB:>l;F:>Z<B>cTTUUSaEBWTSiEB_tUUURWM`^:f
_;J?<Z:N`D>v;^@?jPDJ>MM;HRlDj;QP;HRtHj:R:_::J;J>AM;\\B>:::::::::::::::
ceVQmDG;;B:_;WY:=Z:NCOo>JSjvaj:>:O<O:_;qu=?:c=:_;eS>=:O<SJ:<JSjbQk;JlP
Pt>r:NOQm>JSJm\\k:J_jDjw;<:[V:;B:<Z<B>:Zi;eJ<jPF:C:[Y:F[:>:_;Kg:=J:<jR
<JZG;N@n^@F:;Z:>RZDjSej:jR:N@F`HF:B::d:YY>=J:F@:_;QR?=:cKZ;_;;S@=:];O:
_;Qw@?:]cF;H:;pD>:qQJxI<:[V:B:D:c<::SL<jP>:C:[q:F[:>:_kyG]MF:;B:<Jr><<
jw?Zw]:JBA:<Z<B:GTSWuREuV`K<jPF:C:[Q:F;;JSdjvAj\\@j:>Z:Nf:nIM:_KwAj_Hj
:jmnIJSFY[T<=:?=JSNYgt<=J:vF:_kv_kJFZ:Jk:N@oQCEk:>:I=JSFYuw?=:?=JSNYAX
@=:I=[Z:VY;>;B:>l;Z<b:<kc`a\\wgfS@C:UK:^:>X=j>>:_c<;qHLj:>Z:fAap>JSjcP
j:jP<jc;_;KC==:u;JSJipj:jP:N@nxMF:f?:_;ef@=:u;e:qAB:>l;Z:b:^DP@:jD^:f_
;F:C:[Q:F;;JSJpUj:B:SKj^[:B:_;[u;AB:;JlZ:NOQ=N@vxHF:^?N;n>N;yyyxy:J:><
<jw?<IZ:>:sg:>Z:FZ<NZ\\VdsWhngfgcUC:UK:^:>X=j;<J:<J:::::::::::::::::::
:::::::::::::::::::::::::::::::::::::::::::::::::1:" }}{PARA 263 "" 0 
"" {MPLTEXT 0 21 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 77 "The eigen
values are real and negative and the linearisation is a stable node." 
}}{PARA 0 "" 0 "" {TEXT -1 117 "By the linearisation theorem the nonli
near system mimics this behaviour and the point(2,2) is an attractor a
ndstable." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT 259 26 "To plot the phase portrait" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 24 "A:=diff(x(t),t)=x+y-x^2;" }}{PARA 11 "" 1 "" {TEXT 
-1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "B:=diff(y(t),t)=-
x+y;" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 183 "DEplot(\{A,B\},\{x(t),y(t)\},-50..50,[[x(0)=1,y(0)=1
],[x(0)=.1,y(0)=.1],[x(0)=1.5,y(0)=2],[x(0)=2.6,y(0)=2.1],[x(0)=1.5,y(
0)=1.5],[x(0)=2.2,y(0)=2.2]],x=-.5..2.5,y=-0.5..2.5,stepsize=.1);" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 46 "Th
e phase portrait supports the above analysis" }}{PARA 0 "" 0 "" {TEXT 
-1 0 "" }}}}{MARK "28 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }