{VERSION 2 3 "IBM INTEL NT" "2.3" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{CSTYLE "" -1 257 "" 0 1 74 0 68 0 0 1 0 0 0 0 0 0 0 } {CSTYLE "" -1 259 "" 0 1 189 0 244 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 260 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 261 "" 0 1 181 0 148 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 0 1 182 0 120 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 263 "" 0 1 61 1 176 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 264 "" 0 1 68 1 196 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 265 "" 0 1 60 1 120 0 1 0 0 0 0 0 0 0 0 }{PSTYLE "Normal" -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 0 }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 Output" 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 "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 13 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 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 258 1 {CSTYLE "" -1 -1 "" 0 1 181 0 188 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 259 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 260 1 {CSTYLE "" -1 -1 "" 0 1 182 0 24 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 261 1 {CSTYLE "" -1 -1 "" 0 1 93 0 196 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 262 1 {CSTYLE "" -1 -1 "" 0 1 63 1 240 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 48 "Stability of equilibrium p oints investigation 1." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 9 "Consider " }}{PARA 256 "" 0 "" {OLE 1 3584 1 "[xm]Br= WfoRrB:::wk;nyyI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::::::::fyyyyya:nYf::wyyyqy;:: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ::::::::NDYmq^H;C:ELq^H_mvJ::::::::gjR<:T><::XMxty@iSY Aj;J:L:<:=ja^GE=;:::::::::N;?R:yyyyyyA:yayA:<::::::JDJ:j::F@[KaFFcmnnH EM:>:::::::oJ;Zy=J:B::::::F:;JnJ:j:VBYmp>HYLkNG>::::::::Nj>J?>:QJ:^;f;;JAjA>:wAyA::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::=ZkmN\\=d:GVE?B\\_c>AVX@JSHrvTJ;ZLVjr?:D] Hr;R:[DZGcMB:::::JFIpgJ::::j::yayQ:>:nyAJ@M:Of:F:MZ=f b;J;f<>:EJ:F[Z:^:n_;>=E:]c:=B:;;JZyYjyiAB;Q:G;Sj`@@J:>DEZ:F;NZ:vCS=[LsfFaMR>@> Z::::::::kJ;@:;J;>Z:vYxY:>Z::::::jD_=a=[;;B:::::::JF>:yay =J:B::::::nyyyYB:::::::::::::yay=J:B:::::::::::::::::::jysy:>:<::::::: :Fjx]:JBAJ:UDO;SKkEW_UR\\M<<:US:F[:>Z:N`D^`;>H=B: b;^C=Z:>IKG:M:_;[t:=B:;Jv:N@vb@F:;jv>OJS>j=>E=:m=JSJQEj:>:k=[B:X=j>>:_c<;Lh@j:>Z:fAkp>JSJ^Lj:>:UC:k@N@;Nn@ j:jXjDjw;<:[V:B:D:c\\_;::f_;F:C:[Q:F;;JS>BOf:=J:RM:_;uT;=::[Z:VY[j=>Z:>:sg:B:=b:?bBaTXaEWEUUtP:C:[q:VZ:B:;:_c<;v: gf:=Z:V@1:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 257 43 "To find the equili brium point of the system" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 29 "The equlibrium point is (0,0)" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} {EXCHG {PARA 258 "" 0 "" {TEXT -1 23 "To find the eigenvalues" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 13 "with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 32 "Warning, new definition for norm" }}{PARA 7 "" 1 "" {TEXT -1 33 "W arning, new definition for trace" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "A:=matrix(2,2,[3,-2,1,0]);" }}}{EXCHG {PARA 260 "" 0 "" {TEXT -1 3 "tr(" }{TEXT 262 4 "A)=3" }}{PARA 261 "" 0 "" {TEXT -1 4 "det(" }{TEXT 263 4 "A)=2" }}{PARA 262 "" 0 "" {TEXT -1 8 "discrim( " }{TEXT 264 1 "A" }{TEXT -1 7 ")=9-8=1" }}{PARA 0 "" 0 "" {TEXT -1 49 "Thus there will be two real positive eigenvalues." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "eigenval ues:=eigenvals(A);" }}{PARA 11 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 69 "There are two real eigenvalues so there will be two real directrices." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 259 23 "To find the directrices" }}{PARA 0 "" 0 "" {TEXT -1 48 "To find these you need to find the eigenvectors." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "Eigensystem:=eigenvectors(A); " }} {PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 259 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 28 "The directrices are given by " }}{PARA 257 "" 0 "" {OLE 1 5120 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>: B>N:F:nyyyyy]::yyyyyy::::::::::::::::::::::::::::::::::::::::::::::::: ::::::::::::::::::::::fyyyyya:nYf::G:I:K:wAyA::::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::JcvGYMt>^:fBWMtNH m=;:::::::n:;JZC:bKB::>^RYjxikUpq;V:>ZCZ:j:vCSmlJ::::::: :::OJ;@jyyyyyy;jysy;Z:::::::^<>:F::]KRnC=MtFGgml>:;::::::JGN:ry:>:<::: :::=J:^_;Z:j:VBYmp>HYLkNG>::::::::N;^;N;;j?>:wAU:W:YJ:>\\:B:]:_J:V<^:q:s:u:w:y:;C:< jJ<:?;A;wAyA:::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ::j:b:B::::::::::>_;J<@Z=fZD^Z;FZ;VZ=j:Pj:@j:Dj:R:_rZkm:D[s[AvTAr:R;:NZBH@=b:B:;::Qs:A:E:j:DZjRJ;P;\\RR>Z:NZP<:^cEjdd:Sh;F:MZ=nf=J;D:;>ryYnye^ ;UTR:;JJDj@[;;B:::::::JFNZ;J:N:;B:yayA:;B::::::f:D_=a=[;;B:::::::JF>:yay=J:B::::::nYyA<::::::::::::jysy:>:<: ::::::::::::::::::vYxI:;Z:::::::::eB:klR:?jR B::rDZ?:R:?JRF:B:;b::;B:;::::::::::JDDpg^\\\\K:>::::kMN=j><:_khIJvGZ:> Z:FIK?N@;h:Ay:=J:VGK?N@ud;_q:<:k=JSfmB>Y=:m=JS>jAVv=F:VG><r;>B=B:B:_kYLJZF:nC[n>B:_ kRHJZF:fc:>LMZ:N@;pUAj:>:YLB=j>B:_kYvp;F:Fd:>LMZ:N@]kUAj:Jb>L[Z:VY;><: [V:b:DJ:C:UK:^:>X=j>>Z:N`D>rl;Z:b:<: :jPF:C:[Q:F;;JSR:=H;=J:B:_kYFr>F:nCF;F:fCF;RCY;=J:vC F;D><l;ZZ:V?k@l;Z:b:<::jPF:C:[Q:F;;JSR:=y;=Z:^CF;B:_kY^`CF:FD[n>B:_kR^`CF:>DF;j@>B=B:^C F;B=:WB:_K:wP=J:vCF;jq:f?;JZ:fAZ:VY;:<:[V:B:DZ:::f?=J B:_krGr>F:fCF;J;c;=:YB:_kr?^@F:>D><l;ZX;j>>:_[@f]BF:B:_kyg]BF:nCF;B:_ky?hCF:FDB:_kr?hCF:>Df\\:B:qQ:[:JBA:c\\_;ND:jP> :C:[q:F;;B:_kHAjXPj:>Z:F`:>RM:_khga?F:F@:_kHghBF:>@;H:JSfUeE<=:[;JSfe@ fd?F:B:;@bN[:vyyuy=:;JBB:qQBv: " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 260 27 "To plot the direction field" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 14 "with(DEtools):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 18 "Enter the equation" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "P :=diff(x(t),t)=3*x-2*y;" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "Q:=diff(y(t),t)=x;" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "DEplot ([P,Q],[x(t),y(t)],-2..2,x=-5..5,y=-5..5);" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 13 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 82 "To fill in the initial conditions use the four on the dir ectrices plus some others" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 200 "DEplot([P,Q],[x(t),y(t)],-2..2,[[x(0)=1,y(0)=1],[x(0)=-1,y(0)=-1] ,[x(0)=2,y(0)=1],[x(0)=-2,y(0)=-1],[x(0)=0,y(0)=4],[x(0)=4,y(0)=0],[x( 0)=-4,y(0)=0],[x(0)=0,y(0)=-4]],x=-5..5,y=-5..5,linecolor=black);" }} {PARA 13 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 269 "Look at the direction of the arrows. As the eigenvalues are positive \+ all the trajectories move away from the equilibrium point. There are t wo directrices both with arrows pointing away from the equilibrium poi nt. Thus the equilibrium point is unstable and is called an " }{TEXT 261 13 "unstable node" }}}{EXCHG {PARA 12 "" 1 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 23 "Now consi der the system" }}{PARA 263 "" 0 "" {OLE 1 4096 1 "[xm]Br=WfoRrB:::wk; nyyI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy::::::::::::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::fyyyyya:nYf::G:jy;:::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::JcvGYMt >^:fBWMtNHm=;:::::::n:;JZC:bKB::>^RYjxikUpq;V:>b>B:<:=ja ^GE=;:::::::::N;?R:yyyyyyA:yayA:<::::::JDJ:j::F@[KaFFcmnnHEM:>:::::::o J;Zy=J:B::::::F:;JnJ:j:VBYmp>HYLkNG>::::::::Nj>J?>:QJ:^;f;;JAjA>:[B::yayQ:>:FtBJ@M:sf:F:MZ=nc;J;f<>:EJ:F[Z:^:n_;>=E:]c:=B:ryAvyArAB;Q:G;Sj`@@J:>DEZ:F;NZ:vCS=[LsfFaMR>@>Z::::::::kJ;@:;J;>Z:v YxY:>Z::::::jD_=a=[;;B:::::::JF>:yay=J:B::::::nYyA<:::::: ::::::jysy:>:<:::::::::::::::::::vYxI:;Z:::::::::]:qi:;fyB:>l;>Z:b:DJ: JS^h=nH=:G;eZ:VY;><:[V:=: ;b:^dcgg_WhZnc_whZNdigg[oG[K<>:Uk:^:>x;j>>:_KZ@J`F:o:F;N@>`>F:;Jv: N@v[AF:FIK?N@;v:kl:jv:N@no:UK:^:; JrG:MJ:N`D>BsT:=J:UM:_;_f:=J:f_:>UJSJNPj:jP:N@;Nn@j:jXjDjw;<:[V:B: D:c\\_;::f_;F:C:[Q:F;;JS>BOu:=J:RM:_;uc;=J:V@N;:[Z:V Y[j=>:;JXE: