{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 "" -1 256 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 257 "" 0 1 187 0 152 0 1 0 0 0 0 0 0 0 0 }{PSTYLE "Norma l" -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 "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 1 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 } {PSTYLE "" 0 258 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 259 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 74 0 63 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 63 1 100 0 0 1 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 16 "When p=2 and q=4" }}{PARA 0 "" 0 "" {TEXT -1 20 "The equations become" }}{PARA 256 "" 0 "" {OLE 1 4106 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy:::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::fyyyy ya:nYf::G:jy;::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ::::::::::::::::::::::JcvGYMt>^:fBWMtNHm=;:::::::n:;JZC: bKB::>\\VOpZNgZpq;V:>b>B:<:=ja^GE=;:::::::::N;?R:yyyyyyA:yayA:<::::::J DJ:j::F@[KaFFcmnnHEM:>:::::::oJ;Zy=J:B::::::F:;Jp>:;:=j[vGUMrvC?MoJ::: :::::JCN:yyyxI:;Z::::::j;B:s<;:wA?Z::C:J=j=B:K:M:OJ:V;;J@j@>:W:YJ:>\\: B:wAyA:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: :::::::::F:DZ:B::::::::::>G^Z;r:EbAFX?b[U[lN\\=Hi:>:;r;R:[DZGcMB:::::JvYkg>pSVE;Y>q<;::j::yayQ:>:>]CjQM:sf:F:MZ=nc;J;f<>:EJ:F[Z:^:n_;>=E:]c:=Z:f:V[b>>bBn;n^@v;_jysy;B:KJ:>ryysyEVEZ:F;NZ:vCSml>D_mlVH[KRJ:<:::::::>]:NZ;J:N:;B:yayAZ:>Z:::: ::jD_=a=[;;B:::::::JF>:yay=J:B::::::nYyA<::::::::::::jysy :>:<:::::::::::::::::::vYxI:;Z:::::::::]:qi:;fy>Z:JBAJ:UDO;SKk EW_URdLf<l;F:>Z>C:Uk:^:>x;j>>:_kb@J`JSJQEj:>:mMN;_kBIJfF :FI:_;cD;=J:>I:_;Qt;=:m=[B:C:UK:^:>X=j>>:_c< ]Lh@j:UM:_kB?g;F:fA:_;wf:=J:f_:>UJSJ`Pj:jPjDjw;<:[V:B:D:c\\_;:: f_;F:C:[Q:F;;JSJuHj:RM:_;=e;=J:B:;P?B:Gc;N;yyyxy:J:><:C:[q:VZ:B:;:::::::::::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::j:NM;Bv3:" }}} {EXCHG {PARA 257 "" 0 "" {TEXT -1 30 "To find the equilibrium point." }}{PARA 0 "" 0 "" {TEXT -1 21 "This is clearly (0,0)" }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}}{EXCHG {PARA 258 "" 0 "" {TEXT -1 33 "To classify t he equilibrium point" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "wit h(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 32 "Warning, new definition for norm" }}{PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition for trac e" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "A:=matrix(2,2,[0,1,-4, -2]);" }}}{EXCHG {PARA 259 "" 0 "" {TEXT -1 35 "tr(A)= -2, det(A)=4,di scrim(A)= -12" }}{PARA 0 "" 0 "" {TEXT -1 102 "From the tr-det diagram the equilibrium point is a stable spiral (focus). There will be no d irectrix." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 260 "" 0 " " {TEXT -1 35 "To plot the direction field diagram" }}}{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 18 "P:=diff(x(t),t)=y;" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "Q:=diff(y(t),t)=-4*x-2*y; " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "DEplot([P,Q],[x(t),y(t)],-2. .5,x=-5..5,y=-5..5);" }}{PARA 13 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 261 "" 0 "" {TEXT -1 26 "To draw the phase portrait" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 111 "DEplot([P,Q],[x(t),y(t)],-2..5,[[x(0)=4,y(0)=0],[x(0)=-4,y(0)=0 ]],x=-5..5,y=-5..5,linecolor=black,stepsize=.1);" }}{PARA 13 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 133 "That the solutio ns oscillate with decreasing amplitude as they tend to the rest positi on can be seen by looking at the time series of" }{TEXT 256 2 " x" } {TEXT -1 9 " against " }{TEXT 257 2 "t." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 108 "DEplot([P,Q],[x(t),y(t)],-1..6,[[x(0)=0,y(0)=5]],x=- 1..1.5,y=-5..5,linecolor=black,stepsize=.1,scene=[t,x]);" }}{PARA 13 " " 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} }{MARK "17" 0 }{VIEWOPTS 1 1 0 1 1 1803 }